(Back to Bearclaw Outline)

PURE FUNCTION H2_dust_Recomb(nHneuclei,nH,T,Tdust)

Defined in

contrib/astro/physics/cool/cool.f90

Inputs

• nHneuclei (= 1.d0)
• nH (= 1.d0)
• T
• Tdust (defined in neqcool.data)

Description

• Calculates the cooling from H2 recombination due to condensation on dust grains
• Values are nominally in cm^-3 s^-1, however:
  • When used in initneqCool, nHneuclei & nH are inputted as 1.d0.
  • The actual densities aren't multiplied in until the tables are used in the function H2_dust_Recomb_table (called 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
• Only valid for temperatures up to 10⁵K; this cooling = 0 above that.
• See Hollenbach & McKee (1979ApJS…41..555H):
  • = nH * nHneuclei * (3.d-17 * sqrt(T2) * fa / (1 + 0.4*sqrt(T2+Tdust2) + 0.2*T2 + 0.08*T22))
  • where T2 = T/100; Tdust2=Tdust/100; fa = 0.5

Called in

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

Modules Used

None

Functions Used

None