Changes between Version 7 and Version 8 of u/bliu/coolingLength


Ignore:
Timestamp:
12/07/16 14:14:53 (8 years ago)
Author:
Baowei Liu
Comment:

—

Legend:

Unmodified
Added
Removed
Modified

  • u/bliu/coolingLength

    v7 v8  
    3737||  Cooling Time || 2.0021E+01 ||
    3838|| Cooling Length || 1.0043E+00 ||
     39
     404. Aluminum Cooling
     41
     42   FUNCTION AlCoolingStrength(q,Temp)
     43      REAL(KIND=qPrec) :: q(:), AlCoolingStrength, Temp, ne, nAl, logne, logtemp
     44      ne = (6d0/7d0)*q(1)*nScale  ! number densities based on Z = 6
     45      nAl = (1d0/7d0)*q(1)*nScale
     46      logne = LOG10(ne)
     47      Temp = Temp/1.16045d4       ! K to eV
     48      logtemp = LOG10(Temp)
     49      AlCoolingStrength=ne*nAl*GetTableValue(AlTab, (/logtemp, logne/))/(nScale**2d0)
     50      AlCoolingStrength=AlCoolingStrength*ScaleCool
     51   END FUNCTION AlCoolingStrength