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@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@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org




--
Munier A. Salem // 845.489.6450