Hi Carla,

Interesting.  It's possible, for the call to "temperature" (which
raises NotImplementedError) that the version of yt you're on doesn't
yet support mesh fields in arbitrary grids; that was added reasonably
recently.  For the others, I'm not totally sure why it's breaking --
which version of yt are you using, and what type of data?  It may just
be that they aren't deposit fields.  In that case, try explicitly
stating obj["gas","temperature"] or obj["gas","HI_Density"].  It may
also be that the aliasing is getting in the way, in which case instead
of "HI_Density" you might need to try "Hp0_density".

On Tue, Jul 11, 2017 at 5:42 AM, Carla Bernhardt
<carla.j.bernhardt@gmail.com> wrote:
> Dear yt-users,
>
> I am trying to extract a unigrid data cube within a specific region and
> specified dimensions, thus I am using arbitrary_grid(), however, beyond the
> example given in the docs, my field extractions all fail. I would like to
> extract temperature, HI_Density and HII_Density. Here is a snippet from my
> code:
>
> obj = ds.arbitrary_grid(left_edge, right_edge, dims=[256, 256, 256])
> print(obj["deposit", "all_density"])
> print(obj["temperature"])
> print(obj["HI_Density"])
>
> The line print(obj["deposit", "all_density"]) works, but the two lines after
> both fail. Here is the traceback from the temperature, and I even tried out
> deposited temperature and deposited HI_Density. Am I formatting this
> incorrectly, or is there a different way to extract the fields?
>
> Kind regards,
> Carla Bernhardt
> PhD Student
> ZAH Institut für Theoretische Astrophysik
> Universität Heidelberg
>

> _______________________________________________
> yt-users mailing list
> yt-users@lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org