Grids on Circular Images

[jonhulberg]

Hello,

I've been using radmc3d to create 1D circular images of an AGB star. My plan has been to take the radial intensity profile from the circular images and do a Hankel transformation to get interferometric visibilities. I realized the intensity profile of the images seems to be calculated according to the radial grid the simulations are performed on - this has caused an issue, because the central star cannot lie inside the grid. The intensity profiles I generate do not cover the central star very well, as shown in the figure below.

The image shows intensity vs. radius at two different wavelengths. The vertical blue line marks where the grid begins just outside the stellar radius. Besides the point at the center, there is no coverage of the region interior to the stellar radius. Is there any way to ask radmc3d to create a circular image at a specified set of radial distances? I've noticed this problem doesn't occur if I create a Cartesian image, so I have a backup plan, but it would be nice to use a radial profile that guarantees spherical symmetry.

Thanks for your help/for writing and maintaining this code

Jon Hulberg