DART-Ray
|
Contains the subroutines to specify the stellar emission and dust density distribution for the analytical galaxy geometry. More...
Functions/Subroutines | |
subroutine | read_input_galaxy |
Reads input file for the grid creation for the galaxy 2D model. More... | |
subroutine | initialize_input_galaxy |
Initialize variables specific for the galaxy grid creation program. More... | |
subroutine | check_lambda_list_galaxy |
Checks that lambda_arr contains at least two wavelengths of which one is equal to 0.443 um (B-band). This wavelength is necessary for calc_scaling_factors_dust() to work properly. More... | |
subroutine | check_grid_type |
subroutine | calc_luminosities_pt11 |
Calculates luminosity of the different stellar components depending on wavelength and intensity parameters old and sfr. The unit luminosities are those from Popescu et al. 2011. More... | |
subroutine | calc_luminosities |
Calculates luminosity of the different stellar components depending on wavelength and intensity parameters old and sfr. The unit luminosities are those from the input files file_old_star_sed() and file_young_star_sed(). More... | |
subroutine | calc_scaling_factors_stars |
subroutine | calc_scaling_factors_dust |
subroutine | interpolate_kext |
Interpolates value of kext given wavelength lambda. More... | |
real(kind=real64) function | av_disk (x, y, z, cellsize, disk_comp, disk_type_ID) |
real(kind=real64) function | av_star_bulge (x, y, z, cellsize) |
subroutine | create_grid_arrays_galaxy |
subroutine | fix_dens_stars_arrays |
subroutine | print_3d_grid_file_galaxy |
subroutine | print_profiles_stars |
subroutine | set_galaxy_input |
subroutine | assign_dens_to_parent |
subroutine | read_stellar_sed |
Reads the input stellar emission SEDs in unit luminosities from the input files file_old_star_sed() and file_young_star_sed(). More... | |
subroutine | assign_disk_type_id (disk_type, disk_type_ID) |
Assigns the right value of the disk_type_ID() depending on the input disk_type(). Using numbers instead of character variables is more efficient in IF statement and less error prone while modifying the code. More... | |
Variables | |
real(kind=real64) | lambda_in |
real(kind=real64) | lambda_min |
real(kind=real64) | lambda_max |
real(kind=real64) | old |
real(kind=real64) | lnu_old |
real(kind=real64) | hs_disk_b |
real(kind=real64) | zs_disk |
real(kind=real64) | zs_disk_r1 |
real(kind=real64) | zs_disk_rsun |
real(kind=real64) | hsin |
real(kind=real64) | hs_disk |
real(kind=real64) | hs_disk2 |
real(kind=real64) | eta_disk0 |
real(kind=real64) | chi_disk |
real(kind=real64), dimension(:), allocatable | hs_disk_arr |
real(kind=real64), dimension(:), allocatable | hs_disk2_arr |
integer, dimension(:), allocatable | id_hs_disk_arr |
integer, parameter | maxsize_hs_disk_arr = 1000 |
real(kind=real64) | theta_disk_ellipt |
real(kind=real64) | sfr |
real(kind=real64) | lnu_sf |
real(kind=real64) | hs_tdisk |
real(kind=real64) | hs_tdisk2 |
real(kind=real64) | zs_tdisk |
real(kind=real64) | hs1in |
real(kind=real64) | eta_tdisk0 |
real(kind=real64) | zs_tdisk_r1 |
real(kind=real64) | zs_tdisk_rsun |
real(kind=real64) | chi_tdisk |
real(kind=real64) | theta_tdisk_ellipt |
real(kind=real64) | reff |
real(kind=real64) | lnu_bulge |
real(kind=real64) | acap_bulge |
real(kind=real64) | ellipt |
real(kind=real64) | ellipt_xy |
real(kind=real64) | mtrunc |
real(kind=real64) | bd_ratio |
real(kind=real64) | eta_bulge0 |
integer | nsersic |
real(kind=real64) | theta_bulge |
real(kind=real64) | lnu_tot |
real(kind=real64) | hd_disk |
real(kind=real64) | hd_disk2 |
real(kind=real64) | zd_disk |
real(kind=real64) | zd_disk_r1 |
real(kind=real64) | zd_disk_rsun |
real(kind=real64) | hdin |
real(kind=real64) | tau1 |
real(kind=real64) | kext_disk0 |
real(kind=real64) | chi_dust_disk |
real(kind=real64) | theta_dust_disk_ellipt |
real(kind=real64) | hd_tdisk |
real(kind=real64) | hd_tdisk2 |
real(kind=real64) | zd_tdisk |
real(kind=real64) | hd1in |
real(kind=real64) | tau2 |
real(kind=real64) | kext_tdisk0 |
real(kind=real64) | zd_tdisk_r1 |
real(kind=real64) | zd_tdisk_rsun |
real(kind=real64) | chi_dust_tdisk |
real(kind=real64) | theta_dust_tdisk_ellipt |
real(kind=real64) | rtrun |
real(kind=real64) | rsun |
real(kind=real64) | max_z |
real(kind=real64) | max_rad |
real(kind=real64) | sha |
real(kind=real64) | sha1 |
real(kind=real64) | omega_max |
integer, parameter | lcar_type =30 |
character(len=lcar_type) | old_disk_type |
integer | old_disk_type_id |
character(len=lcar_type) | young_disk_type |
integer | young_disk_type_id |
character(len=lcar_type) | thick_disk_type |
integer | thick_disk_type_id |
character(len=lcar_type) | thin_disk_type |
integer | thin_disk_type_id |
character(len=lcar_type) | grid_type |
integer, parameter | expr_expz_id = 0 |
integer, parameter | expr_sech2z_id = 1 |
integer, parameter | flared_expz_id = 2 |
integer, parameter | flared_sech2z_id = 3 |
integer, parameter | ellipt_expr_expz_id = 4 |
integer, parameter | ellipt_expr_sech2z_id = 5 |
real(kind=real64), dimension(:), allocatable | dens_disk |
real(kind=real64), dimension(:), allocatable | dens_tdisk |
real(kind=real64), dimension(:), allocatable | dens_bulge |
real(kind=real64), dimension(:), allocatable | dens_dust_disk |
real(kind=real64), dimension(:), allocatable | dens_dust_tdisk |
real(kind=real64) | omega_cell |
integer, parameter | step_int =10 |
character(len=lcar) | file_old_star_sed |
character(len=lcar) | file_young_star_sed |
real(kind=real64), dimension(:), allocatable | lambda_arr_old_sed |
integer | nlambda_old_sed |
real(kind=real64), dimension(:), allocatable | lambda_arr_sf_sed |
integer | nlambda_sf_sed |
real(kind=real64), dimension(:), allocatable | lnu_old_unit |
real(kind=real64), dimension(:), allocatable | lnu_sf_unit |
real(kind=real64) | z_subd_lim |
real(kind=real64) | z_subd_lim2 |
real(kind=real64) | r_subd_lim |
character(len=lcar_type) | subdivision_criteria |
Contains the subroutines to specify the stellar emission and dust density distribution for the analytical galaxy geometry.
subroutine user_routines_galaxy::assign_dens_to_parent | ( | ) |
subroutine user_routines_galaxy::assign_disk_type_id | ( | character(len=lcar_type) | disk_type, |
integer | disk_type_ID | ||
) |
Assigns the right value of the disk_type_ID() depending on the input disk_type(). Using numbers instead of character variables is more efficient in IF statement and less error prone while modifying the code.
real(kind=real64) function user_routines_galaxy::av_disk | ( | real(kind=real64) | x, |
real(kind=real64) | y, | ||
real(kind=real64) | z, | ||
real(kind=real64) | cellsize, | ||
character(len=lcar_type) | disk_comp, | ||
integer | disk_type_ID | ||
) |
real(kind=real64) function user_routines_galaxy::av_star_bulge | ( | real(kind=real64) | x, |
real(kind=real64) | y, | ||
real(kind=real64) | z, | ||
real(kind=real64) | cellsize | ||
) |
subroutine user_routines_galaxy::calc_luminosities | ( | ) |
Calculates luminosity of the different stellar components depending on wavelength and intensity parameters old and sfr. The unit luminosities are those from the input files file_old_star_sed() and file_young_star_sed().
subroutine user_routines_galaxy::calc_luminosities_pt11 | ( | ) |
Calculates luminosity of the different stellar components depending on wavelength and intensity parameters old and sfr. The unit luminosities are those from Popescu et al. 2011.
subroutine user_routines_galaxy::calc_scaling_factors_dust | ( | ) |
subroutine user_routines_galaxy::calc_scaling_factors_stars | ( | ) |
subroutine user_routines_galaxy::check_grid_type | ( | ) |
subroutine user_routines_galaxy::check_lambda_list_galaxy | ( | ) |
Checks that lambda_arr contains at least two wavelengths of which one is equal to 0.443 um (B-band). This wavelength is necessary for calc_scaling_factors_dust() to work properly.
subroutine user_routines_galaxy::create_grid_arrays_galaxy | ( | ) |
subroutine user_routines_galaxy::fix_dens_stars_arrays | ( | ) |
subroutine user_routines_galaxy::initialize_input_galaxy | ( | ) |
Initialize variables specific for the galaxy grid creation program.
subroutine user_routines_galaxy::interpolate_kext | ( | ) |
Interpolates value of kext given wavelength lambda.
subroutine user_routines_galaxy::print_3d_grid_file_galaxy | ( | ) |
subroutine user_routines_galaxy::print_profiles_stars | ( | ) |
subroutine user_routines_galaxy::read_input_galaxy | ( | ) |
Reads input file for the grid creation for the galaxy 2D model.
subroutine user_routines_galaxy::read_stellar_sed | ( | ) |
Reads the input stellar emission SEDs in unit luminosities from the input files file_old_star_sed() and file_young_star_sed().
subroutine user_routines_galaxy::set_galaxy_input | ( | ) |
real(kind=real64) user_routines_galaxy::acap_bulge |
acap_bulge | Inner truncation radius for the bulge profile |
real(kind=real64) user_routines_galaxy::bd_ratio |
bd_ratio | Bulge-to-Disk luminosity ratio. Used to define luminosity scaling for the bulge, that is, bd_ratio*old*lnu_old. |
real(kind=real64) user_routines_galaxy::chi_disk |
chi_disk | Parameter used in the "flared" profile to specify the steepness of the radial profile within the inner radius hsin(). |
real(kind=real64) user_routines_galaxy::chi_dust_disk |
chi_dust_disk | Parameter used in the "flared" profile to specify the steepness of the radial profile within the inner radius hdin(). |
real(kind=real64) user_routines_galaxy::chi_dust_tdisk |
chi_dust_tdisk | Parameter used in the "flared" profile to specify the steepness of the radial profile within the inner radius hd1in(). |
real(kind=real64) user_routines_galaxy::chi_tdisk |
chi_tdisk | Parameter used in the "flared" profile to specify the steepness of the radial profile within the inner radius hs1in(). |
real(kind=real64), dimension(:), allocatable user_routines_galaxy::dens_bulge |
dens_bulge | Volume emissivity for the bulge. |
real(kind=real64), dimension(:), allocatable user_routines_galaxy::dens_disk |
dens_disk | Volume emissivity for the old stellar disk. |
real(kind=real64), dimension(:), allocatable user_routines_galaxy::dens_dust_disk |
dens_dust_disk | Dust density for the thick dust disk. |
real(kind=real64), dimension(:), allocatable user_routines_galaxy::dens_dust_tdisk |
dens_dust_tdisk | Dust density for the thin dust disk. |
real(kind=real64), dimension(:), allocatable user_routines_galaxy::dens_tdisk |
dens_tdisk | Volume emissivity for the young stellar disk. |
real(kind=real64) user_routines_galaxy::ellipt |
ellipt | Bulge ellipticity. The normalized radial coordinate m for the bulge volume emissivity is defined as
|
integer, parameter user_routines_galaxy::ellipt_expr_expz_id = 4 |
integer, parameter user_routines_galaxy::ellipt_expr_sech2z_id = 5 |
real(kind=real64) user_routines_galaxy::ellipt_xy |
ellipt_xy | Bulge XY ellipticity. The radial coordinate R for the bulge volume emissivity is defined as
|
real(kind=real64) user_routines_galaxy::eta_bulge0 |
eta_bulge0 | Central value of the bulge volume emissivity. |
real(kind=real64) user_routines_galaxy::eta_disk0 |
eta_disk0 | Central value of the old stellar disk volume emissivity. |
real(kind=real64) user_routines_galaxy::eta_tdisk0 |
eta_tdisk0 | Central value of the young stellar disk volume emissivity. |
integer, parameter user_routines_galaxy::expr_expz_id = 0 |
integer, parameter user_routines_galaxy::expr_sech2z_id = 1 |
character(len=lcar) user_routines_galaxy::file_old_star_sed |
file_old_star_sed | File containing the SED to be used by the old stellar components in units luminosities to be multiplied by the OLD parameter. |
character(len=lcar) user_routines_galaxy::file_young_star_sed |
file_young_star_sed | File containing the SED to be used by the young stellar components in units luminosities to be multiplied by the SFR parameter. |
integer, parameter user_routines_galaxy::flared_expz_id = 2 |
integer, parameter user_routines_galaxy::flared_sech2z_id = 3 |
character(len=lcar_type) user_routines_galaxy::grid_type |
Specify | which stellar emission components have to be included: 'all'= all components; 'disk' = only old stellar disk; 'tdisk' = only thin stellar disk; 'bulge' = only bulge component. |
real(kind=real64) user_routines_galaxy::hd1in |
hd1in | Inner radius for the thin dust disk. |
real(kind=real64) user_routines_galaxy::hd_disk |
hd_disk | Scale length of the thick dust disk or X'-axis scale length in the case of the ellipsoidal exponential profile |
real(kind=real64) user_routines_galaxy::hd_disk2 |
hd_disk2 | Y'-axis thick dust disk scale length for the ellipsoidal exponential profile. |
real(kind=real64) user_routines_galaxy::hd_tdisk |
hd_tdisk | Scale length of the thin dust disk or X'-axis scale length for the ellipsoidal exponential profile. |
real(kind=real64) user_routines_galaxy::hd_tdisk2 |
hd_tdisk2 | Y'-axis thin dust disk scale length for the ellipsoidal exponential profile. |
real(kind=real64) user_routines_galaxy::hdin |
hdin | Inner radius for the thick dust disk. |
real(kind=real64) user_routines_galaxy::hs1in |
hs1in | Inner radius of the young stellar disk. |
real(kind=real64) user_routines_galaxy::hs_disk |
hs_disk | Old disk scale length or X'-axis scale length if elliptical exponential model is used. |
real(kind=real64) user_routines_galaxy::hs_disk2 |
hs_disk2 | Y'-axis old disk scale length for the elliptical exponential model. |
real(kind=real64), dimension(:), allocatable user_routines_galaxy::hs_disk2_arr |
hs_disk2_arr | Same as hs_disk_arr() but for the Y'-axis scale length required in the elliptical exponential disk profile type. |
real(kind=real64), dimension(:), allocatable user_routines_galaxy::hs_disk_arr |
hs_disk_arr | Scale lenghts values for the old stellar disk at each wavelength of the input wavelength grid. In the case of the elliptical exponential profile, this is the X'-axis scale length. By default, all values are set equal to hs_disk_b. If different values are required, hs_disk_arr can be specified in the input. The wavelength indeces of the input hs_disk values have to be specified in id_hs_disk_arr() |
real(kind=real64) user_routines_galaxy::hs_disk_b |
hs_disk_b | B-band (4430 A) old stellar disk scale length. |
real(kind=real64) user_routines_galaxy::hs_tdisk |
hs_tdisk | Scale length of the young stellar disk or, in the case of the ellipsoidal exponential profile, X'-axis scale length |
real(kind=real64) user_routines_galaxy::hs_tdisk2 |
hs_tdisk2 | Y'-axis young stellar disk scale length for the ellipsoidal esponential profile. |
real(kind=real64) user_routines_galaxy::hsin |
hsin | Inner radius of the old stellar disk. |
integer, dimension(:), allocatable user_routines_galaxy::id_hs_disk_arr |
id_hs_disk_arr | Wavelength grid indeces of the values of hs_disk specified in the input hs_disk_arr() and hs_disk2_arr(). |
real(kind=real64) user_routines_galaxy::kext_disk0 |
kext_disk0 | Central value for the thick dust disk density profile. |
real(kind=real64) user_routines_galaxy::kext_tdisk0 |
kext_tdisk0 | Central value for the thin dust disk density profile. |
real(kind=real64), dimension(:), allocatable user_routines_galaxy::lambda_arr_old_sed |
lambda_arr_old_sed | Wavelength array in [um] for the old stellar luminosity sed |
real(kind=real64), dimension(:), allocatable user_routines_galaxy::lambda_arr_sf_sed |
lambda_arr_sf_sed | Wavelength array in [um] for the young stellar luminosity sed |
real(kind=real64) user_routines_galaxy::lambda_in |
lambda_in | Wavelength variable |
real(kind=real64) user_routines_galaxy::lambda_max |
lambda_max | Maximum wavelength used in the calculation of the lambda grids. |
real(kind=real64) user_routines_galaxy::lambda_min |
lambda_min | Minimum wavelength used in the calculation of the lambda grids. |
integer, parameter user_routines_galaxy::lcar_type =30 |
lcar_type | String variable length parameter. |
real(kind=real64) user_routines_galaxy::lnu_bulge |
lnu_bulge | Bulge volume emissivity. |
real(kind=real64) user_routines_galaxy::lnu_old |
lnu_old | Old stellar luminosity spectra. |
real(kind=real64), dimension(:), allocatable user_routines_galaxy::lnu_old_unit |
lnu_old_unit | Old stellar luminosity sed in unit luminosities for the old() parameter. |
real(kind=real64) user_routines_galaxy::lnu_sf |
lnu_sf | Young stellar disk luminosity. |
real(kind=real64), dimension(:), allocatable user_routines_galaxy::lnu_sf_unit |
lnu_sf_unit | Young stellar luminosity sed in unit luminosities for the young() parameter. |
real(kind=real64) user_routines_galaxy::lnu_tot |
lnu_tot | Model total luminosity. |
real(kind=real64) user_routines_galaxy::max_rad |
max_rad | Maximum R allowed for non zero stellar emissivity and dust density. |
real(kind=real64) user_routines_galaxy::max_z |
max_z | Maximum |z| allowed for non zero stellar emissivity and dust density. |
integer, parameter user_routines_galaxy::maxsize_hs_disk_arr = 1000 |
maxsize_hs_disk_arr | Maximum number of input values in hs_disk_arr() and hs_disk2_arr() |
real(kind=real64) user_routines_galaxy::mtrunc |
mtrunc | Outern truncation radius for the bulge in units of effective radii. |
integer user_routines_galaxy::nlambda_old_sed |
nlambda_old_sed | Number of elements in lambda_arr_old_sed() |
integer user_routines_galaxy::nlambda_sf_sed |
nlambda_sf_sed | Number of elements in lambda_arr_sf_sed() |
integer user_routines_galaxy::nsersic |
nsersic | Bulge sersic index (see inside subroutine av_star_bulge() for the volume emissivity formula used in each case). |
real(kind=real64) user_routines_galaxy::old |
old | Scaling factor applied to the input old stellar luminosity spectra. |
character(len=lcar_type) user_routines_galaxy::old_disk_type |
old_disk_type | Old stellar disk profile type. Possible choices:
|
integer user_routines_galaxy::old_disk_type_id |
old_disk_type_ID | ID number for old_disk_type(). |
real(kind=real64) user_routines_galaxy::omega_cell |
Solid | angle subtended by a cell |
real(kind=real64) user_routines_galaxy::omega_max |
omega_max | Maximum subtended solid angle of a cell seen from the position (rsun(), 0,0). To use this variable you have to make sure that the corresponding line in the function subdivision is used while compiling the program. |
real(kind=real64) user_routines_galaxy::r_subd_lim |
R_subd_lim | R coordinate below which cells have to be subdivided if their z is less than z_subd_lim() |
real(kind=real64) user_routines_galaxy::reff |
reff | Bulge effective radius. |
real(kind=real64) user_routines_galaxy::rsun |
rsun | Third radial position used in the "flared" profile. Called rsun because it has been used to set the sun position in Milky Way model. |
real(kind=real64) user_routines_galaxy::rtrun |
rtrun | Radial truncation radius for all disks. |
real(kind=real64) user_routines_galaxy::sfr |
sfr | Scaling factor applied to the young stellar disk luminosity spectra. |
real(kind=real64) user_routines_galaxy::sha |
real(kind=real64) user_routines_galaxy::sha1 |
integer, parameter user_routines_galaxy::step_int =10 |
character(len=lcar_type) user_routines_galaxy::subdivision_criteria |
subdivision_criteria | Subdivision criteria to be used. Choices are:
|
real(kind=real64) user_routines_galaxy::tau1 |
tau1 | B-band (4430 A) central face-on optical depth of the thick dust disk in the case of the pure exponential profile. In the case of the flared profiles, this is the optical depth one would have if the exponential radial profile for R>R1 would continue until R=0. |
real(kind=real64) user_routines_galaxy::tau2 |
tau2 | B-band (4430 A) central face-on optical depth of the thin dust disk in the case of the pure exponential profile. Other profiles still use this factor to set the dust density although its meaning is different. |
real(kind=real64) user_routines_galaxy::theta_bulge |
theta_bulge | Rotation angle around the Z-axis for the bulge profile. |
real(kind=real64) user_routines_galaxy::theta_disk_ellipt |
theta_disk_ellipt | Rotation angle around the Z-axis for the elliptical exponential profile for the old stellar component [deg] |
real(kind=real64) user_routines_galaxy::theta_dust_disk_ellipt |
theta_dust_disk_ellipt | Rotation angle around the Z-axis for the elliptical exponential profile for the thick dust disk component [deg]. |
real(kind=real64) user_routines_galaxy::theta_dust_tdisk_ellipt |
theta_dust_tdisk_ellipt | Rotation angle around the Z-axis for the elliptical exponential profile for the thin dust disk component [deg]. |
real(kind=real64) user_routines_galaxy::theta_tdisk_ellipt |
theta_tdisk_ellipt | Rotation angle around the Z-axis for the elliptical exponential profile for the young stellar component [deg]. |
character(len=lcar_type) user_routines_galaxy::thick_disk_type |
thick_disk_type | Thick dust disk profile type. Same choices as for old_disk_type(). |
integer user_routines_galaxy::thick_disk_type_id |
thick_disk_type_ID | ID number for thick_disk_type(). |
character(len=lcar_type) user_routines_galaxy::thin_disk_type |
thin_disk_type | Thin dust disk profile type. Same choices as for old_disk_type(). |
integer user_routines_galaxy::thin_disk_type_id |
thin_disk_type_ID | ID number for thin_disk_type(). |
character(len=lcar_type) user_routines_galaxy::young_disk_type |
young_disk_type | Young stellar disk profile type. Same choices as for old_disk_type(). |
integer user_routines_galaxy::young_disk_type_id |
young_disk_type_ID | ID number for young_disk_type(). |
real(kind=real64) user_routines_galaxy::z_subd_lim |
z_subd_lim | Z coordinate below which cells have to be subdivided if their R is less than R_subd_lim() |
real(kind=real64) user_routines_galaxy::z_subd_lim2 |
z_subd_lim2 | Z coordinate used in the 'milky_way' subdivision criteria to limit the subdivisions of cells with small subtended solid angle from the position (rsun(), 0,0). See subdivision_criteria(). |
real(kind=real64) user_routines_galaxy::zd_disk |
zd_disk | Scale height for the thick dust disk (basic value). |
real(kind=real64) user_routines_galaxy::zd_disk_r1 |
zd_disk_r1 | Scale height for the thick dust disk at the inner radius hdin(). |
real(kind=real64) user_routines_galaxy::zd_disk_rsun |
zd_disk_rsun | Scale height for the thick dust disk at the radius rsun(). |
real(kind=real64) user_routines_galaxy::zd_tdisk |
zd_tdisk | Scale height for the thin dust disk (basic value). |
real(kind=real64) user_routines_galaxy::zd_tdisk_r1 |
zd_tdisk_r1 | Scale height for the thin dust disk at the inner radius hd1in(). |
real(kind=real64) user_routines_galaxy::zd_tdisk_rsun |
zd_tdisk_rsun | Scale height for the thin dust disk at the radius rsun(). |
real(kind=real64) user_routines_galaxy::zs_disk |
zs_disk | Scale height for the old stellar disk (basic value). |
real(kind=real64) user_routines_galaxy::zs_disk_r1 |
zs_disk_r1 | Scale height for the old stellar disk at the inner radius hsin(). |
real(kind=real64) user_routines_galaxy::zs_disk_rsun |
zs_disk_rsun | Scale height for the old stellar disk at the radius rsun(). |
real(kind=real64) user_routines_galaxy::zs_tdisk |
zs_tdisk | Scale height of the young stellar disk (basic value). |
real(kind=real64) user_routines_galaxy::zs_tdisk_r1 |
zs_tdisk_r1 | Scale height of the young stellar disk at the inner radius hs1in(). |
real(kind=real64) user_routines_galaxy::zs_tdisk_rsun |
zs_tdisk_rsun | Scale height of the young stellar disk at the radius rsun(). |