<div dir="ltr"><div style>Hey Guys,</div><div style><br></div><div style>      I've been using yt's star analysis module to produce plots of SFR vs. time to compare the star formation histories across multiple runs. It works well except when it doesn't. </div>
<div style><br></div><div style>     I've been working with enzo and for late time files (I believe produced after a run was restarted) the star formation history cuts off and drops to 0 for the earlier times. Weirder still, when I open earlier data files and plot them atop, the earlier results disagree! I attached a plot of this below.</div>
<div style><br></div><div style>    Does anyone know what's happening here? Is there a better way to do this?</div><div style><br></div><div style>     best,</div><div style>        Munier</div><div style><br></div><div style>
<br></div><div style>The Code:</div><div><br></div><div><br></div><div><div><font face="courier new, monospace" size="1"><b>import matplotlib as mpl</b></font></div><div><font face="courier new, monospace" size="1"><b>mpl.use('agg')</b></font></div>
<div><font face="courier new, monospace" size="1"><b>import matplotlib.pyplot as plt</b></font></div><div><font face="courier new, monospace" size="1"><b>from yt.mods import *</b></font></div><div><font face="courier new, monospace" size="1"><b>from yt.analysis_modules.star_analysis.api import *</b></font></div>
<div><font face="courier new, monospace" size="1"><b><br></b></font></div><div><font face="courier new, monospace" size="1"><b>def SFR(fName,stroke):</b></font></div><div><font face="courier new, monospace" size="1"><b><br>
</b></font></div><div><font face="courier new, monospace" size="1"><b>  # open file</b></font></div><div><font face="courier new, monospace" size="1"><b>  try:</b></font></div><div><font face="courier new, monospace" size="1"><b>    pf=load(fName)</b></font></div>
<div><font face="courier new, monospace" size="1"><b>  except YTOutputNotIdentified:</b></font></div><div><font face="courier new, monospace" size="1"><b>    return False</b></font></div><div><font face="courier new, monospace" size="1"><b><br>
</b></font></div><div><font face="courier new, monospace" size="1"><b>  # calculate SFR</b></font></div><div><font face="courier new, monospace" size="1"><b>  dd = pf.h.all_data()</b></font></div><div><font face="courier new, monospace" size="1"><b>  sfr = StarFormationRate(pf,data_source=dd)</b></font></div>
<div><font face="courier new, monospace" size="1"><b><br></b></font></div><div><font face="courier new, monospace" size="1"><b>  plt.plot(sfr.time,sfr.Msol_yr,stroke,linewidth=1)</b></font></div><div><font face="courier new, monospace" size="1"><b>  return ( sfr.time.min() , sfr.time.max() )</b></font></div>
<div><font face="courier new, monospace" size="1"><b><br></b></font></div><div><font face="courier new, monospace" size="1"><b>MaxVals = []</b></font></div><div><font face="courier new, monospace" size="1"><b>SFR("../FIDUCIAL/DD0005/test_sim_0005","k-")</b></font></div>
<div><font face="courier new, monospace" size="1"><b>SFR("../FIDUCIAL/DD0010/test_sim_0010","k-")</b></font></div><div><font face="courier new, monospace" size="1"><b>SFR("../FIDUCIAL/DD0017/test_sim_0017","k-")</b></font></div>
<div><font face="courier new, monospace" size="1"><b><br></b></font></div><div><font face="courier new, monospace" size="1"><b>plt.xlim((1.27E10,1.30E10))</b></font></div><div><font face="courier new, monospace" size="1"><b>plt.savefig("SFR")</b></font></div>
</div><div><br></div><br clear="all"><div><br></div>-- <br>Munier A. Salem // 845.489.6450
</div>