h2chem

Module

Description

Contains routines for Hydrogen Chemistry

Routines are originally by Clare Dobbs Modified extensively by Alex Pettitt Translated to Fortran 90 and adapted for use in Phantom by Daniel Price (2011)

References:

Nelson & Langer (1997) Glover et al. (2010) Bergin et al. (2004)

Owner:

Lionel Siess

Runtime parameters:

None

Dependencies:

cooling_ism, part, physcon, units

Quick access

Variables:

evolve_abundances, get_extra_abundances, h2fd_rate, hchem, k0, k1, phrates, rate_diss0, rconst

Routines:

get_dphot(), init_chem(), update_abundances()

Variables

  • chem/evolve_abundances [private]
  • chem/get_extra_abundances [private]
  • chem/h2fd_rate [private]
  • chem/hchem [private]
  • chem/k0 [real,private/parameter/optional/default=5.0d-16]
  • chem/k1 [real,private/parameter/optional/default=5.0d-10]
  • chem/phrates [real,private/parameter/optional/default=1.56d-10]
  • chem/rate_diss0 [real,private/parameter/optional/default=4.17d-11]
  • chem/rconst [real,private/parameter/optional/default=3.d-1*2.2d-18]

Subroutines and functions

subroutine  chem/init_chem()
Called from:

test_cooling_rate()

function  chem/get_dphot(dphotflag, dphot0, xi, yi, zi)
Parameters:

  • dphotflag [integer,in]

  • dphot0 [real,in]

  • xi [real,in]

  • yi [real,in]

  • zi [real,in]

Return:

get_dphot [real]

Use :

units (udist(), umass()), physcon (solarm(), kpc(), pi())

Called from:

test_cooling_rate()

subroutine  chem/update_abundances(ui, rhoi, chemarrays, nchem, dphot, dt, abund, nabn, gmwvar, abundc, abunde, abundo, abundsi)
Parameters:

  • ui [real,in]

  • rhoi [real,in]

  • chemarrays (nchem) [real,inout]

  • nchem [integer,in,]

  • dphot [real,in]

  • dt [real,in]

  • abund (nabn) [real,out]

  • nabn [integer,in]

  • gmwvar [real,out]

  • abundc [real,in]

  • abunde [real,in]

  • abundo [real,in]

  • abundsi [real,in]

Called from:

test_cooling_rate()