prop = {
'dens_H2': density.total,
'temp_gas': temperature.total,
'vel_x': vel.x,
'vel_y': vel.y,
'vel_z': vel.z,
'abundance': abundance,
'gtdratio': gtdratio
}
lime = rt.Lime(GRID)
lime.finalmodel(prop)
#-----------------------------
#PRINTING resultant PROPERTIES
#-----------------------------
Model.PrintProperties(density, temperature, GRID)
#-------
#TIMING
#-------
print('Ellapsed time: %.3fs' % (time.time() - t0))
print(
'-------------------------------------------------\n-------------------------------------------------\n'
)
#----------------------------------------
#3D PLOTTING (weighting with temperature)
#----------------------------------------
tag = 'Burger_Tapering'
dens_plot = density.total / 1e6
#---------------
#GRID Definition
#---------------
sizex = sizey = sizez = 2600 * u.au
Nx = Ny = Nz = 63 #Number of divisions for each axis
GRID = Model.grid([sizex, sizey, sizez], [Nx, Ny, Nz], rt_code='radmc3d')
NPoints = GRID.NPoints #Final number of nodes in the grid
#-------------------
#PHYSICAL PROPERTIES
#-------------------
density = Model.density_Constant(r_max, GRID, envDens = dens_e)
temperature = Model.temperature_Constant(density, GRID, envTemp = t_e, backTemp = 2.725)
Model.PrintProperties(density, temperature, GRID) #Printing resultant properties (mass, mean temperature, etc)
#----------------------
#WRITING RADMC-3D FILES
#----------------------
prop = {'dens_e': density.total,
'dens_ion': density.total,
'temp_gas': temperature.total}
Rad = rt.Radmc3dDefaults(GRID)
Rad.freefree(prop)
#------------------------------------
#3D PLOTTING (weighting with density)
#------------------------------------
tag = 'ctsphere_HII'