temperature_bcs

Module

Quick access

Variables:

bc_stellar_surface, cp, density_scale, density_scale1, hcondxbot, hcondxtop, hcondzbot, hcondztop, ilnrho_cs2, ilnrho_ee, ilnrho_eth, ilnrho_lntt, ilnrho_pp, ilnrho_ss, ilnrho_tt, ipp_cs2, ipp_ss, irho_cs2, irho_ee, irho_eth, irho_lntt, irho_pp, irho_ss, irho_tt, reference_state, ttbot, tttop

Routines:

bc_ism_energ(), bc_lnrho_hds_z_iso_energ(), bdry_magnetic(), initialize_energy_bcs()

Needed modules

Variables

  • energybcs/cp [real,private]
  • energybcs/density_scale [real,private]
  • energybcs/density_scale1 [real,private]
  • energybcs/gamma [real,private]
  • energybcs/gamma_m1 [real,private]
  • energybcs/ilnrho_cs2 [integer,private/parameter/optional/default=5]
  • energybcs/ilnrho_ee [integer,private/parameter/optional/default=2]
  • energybcs/ilnrho_eth [integer,private/parameter/optional/default=14]
  • energybcs/ilnrho_lntt [integer,private/parameter/optional/default=4]
  • energybcs/ilnrho_pp [integer,private/parameter/optional/default=3]
  • energybcs/ilnrho_ss [integer,private/parameter/optional/default=1]
  • energybcs/ilnrho_tt [integer,private/parameter/optional/default=9]
  • energybcs/ipp_cs2 [integer,private/parameter/optional/default=12]
  • energybcs/ipp_ss [integer,private/parameter/optional/default=11]
  • energybcs/irho_cs2 [integer,private/parameter/optional/default=6]
  • energybcs/irho_ee [integer,private/parameter/optional/default=15]
  • energybcs/irho_eth [integer,private/parameter/optional/default=13]
  • energybcs/irho_lntt [integer,private/parameter/optional/default=8]
  • energybcs/irho_pp [integer,private/parameter/optional/default=16]
  • energybcs/irho_ss [integer,private/parameter/optional/default=7]
  • energybcs/irho_tt [integer,private/parameter/optional/default=10]
  • energybcs/ttbot [real,private/pointer]
  • energybcs/tttop [real,private/pointer]

Subroutines and functions

subroutine  energybcs/initialize_energy_bcs()
Use :

equationofstate (get_gamma_etc()), sharedvariables (get_shared_variable())

Call to:

get_gamma_etc(), fatal_error(), set_ghosts_for_onesided_ders()

subroutine  energybcs/bc_lnrho_hds_z_iso_energ(f, topbot)

Boundary condition for density and entropy.

This sets

partial_{z} lnrho

such that

partial_{z} p = rho g_{z},

i.e. it enforces hydrostatic equlibrium at the boundary.

Currently this is only correct if

partial_{z} lnT = 0

at the boundary.

12-Juil-2006/dintrans: coded

Parameters:

  • f (,,*,*) [real,inout] :: factor 1/cp corrected

  • topbot [integer,in]

Use :

equationofstate (eoscalc()), gravity (gravz())

Call to:

fatal_error(), set_ghosts_for_onesided_ders()

subroutine  energybcs/bc_ism_energ(f, topbot, j)

30-nov-15/fred: Replaced bc_ctz and bc_cdz. Apply observed scale height locally from Reynolds 1991, Manchester & Taylor 1981 for warm ionized gas - dominant scale height above 500 parsecs. Apply constant local temperature across boundary for entropy. Motivation to prevent numerical spikes in shock fronts, which cannot be absorbed in only three ghost cells, but boundary thermodynamics still responsive to interior dynamics.

Parameters:

  • f (,,*,*) [real]

  • topbot [integer,in]

  • j [integer]

Call to:

fatal_error()