[yt-users] Getting yt to ignore empty bins in Radial profile

[Matthew Turk]

Hi Kearn,

On Thu, Nov 21, 2013 at 7:41 AM,  <k.grisdale at surrey.ac.uk> wrote:
> Dear yt Users
>
> I am plotting some radial profiles from Ramses data. However some of the bins that yt is using are too small for the course grid from the Ramses data and seem to be missing cell centres from the data.
>
> Is there anyway I can either change the bin size between certain radii or get yt to ignore data points with values of 0 and just contact the bins with data?
>

I think right now, it will just do either linearspace or logspace
spacing for the bins.

But, there's also the additional field "UsedBins" which is a boolean
of the bins that contain items and the ones that don't.  So for
instance, this:

m = rad["UsedBins"]
rad["Radius"][m]
rad["Density"][m]

will just give back the bins that have things in them.

Does that help?

-Matt

> Thanks in advance.
>
> Kearn
>
> PS: The commands that I am using are below in case they are needed:
>
> p6 = load("Ramses_sedov_3d/HD_Runs/l9/output_00006/info_00006.txt", fields = ["Density","x-velocity", "y-velocity", "z-velocity", "Pressure"])
> cen = [0, 0, 0]
>
> sphere = p6.h.sphere(cen, 0.5/p6['cm'])
> r_max = sphere['Radius'][sphere['Density'] == sphere['Density'].max()]
> sphere.set_field_parameter('r_max', r_max)
>
> rsh = r_max[0]
> rad_profile = BinnedProfile1D(sphere, 1500, 'Radius', 0.0, 1.4, log_space=False)
> rad_profile.add_fields("VelocityMagnitude")
> rad_profile.add_fields("Density")
> rad_profile.add_fields("Pressure”)
>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org