MPI-AMRVAC  3.1
The MPI - Adaptive Mesh Refinement - Versatile Advection Code (development version)
Modules | Functions/Subroutines | Variables
mod_cak_opacity.t File Reference

Go to the source code of this file.

Modules

module  mod_cak_opacity
 This module reads in CAK line opacities in the Gayley (1995) notation (alpha, Qbar, Q0, kappae) from corresponding tables. Tabulated values assume LTE conditions and are a function of mass density (D) and temperature (T), which are both given in base 10 logarithm. The construction of the tables is outlined in Poniatowski+ (2021), A&A, 667.
 

Functions/Subroutines

subroutine, public mod_cak_opacity::init_cak_table (tabledir)
 This routine is called when the FLD radiation module is initialised. More...
 
subroutine, public mod_cak_opacity::set_cak_opacity (rho, temp, alpha_output, Qbar_output, Q0_output, kappae_output)
 This subroutine calculates the opacity for a given temperature-density structure. Opacities are read from a table with given metalicity. More...
 
subroutine mod_cak_opacity::read_table (D, T, K, filename)
 This routine reads in 1-D (T,rho) values from a line opacity table and gives back as output the 1-D (T,rho) values and 2-D line opacity value. More...
 
subroutine mod_cak_opacity::get_val_comb (K1_vals, K2_vals, K3_vals, K4_vals, logD_list, logT_list, D, T, K1, K2, K3, K4)
 This subroutine looks in the table for the four couples (T,rho) surrounding a given input T and rho. More...
 
subroutine mod_cak_opacity::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 mod_cak_opacity::interpolate_krt (low_r, up_r, low_t, up_t, logD_list, logT_list, Kappa_vals, D, T, kappa_interp)
 This subroutine does a bilinear interpolation in the R,T-plane. More...
 
subroutine mod_cak_opacity::interpolate1d (x1, x2, x, y1, y2, y)
 Interpolation in one dimension. More...
 

Variables

integer, parameter mod_cak_opacity::idmin = 2
 min and max indices for R,T-range in opacity table More...
 
integer, parameter mod_cak_opacity::idmax = 16
 
integer, parameter mod_cak_opacity::itmin = 2
 
integer, parameter mod_cak_opacity::itmax = 51
 
double precision, dimension(idmin:idmax, itmin:itmax), public mod_cak_opacity::alpha_vals
 The opacity tables are read once and stored globally. More...
 
double precision, dimension(idmin:idmax, itmin:itmax), public mod_cak_opacity::qbar_vals
 
double precision, dimension(idmin:idmax, itmin:itmax), public mod_cak_opacity::q0_vals
 
double precision, dimension(idmin:idmax, itmin:itmax), public mod_cak_opacity::kappae_vals
 
double precision, dimension(idmin:idmax), public mod_cak_opacity::logd_list
 
double precision, dimension(itmin:itmax), public mod_cak_opacity::logt_list