<div dir="ltr">Hi Kacper,<div><br></div><div>Can I add the gradient fields globally instead of to a specific dataset? </div><div><br></div><div>Yuan</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jul 20, 2015 at 6:08 PM, Kacper Kowalik <span dir="ltr"><<a href="mailto:xarthisius.kk@gmail.com" target="_blank">xarthisius.kk@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">On 07/20/2015 05:03 PM, Suoqing Ji wrote:<br>
> Hi Yuan,<br>
><br>
> Unfortunately I can not reproduce this error with either yt 2.6 or yt<br>
> 3… If you did not mind, could you provide your script and dataset (or<br>
> send me offline) so I can take a quick look at it?<br>
><br>
> Best wishes, -- Suoqing JI Ph.D Candidate Department of Physics<br>
> University of California, Santa Barbara<br>
> <a href="http://web.physics.ucsb.edu/~suoqing" rel="noreferrer" target="_blank">http://web.physics.ucsb.edu/~suoqing</a><br>
<br>
</span>Hi,<br>
could you show us the method you're using for your fields?<br>
<br>
Maybe, it's also a good opportunity to advertise yt-3.x a bit ;)<br>
In latest version you can do:<br>
<br>
ds.add_gradient_fields(("gas", "temperature"))<br>
<br>
where ds is your loaded dataset. It should create following fields:<br>
<br>
[('gas', 'temperature_gradient_x'),<br>
 ('gas', 'temperature_gradient_y'),<br>
 ('gas', 'temperature_gradient_z'),<br>
 ('gas', 'temperature_gradient_magnitude')]<br>
<br>
Cheers,<br>
Kacper<br>
<span class=""><br>
>> On Jul 20, 2015, at 2:29 PM, Yuan Li <<a href="mailto:bear0980@gmail.com">bear0980@gmail.com</a>> wrote:<br>
>><br>
>> Hi Suoqing,<br>
>><br>
>> I added validators=[ValidateSpatial(1, ["Pressure"])]<br>
>><br>
>> Now the error message is ValueError: could not broadcast input<br>
>> array from shape (64,64,66) into shape (64,64,64)<br>
>><br>
>> Yuan<br>
>><br>
>> On Sun, Jul 19, 2015 at 6:25 PM, Suoqing Ji<br>
</span>>> <<a href="mailto:suoqing@physics.ucsb.edu">suoqing@physics.ucsb.edu</a> <mailto:<a href="mailto:suoqing@physics.ucsb.edu">suoqing@physics.ucsb.edu</a>>><br>
<span class="">>> wrote: Hi Yuan,<br>
>><br>
>> Just to check, have you included the “ValidateSpatial" object when<br>
>> adding the field? See the definition of DivV field:<br>
>> <a href="http://yt-project.org/docs/2.6/analyzing/creating_derived_fields.html#some-more-complicated-examples" rel="noreferrer" target="_blank">http://yt-project.org/docs/2.6/analyzing/creating_derived_fields.html#some-more-complicated-examples</a><br>
>> <<a href="http://yt-project.org/docs/2.6/analyzing/creating_derived_fields.html#some-more-complicated-examples" rel="noreferrer" target="_blank">http://yt-project.org/docs/2.6/analyzing/creating_derived_fields.html#some-more-complicated-examples</a>><br>
>><br>
>><br>
>><br>
Best wishes,<br>
>> -- Suoqing JI Ph.D Candidate Department of Physics University of<br>
>> California, Santa Barbara <a href="http://web.physics.ucsb.edu/~suoqing" rel="noreferrer" target="_blank">http://web.physics.ucsb.edu/~suoqing</a><br>
>> <<a href="http://web.physics.ucsb.edu/~suoqing" rel="noreferrer" target="_blank">http://web.physics.ucsb.edu/~suoqing</a>><br>
</span>>>> On Jul 19, 2015, at 2:52 PM, Yuan Li <<a href="mailto:yuan@astro.columbia.edu">yuan@astro.columbia.edu</a><br>
<span class="">>>> <mailto:<a href="mailto:yuan@astro.columbia.edu">yuan@astro.columbia.edu</a>>> wrote:<br>
>>><br>
>>> Dear all,<br>
>>><br>
>>> I am using yt2.6. I was trying to define a field of the<br>
>>> temperature gradient like gradPressureX, but I got an error<br>
>>> message: ... new_field[1:-1,1:-1,1:-1]  =<br>
>>> data["Temperature"][sl_right,1:-1,1:-1]/ds ValueError: could not<br>
>>> broadcast input array from shape (62,62,64) into shape<br>
>>> (62,62,62)<br>
>>><br>
>>> I then tried to re-define gradPressureX by copying the original<br>
>>> definition<br>
>>> (<a href="http://yt-project.org/docs/2.6/reference/field_list.html#gradpressurex" rel="noreferrer" target="_blank">http://yt-project.org/docs/2.6/reference/field_list.html#gradpressurex</a><br>
</span>>>> <<a href="http://yt-project.org/docs/2.6/reference/field_list.html#gradpressurex" rel="noreferrer" target="_blank">http://yt-project.org/docs/2.6/reference/field_list.html#gradpressurex</a>>)<br>
<span class="">>>> and giving it a different name, and I got the same error. The<br>
>>> original gradPressureX works fine.<br>
>>><br>
>>> Is it because I did not load some specific module or something?<br>
>>><br>
>>><br>
>>> Thank you! Yuan _______________________________________________<br>
</span>>>> yt-users mailing list <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
>>> <mailto:<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a>><br>
<span class="">>>> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
>>> <<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a>><br>
>><br>
>><br>
>> _______________________________________________ yt-users mailing<br>
</span>>> list <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
>> <mailto:<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a>><br>
<span class="">>> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
>> <<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a>><br>
>><br>
>><br>
>> _______________________________________________ yt-users mailing<br>
</span>>> list <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<span class="">>> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
><br>
><br>
><br>
><br>
> _______________________________________________ yt-users mailing<br>
</span>> list <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/listinfo.cgi/yt-users-spacepope.org</a><br>
><br>
<br>
<br>
<br>_______________________________________________<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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>