Alternative small scale meteorology input to a chemical transport ...

C
REAL DTEMPMAX ! maximum sfc_T - air_T
REAL DTEMPMIN ! minimum sfc_T - air_T
DATA DTEMPMAX / 5.0 /, DTEMPMIN / -3.0 /
INTEGER ITER, IWAT
DATA ITER / 6 /, IWAT / 6 /
C****************************************
C WESELY RADM LAND-USE TYPES
c 1. urban land
c 2. agriculture
c 3. range
c 4. deciduous forest
c 5. coniferous forest
c 6. mixed forest wetland
c 7. water
c 8. barren land
c 9. non-forested wetland
c 10. mixed agriculture/rangeland
c 11. rocky open areas with low shrubs
C***********************************************************************
C **************** STATEMENT FUNCTIONS *****************
C** FOR PBL
REAL RASP ! aerodynamic resistance function for stable PBL
REAL RASS ! aerodynamic resistance function for stable sfc layer
REAL RAUS ! aerodynamic resistance function - unstable sfc layer
REAL RAUM ! aerodynamic resistance function for mixed layer
C
C** for moisture & virtual temperature
REAL E_AERK
REAL F_VTEMP
REAL F_PTEMP
C ------------------------------------------------------------------
C variables used in the function definitions
REAL DD ! displacement height
REAL ZU ! upper limit of integration
REAL AMOL ! Monin-Obukhov length
REAL CKUST ! VKAR*USTAR
REAL A1 ! 9.*ZU/ABS(MOL)
REAL B1 ! 9.*D /ABS(MOL)
REAL X1,X2 ! Z1/PBLHT, Z2/PBLHT
REAL BETAB ! (4.7/0.74)*PBLHT/MOL
REAL CKWST ! VKAR*WSTAR
REAL ZETA
C
REAL D, TEMP, TEMPC, QMIX, PRESS
C ------------------------------------------------------------------
C aerodynamic resistance for a stable surface layer
RASS( D, ZU, AMOL, CKUST ) = ( PRO * LOG( ZU / D )
& + BETAM * (ZU-D) / AMOL ) / CKUST
C
C aerodynamic resistance for an unstable surface layer
102