Profile each step of model generation to see where things should be optimized
In [1]:
using model
In [3]:
model.plot_vel(params)
In [2]:
write_grid()
In [3]:
@time write_model(params)
In addition to the gridding chosen by the model, there are many options here to tweak
Pg 25: 2D axially-symmetric with mirror symmetry. Do we just specify it and see what happens? Or does this do it automatically? If the largest value is = pi/2, then it is switched on automatically (should say in the output).
try running a single channel vs. multiple simultaneously
How many npix should we really be using? Sean says we want to oversample the SMA beam by a factor of 5x10, so let's try pixels that are 0.05"
At this point, I think any tweaking is really going to be either in the number of pixels or number of grid cells.
In [13]:
12/0.05
Out[13]:
In [4]:
using constants
incl = 33. #33. # deg. 0 deg = face on, 90 = edge on.
vel = 0.0 # km/s
PA = 90 - 73. # 73 deg. Position angle, runs counter clockwise, due to looking at sky.
npix = 240 # number of pixels, can alternatively specify x and y separately
lam0 = cc/230.538e9 * 1e6 # [microns]
Out[4]:
In [5]:
#loads the camera_wavelength_micron.inp file
@time run(`radmc3d image incl $incl posang $PA vkms $vel npix $npix loadlambda`)
In [6]:
@time run(`radmc3d image incl $incl posang $PA vkms $vel npix 240 lambda $lam0`)
In [ ]:
using read_image
In [ ]:
In [ ]: