import numpy as np

currenttime = 10.0

volume=138.87
tunit = 6.88027e+14/3.15569e7 
covertion = 0.08*tunit*volume

fi=open('PopII_evolution.txt','r')
time = []
SFR = []
line=fi.readline()
while line:
   time.append(float(line.split()[1]))
   SFR.append(float(line.split()[4])/covertion)
   print float(line.split()[1]),float(line.split()[4])/covertion
   line=fi.readline()

for t,s in zip(time,SFR):
  if t < currenttime:
     print s+(old_s-s)/(old_t-t)*(currenttime-t)
     break
  old_t = t; old_s=s