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