import numpy as np

f=open('PopII_list_0038.txt','r')

fo = open('PopII_evolution.txt','w')

currenttime = 11.65522741145

t_min = 0
t_max = 8.0

bins = 100
dt = (t_max-t_min) /bins

line = f.readline()
line = f.readline()

time_bins = np.zeros(bins)
mass_bins = np.zeros(bins)

while line !='':
   mass = float(line.split()[1])
   ct = float(line.split()[2])
   index = int((currenttime-ct)/dt) 
   time_bins[index] = time_bins[index]+1
   mass_bins[index] = mass_bins[index]+mass
   line = f.readline()

total_stars = time_bins.sum()
total_mass = mass_bins.sum()
for i in xrange(bins):
   fo.write(str(i*dt)+" "+str(currenttime-i*dt)+" "+str(time_bins[i])+" "+str(total_stars)+" " \
    +str(mass_bins[i])+" "+str(total_mass)+"\n")
   total_stars = total_stars - time_bins[i]
   total_mass = total_mass - mass_bins[i]
   if total_mass < 1: total_mass = 0.0

f.close()
fo.close()

