[Yt-svn] commit/yt: Michael Kuhlen: Fixed comments h^3 multiplication in analytical mass function. Changed to outside-in determination of the virial radius for the halo profiler.

[Bitbucket]

1 new changeset in yt:

http://bitbucket.org/yt_analysis/yt/changeset/173a08b2e892/
changeset:   r4280:173a08b2e892
branch:      yt
user:        Michael Kuhlen
date:        2011-05-24 23:13:25
summary:     Fixed comments h^3 multiplication in analytical mass function. Changed to outside-in determination of the virial radius for the halo profiler.
affected #:  2 files (47 bytes)

--- a/yt/analysis_modules/halo_mass_function/halo_mass_function.py	Tue May 24 15:04:50 2011 -0600
+++ b/yt/analysis_modules/halo_mass_function/halo_mass_function.py	Tue May 24 14:13:25 2011 -0700
@@ -323,7 +323,7 @@
             dn_M_z = -1.0 / thissigma * dsigmadm * rho0 / self.massarray[i] * \
             self.multiplicityfunction(thissigma)*(self.massarray[i+1] - self.massarray[i]);
 
-            # scale by h^4 to get rid of all factors of h
+            # scale by h^3 to get rid of all factors of h
             dn_M_z *= math.pow(self.hubble0, 3.0);
             
             # keep track of cumulative number density


--- a/yt/analysis_modules/halo_profiler/halo_filters.py	Tue May 24 15:04:50 2011 -0600
+++ b/yt/analysis_modules/halo_profiler/halo_filters.py	Tue May 24 14:13:25 2011 -0700
@@ -102,8 +102,8 @@
     elif (overDensity[-1] >= virial_overdensity):
         index = -2
     else:
-        for q in (na.arange(len(overDensity)-2))+2:
-            if (overDensity[q] < virial_overdensity):
+        for q in (na.arange(len(overDensity),0,-1)-1):
+            if (overDensity[q] < virial_overdensity) and (overDensity[q-1] > virial_overdensity):
                 index = q - 1
                 break

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.