(Back to Bearclaw Outline)

PURE FUNCTION H2_diss(nH2, nH, nHe, ne, T)

Defined in

contrib/astro/physics/cool/cool.f90

Inputs

• nH2
• nH
• nHe
• ne
• T

Description

• Computes H2 dissociation due to collisions with H, He, H2, and e
• Values are in cm^-3 s^-1

Specifics

• Sets H2_diss to zero
  • If T < 50 (or less than 1000, since that can cause a divide by zero in function kd), H2_diss = 0.
  • Also, if nH2 < 10^-6 cm^-3, set H2_diss = 0.
• Limit T ⇐ 10⁵ (although comment says limit to 10k, = 10⁴?), to avoid divergence issues
• Add in H2_diss due to collisions with the different species
• Calculation is done by PURE FUNCTION kd

Called in

contrib/astro/physics/cool/i_evolove.f90

Functions Used

• PURE FUNCTION kd
• PURE FUNCTION kdH_LS
• PURE FUNCTION kdL_MK
• PURE FUNCTION kdL_DM
• PURE FUNCTION kdL_LS (obsolete)