<div dir="ltr">Hi Pedro,<div style>Does the time look so bad? That number of seconds corresponds to 13.88 Gyrs. As for the KeyError, most codes use pf["Time"] to indicate the conversion factor between code units and time units. I've never filled this out, since as you point out, there's a complicated nonlinear relation between code time units and seconds. Alternatively we could write pf.conversion_factors["Time"]=1.0 and avoid this error.</div>
<div style><br></div><div style>Nevertheless, you can check out how the code time -> seconds time works here in the b2t function:</div><div style><a href="https://bitbucket.org/juxtaposicion/yt/src/b529d44c0de571a1d5afb5695b9a1501052428d6/yt/frontends/art/io.py?at=yt#cl-550">https://bitbucket.org/juxtaposicion/yt/src/b529d44c0de571a1d5afb5695b9a1501052428d6/yt/frontends/art/io.py?at=yt#cl-550</a> The above function (and the set of functions it calls) is translated from some of Daniel Ceverino's ART-analyzing Fortran code. Because they involve finding the zero-point of a function that in turn requires a numerical integration, the routine can be quite slow, despite using interpolation. Nevertheless, I think it still gives correct answers.</div>
<div style><br></div><div style>Thanks!</div><div style>chris</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Feb 19, 2013 at 2:22 PM, Pedro Colin <span dir="ltr"><<a href="mailto:p.colin@crya.unam.mx" target="_blank">p.colin@crya.unam.mx</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
I have downloaded several ART files and parameters like current redshift or<br>
omega lambda are okay but the current time = 4.38272802952e+17 looks very<br>
suspicious. This may be the reason why I got an error when I try to print_stats<br>
<br>
raise KeyError(key)<br>
KeyError: 'Time'<br>
<br>
I also think this error is the reason why I got wrong units. By the way, I<br>
know ART time code unit is negative and reached zero at the present time. Any<br>
help?<br>
<br>
Pedro Colin<br>
<br>
<br>
_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</blockquote></div><br></div>