[Yt-svn] commit/yt: MatthewTurk: Moving toward a more full field list

[Bitbucket]

1 new changeset in yt:

http://bitbucket.org/yt_analysis/yt/changeset/f9d3ff8fd862/
changeset:   r4255:f9d3ff8fd862
branch:      yt
user:        MatthewTurk
date:        2011-05-13 17:18:39
summary:     Moving toward a more full field list
affected #:  1 file (906 bytes)

--- a/yt/analysis_modules/halo_profiler/standard_analysis.py	Fri May 13 10:25:20 2011 -0400
+++ b/yt/analysis_modules/halo_profiler/standard_analysis.py	Fri May 13 11:18:39 2011 -0400
@@ -30,9 +30,35 @@
 
 analysis_field_list = [
     "Density",
-    "Temperature"
+    "Temperature",
+    "VelocityMagnitude",
+    ("Ones", None),
+    "Entropy",
+    "RadialVelocity",
+    "SpecificAngularMomnetumX",
+    "SpecificAngularMomnetumY",
+    "SpecificAngularMomnetumZ",
+    "CoolingTime",
+    "DynamicalTime",
+    ("CellMassMsun", None),
+    "Dark_Matter_Density",
+    #("ParticleSpecificAngularMomentumX", "ParticleMassMsun"),
+    #("ParticleSpecificAngularMomentumY", "ParticleMassMsun"),
+    #("ParticleSpecificAngularMomentumZ", "ParticleMassMsun"),
+    ("TotalMass", None),
+    "OverDensity",
+    #("ParticleMassMsun", None),
+    ("StarParticleDensity", "StarParticleMassMsun"), # How do we weight this?
+    #("StarParticleMassMsun", None), 
+    ("StarParticleDensity", "StarParticleMassMsun"), # How do we weight this?
 ]
 
+analysis_field_list += ["%s_Fraction" % (s) for s in
+    ["HI","HII","HeI","HeII","HeIII","H2I","H2II","HM","Electron",
+    "DI","DII","HDI","Metal"]
+]
+    
+
 class StandardRadialAnalysis(object):
     def __init__(self, pf, center, radius, n_bins = 128, inner_radius = None):
         self.pf = pf
@@ -78,7 +104,7 @@
         import matplotlib; matplotlib.use("Agg")
         import pylab
         for field in self.prof.keys():
-            if field == "Radius": continue
+            if field in ("UsedBins", "Radius"): continue
             pylab.clf()
             pylab.loglog(self.prof["Radius"], self.prof[field], '-x')
             pylab.xlabel("Radius [cm]")

Repository URL: https://bitbucket.org/yt_analysis/yt/

--

This is a commit notification from bitbucket.org. You are receiving
this because you have the service enabled, addressing the recipient of
this email.