temperature_TT

Module

Description

$Id$

This module can replace the entropy module by using _T_ as dependent variable. For a perfect gas with constant coefficients (no ionization) we have (1-1/gamma) * cp*T = cs02 * exp( (gamma-1)*ln(rho/rho0)-gamma*s/cp )

At a later point we may want to rename the module Entropy into Energy

Quick access

Variables:

adi_kprof, ampl_ss, calc_heatcond_simple, calc_tau_ss_exterior, chi_t, cooltype, cs2cool, dss_dt, epsilon_ss, fbotkbot, ftop, ftopktop, get_slices_entropy, gradloghcond, heatcond, idiag_ugradpm, iglobal_glhc, iglobal_hcond, ilntt, init_ss, initialize_entropy, kbot, khor_ss, lcalc_heatcond_constchi, lcalc_heatcond_simple, radius_ss, read_entropy_init_pars, read_entropy_run_pars, register_entropy, rprint_entropy, ss0, ss_left, ss_right, tau_ss_exterior, thermal_background, thermal_peak, thermal_scaling, widthss, write_entropy_init_pars, write_entropy_run_pars

Routines:

adi_constk(), bc_ct(), bc_flux(), calc_heatcond_arctan(), calc_pencils_entropy(), cyclic(), hcond_adi(), pencil_criteria_entropy(), pencil_interdep_entropy(), rad_equil(), tridag()

Needed modules

Variables

  • entropy/adi_kprof [public]
  • entropy/ampl_lntt [real,optional/default=0.0]
  • entropy/center1_x [real,optional/default=0.0]
  • entropy/center1_y [real,optional/default=0.0]
  • entropy/center1_z [real,optional/default=0.0]
  • entropy/chi [real,optional/default=impossible]

    ‘)

  • entropy/fbot [real]
  • entropy/fbotkbot [real]
  • entropy/ftop [real]
  • entropy/ftopktop [real]
  • entropy/hcond0 [real,optional/default=impossible]

    ‘)

  • entropy/hcond1 [real,optional/default=1.0]
  • entropy/heat_uniform [real,optional/default=0.0]
  • entropy/idiag_csm [integer,optional/default=0]
  • entropy/idiag_dtc [integer,optional/default=0]
  • entropy/idiag_dtchi [integer,optional/default=0]
  • entropy/idiag_eem [integer,optional/default=0]
  • entropy/idiag_ethm [integer,optional/default=0]
  • entropy/idiag_ppm [integer,optional/default=0]
  • entropy/idiag_ssm [integer,optional/default=0]
  • entropy/idiag_ttm [integer,optional/default=0]
  • entropy/idiag_ttmax [integer,optional/default=0]
  • entropy/idiag_ttmin [integer,optional/default=0]
  • entropy/iglobal_glhc [integer,optional/default=0]
  • entropy/iglobal_hcond [integer,optional/default=0]
  • entropy/iheatcond [character,optional/default='nothing']
  • entropy/iinit_str [character]
  • entropy/initlntt (ninit) [character,optional/default='nothing']
  • entropy/kbot [real]
  • entropy/kmax [real,optional/default=0.0]
  • entropy/kmin [real,optional/default=0.0]
  • entropy/kx_lntt [real,optional/default=1.0]
  • entropy/ky_lntt [real,optional/default=1.0]
  • entropy/kz_lntt [real,optional/default=1.0]
  • entropy/ladvection_temperature [logical,optional/default=.true.]
  • entropy/lhcond_global [logical,optional/default=.false.]
  • entropy/lheatc_chiconst [logical,optional/default=.false.]
  • entropy/lheatc_chiconst_accurate [logical,optional/default=.false.]
  • entropy/lheatc_karctan [logical,optional/default=.false.]
  • entropy/lheatc_kconst [logical,optional/default=.false.]
  • entropy/lheatc_kprof [logical,optional/default=.false.]
  • entropy/lmultilayer [logical,optional/default=.false.]
  • entropy/lntt_const [real,optional/default=0.0]
  • entropy/lntt_left [real,optional/default=1.0]
  • entropy/lntt_right [real,optional/default=1.0]
  • entropy/lpressuregradient_gas [logical,optional/default=.true.]
  • entropy/lupw_lntt [logical,optional/default=.false.]
  • entropy/lviscosity_heat [logical,optional/default=.true.]
  • entropy/nheatc_max [integer,parameter=1]
  • entropy/r_bcz [real,optional/default=0.0]
  • entropy/radius_lntt [real,optional/default=0.1]
  • entropy/tau_damp_cor [real,optional/default=-1.0]
  • entropy/tau_heat_cor [real,optional/default=-1.0]
  • entropy/tbump [real,optional/default=0.0]
  • entropy/tt_const [real,optional/default=1.0]
  • entropy/tt_cor [real,optional/default=0.0]
  • entropy/widthlntt [real,optional/default=2*epsi]
  • entropy/zbot [real,optional/default=0.0]
  • entropy/zcor [real,optional/default=0.0]
  • entropy/ztop [real,optional/default=0.0]

Subroutines and functions

subroutine  entropy/pencil_criteria_entropy()
Use :

cdata

Call to:

identify_bcs(), calc_viscous_heat(), calc_heatcond_arctan(), max_mn_name(), der_step(), cyclic(), bc_ct(), tridag(), hcond_adi(), bc_flux()

subroutine  entropy/pencil_interdep_entropy(lpencil_in)

Interdependency among pencils from the Entropy module is specified here.

20-11-04/anders: coded

Parameters:

lpencil_in (npencils) [logical]

Call to:

identify_bcs(), calc_viscous_heat(), calc_heatcond_arctan(), max_mn_name(), der_step(), cyclic(), bc_ct(), tridag(), hcond_adi(), bc_flux()

subroutine  entropy/calc_pencils_entropy(f, p)

Calculate Entropy pencils. Most basic pencils should come first, as others may depend on them.

20-11-04/anders: coded

Parameters:
Use :

equationofstate, sub

Call to:

identify_bcs(), calc_viscous_heat(), calc_heatcond_arctan(), max_mn_name(), der_step(), cyclic(), bc_ct(), tridag(), hcond_adi(), bc_flux()

subroutine  entropy/rad_equil(f)

16-mai-2007/tgastine+dintrans: compute the radiative and hydrostatic equilibria for a given radiative profile (here a hole for the moment)

Parameters:

f (mx,my,mz,mfarray) [real,inout]

Use :

cdata, gravity (gravz()), equationofstate (cs20(), lnrho0(), cs2top(), cs2bot(), gamma(), gamma_m1())

Call to:

max_mn_name(), der_step(), cyclic(), bc_ct(), tridag(), hcond_adi(), bc_flux()

subroutine  entropy/calc_heatcond_arctan(df, p)
16-mai-2007/tgastine+dintrans: radiative diffusion with an arctan

profile for the conductivity calculate gamma/rho*cp*div(K T*grad lnTT)=

gamma*K/rho*cp*(gradlnTT.gradlnTT + del2ln TT + gradlnTT.gradlnK)

Parameters:
  • df (mx,my,mz,mvar) [real]

  • p [pencil_case]

Use :

diagnostics, equationofstate (gamma()), hdf5io (output_profile()), sub

Called from:

pencil_criteria_entropy(), pencil_interdep_entropy(), calc_pencils_entropy()

Call to:

max_mn_name(), der_step(), cyclic(), bc_ct(), tridag(), hcond_adi(), bc_flux()

subroutine  entropy/adi_constk(finit, f)
Parameters:
Use :

cdata, cparam, equationofstate (gamma(), gamma_m1(), cs2bot(), cs2top())

Call to:

cyclic(), bc_ct(), tridag(), hcond_adi(), bc_flux()

subroutine  entropy/bc_ct(f_2d)
Parameters:

f_2d (mx,mz) [real]

Called from:

pencil_criteria_entropy(), pencil_interdep_entropy(), calc_pencils_entropy(), rad_equil(), calc_heatcond_arctan(), adi_constk()

Call to:

tridag()

subroutine  entropy/bc_flux(f, chiprof)
Parameters:
Called from:

pencil_criteria_entropy(), pencil_interdep_entropy(), calc_pencils_entropy(), rad_equil(), calc_heatcond_arctan(), adi_constk()

Call to:

tridag()

subroutine  entropy/hcond_adi(f, chiprof, dchi)
Parameters:
Use :

sub (write_zprof())

Called from:

pencil_criteria_entropy(), pencil_interdep_entropy(), calc_pencils_entropy(), rad_equil(), calc_heatcond_arctan(), adi_constk()

Call to:

tridag()

subroutine  entropy/tridag(a, b, c, r, u, n)
Parameters:
  • a (n) [real]

  • b (n) [real]

  • c (n) [real]

  • r (n) [real]

  • u (n) [real]

  • n [integer]

Called from:

register_magnetic(), initialize_magnetic(), init_aa(), pencil_criteria_magnetic(), pencil_interdep_magnetic(), magnetic_after_boundary(), daa_dt(), time_integrals_magnetic(), df_diagnos_magnetic(), rescaling_magnetic(), register_density(), initialize_density(), init_lnrho(), density_before_boundary(), pencil_criteria_density(), pencil_interdep_density(), density_after_boundary(), dlnrho_dt(), impose_density_floor(), read_density_init_pars(), write_density_init_pars(), read_density_run_pars(), write_density_run_pars(), rprint_density(), get_slices_density(), get_slices_pressure(), get_init_average_density(), calc_pencils_density(), in_solid_cell(), wsnap_ogrid(), map_nearest_grid_ogrid(), interpolate_particles_ogrid(), anelastic_after_mn(), calc_pencils_magnetic(), sine_avoid_boundary(), pencil_criteria_entropy(), pencil_interdep_entropy(), calc_pencils_entropy(), rad_equil(), calc_heatcond_arctan(), adi_constk(), bc_ct(), bc_flux(), hcond_adi(), tridag(), cyclic()

Call to:

tridag()

subroutine  entropy/cyclic(a, b, c, alpha, beta, r, x, n)
Parameters:
  • a (n) [real]

  • b (n) [real]

  • c (n) [real]

  • alpha [real]

  • beta [real]

  • r (n) [real]

  • x (n) [real]

  • n [integer]

Called from:

in_solid_cell(), wsnap_ogrid(), map_nearest_grid_ogrid(), interpolate_particles_ogrid(), pencil_criteria_entropy(), pencil_interdep_entropy(), calc_pencils_entropy(), rad_equil(), calc_heatcond_arctan(), adi_constk()

Call to:

tridag()