[yt-users] Subtracting datasets

[Jonah Miller]

Hi Matt, Boyan,

Although I'm not currently struggling with this, I will eventually 
encounter a similar situation... What would one do if the grid was /not/ 
uniform? Would a derived field still work if the grids were identical?

Best,
Jonah

On 15-06-16 09:34 AM, Matthew Turk wrote:
> Hi Boyan,
>
> Since the grids are uniform, you can do one of a few different things.
> The first would be if you want to set up a derived field, you could do
> something like:
>
> ds1 = ...
> ds2 = ...
>
> @yt.derived_field(name = "my_field", validators = [yt.ValidateSpatial()])
> def my_field(field, data):
>      f1 = ds1.index.grids[data.id - data._id_offset][field.name]
>      f2 = data[field.name]
>      return f1 - f2
>
> If you want to simply inspect, you can do
> ds1.index.grids[whatever][field] - ds2.index.grids[whatever][field].
>
> -Matt
>
> On Mon, Jun 15, 2015 at 11:05 PM, bh11e <bh11e at my.fsu.edu> wrote:
>> Hi,
>>
>> I was wondering if someone could advise on the following problem. Consider two (Enzo) datasets with identical grids. What is the way to plot the difference between a field in one set and a field in the other. Alternatively is there a way to create a new dataset with the calculated difference, after which plotting will be as usual.
>>
>> Thanks,
>>
>> --Boyan Hristov
>> _______________________________________________
>> yt-users mailing list
>> yt-users at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> _______________________________________________
> 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/20150616/e707a86f/attachment.htm>