[Yt-svn] yt-commit r696 - trunk/yt/lagos
britton at wrangler.dreamhost.com
britton at wrangler.dreamhost.com
Fri Jul 25 16:19:57 PDT 2008
Author: britton
Date: Fri Jul 25 16:19:56 2008
New Revision: 696
URL: http://yt.spacepope.org/changeset/696
Log:
Added JeansMassMsun field.
Finished the last of the root beer.
25 July, 2008
Britton Smith
Modified:
trunk/yt/lagos/DerivedFields.py
Modified: trunk/yt/lagos/DerivedFields.py
==============================================================================
--- trunk/yt/lagos/DerivedFields.py (original)
+++ trunk/yt/lagos/DerivedFields.py Fri Jul 25 16:19:56 2008
@@ -44,6 +44,9 @@
clight = 3.0e10 # cm/s
kboltz = 1.38e-16 # erg K^-1
G = 6.67e-8 # cm^3 g^-1 s^-2
+Msun2g = 1.989e33
+MJ_constant = (((5*kboltz)/(G*mh))**(1.5)) * (3/(4*pi))**(0.5) / Msun2g
+
class FieldInfoContainer: # We are all Borg.
_shared_state = {}
@@ -575,6 +578,13 @@
add_field("NumberDensity", units=r"\rm{cm}^{-3}",
convert_function=_ConvertNumberDensity)
+def JeansMassMsun(field,data):
+ return (MJ_constant *
+ ((data["Temperature"]/data["MeanMolecularWeight"])**(1.5)) *
+ (data["Density"]**(-0.5)))
+add_field("JeansMassMsun",function=JeansMassMsun,units=r"\rm{Msun}")
+
+
def _CellMass(field, data):
return data["Density"] * data["CellVolume"]
def _convertCellMassMsun(data):
More information about the yt-svn
mailing list