Sources
Includes all source terms. More...
Modules | |
module | sources_base_mod |
generic source terms module providing functionaly common to all source terms | |
module | sources_c_accel_mod |
source terms module for constant acceleration | |
module | sources_diskcooling_mod |
source terms module for cooling of geometrically thin accretion disks | |
module | sources_gravity_mod |
generic gravity terms module providing functionaly common to all gravity terms | |
module | sources_shearbox_mod |
Source terms module for fictious forces in a shearingsheet. | |
module | sources_viscosity_mod |
computes momentum and energy sources due to shear stresses | |
Detailed Description
Includes all source terms.
source term module for shearing forces in a shearingsheet
- general parameters of sources group as key-values
- INTEGER stype of source
- parameters of sources_diskcooling as key-values
- INTEGER method model (GRAY or GAMMIE)
- REAL cvis = 0.1 factor for numerical stability
- REAL Tmin = 1.0E-30 a minimum temperature
- REAL rhomin = 1.0E-30 a minimum density
- REAL b_cool = 1.0 parameter
- REAL b_start = 1.0 cooling parameter
- REAL b_final = 1.0 cooling parameter
- REAL t_start = 0.0 to start decreasing value of beta from b_start
- REAL dt_bdec = 0.0 over which b_cool should reach its final value
- beginning at t_start
- INTEGER output/Qcool = 0(=1) output of cooling function
- general parameters of gravity group as key-values
- INTEGER energy source terms to energy equation?
- parameters of sources_rotframe as key-values
- REAL gparam parameter (needed for bianglspherical
- geometry - corresponds to the radius)
- parameters of sources_shearbox as key-values
- Todo:
- This module can only be used with fixed 2D plane in x-y direction. It should be no problem to make it also run in 3D.
- parameters of sources_viscosity as key-values
- INTEGER vismodel model
- REAL dynconst = 0.1 viscosity constant
- REAL bulkconst = -2/3*dynconst viscosity constant
- REAL exponent = 0.0 law exponent
- REAL cvis = 0.5 courant number
- INTEGER output/stress = 0(=1) output of the stress tensor
- INTEGER output/dynvis = 0(=1) output of dynamic viscosity
- INTEGER output/kinvis = 0(=1) output of kinematic viscosity
- INTEGER output/bulkvis = 0(=1) output of bulk viscosity