[yt-users] FLASH 1D Data Radial Profiles

Jill Naiman jnaiman at gmail.com
Wed Jul 9 15:25:25 PDT 2014


The following has worked for me:

pf = load(fname)
sphere = pf.h.sphere([0.,0.,0.], (pf.domain_right_edge[0]/3.1e18, "pc"))
rad_profile = BinnedProfile1D(sphere, ninterp, "Radiuspc", 0.0, \
                                  pf.domain_right_edge[0]/3.1e18,
log_space=False)
radius = rad_profile['Radiuspc']

rad_profile.add_fields("Density")
rad_profile.add_fields("Temperature")
rad_profile.add_fields("x-velocity")


And then I just plot using "fig"


On Wed, Jul 9, 2014 at 2:47 PM, Melinda Soares-Furtado <
msoares.physics at gmail.com> wrote:

> Is there a simple way to create density/temp/velocity profiles using 1D
> FLASH data? I've done this in the past for 3-dimensional data, which is
> very user friendly.
>
> Is the best method of attack is dumping the data to np arrays and then
> plotting from these arrays?
>
> Regards,
> Melinda Soares-Furtado
> msoares at princeton.edu
> http://cargo.ucsc.edu/msoares
> 415-860-0438
>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20140709/9d501324/attachment.htm>


More information about the yt-users mailing list