[yt-users] derived field and TotalQuantity from two sources possible?

[Britton Smith]

Hi G.S.,

This is pretty far outside of how fields are designed to work in yt.  The
access of data in yt isn't currently designed to work with more than one
pf.  I don't actually see a way to do this without combining the data so
that it can be seen as a single pf.  If parallel wasn't an issue and this
were only a unigrid run, you might be able to read in all the data at once
from both pfs and do some of this by hand.  If this is AMR, then I think
this may not be possible at all, considering the hierarchy is unlikely to be
exactly the same between two different datasets.

Britton

On Mon, Jul 18, 2011 at 5:27 PM, <gso at physics.ucsd.edu> wrote:

> Hi,
>
> I was wondering if it is possible to specify two different sources of
> input to create a derived field.  Currently with a single source, I have a
> derived field defined as:
>
> def _HII_HFraction(field, data):
>   return data["HII_Density"]/(data["HI_Density"]+data["HII_Density"])
>
> add_field("HII_HFraction", function=_HII_HFraction,
> units=r"\frac{\rho_{HII}}{\rho_H}")
>
> And I can access the total quantity in parallel
>
> pf=load(file)
>
> dd = pf.h.all_data()
>
> dd.quantities["TotalQuantity"]("HII_HFraction")[0]
>
> But I'm dealing with
>
> pf1=load(file1)
> pf2=load(file2)
>
> pf1 has the HI_Density field data and pf2 has the HII_Density field data.
> Is there a way to create a derived field and use TotalQuantity to operate
> on the data in parallel?
>
> From
> G.S.
>
> PS. I guess my alternative is to glue the two HDF5 files into one, but
> want to avoid it if possible.
>
>
> _______________________________________________
> 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/20110718/7ee7e1ac/attachment.html>