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

[Kacper Kowalik]

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


-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 901 bytes
Desc: OpenPGP digital signature
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20130807/437eceff/attachment.sig>