This module reads in Rosseland-mean opacities from OPAL tables. Table opacity values are given in base 10 logarithm and are a function of mass density (R) and temperature (T), which are also both given in base 10 logarithm. More...
Functions/Subroutines | |
| subroutine, public | init_opal_table (tablename) |
| This subroutine is called when the FLD radiation module is initialised. The OPAL tables for different helium abundances are read and interpolated. More... | |
| subroutine, public | set_opal_opacity (rho, temp, kappa) |
| This subroutine calculates the opacity for a given temperature-density structure. Opacities are read from a table with given metalicity. More... | |
| subroutine | read_table (R, T, K, filename) |
| This routine reads in 1-D (rho,T) values from an opacity table and gives back as output the 1-D (rho,T) values and 2-D opacity. More... | |
| subroutine | get_kappa (Kappa_vals, logR_list, logT_list, R, T, K) |
| This subroutine looks in the table for the four couples (T,rho) surrounding a given input T and rho. More... | |
| subroutine | get_low_up_index (var_in, var_list, imin, imax, low_i, up_i) |
| This subroutine finds the indices in rho and T arrays of the two values surrounding the input rho and T. More... | |
| subroutine | interpolate_krt (low_r, up_r, low_t, up_t, logR_list, logT_list, Kappa_vals, R, T, kappa_interp) |
| This subroutine does a bilinear interpolation in the R,T-plane. More... | |
| subroutine | interpolate1d (x1, x2, x, y1, y2, y) |
| Interpolation in one dimension. More... | |
Variables | |
| integer, parameter | irmin = 2 |
| min and max indices for R,T-range in opacity table More... | |
| integer, parameter | irmax = 20 |
| integer, parameter | itmin = 7 |
| integer, parameter | itmax = 76 |
| double precision, dimension(itmin:itmax, irmin:irmax), public | kappa_vals |
| The opacity tables are read once and stored globally in Kappa_vals. More... | |
| double precision, dimension(irmin:irmax), public | logr_list |
| double precision, dimension(itmin:itmax), public | logt_list |
Detailed Description
This module reads in Rosseland-mean opacities from OPAL tables. Table opacity values are given in base 10 logarithm and are a function of mass density (R) and temperature (T), which are also both given in base 10 logarithm.
Function/Subroutine Documentation
◆ get_kappa()
| subroutine mod_opal_opacity::get_kappa | ( | double precision, dimension(itmin:itmax,irmin:irmax), intent(in) | Kappa_vals, |
| double precision, dimension(irmin:irmax), intent(in) | logR_list, | ||
| double precision, dimension(itmin:itmax), intent(in) | logT_list, | ||
| double precision, intent(in) | R, | ||
| double precision, intent(in) | T, | ||
| double precision, intent(out) | K | ||
| ) |
This subroutine looks in the table for the four couples (T,rho) surrounding a given input T and rho.
Definition at line 107 of file mod_opal_opacity.t.
◆ get_low_up_index()
| subroutine mod_opal_opacity::get_low_up_index | ( | double precision, intent(in) | var_in, |
| double precision, dimension(imin:imax), intent(in) | var_list, | ||
| integer, intent(in) | imin, | ||
| integer, intent(in) | imax, | ||
| integer, intent(out) | low_i, | ||
| integer, intent(out) | up_i | ||
| ) |
This subroutine finds the indices in rho and T arrays of the two values surrounding the input rho and T.
Definition at line 148 of file mod_opal_opacity.t.
◆ init_opal_table()
| subroutine, public mod_opal_opacity::init_opal_table | ( | character(len=*), intent(in) | tablename | ) |
This subroutine is called when the FLD radiation module is initialised. The OPAL tables for different helium abundances are read and interpolated.
Definition at line 24 of file mod_opal_opacity.t.
◆ interpolate1d()
| subroutine mod_opal_opacity::interpolate1d | ( | double precision, intent(in) | x1, |
| double precision, intent(in) | x2, | ||
| double precision, intent(in) | x, | ||
| double precision, intent(in) | y1, | ||
| double precision, intent(in) | y2, | ||
| double precision, intent(out) | y | ||
| ) |
Interpolation in one dimension.
Definition at line 217 of file mod_opal_opacity.t.
◆ interpolate_krt()
| subroutine mod_opal_opacity::interpolate_krt | ( | integer, intent(in) | low_r, |
| integer, intent(in) | up_r, | ||
| integer, intent(in) | low_t, | ||
| integer, intent(in) | up_t, | ||
| double precision, dimension(irmin:irmax), intent(in) | logR_list, | ||
| double precision, dimension(itmin:itmax), intent(in) | logT_list, | ||
| double precision, dimension(itmin:itmax,irmin:irmax), intent(in) | Kappa_vals, | ||
| double precision, intent(in) | R, | ||
| double precision, intent(in) | T, | ||
| double precision, intent(out) | kappa_interp | ||
| ) |
This subroutine does a bilinear interpolation in the R,T-plane.
Definition at line 172 of file mod_opal_opacity.t.
◆ read_table()
| subroutine mod_opal_opacity::read_table | ( | double precision, dimension(irmin:irmax), intent(out) | R, |
| double precision, dimension(itmin:itmax), intent(out) | T, | ||
| double precision, dimension(itmin:itmax,irmin:irmax), intent(out) | K, | ||
| character(*), intent(in) | filename | ||
| ) |
This routine reads in 1-D (rho,T) values from an opacity table and gives back as output the 1-D (rho,T) values and 2-D opacity.
Definition at line 76 of file mod_opal_opacity.t.
◆ set_opal_opacity()
| subroutine, public mod_opal_opacity::set_opal_opacity | ( | double precision, intent(in) | rho, |
| double precision, intent(in) | temp, | ||
| double precision, intent(out) | kappa | ||
| ) |
This subroutine calculates the opacity for a given temperature-density structure. Opacities are read from a table with given metalicity.
Definition at line 43 of file mod_opal_opacity.t.
Variable Documentation
◆ irmax
| integer, parameter mod_opal_opacity::irmax = 20 |
Definition at line 10 of file mod_opal_opacity.t.
◆ irmin
| integer, parameter mod_opal_opacity::irmin = 2 |
min and max indices for R,T-range in opacity table
Definition at line 10 of file mod_opal_opacity.t.
◆ itmax
| integer, parameter mod_opal_opacity::itmax = 76 |
Definition at line 11 of file mod_opal_opacity.t.
◆ itmin
| integer, parameter mod_opal_opacity::itmin = 7 |
Definition at line 11 of file mod_opal_opacity.t.
◆ kappa_vals
| double precision, dimension(itmin:itmax,irmin:irmax), public mod_opal_opacity::kappa_vals |
The opacity tables are read once and stored globally in Kappa_vals.
Definition at line 14 of file mod_opal_opacity.t.
◆ logr_list
| double precision, dimension(irmin:irmax), public mod_opal_opacity::logr_list |
Definition at line 15 of file mod_opal_opacity.t.
◆ logt_list
| double precision, dimension(itmin:itmax), public mod_opal_opacity::logt_list |
Definition at line 15 of file mod_opal_opacity.t.
