[yt-users] FLASH 1D Data Radial Profiles
John Zuhone
jzuhone at gmail.com
Wed Jul 9 16:06:37 PDT 2014
Hi Jill,
No, not your bad at all! 2.x is still the most stable version. Just
wanted to make sure that code was out there.
Best,
John
On Wed, Jul 9, 2014 at 7:01 PM, Jill Naiman <jnaiman at gmail.com> wrote:
> My bad, yeah, this is yt-2.x stuff. Also, ninterp is just some number to
> bin things in. In my code I have the note:
>
> ninterp = 1400 # for example: nblockx =1500 => ninterp ~ 500-1400
>
>
>
> On Wed, Jul 9, 2014 at 3:58 PM, John ZuHone <jzuhone at gmail.com> wrote:
>>
>> Hi Melinda and all,
>>
>> Just for clarification, Jill's code works in yt 2.x. The 3.0 beta version
>> (to be released soon) of this code would be:
>>
>> import yt
>> pf = yt.load(fname)
>> sphere = pf.sphere([0.,0.,0.], (pf.domain_right_edge[0]/3.1e18, "pc"))
>> rad_profile = yt.create_profile(sphere, "radius",
>> ["density","temperature","velocity_x"],
>> extrema={"radius":(0.0,pf.domain_right_edge[0]/3.1e18*yt.units*pc)},
>> n_bins=ninterp,
>> units={"radius":"pc"}, logs={"radius":False})
>>
>> radius = rad_profile.x
>>
>> Then you can access the fields like so:
>>
>> rad_profile["density"]
>> rad_profile["temperature"]
>> rad_profile["velocity_x"]
>>
>> Also, yt 3.0 is aware of spherical geometries, so you could also just use
>> the unbinned quantities:
>>
>> sp = pf.sphere([0.0,0.0,0.0], (pf.domain_right_edge[0]/3.1e18, "pc"))
>> radius = sp["radius"].in_units("pc")
>> density = sp["density"]
>>
>> etc...
>>
>> Best,
>>
>> John ZuHone
>>
>> On Jul 9, 2014, at 5:25 PM, Jill Naiman <jnaiman at gmail.com> wrote:
>>
>> 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
>>>
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>>
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
--
John ZuHone
Postdoctoral Researcher
NASA/Goddard Space Flight Center
jzuhone at gmail.com
john.zuhone at nasa.gov
More information about the yt-users
mailing list