noneutraldensity

Module

Description

$Id$

This module is used both for the initial condition and during run time. It contains dlnrhon_dt and init_lnrhon, among other auxiliary routines.

** AUTOMATIC CPARAM.INC GENERATION ************************ Declare (for generation of cparam.inc) the number of f array variables and auxiliary variables added by this module

CPARAM logical, parameter :: lneutraldensity = .false.

MVAR CONTRIBUTION 0 MAUX CONTRIBUTION 0

Quick access

Variables:

alpha_prescription, alpha_time, ampllnrhon, borderlnrhon, diffrhon, diffrhon_hyper3, diffrhon_hyper3_aniso, diffrhon_shock, diffus_diffrhon, enum_alpha_prescription, enum_borderlnrhon, idiag_alprec, idiag_dtnd, idiag_lnrhon2m, idiag_lnrhonmphi, idiag_neutralmass, idiag_rhon2m, idiag_rhonm, idiag_rhonmax, idiag_rhonmin, idiag_rhonmphi, idiag_rhonmr, idiag_rhonmx, idiag_rhonmxy, idiag_rhonmy, idiag_rhonmz, idiag_unglnrhonm, idiffn, initlnrhon, kx_lnrhon, ky_lnrhon, kz_lnrhon, lcontinuity_neutral, ldiffn_hyper3, ldiffn_hyper3_aniso, ldiffn_hyper3_polar, ldiffn_hyper3lnrhon, ldiffn_normal, ldiffn_shock, lfreeze_lnrhonext, lfreeze_lnrhonint, lnrhon0, lnrhon_const, lnrhon_ext, lnrhon_int, lnrhon_left, lnrhon_right, lpretend_star, lramp_up, lupw_lnrhon, lupw_rhon, luse_as_ionization, rhon_const, rhon_left, rhon_right, set_border_neutraldensity, star_form_exponent, star_form_threshold

Routines:

calc_diagnostics_neutraldens(), calc_pencils_neutraldensity(), dlnrhon_dt(), init_lnrhon(), initialize_neutraldensity(), neutraldensity_after_boundary(), pencil_criteria_neutraldensity(), pencil_interdep_neutraldensity(), read_neutraldensity_init_pars(), read_neutraldensity_run_pars(), register_neutraldensity(), rprint_neutraldensity(), write_neutraldensity_init_pars(), write_neutraldensity_run_pars()

Needed modules

Variables

Subroutines and functions

subroutine  neutraldensity/register_neutraldensity()

Initialise variables which should know that we solve the compressible hydro equations: ilnrhon; increase nvar accordingly.

18-mar-03/axel: adapted from neutraldensity

Identify version number (generated automatically by SVN).

Called from:

register_modules()

Call to:

svn_id()

subroutine  neutraldensity/initialize_neutraldensity()

Perform any post-parameter-read initialization i.e. calculate derived parameters.

18-mar-03/axel: adapted from neutraldensity

Called from:

register_modules(), initialize_modules()

subroutine  neutraldensity/init_lnrhon(f)

Initialise lnrhon; called from start.f90.

18-mar-03/axel: adapted from neutraldensity

Parameters:

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

Called from:

start

subroutine  neutraldensity/pencil_criteria_neutraldensity()

All pencils that the Neutraldensity module depends on are specified here.

20-11-04/anders: coded

Called from:

register_modules(), initialize_modules(), finalize_modules(), choose_pencils()

subroutine  neutraldensity/pencil_interdep_neutraldensity(lpencil_in)

Interdependency among pencils provided by the Neutraldensity module is specified here.

20-11-04/anders: coded

Parameters:

lpencil_in (npencils) [logical]

Called from:

register_modules(), initialize_modules(), finalize_modules(), choose_pencils()

subroutine  neutraldensity/calc_pencils_neutraldensity(f, p)

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

13-nov-04/anders: coded

Parameters:

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

  • p [pencil_case,in]

Called from:

pde(), write_diagnostics(), calc_all_module_diagnostic_auxiliaries(), perform_diagnostics(), finalize_diagnostics()

subroutine  neutraldensity/calc_diagnostics_neutraldens(p)
Parameters:

p [pencil_case]

Called from:

pde(), write_diagnostics(), calc_all_module_diagnostic_auxiliaries()

subroutine  neutraldensity/neutraldensity_after_boundary(f)
Parameters:

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

Called from:

pde(), write_diagnostics(), calc_all_module_diagnostic_auxiliaries(), perform_diagnostics(), finalize_diagnostics()

subroutine  neutraldensity/dlnrhon_dt(f, df, p)

Continuity equation. Calculate dlnrhon/dt = - u.gradlnrhon - divud

18-mar-03/axel: adapted from neutraldensity

Parameters:

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

  • df (,,*,*) [real,contiguous]

  • p [pencil_case]

Called from:

pde(), write_diagnostics(), calc_all_module_diagnostic_auxiliaries(), perform_diagnostics(), finalize_diagnostics()

subroutine  neutraldensity/read_neutraldensity_init_pars(iomsg)
Parameters:

iomsg [character,out]

subroutine  neutraldensity/write_neutraldensity_init_pars(unit)
Parameters:

unit [integer,in]

Called from:

get_datadir(), get_snapdir(), read_all_init_pars(), read_all_run_pars(), get_downpars()

subroutine  neutraldensity/read_neutraldensity_run_pars(iomsg)
Parameters:

iomsg [character,out]

subroutine  neutraldensity/write_neutraldensity_run_pars(unit)
Parameters:

unit [integer,in]

Called from:

get_datadir(), get_snapdir(), read_all_init_pars(), read_all_run_pars(), get_downpars()

subroutine  neutraldensity/rprint_neutraldensity(lreset[, lwrite])

Reads and registers print parameters relevant for neutral density.

7-jun-02/axel: incoporated from subroutine pde

Parameters:

  • lreset [logical]

  • lwrite [logical]

Called from:

register_modules(), initialize_modules(), finalize_modules(), choose_pencils(), rprint_list()