(Back to Bearclaw Outline)

PURE FUNCTION H_recomb(nHII,ne,T)

Defined in

contrib/astro/physics/cool/cool.f90

Inputs

• nHII (=1.d0)
• ne (=1.d0)
• T (in K)

Description

• Calculates rates for HII + e^- recombination using the expression of Verner & Ferland (1996ApJS..103..467V):
  • H_recomb = nHII*ne* (a / [sqrt(T/T0)*(1+sqrt(T/T0))^(1-b)*(1+sqrt(Tin/T1))^(1+b)])
    • where a,b,T0,T1 are constants, and Tin = max(200.,T) (I'm not sure why)
• Values are nominally in cm^-3 s^-1, however:
  • When used in initneqCool, nHII & ne are inputted as 1.d0.
  • The actual densities aren't multiplied in until the tables are used in i_evolove.f90 to calculate the actual cooling rates.
• Is used in Subroutine initneqCool (in cool.f90) to build table of cooling values for this process

Called in

contrib/astro/physics/cool/cool.f90, in SUBROUTINE initneqCool

Modules Used

None

Functions Used

None