colorlabel=r'$\rho [cm^{-3}]$',
output='3Dpoints%s.png' % tag,
show=False)
#----------------------------------------
#2D PLOTTING (Density and Temperature)
#----------------------------------------
vmin, vmax = np.array([9e15, 5e19]) / 1e6
norm = colors.LogNorm(vmin=vmin, vmax=vmax)
Plot_model.plane2D(GRID,
dens_plot,
axisunit=u.au,
cmap='ocean_r',
plane={'z': 0 * u.au},
norm=norm,
colorlabel=r'$[\rm cm^{-3}]$',
output='DensMidplane_%s.png' % tag,
show=False)
vmin, vmax = np.array([5e14, 1e18]) / 1e6
norm = colors.LogNorm(vmin=vmin, vmax=vmax)
Plot_model.plane2D(GRID,
dens_plot,
axisunit=u.au,
cmap='ocean_r',
plane={'y': 0 * u.au},
norm=norm,
colorlabel=r'$[\rm cm^{-3}]$',
}
Model.PrintProperties(density, temperature, GRID)
print('Ellapsed time: %.3fs' % (time.time() - t0))
print(
'-------------------------------------------------\n-------------------------------------------------\n'
)
#---------------
#PLOTTING 2D z=0
#---------------
Plot_model.plane2D(GRID,
prop['dens_H'],
axisunit=u.au,
plane={'z': 0},
norm=LogNorm(vmin=1e9, vmax=1e14),
colorlabel=r'$n_{H}$ $\rm{[m^{-3}]}$',
cmap='cubehelix',
output='densH_2D.png',
show=True)
Plot_model.plane2D(GRID,
prop['temp_gas'],
axisunit=u.au,
plane={'z': 0},
norm=LogNorm(vmin=1e1, vmax=1e3),
colorlabel=r'T [K]',
cmap='hot',
output='tempgas_2D.png',
show=True)
[-u.TSun]
], #flux --> if negative, radmc assumes the input number as the blackbody temperature of the star
lam=wavelength_intervals,
nxx=wavelength_divisions)
radmc.write_wavelength_micron(
lam=wavelength_intervals, nxx=wavelength_divisions
) #lam --> wavelengths in microns, nxx --> number of divisions in between wavelengths
Model.PrintProperties(density, temperature, GRID)
#---------------
#PLOTTING 2D z=0
#---------------
Plot_model.plane2D(GRID,
prop['dens_H2'],
axisunit=u.au,
plane={'z': 0},
norm=LogNorm(vmin=1e9, vmax=1e13),
output='densH2_2D.png')
Plot_model.plane2D(GRID,
prop['temp_gas'],
axisunit=u.au,
plane={'z': 0},
norm=LogNorm(vmin=1e1, vmax=1e3),
output='tempgas_2D.png')
#-------
#TIMING
#-------
print('Ellapsed time: %.3fs' % (time.time() - t0))
print(
'-------------------------------------------------\n-------------------------------------------------\n'
colorlabel=r'${\rm log}_{10}(r [au])$',
output='3Dpoints%s.png' % tag,
show=False)
#---------------------
#2D PLOTTING (Density)
#---------------------
vmin, vmax = np.array([2e13, 1e19]) / 1e6
norm = colors.LogNorm(vmin=vmin, vmax=vmax)
Plot_model.plane2D(GRID,
dens_plot,
axisunit=u.au,
cmap='jet',
plane={'z': 0 * u.au},
norm=norm,
colorlabel=r'$[\rm cm^{-3}]$',
output='DensMidplane_%s.png' % tag,
show=False)
vmin, vmax = np.array([1e13, 3e17]) / 1e6
norm = colors.LogNorm(vmin=vmin, vmax=vmax)
Plot_model.plane2D(GRID,
dens_plot,
axisunit=u.au,
cmap='jet',
plane={'y': 0 * u.au},
norm=norm,
colorlabel=r'$[\rm cm^{-3}]$',