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

Go to the source code of this file.

Modules

module  mod_cak_force
 Module to include CAK radiation line force in (magneto)hydrodynamic models Computes both the force from free electrons and the force from an ensemble of lines (various possibilities for the latter). There is an option to only simulate the pure radial CAK force (with various corrections applied) as well as the full vector CAK force. Depending on the chosen option additional output are the CAK line force component(s) and, when doing a 1-D radial force, the finite disc factor.
 

Functions/Subroutines

subroutine mod_cak_force::cak_params_read (files)
 Public method. More...
 
subroutine mod_cak_force::cak_init (phys_gamma)
 Initialize the module. More...
 
subroutine, public mod_cak_force::set_cak_force_norm (rstar, twind)
 Compute some (unitless) variables for CAK force normalisation. More...
 
subroutine mod_cak_force::cak_add_source (qdt, ixIL, ixOL, wCT, w, x, energy, qsourcesplit, active)
 w[iw]=w[iw]+qdt*S[wCT,qtC,x] where S is the source based on wCT within ixO More...
 
subroutine mod_cak_force::get_cak_force_radial (ixIL, ixOL, wCT, w, x, gcak)
 1-D CAK line force in the Gayley line-ensemble distribution parametrisation More...
 
subroutine mod_cak_force::get_cak_force_vector (ixIL, ixOL, wCT, w, x, gcak)
 Vector CAK line force in the Gayley line-ensemble distribution parametrisation. More...
 
subroutine mod_cak_force::get_gelectron (ixIL, ixOL, w, x, ge)
 Compute continuum radiation force from Thomson scattering. More...
 
subroutine mod_cak_force::cak_get_dt (w, ixIL, ixOL, dtnew, dxD, x)
 Check time step for total radiation contribution. More...
 
subroutine mod_cak_force::get_velocity_gradient (ixIL, ixOL, vfield, x, idir, grad_vn)
 Compute velocity gradient in direction 'idir' on a non-uniform grid. More...
 
subroutine mod_cak_force::rays_init (ntheta_point, nphi_point)
 Initialise (theta',phi') radiation angles coming from stellar disc. More...
 
subroutine mod_cak_force::gauss_legendre_quadrature (xlow, xhi, n, x, w)
 Fast Gauss-Legendre N-point quadrature algorithm by G. Rybicki. More...
 

Variables

real(8), public mod_cak_force::cak_alpha
 Line-ensemble parameters in the Gayley (1995) formalism. More...
 
real(8), public mod_cak_force::gayley_qbar
 
real(8), public mod_cak_force::gayley_q0
 
integer mod_cak_force::cak_1d_opt
 Switch to choose between the 1-D CAK line force options. More...
 
logical mod_cak_force::cak_split =.false.
 To treat source term in split or unsplit (default) fashion. More...
 
logical mod_cak_force::cak_1d_force =.false.
 To activate the original CAK 1-D line force computation. More...
 
logical mod_cak_force::cak_vector_force =.false.
 To activate the vector CAK line force computation. More...
 
logical mod_cak_force::fix_vector_force_1d =.false.
 To activate the pure radial vector CAK line force computation. More...
 
integer mod_cak_force::nthetaray
 Amount of rays in radiation polar and radiation azimuthal direction. More...
 
integer mod_cak_force::nphiray
 
integer mod_cak_force::gcak1_
 Extra slots to store quantities in w-array. More...
 
integer mod_cak_force::gcak2_
 
integer mod_cak_force::gcak3_
 
integer mod_cak_force::fdf_