<div dir="ltr">I forgot to mention, the X-axis is time and the Y-axis is SFR (M_solar / yr)</div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Mar 5, 2013 at 1:41 PM, Munier Azzam Salem <span dir="ltr"><<a href="mailto:msalem@astro.columbia.edu" target="_blank">msalem@astro.columbia.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Hey Guys,</div><div><br></div><div>      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><br></div><div>     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><br></div><div>    Does anyone know what's happening here? Is there a better way to do this?</div><div><br></div><div>     best,</div><div>        Munier</div><div><br></div><div>
<br></div><div>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>
<span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888"><div><br></div><br clear="all"><div><br></div>-- <br>Munier A. Salem // <a href="tel:845.489.6450" value="+18454896450" target="_blank">845.489.6450</a>
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Munier A. Salem // 845.489.6450
</div>