[yt-users] Cold Gas

[Chen Jingjing]

Dear yt-users,
I was trying to calculate only the total mass of cold gas(<1e4.5) within a
sphere and I tried the following script:

ds=yt.load('RedshiftOutputxxxx')
sp = ds.sphere(center,radius)
*spcut = sp.cut_region(["obj['temperature']<1e4.5"])*
*spcut.quantities.total_mass()*

But  I got an error that says:

yt.utilities.exceptions.YTMixedCutRegion: Can't mix particle/discrete and
fluid/mesh conditions or

quantities.  Field: ('all', 'particle_mass') and Conditions specified:

            obj['temperature']<1e4.5


Does anyone know how to fix it?

Thanks,

Jingjing
-- 

Jingjing Chen [ chenjj235 at gmail.com ]
Department of Astronomy
Columbia University
550 West 120th Street
New York, New York 10027
Cell Phone: +1(347)5746709
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20150407/11e0825b/attachment.htm>
-------------- next part --------------
_______________________________________________
yt-users mailing list
yt-users at lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org