<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Dec 19, 2014 at 10:14 AM, José Utreras <span dir="ltr"><<a href="mailto:jutreras@das.uchile.cl" target="_blank">jutreras@das.uchile.cl</a>></span> wrote:<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr">Thanks for your answer John<br>I have to define  sl_right, <span>sl_center</span> :global name 'sl_right' is not defined<br><span></span><div><br></div><div>I'm not really sure how to do it. I want to do something like these:<br></div><div><br>vx = data["gas", "velocity_x"][sl_right,sl_center,sl_center]<br>vx+= data["gas", "velocity_x"][sl_center,sl_right,sl_center]<br></div></div></blockquote><div><br></div><div>Hi Jose,</div><div><br></div><div>This code looks like came from the definition of the vorticity field.  You're correct that ftype, sl_right, sl_center, and sl_left aren't defined inside the voriticity field definition.  That's because the voriticity field is a closure (a special type of python function that is able to use varibles defined in the enclosing scope around the function definition).  I'd encourage you to take a look at yt/fields/fluid_vector_fields.py to see how this works inside yt:</div><div><br></div><div><a href="https://bitbucket.org/yt_analysis/yt/src/5007c9a46622118abab89dc153dbfb6591cc5177/yt/fields/fluid_vector_fields.py?at=yt">https://bitbucket.org/yt_analysis/yt/src/5007c9a46622118abab89dc153dbfb6591cc5177/yt/fields/fluid_vector_fields.py?at=yt</a><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div>...<br><br></div><div>add the x velocity of 6 neighbour cells.<br><br></div></div></blockquote><div><br></div><div>Your field definition will need a ValidateSpatial field validator (like vorticity). ValidateSpatial will make it so data is passed into a field definition along with a user-selectable number of ghost zones.</div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div></div><div>Thanks<br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2014-12-19 13:31 GMT-03:00 John ZuHone <span dir="ltr"><<a href="mailto:jzuhone@gmail.com" target="_blank">jzuhone@gmail.com</a>></span>:<div><div class="h5"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Jose,<br>
<br>
Presuming that you have sl_right and sl_center defined in your field already, you should be able to just substitute “gas” for ftype, since the velocity fields are “gas" fields.<br>
<br>
Best,<br>
<br>
John<br>
<div><div><br>
> On Dec 19, 2014, at 11:04 AM, José Utreras <<a href="mailto:jutreras@das.uchile.cl" target="_blank">jutreras@das.uchile.cl</a>> wrote:<br>
><br>
> Hi everyone,<br>
><br>
> I want to write a new field based on the neighbour cells like this:<br>
><br>
> vx = data[ftype, "velocity_x"][sl_right,sl_center,sl_center]<br>
><br>
> but I got " global name 'ftype' is not defined". I suppose that I have the same problem with sl_left, sl_right , sl_center.<br>
><br>
> Thanks,<br>
><br>
> Jose<br>
</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>
<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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</blockquote></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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div></div></div>