(Back to Bearclaw Outline)

PURE FUNCTION dust_cool(n,T)

Defined in

contrib/astro/physics/cool/cool.f90

Inputs

• n (n = 1.d0 apparently, from calling arguments)
• T

Description

• Cooling due to H atom sticking to dust (ergs/cm³/s)
• Used to build tabulated form of cooling, dust_cool_TAB
• dust_cool = 0. for T < 10.
• Otherwise, dust_cool is calculated using the expression from Smith & Rosen (2003MNRAS.339..133S) (Appendix, A3):
  • n2 * 3.8d-33 * sqrt(T) * (T-Tdust) * (1 - 0.8*exp(-75/T))
• Uses Tdust = 20., defined in dev/contrib/astro/sample_data_files/neqcool.data

Called in

• contrib/astro/physics/cool/cool.f90 in SUBROUTINE initneqCool
• However, in i_evolove.f90, where the actual cooling is calculated from dust_cool_table, this is commented out, so this form of dust cooling is not currently being used.

Functions Used

None