(Back to Bearclaw Outline)

PURE FUNCTION He_recomb(nHeII,ne,T)

Defined in

contrib/astro/physics/cool/cool.f90

Inputs

• nHeII (=1.d0) (in cm^-3)
• ne (=1.d0) (in cm^-3)
• T (in K)

Description

• Calculates rates for HeII radiative and dielectric recombination:
  • He_recomb = nHeII * ne * (He_rad_recomb(Tin) + He_dielec_recomb(Tin)), with
    • radiative recombination from Verner & Ferland (1996ApJS..103..467V) and
    • dielectric recombination from Mazzotta et al. (1998A&AS..133..403M)
    • where Tin = max(200.,T) (I'm not sure why)
• Values are nominally in cm^-3 s^-1, however:
  • When used in initneqCool, nHeII & 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

Functions Used

• He_rad_recomb(Tin) in this function (He_recomb)
• He_dielec_recomb(Tin) in this function (He_recomb)

Modules Used

None