<div dir="ltr">Thanks Everybody!!!<div><br></div><div>There was no problem after all :P</div><div>I hope this to be useful for someone else.</div><div><br></div><div>Matías Suazo,</div><div>MSc student at Universidad de Chile</div></div><div class="gmail_extra"><br><div class="gmail_quote">2017-11-14 17:42 GMT-03:00 <span dir="ltr"><<a href="mailto:yt-users-request@lists.spacepope.org" target="_blank">yt-users-request@lists.spacepope.org</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send yt-users mailing list submissions to<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:yt-users-request@lists.spacepope.org">yt-users-request@lists.<wbr>spacepope.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:yt-users-owner@lists.spacepope.org">yt-users-owner@lists.<wbr>spacepope.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of yt-users digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Re: Projection (Jos? Mauricio Utreras)<br>
2. Re: Projection (Nathan Goldbaum)<br>
3. Re: Projection (Matthew Turk)<br>
<br>
<br>
------------------------------<wbr>------------------------------<wbr>----------<br>
<br>
Message: 1<br>
Date: Tue, 14 Nov 2017 15:45:23 -0300<br>
From: Jos? Mauricio Utreras <<a href="mailto:jutreras@ug.uchile.cl">jutreras@ug.uchile.cl</a>><br>
To: Discussion of the yt analysis package<br>
<<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a>><br>
Subject: Re: [yt-users] Projection<br>
Message-ID:<br>
<<wbr>CANhQbF5OzXEhNzguuXR8iNRrAryy7<wbr>4iQ=<a href="mailto:sVg8zZg5QCzyzqvJg@mail.gmail.com">sVg8zZg5QCzyzqvJg@mail.<wbr>gmail.com</a>><br>
Content-Type: text/plain; charset="UTF-8"<br>
<br>
Hi everyone,<br>
<br>
I got the same result using the IsolatedGalaxy dataset<br>
<br>
ds=yt.load('IsolatedGalaxy/<wbr>galaxy0030/galaxy0030')<br>
disk = ds.disk('c',[0,0,1],(1000,'pc'<wbr>),(100,'pc'))<br>
proj=ds.proj('density',axis='<wbr>z',weight_field='density',<wbr>data_source=disk,method='<wbr>integrate')<br>
proj['density']<br>
>>> YTArray([ nan, nan, nan, ..., nan, nan, nan]) g/cm**3<br>
<br>
2017-11-14 15:38 GMT-03:00, Nathan Goldbaum <<a href="mailto:nathan12343@gmail.com">nathan12343@gmail.com</a>>:<br>
> Can I get one of you to open an issue?<br>
><br>
> <a href="https://github.com/yt-project/yt/issues/new" rel="noreferrer" target="_blank">https://github.com/yt-project/<wbr>yt/issues/new</a><br>
><br>
> If you can, please include a script to reproduce the problem, making use of<br>
> one of the public datasets on <a href="http://yt-project.org/data" rel="noreferrer" target="_blank">yt-project.org/data</a>.<br>
><br>
> -Nathan<br>
><br>
> On Tue, Nov 14, 2017 at 12:34 PM, Andrew James Emerick<br>
> <<a href="mailto:aje2123@columbia.edu">aje2123@columbia.edu</a><br>
>> wrote:<br>
><br>
>> Hello Matias and Nathan,<br>
>><br>
>> I was able to reproduce this error on my own Enzo data set(s) using<br>
>> similar code. I did a little more playing around and it looks like it<br>
>> only<br>
>> breaks when having a weight_field. weight_field = None seems to be O.K. .<br>
>> It also occurs for the "region" objects, in addition to the disk.<br>
>> Interestingly, as Matias pointed out, doing something like this works:<br>
>><br>
>> region = ds.region([0.5,0.5,0.5], [0.0, 0.0, 0.0], [1.0, 1.0, 1.0])<br>
>> proj=ds.proj('density',axis='<wbr>z',weight_field='density',<wbr>data_<br>
>> source=region,method='<wbr>integrate')<br>
>><br>
>> but this does not:<br>
>><br>
>> region = ds.region([0.5,0.5,0.5], [0.1,0.1,0.1], [0.9,0.9,0.9] )<br>
>> proj =ds.proj('density',axis='z',<wbr>weight_field='density',data_<br>
>> source=region,method='<wbr>integrate')<br>
>><br>
>> In short, this indeed appears to be a bug, and possibly not specific to a<br>
>> single dataset.<br>
>><br>
>><br>
>><br>
>><br>
>> Andrew<br>
>><br>
>><br>
>> On Tue, Nov 14, 2017 at 12:40 PM, Nathan Goldbaum <<a href="mailto:nathan12343@gmail.com">nathan12343@gmail.com</a>><br>
>> wrote:<br>
>><br>
>>> This sounds like a bug. I don't think you should ever get back NaNs from<br>
>>> a projection operation.<br>
>>><br>
>>> Is there any chance you can share the dataset that's triggering this<br>
>>> behavior?<br>
>>><br>
>>> On Tue, Nov 14, 2017 at 11:34 AM, Matias Enrique Suazo Campos <<br>
>>> <a href="mailto:matias.suazo@ug.uchile.cl">matias.suazo@ug.uchile.cl</a>> wrote:<br>
>>><br>
>>>> Dear YT users,<br>
>>>><br>
>>>> I'm trying to make a density projection, in the following lines I show<br>
>>>> you the lines that I have.<br>
>>>><br>
>>>> #Ramses output<br>
>>>> ds = yt.load(path+"output_%05d"%<wbr>snap+"/info_%05d.txt" %snap)<br>
>>>> disk = ds.disk('max',[0,0,1],(20,'pc'<wbr>),(5,'pc'))<br>
>>>> proj=ds.proj('density',axis='<wbr>z',weight_field='density',<wbr>data_<br>
>>>> source=disk,method='integrate'<wbr>)<br>
>>>><br>
>>>> However when I print proj["density"] I get a list of nan, however when<br>
>>>> I<br>
>>>> use the whole domain as my data_source I get reasonable numbers.<br>
>>>><br>
>>>> Thanks,<br>
>>>> Mat?as Suazo<br>
>>>> MSc student at Universidad de Chile<br>
>>>><br>
>>>> ______________________________<wbr>_________________<br>
>>>> yt-users mailing list<br>
>>>> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
>>>> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
>>>><br>
>>>><br>
>>><br>
>>> ______________________________<wbr>_________________<br>
>>> yt-users mailing list<br>
>>> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
>>> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
>>><br>
>>><br>
>><br>
>><br>
>> --<br>
>> NSF Graduate Fellow<br>
>> Columbia University<br>
>> Department of Astronomy<br>
>><br>
>> ______________________________<wbr>_________________<br>
>> yt-users mailing list<br>
>> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
>> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
>><br>
>><br>
><br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Tue, 14 Nov 2017 13:04:04 -0600<br>
From: Nathan Goldbaum <<a href="mailto:nathan12343@gmail.com">nathan12343@gmail.com</a>><br>
To: Discussion of the yt analysis package<br>
<<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a>><br>
Subject: Re: [yt-users] Projection<br>
Message-ID:<br>
<CAJXewOm3f91Z66L8WPYZ3ztnVx=<a href="mailto:BPaKAO9FuZ1Lq28qJZHNVZg@mail.gmail.com">B<wbr>PaKAO9FuZ1Lq28qJZHNVZg@mail.<wbr>gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Ah, so this isn't actually a bug. There's even a comment in the code<br>
indicating that this is the intended behavior:<br>
<br>
<a href="https://github.com/yt-project/yt/blob/master/yt/data_objects/construction_data_containers.py#L373" rel="noreferrer" target="_blank">https://github.com/yt-project/<wbr>yt/blob/master/yt/data_<wbr>objects/construction_data_<wbr>containers.py#L373</a><br>
<br>
So the choices about what to do here are, in my opinion:<br>
<br>
* Change this behavior to fill the array with zeros instead of NaNs<br>
* Improve the documentation to make this choice clearer<br>
<br>
I'm not sure offhand what is the best approach, especially given that<br>
changing this behavior would be a backward compatibility break.<br>
<br>
For the original question: your data aren't all NaNs I suspect. E.g. take a<br>
look at this example script:<br>
<br>
import yt<br>
import numpy as np<br>
ds = yt.load('IsolatedGalaxy/<wbr>galaxy0030/galaxy0030')<br>
disk = ds.disk('c', [0, 0, 1], (1000, 'pc'), (100, 'pc'))<br>
proj = ds.proj('density', axis='z', weight_field='density',<br>
data_source=disk,<br>
method='integrate')<br>
dens_proj = proj['density']<br>
print(dens_proj[np.isfinite(<wbr>dens_proj)])<br>
print(np.nanmax(dens_proj))<br>
print(np.nanmin(dens_proj))<br>
<br>
Which prints out values in regions where the projection isn't filled with<br>
NaNs.<br>
<br>
<br>
On Tue, Nov 14, 2017 at 12:45 PM, Jos? Mauricio Utreras <<br>
<a href="mailto:jutreras@ug.uchile.cl">jutreras@ug.uchile.cl</a>> wrote:<br>
<br>
> Hi everyone,<br>
><br>
> I got the same result using the IsolatedGalaxy dataset<br>
><br>
> ds=yt.load('IsolatedGalaxy/<wbr>galaxy0030/galaxy0030')<br>
> disk = ds.disk('c',[0,0,1],(1000,'pc'<wbr>),(100,'pc'))<br>
> proj=ds.proj('density',axis='<wbr>z',weight_field='density',<br>
> data_source=disk,method='<wbr>integrate')<br>
> proj['density']<br>
> >>> YTArray([ nan, nan, nan, ..., nan, nan, nan]) g/cm**3<br>
><br>
> 2017-11-14 15:38 GMT-03:00, Nathan Goldbaum <<a href="mailto:nathan12343@gmail.com">nathan12343@gmail.com</a>>:<br>
> > Can I get one of you to open an issue?<br>
> ><br>
> > <a href="https://github.com/yt-project/yt/issues/new" rel="noreferrer" target="_blank">https://github.com/yt-project/<wbr>yt/issues/new</a><br>
> ><br>
> > If you can, please include a script to reproduce the problem, making use<br>
> of<br>
> > one of the public datasets on <a href="http://yt-project.org/data" rel="noreferrer" target="_blank">yt-project.org/data</a>.<br>
> ><br>
> > -Nathan<br>
> ><br>
> > On Tue, Nov 14, 2017 at 12:34 PM, Andrew James Emerick<br>
> > <<a href="mailto:aje2123@columbia.edu">aje2123@columbia.edu</a><br>
> >> wrote:<br>
> ><br>
> >> Hello Matias and Nathan,<br>
> >><br>
> >> I was able to reproduce this error on my own Enzo data set(s) using<br>
> >> similar code. I did a little more playing around and it looks like it<br>
> >> only<br>
> >> breaks when having a weight_field. weight_field = None seems to be O.K.<br>
> .<br>
> >> It also occurs for the "region" objects, in addition to the disk.<br>
> >> Interestingly, as Matias pointed out, doing something like this works:<br>
> >><br>
> >> region = ds.region([0.5,0.5,0.5], [0.0, 0.0, 0.0], [1.0, 1.0, 1.0])<br>
> >> proj=ds.proj('density',axis='<wbr>z',weight_field='density',<wbr>data_<br>
> >> source=region,method='<wbr>integrate')<br>
> >><br>
> >> but this does not:<br>
> >><br>
> >> region = ds.region([0.5,0.5,0.5], [0.1,0.1,0.1], [0.9,0.9,0.9] )<br>
> >> proj =ds.proj('density',axis='z',<wbr>weight_field='density',data_<br>
> >> source=region,method='<wbr>integrate')<br>
> >><br>
> >> In short, this indeed appears to be a bug, and possibly not specific to<br>
> a<br>
> >> single dataset.<br>
> >><br>
> >><br>
> >><br>
> >><br>
> >> Andrew<br>
> >><br>
> >><br>
> >> On Tue, Nov 14, 2017 at 12:40 PM, Nathan Goldbaum <<br>
> <a href="mailto:nathan12343@gmail.com">nathan12343@gmail.com</a>><br>
> >> wrote:<br>
> >><br>
> >>> This sounds like a bug. I don't think you should ever get back NaNs<br>
> from<br>
> >>> a projection operation.<br>
> >>><br>
> >>> Is there any chance you can share the dataset that's triggering this<br>
> >>> behavior?<br>
> >>><br>
> >>> On Tue, Nov 14, 2017 at 11:34 AM, Matias Enrique Suazo Campos <<br>
> >>> <a href="mailto:matias.suazo@ug.uchile.cl">matias.suazo@ug.uchile.cl</a>> wrote:<br>
> >>><br>
> >>>> Dear YT users,<br>
> >>>><br>
> >>>> I'm trying to make a density projection, in the following lines I show<br>
> >>>> you the lines that I have.<br>
> >>>><br>
> >>>> #Ramses output<br>
> >>>> ds = yt.load(path+"output_%05d"%<wbr>snap+"/info_%05d.txt" %snap)<br>
> >>>> disk = ds.disk('max',[0,0,1],(20,'pc'<wbr>),(5,'pc'))<br>
> >>>> proj=ds.proj('density',axis='<wbr>z',weight_field='density',<wbr>data_<br>
> >>>> source=disk,method='integrate'<wbr>)<br>
> >>>><br>
> >>>> However when I print proj["density"] I get a list of nan, however when<br>
> >>>> I<br>
> >>>> use the whole domain as my data_source I get reasonable numbers.<br>
> >>>><br>
> >>>> Thanks,<br>
> >>>> Mat?as Suazo<br>
> >>>> MSc student at Universidad de Chile<br>
> >>>><br>
> >>>> ______________________________<wbr>_________________<br>
> >>>> yt-users mailing list<br>
> >>>> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
> >>>> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
> >>>><br>
> >>>><br>
> >>><br>
> >>> ______________________________<wbr>_________________<br>
> >>> yt-users mailing list<br>
> >>> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
> >>> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
> >>><br>
> >>><br>
> >><br>
> >><br>
> >> --<br>
> >> NSF Graduate Fellow<br>
> >> Columbia University<br>
> >> Department of Astronomy<br>
> >><br>
> >> ______________________________<wbr>_________________<br>
> >> yt-users mailing list<br>
> >> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
> >> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
> >><br>
> >><br>
> ><br>
> ______________________________<wbr>_________________<br>
> yt-users mailing list<br>
> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20171114/0e8395d8/attachment-0001.html" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>pipermail/yt-users-spacepope.<wbr>org/attachments/20171114/<wbr>0e8395d8/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Tue, 14 Nov 2017 14:42:01 -0600<br>
From: Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>><br>
To: Discussion of the yt analysis package<br>
<<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a>><br>
Subject: Re: [yt-users] Projection<br>
Message-ID:<br>
<CALO3=<a href="mailto:5GimZpK9u2i40ZWfECXx3btj1GFT1jLO%2B6OhVn6_DqQGg@mail.gmail.com">5GimZpK9u2i40ZWfECXx3bt<wbr>j1GFT1jLO+6OhVn6_DqQGg@mail.<wbr>gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
On Tue, Nov 14, 2017 at 1:04 PM, Nathan Goldbaum <<a href="mailto:nathan12343@gmail.com">nathan12343@gmail.com</a>><br>
wrote:<br>
<br>
> Ah, so this isn't actually a bug. There's even a comment in the code<br>
> indicating that this is the intended behavior:<br>
><br>
> <a href="https://github.com/yt-project/yt/blob/master/yt/data_" rel="noreferrer" target="_blank">https://github.com/yt-project/<wbr>yt/blob/master/yt/data_</a><br>
> objects/construction_data_<wbr>containers.py#L373<br>
><br>
> So the choices about what to do here are, in my opinion:<br>
><br>
> * Change this behavior to fill the array with zeros instead of NaNs<br>
><br>
<br>
I think this would give the wrong impression, especially for things like<br>
weighted average velocity fields.<br>
<br>
<br>
> * Improve the documentation to make this choice clearer<br>
><br>
<br>
Probably a good thing.<br>
<br>
I think I would suggest another option, which would be to filter it based<br>
on which values are non-NaN. In the cases identified in the emails, I<br>
would *suspect* that what's happening is that there are elements in the<br>
quad tree that are initialized from the tree initialization step (which<br>
just looks at the mesh from the container, not the selected elements in the<br>
mesh) and left at zero following initialization. We could either have a<br>
"touched" flag on all of the mesh elements (i.e., "has this been touched by<br>
data, or is it left at initialized value) and filter in the conversion of<br>
the quad tree to the flat arrays, or we could filter in the python code<br>
based on NaNs. My preference would be to filter it in the creation of the<br>
flat arrays, as then it will be uniformly applied even without weight<br>
fields.<br>
<br>
<br>
><br>
> I'm not sure offhand what is the best approach, especially given that<br>
> changing this behavior would be a backward compatibility break.<br>
><br>
<br>
I think filtering NaNs is not a big change in compatibility. I think we'd<br>
just need to make sure that the cmap normalization has the same behavior<br>
for missing values as for NaN values, which I think it does.<br>
<br>
<br>
><br>
> For the original question: your data aren't all NaNs I suspect. E.g. take<br>
> a look at this example script:<br>
><br>
> import yt<br>
> import numpy as np<br>
> ds = yt.load('IsolatedGalaxy/<wbr>galaxy0030/galaxy0030')<br>
> disk = ds.disk('c', [0, 0, 1], (1000, 'pc'), (100, 'pc'))<br>
> proj = ds.proj('density', axis='z', weight_field='density',<br>
> data_source=disk,<br>
> method='integrate')<br>
> dens_proj = proj['density']<br>
> print(dens_proj[np.isfinite(<wbr>dens_proj)])<br>
> print(np.nanmax(dens_proj))<br>
> print(np.nanmin(dens_proj))<br>
><br>
> Which prints out values in regions where the projection isn't filled with<br>
> NaNs.<br>
><br>
><br>
> On Tue, Nov 14, 2017 at 12:45 PM, Jos? Mauricio Utreras <<br>
> <a href="mailto:jutreras@ug.uchile.cl">jutreras@ug.uchile.cl</a>> wrote:<br>
><br>
>> Hi everyone,<br>
>><br>
>> I got the same result using the IsolatedGalaxy dataset<br>
>><br>
>> ds=yt.load('IsolatedGalaxy/<wbr>galaxy0030/galaxy0030')<br>
>> disk = ds.disk('c',[0,0,1],(1000,'pc'<wbr>),(100,'pc'))<br>
>> proj=ds.proj('density',axis='<wbr>z',weight_field='density',<wbr>data_<br>
>> source=disk,method='integrate'<wbr>)<br>
>> proj['density']<br>
>> >>> YTArray([ nan, nan, nan, ..., nan, nan, nan]) g/cm**3<br>
>><br>
>> 2017-11-14 15:38 GMT-03:00, Nathan Goldbaum <<a href="mailto:nathan12343@gmail.com">nathan12343@gmail.com</a>>:<br>
>> > Can I get one of you to open an issue?<br>
>> ><br>
>> > <a href="https://github.com/yt-project/yt/issues/new" rel="noreferrer" target="_blank">https://github.com/yt-project/<wbr>yt/issues/new</a><br>
>> ><br>
>> > If you can, please include a script to reproduce the problem, making<br>
>> use of<br>
>> > one of the public datasets on <a href="http://yt-project.org/data" rel="noreferrer" target="_blank">yt-project.org/data</a>.<br>
>> ><br>
>> > -Nathan<br>
>> ><br>
>> > On Tue, Nov 14, 2017 at 12:34 PM, Andrew James Emerick<br>
>> > <<a href="mailto:aje2123@columbia.edu">aje2123@columbia.edu</a><br>
>> >> wrote:<br>
>> ><br>
>> >> Hello Matias and Nathan,<br>
>> >><br>
>> >> I was able to reproduce this error on my own Enzo data set(s) using<br>
>> >> similar code. I did a little more playing around and it looks like it<br>
>> >> only<br>
>> >> breaks when having a weight_field. weight_field = None seems to be<br>
>> O.K. .<br>
>> >> It also occurs for the "region" objects, in addition to the disk.<br>
>> >> Interestingly, as Matias pointed out, doing something like this works:<br>
>> >><br>
>> >> region = ds.region([0.5,0.5,0.5], [0.0, 0.0, 0.0], [1.0, 1.0, 1.0])<br>
>> >> proj=ds.proj('density',axis='<wbr>z',weight_field='density',<wbr>data_<br>
>> >> source=region,method='<wbr>integrate')<br>
>> >><br>
>> >> but this does not:<br>
>> >><br>
>> >> region = ds.region([0.5,0.5,0.5], [0.1,0.1,0.1], [0.9,0.9,0.9] )<br>
>> >> proj =ds.proj('density',axis='z',<wbr>weight_field='density',data_<br>
>> >> source=region,method='<wbr>integrate')<br>
>> >><br>
>> >> In short, this indeed appears to be a bug, and possibly not specific<br>
>> to a<br>
>> >> single dataset.<br>
>> >><br>
>> >><br>
>> >><br>
>> >><br>
>> >> Andrew<br>
>> >><br>
>> >><br>
>> >> On Tue, Nov 14, 2017 at 12:40 PM, Nathan Goldbaum <<br>
>> <a href="mailto:nathan12343@gmail.com">nathan12343@gmail.com</a>><br>
>> >> wrote:<br>
>> >><br>
>> >>> This sounds like a bug. I don't think you should ever get back NaNs<br>
>> from<br>
>> >>> a projection operation.<br>
>> >>><br>
>> >>> Is there any chance you can share the dataset that's triggering this<br>
>> >>> behavior?<br>
>> >>><br>
>> >>> On Tue, Nov 14, 2017 at 11:34 AM, Matias Enrique Suazo Campos <<br>
>> >>> <a href="mailto:matias.suazo@ug.uchile.cl">matias.suazo@ug.uchile.cl</a>> wrote:<br>
>> >>><br>
>> >>>> Dear YT users,<br>
>> >>>><br>
>> >>>> I'm trying to make a density projection, in the following lines I<br>
>> show<br>
>> >>>> you the lines that I have.<br>
>> >>>><br>
>> >>>> #Ramses output<br>
>> >>>> ds = yt.load(path+"output_%05d"%<wbr>snap+"/info_%05d.txt" %snap)<br>
>> >>>> disk = ds.disk('max',[0,0,1],(20,'pc'<wbr>),(5,'pc'))<br>
>> >>>> proj=ds.proj('density',axis='<wbr>z',weight_field='density',<wbr>data_<br>
>> >>>> source=disk,method='integrate'<wbr>)<br>
>> >>>><br>
>> >>>> However when I print proj["density"] I get a list of nan, however<br>
>> when<br>
>> >>>> I<br>
>> >>>> use the whole domain as my data_source I get reasonable numbers.<br>
>> >>>><br>
>> >>>> Thanks,<br>
>> >>>> Mat?as Suazo<br>
>> >>>> MSc student at Universidad de Chile<br>
>> >>>><br>
>> >>>> ______________________________<wbr>_________________<br>
>> >>>> yt-users mailing list<br>
>> >>>> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
>> >>>> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
>> >>>><br>
>> >>>><br>
>> >>><br>
>> >>> ______________________________<wbr>_________________<br>
>> >>> yt-users mailing list<br>
>> >>> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
>> >>> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
>> >>><br>
>> >>><br>
>> >><br>
>> >><br>
>> >> --<br>
>> >> NSF Graduate Fellow<br>
>> >> Columbia University<br>
>> >> Department of Astronomy<br>
>> >><br>
>> >> ______________________________<wbr>_________________<br>
>> >> yt-users mailing list<br>
>> >> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
>> >> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
>> >><br>
>> >><br>
>> ><br>
>> ______________________________<wbr>_________________<br>
>> yt-users mailing list<br>
>> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
>> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
>><br>
><br>
><br>
> ______________________________<wbr>_________________<br>
> yt-users mailing list<br>
> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
><br>
><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20171114/49378c57/attachment.html" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>pipermail/yt-users-spacepope.<wbr>org/attachments/20171114/<wbr>49378c57/attachment.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
______________________________<wbr>_________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
<br>
<br>
------------------------------<br>
<br>
End of yt-users Digest, Vol 117, Issue 22<br>
******************************<wbr>***********<br>
</blockquote></div><br></div>