|
| subroutine | mod_fld::fld_params_read (files) |
| | public methods these are called in mod_rhd_phys
|
| |
| subroutine, public | mod_fld::fld_init (he_abundance, radiation_diffusion, energy_interact, r_gamma) |
| | Initialising FLD-module: Read opacities Initialise Multigrid adimensionalise kappa Add extra variables to w-array, flux, kappa, eddington Tensor Lambda and R ...
|
| |
| subroutine, public | mod_fld::get_fld_rad_force (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 This subroutine handles the radiation force
|
| |
| subroutine, public | mod_fld::fld_radforce_get_dt (w, ixil, ixol, dtnew, dxd, x) |
| |
| subroutine, public | mod_fld::fld_get_opacity (w, x, ixil, ixol, fld_kappa) |
| | Sets the opacity in the w-array by calling mod_opal_opacity.
|
| |
| subroutine, public | mod_fld::fld_get_fluxlimiter (w, x, ixil, ixol, fld_lambda, fld_r, nth) |
| | Calculate fld flux limiter This subroutine calculates flux limiter lambda using the prescription stored in fld_fluxlimiter. It also calculates the ratio of radiation scaleheight and mean free path.
|
| |
| subroutine, public | mod_fld::fld_get_radflux (w, x, ixil, ixol, rad_flux, nth) |
| | Calculate Radiation Flux stores radiation flux in w-array.
|
| |
| subroutine | mod_fld::fld_get_eddington (w, x, ixil, ixol, eddington_tensor, nth) |
| | Calculate Eddington-tensor Stores Eddington-tensor in w-array.
|
| |
| subroutine, public | mod_fld::fld_get_radpress (w, x, ixil, ixol, rad_pressure, nth) |
| | Calculate Radiation Pressure Returns Radiation Pressure as tensor.
|
| |
| subroutine | mod_fld::fld_implicit_update (dtfactor, qdt, qtc, psa, psb) |
| |
| subroutine | mod_fld::energy_interaction (w, x, ixil, ixol, dtfactor, qdt) |
| | Energy interaction and photon tiring.
|
| |
| subroutine | mod_fld::diffuse_e_rad_mg (dtfactor, qdt, qtc, psa, psb) |
| | Calling all subroutines to perform the multigrid method Communicates rad_e and diff_coeff to multigrid library.
|
| |
| subroutine | mod_fld::evaluate_e_rad_mg (qtc, psa) |
| | inplace update of psa==>F_im(psa)
|
| |
| subroutine | mod_fld::put_diffterm_onegrid (ixil, ixol, w) |
| | inplace update of psa==>F_im(psa)
|
| |
| subroutine | mod_fld::fld_get_diffcoef_central (w, wct, x, ixil, ixol) |
| | Calculates cell-centered diffusion coefficient to be used in multigrid.
|
| |
| subroutine | mod_fld::update_diffcoeff (psa) |
| |
| subroutine | mod_fld::fld_smooth_diffcoef (w, ixil, ixol) |
| | Use running average on Diffusion coefficient.
|
| |
| subroutine | mod_fld::bisection_method (e_gas, e_rad, c0, c1) |
| | Find the root of the 4th degree polynomial using the bisection method.
|
| |
| subroutine | mod_fld::newton_method (e_gas, e_rad, c0, c1) |
| | Find the root of the 4th degree polynomial using the Newton-Ralphson method.
|
| |
| subroutine | mod_fld::halley_method (e_gas, e_rad, c0, c1) |
| | Find the root of the 4th degree polynomial using the Halley method.
|
| |
| double precision function | mod_fld::polynomial_bisection (e_gas, c0, c1) |
| | Evaluate polynomial at argument e_gas.
|
| |
| double precision function | mod_fld::dpolynomial_bisection (e_gas, c0, c1) |
| | Evaluate first derivative of polynomial at argument e_gas.
|
| |
| double precision function | mod_fld::ddpolynomial_bisection (e_gas, c0, c1) |
| | Evaluate second derivative of polynomial at argument e_gas.
|
| |
