[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