!
! Cooling terms are calculated at dz=0.25 regardless of model vertical 
!   resolution, see radcool.F.
! This file included in: radcool.F, recur.F, and vicool.F
!
      real,parameter :: RKO=3.0E-12 ! O-CO2 RATE
!     real,parameter :: RKO=4.0E-12 ! O-CO2 RATE
      real,parameter ::
     |  dzcool = 0.25,    ! delta zp resolution for cooling (includes tndown)
     |  zpsrf  = -21.5,   ! zp at ground surface
     |  zpm9   = -9.,     ! zp -9.0 level
     |  zpm5   = -5.,     ! zp -5.0 level
     |  zpm1925 = -19.25  ! zp -19.25 level
      integer,parameter ::
     |  lmax = 59, ! number of wavelengths (for hco2,ho3, see vicool)
     |  nzpsrf_zpm5 = int((zpm5-zpsrf)/dzcool)+1,  ! 67 for tv,o3
     |  nzpm9_zpm5  = int((zpm5-zpm9 )/dzcool)+1   ! 17 for sn2,so2,etc
!
      real :: tv,sn2,so2,o,co2,o3,flux,xpz,hco2,ho3,uco2
      COMMON/RADTRF/
     |  TV  (zimxp,nzpsrf_zpm5),  ! surface to zp -5  (76,67) 
     |  O3  (zimxp,nzpsrf_zpm5),  ! surface to zp -5  (76,67)
     |  XPZ (nzpsrf_zpm5),        ! surface to zp -5  (67)
     |  SN2 (zimxp,nzpm9_zpm5),   ! zp -9 to zp -5    (76,17)
     |  SO2 (zimxp,nzpm9_zpm5),   ! zp -9 to zp -5    (76,17)
     |  O   (zimxp,nzpm9_zpm5),   ! zp -9 to zp -5    (76,17)
     |  CO2 (zimxp,nzpm9_zpm5),   ! zp -9 to zp -5    (76,17)
     |  UCO2(zimxp,nzpm9_zpm5),   ! zp -9 to zp -5    (76,17)
     |  HCO2(zimxp,lmax),         ! set by vicool     (76,59)
     |  HO3 (zimxp,lmax),         ! set by vicool     (76,59)
     |  FLUX(zimxp)               ! set by vicool     (76)
!
      real :: am,den
      COMMON/RADTRE/
     |  AM (zimxp,nzpm9_zpm5),    ! zp -9 to zp -5 (76,17)
     |  DEN(zimxp,nzpm9_zpm5)     ! zp -9 to zp -5 (76,17)
!
!$OMP THREADPRIVATE (/radtrf/,/radtre/)
!DIR$ TASKCOMMON RADTRF,RADTRE