[yt-users] Cold Gas

[Nathan Goldbaum]

On Tue, Apr 7, 2015 at 2:03 PM, Chen Jingjing <chenjj235 at gmail.com> wrote:

> 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
>
>
>
The total_mass() quantity includes both gas and particle mass. As you've
discovered, cut_regions don't support selecting particle fields.

Instead you should use the total_quantity() derived quantity to find the
total of the 'cell_mass' field, which will only include the gas mass.


> 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
>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20150407/da354cba/attachment.htm>
-------------- next part --------------
_______________________________________________
yt-users mailing list
yt-users at lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org