sorry, that should be "np.sum" instead of "sum"

Andrew E.

On Fri, Apr 10, 2015 at 4:32 PM, Andrew James Emerick <aje2123@columbia.edu> wrote:

Hi Jingjing,There may be a prettier way to do this, but I would go about this like so:1) calculate the HI mass in each cell as HI_density * cell_volume:HI_mass = sp['enzo','HI_Density'] * sp['dx']**32) then summing:total_mass = sum(HI_mass)Andrew E.On Fri, Apr 10, 2015 at 4:20 PM, Chen Jingjing <chenjj235@gmail.com> wrote:_______________________________________________Dear yt-users,I am trying to calculate the total HI mass within a sphere. I have a field called 'HI_Density'. And I tried to use the scriptHI_Mass = sphere['enzo','HI_Density'].in_units('Msun').sum()But I got an error that says the unit code_mass/code_length**3 cannot be converted to code_mass.So I thought I might need to create a derived_field called 'HI_Mass' which mutiplies the HI_Density with code_length**3. But I don't see that there is any existing field of cell length/volume. So how should I do this? Or is there an easier way to compute the total HI mass without deriving a field?Thanks,Jingjing--Jingjing Chen [ chenjj235@gmail.com ] Department of Astronomy Columbia University 550 West 120th Street New York, New York 10027

Cell Phone: +1(347)5746709

yt-users mailing list

yt-users@lists.spacepope.org

http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

--Graduate StudentColumbia UniversityDepartment of Astronomy

Graduate Student

Columbia University

Department of Astronomy