<div dir="ltr">Thanks Guys,<div><br></div><div>     Kacper, that looks like a work-around for now, thanks.</div><div><br></div><div>         Munier</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Aug 7, 2013 at 2:53 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"><div class="im">On 07.08.2013 20:37, Munier Azzam Salem wrote:<br>
> Hey guys,<br>
><br>
>        I've been trying to make slices of density only looking at cooler<br>
> gas (i.e. with a cut in Temperature). I've found a nice way to do this for<br>
> simple projections. But my strategy breaks for Off-Axis projections.<br>
><br>
> Specifically, this works fine:<br>
</div>> *<br>
> *<br>
> *field_parameters={'field_cuts':["grid['Temperature'] < 1e6 "]})*<br>
> *cold_proj =<br>
> ProjectionPlot(pf,0,'Density',field_parameters=field_parameters)*<br>
> *               *<br>
><br>
> but this does not:<br>
><br>
> *cold_proj = OffAxisProjectionPlot(pf,axis,'Density',north_vector=north,*<br>
> *                                  field_parameters= field_parameters)*<br>
<div class="im">><br>
> sensibly enough, since field_parameters is not a keyword argument for off<br>
> axis projections.<br>
><br>
> Is there a decent workaround for this? Or would this be simple to add into<br>
> yt's functionality?<br>
<br>
</div>Hi,<br>
I think that you could define auxiliary derived field like<br>
<br>
@derived_field(name = "cut_dens",<br>
   validators = [ValidateSpatial(1)],<br>
   units = r"g/cm^3")<br>
def _cut_dens(field, data):<br>
   cut = data['Density'].copy()<br>
   cut[data['Temperature'] < 1.0e6] = 0.0<br>
   return cut<br>
<br>
and make OffAxisProjection of that.<br>
Cheers,<br>
Kacper<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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>Munier A. Salem // 845.489.6450
</div>