<div dir="ltr">Hi Suoqing,<div><br></div><div>I basically copy-pasted from universal_fields.py:</div><div><br></div><div><div>def _gradPX(field, data):</div><div> # We need to set up stencils</div><div> if <a href="http://data.pf">data.pf</a>["HydroMethod"] == 2:</div><div> sl_left = slice(None,-2,None)</div><div> sl_right = slice(1,-1,None)</div><div> div_fac = 1.0</div><div> else:</div><div> sl_left = slice(None,-2,None)</div><div> sl_right = slice(2,None,None)</div><div> div_fac = 2.0</div><div> new_field = np.zeros(data["Pressure"].shape, dtype='float64')</div><div> ds = div_fac * data['dx'].flat[0]</div><div> new_field[1:-1,1:-1,1:-1] = data["Pressure"][sl_right,1:-1,1:-1]/ds</div><div> new_field[1:-1,1:-1,1:-1] -= data["Pressure"][sl_left ,1:-1,1:-1]/ds</div><div> return new_field</div><div>def _convertgradPressure(data):</div><div> return 1.0/data.convert("cm")</div><div>add_field("gradPX",function=_gradPX,convert_function=_convertgradPressure,validators=[ValidateSpatial(1, ["Pressure"])])</div></div><div><br></div><div><br></div><div>It gives me the array shape error when I use HydroMethod=2 in the simulation. Interestingly, if I use a different HydroMethod, the script will always try to make a slice plot even if I want a projection of gradPX. </div><div><br></div><div>Yuan</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jul 20, 2015 at 6:03 PM, Suoqing Ji <span dir="ltr"><<a href="mailto:suoqing@physics.ucsb.edu" target="_blank">suoqing@physics.ucsb.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Hi Yuan,<div><br></div><div>Unfortunately I can not reproduce this error with either yt 2.6 or yt 3… If you did not mind, could you provide your script and dataset (or send me offline) so I can take a quick look at it?</div><div><span class=""><br><div>
<div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div>Best wishes,</div><div>--</div><div>Suoqing JI</div><div>Ph.D Candidate</div><div>Department of Physics</div><div>University of California, Santa Barbara</div><div><a href="http://web.physics.ucsb.edu/~suoqing" target="_blank">http://web.physics.ucsb.edu/~suoqing</a></div></div></div></div></div></div></div></div>
</div>
<br></span><div><div class="h5"><div><blockquote type="cite"><div>On Jul 20, 2015, at 2:29 PM, Yuan Li <<a href="mailto:bear0980@gmail.com" target="_blank">bear0980@gmail.com</a>> wrote:</div><br><div><div dir="ltr">Hi Suoqing,<br><br>I added validators=[ValidateSpatial(1, ["Pressure"])]<br><br>Now the error message is<br>ValueError: could not broadcast input array from shape (64,64,66) into shape (64,64,64)<div><br></div><div>Yuan</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Jul 19, 2015 at 6:25 PM, Suoqing Ji <span dir="ltr"><<a href="mailto:suoqing@physics.ucsb.edu" target="_blank">suoqing@physics.ucsb.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Hi Yuan,<div><br></div><div>Just to check, have you included the “ValidateSpatial" object when adding the field? See the definition of DivV field: <a href="http://yt-project.org/docs/2.6/analyzing/creating_derived_fields.html#some-more-complicated-examples" target="_blank">http://yt-project.org/docs/2.6/analyzing/creating_derived_fields.html#some-more-complicated-examples</a></div><div><br></div><div><div>
<div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div>Best wishes,</div><div>--</div><div>Suoqing JI</div><div>Ph.D Candidate</div><div>Department of Physics</div><div>University of California, Santa Barbara</div><div><a href="http://web.physics.ucsb.edu/~suoqing" target="_blank">http://web.physics.ucsb.edu/~suoqing</a></div></div></div></div></div></div></div></div>
</div>
<br><div><blockquote type="cite"><div><div><div>On Jul 19, 2015, at 2:52 PM, Yuan Li <<a href="mailto:yuan@astro.columbia.edu" target="_blank">yuan@astro.columbia.edu</a>> wrote:</div><br></div></div><div><div><div><div dir="ltr">Dear all,<div><br></div><div><span style="font-size:13px">I am using yt2.6. I was trying to define a field of the temperature gradient like gradPressureX, but I got an error message:</span><br style="font-size:13px"><span style="font-size:13px">...</span><br style="font-size:13px"><span style="font-size:13px">new_field[1:-1,1:-1,1:-1] = data["Temperature"][sl_right,</span><span style="font-size:13px">1:-1,1:-1]/ds</span><br style="font-size:13px"><span style="font-size:13px">ValueError: could not broadcast input array from shape (62,62,64) into shape (62,62,62)</span><br style="font-size:13px"><br style="font-size:13px"><span style="font-size:13px">I then tried to re-define gradPressureX by copying the original definition (</span><a href="http://yt-project.org/docs/2.6/reference/field_list.html#gradpressurex" target="_blank">http://yt-project.org/docs/2.6/reference/field_list.html#gradpressurex</a>) <span style="font-size:13px">and giving it a different name, and I got the same error. The original gradPressureX works fine.</span><span style="font-size:13px"> </span></div><div><br style="font-size:13px"><span style="font-size:13px">Is it because I did not load some specific module or something?</span><div style="font-size:13px"><br></div><div style="font-size:13px"><br></div><div style="font-size:13px">Thank you!</div></div><div style="font-size:13px">Yuan</div></div></div></div>
_______________________________________________<br>yt-users mailing list<br><a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br><a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br></div></blockquote></div><br></div></div><br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">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>
_______________________________________________<br>yt-users mailing list<br><a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br><a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br></div></blockquote></div><br></div></div></div></div><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>