[yt-users] Projection of 2D cylindrical FLASH data

[Suoqing JI]

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.html>
-------------- next part --------------
_______________________________________________
yt-users mailing list
yt-users at lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org