[yt-users] Field Cuts in Off-Axis Projections

[Munier Azzam Salem]

Thanks Guys,

     Kacper, that looks like a work-around for now, thanks.

         Munier


On Wed, Aug 7, 2013 at 2:53 PM, Kacper Kowalik <xarthisius.kk at gmail.com>wrote:

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


-- 
Munier A. Salem // 845.489.6450
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20130807/56b19ecb/attachment.html>