[Yt-svn] yt-commit r1150 - trunk/yt/lagos

mturk at wrangler.dreamhost.com mturk at wrangler.dreamhost.com
Thu Jan 29 12:28:25 PST 2009


Author: mturk
Date: Thu Jan 29 12:28:23 2009
New Revision: 1150
URL: http://yt.spacepope.org/changeset/1150

Log:
Fix for dynamical time calculation



Modified:
   trunk/yt/lagos/UniversalFields.py

Modified: trunk/yt/lagos/UniversalFields.py
==============================================================================
--- trunk/yt/lagos/UniversalFields.py	(original)
+++ trunk/yt/lagos/UniversalFields.py	Thu Jan 29 12:28:23 2009
@@ -313,15 +313,9 @@
     M{sqrt(3pi/(16*G*rho))} or M{sqrt(3pi/(16G))*rho^-(1/2)}
     Note that we return in our natural units already
     """
-    return data["Density"]**(-1./2.)
-def _ConvertDynamicalTime(data):
-    G = data.pf["GravitationalConstant"]
-    t_dyn_coeff = (3*pi/(16*G))**0.5 \
-                * data.convert("Time")
-    return t_dyn_coeff
+    return (3.0*pi/(16*G*data["Density"]))**(-1./2.)
 add_field("DynamicalTime", function=_DynamicalTime,
-           units=r"\rm{s}",
-          convert_function=_ConvertDynamicalTime)
+           units=r"\rm{s}")
 
 def JeansMassMsun(field,data):
     return (MJ_constant * 



More information about the yt-svn mailing list