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_