cool.f90

Description

Creates cooling tables for the different cooling factiors (H recombination, etc.)
Includes cooling due to (and makes tables for):

• H ionization and recombination for T < 10⁵
• H2 cooling for T < 10⁶
• H2 dissociation recombination for T < 10⁵
• OI cooling (using H2 as a tracer for OI) for T < 1000
• dust cooling (not used)
• HeII ion recombination for T < 10⁵

USEd by

• i_setprob.f90
  • contrib/astro/i_setprob.f90.
• i_src.f90
  • contrib/astro/physics/i_src.f90 (ONLY EOS_vars, lProtectNeq).
• i_protect.f90
  • contrib/astro/physics/i_protect.f90 (ONLY EOS_vars, lProtectNeq).
• i_i_ODEs.f90
  • contrib/astro/physics/i_i_ODEs.f90.
• contrib/astro/physics/i_i_ExactRS.f90.
  

(Also USEd by CalcComp (ONLY EOS_vars) in set-up files)

Location

contrib/astro/physics/cool/cool.f90

MODULE cool

USES

GlobalDeclarations, ONLY:

• cellinfo, srcPrecision,
• muH, muH2, muHe, mue, Xmu, boltzmann, ev, amu,
• BindH2,IonH,IonHe,gammaH,gammaH2,
• rScale,nScale,TempScale,pScale,lScale,VelScale,
• Lumin,mCentral,alpha,MinTemp,RunTimesc,ViscCD,ScaleGrav,ScaleCool,EOSConstants,
• verbosity,
• NrVars, nDim,
• iH2,iH,iHII,iHaux,iHe,iHeII,iHeIII,iE,iEntropy,nSpecies,nSpeciesLO,nSpeciesHI,iSpeedHI,
• lH,lHe,lHII,lHe,lHeII,lHeIII

Contains:

Subroutines:

• initneqCool
  • CALLs InitDMCool
• EOS_vars (PURE) (q,nvec,T,mu,gamma,gammais,eth,emicro,ke)
• nparticle (PURE) (nvec,npart,nneuc,ne)
• InitDMCool

Functions:

• get_gamma (PURE) (nvec)
  • CALLs SUBROUTINE nparticle
• get_mu (PURE) (nvec)
  • CALLs SUBROUTINE nparticle
• H2_cool (PURE) (nH2, nH, T)
  • Contains:
    • Functions:
      • Lvh (PURE) (T)
      • Lvl (PURE) (nH2,nH,T)
      • Qfunc (PURE) (nH2, nH)
      • kh2 (PURE) (T)
      • kh (PURE) (T)
      • Lrh (PURE) (T)
      • Lrl (PURE) (nH2,nH,T)
• H2_critical_cool (PURE) (nH2, T) (Commented out…obsolete)
  • Contains:
    • Functions:
      • fj (PURE) (j)
      • fv (PURE) (v)
      • gs (PURE) (j)
      • Bv (PURE) (v)
• H2_critical_cool_table (PURE) (nH2, T) (Commented out)
• CriticalDenH2 (PURE) (nH, T) (Commented out (had been used in H2_critical_cool))
• dust_cool (PURE) (n,T)
• dust_cool_table (PURE) (n,T)
• DMCoolingRate (PURE) (Temp)
• OI_cool (PURE) (nH2,nH,T)
• H2_diss (PURE) (nH2,nH,nHe,ne,T)
  • Contains:
    • FUNCTIONS:
      • kd (PURE) (i,T)
      • kdH_LS (PURE (i,T)
      • kdL_MK (PURE) (i,T)
      • kdL_DM (PURE) (i,T)
      • kdL_LS (PURE) (i,T)
• H2_dust_Recomb (PURE) (nHneuclei,nH,T,Tdust)
• H2_dust_Recomb_table (PURE) (nHneuclei,nH,T)
• H_recomb (PURE) (nHII,ne,T)
• H_recomb_table (PURE) (nHII,ne,T)
• H_ioniz (nH,ne,T) (why not pure?)
• H_ioniz_table (PURE) (nH,ne,T)
• He_ioniz (nHe,ne,T)
• He_ioniz_table (PURE) (nHeII,ne,T)
• He_recomb (PURE) (nHeII,ne,T)
  • Contains:
    • FUNCTIONS:
      • He_rad_recomb(PURE) (T)
      • He_dielec_recomb(PURE) (T)
• He_recomb_table (PURE) (nHeII,ne,T)
• df2 (PURE) (x)
• df1 (PURE) (x) - NOT USED (alternate for df2 but apparently doesn't work too well)

INCLUDES:

• i_evolove.f90
  • contrib/astro/physics/cool/i_evolove.f90
• (maybe should be i_evolve.f90 ?… i_evolove.f90 is the actual name, though)