[yt-users] Projection of 2D cylindrical FLASH data
Suoqing JI
suoqing at physics.ucsb.edu
Thu May 14 16:47:43 PDT 2015
Hi Ken,
Just to clarify: when you referred to the column density plot, do you mean the 1D plot of column density which is accumulated along z direction? If so, you could something like the following:
profile = BinnedProfile1D(mydisk, Nbin, 'cylindrical_r', rmin, rmax, log_space=True, lazy_reader=True, end_collect=False)
profile.add_fields('cell_mass', weight=None)
R = profile['cylindrical_r’]
Sigma = profile[‘cell_mass']
R_edge = np.logspace(np.log10(rmin), np.log10(rmax), num=Nbin+1)
for i in range(0, Nbin): Sigma[i] = Sigma[i] / (np.pi*(R_edge[i+1]**2 - R_edge[i]**2.))
Best wishes,
--
Suoqing JI
Ph.D Student
Department of Physics
University of California, Santa Barbara
http://web.physics.ucsb.edu/~suoqing
> On May 14, 2015, at 4:31 PM, Ken Shen <kenshen at astro.berkeley.edu> wrote:
>
> Hi all, I'd like to make column density plots of my 2D cylindrical FLASH data. Ideally, these would be at arbitrary angles to the z-axis. I'm new to yt, so I could very well have missed it, but I couldn't seem to find this in the documentation. Any hints?
>
>
> Thanks,
> Ken
> _______________________________________________
> 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/20150514/af4445b4/attachment.htm>
-------------- next part --------------
_______________________________________________
yt-users mailing list
yt-users at lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
More information about the yt-users
mailing list