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

C 3: winter dec, jan, feb
C 5: spring mar, apr, may ! season index 4 not used in PBLPKG
C 1: summer jun, jul, aug
C 2: autumn sep, oct, nov
C
C** local indices
INTEGER ISESN ! season #
INTEGER IP1, JP1, L ! loop indices
INTEGER C, R, K
INTEGER ILU, IZ, IW
C
C ------------------------------------------------------------------
C For "MAXK" layer meteorology
C new method for the estimation of surface fluxes
C a la Byun (Nov 1990, Jan 1997)
C-- NDEP: number of mm5 layers averaged into the lowest ctm layer
C-- to compute deposition velocities.
C
C** local arrays
C 2-D/3-D
REAL DTMPV (NCOLS_X,NROWS_X) ! virt pot. temp. difference
REAL FCORIS (NCOLS_X,NROWS_X) ! Coriolis factor
REAL TSURF (NCOLS_X,NROWS_X) ! surface temperature used
REAL QSURF (NCOLS_X,NROWS_X) ! mixing ratio at surface
REAL QRATIO (NCOLS_X,NROWS_X) ! mixing ratio at layer 1
c REAL RADIAT (NCOLS_X,NROWS_X) ! radiation
REAL WSPEED (NCOLS_X,NROWS_X,MAXK) ! wind speed
REAL UUSTC (NCOLS_X,NROWS_X) ! SPEED * USTAR = constant
REAL VPTEMP (NCOLS_X,NROWS_X) ! virt. pot. temperature
C
REAL PPSAVE (NCOLS_X,NROWS_X) ! array for saving PBL power
REAL RRSAVE (NCOLS_X,NROWS_X) ! array for saving PBL power
REAL Z0B1 (NCOLS_X,NROWS_X) ! temp. array for roughness
REAL Z0 (NCOLS_X,NROWS_X,LURADM) ! landuse dependent roughness
REAL BRATIO (NCOLS_X,NROWS_X) ! Estimate of Bowen Ratio
C 1-D
C
INTEGER ILANDU (MAXI,MAXJ) ! landuse from calmet
REAL WSTR (MAXI,MAXJ) ! wstar from calmet
COMMON ILANDU ! directly from calmet
COMMON WSTR ! directly from calmet
REAL RIBK ( MAXK )
REAL VPTEMPK( MAXK )
REAL ZL(MAXK), ZR(MAXK+1), UL(MAXK), DTHDZ(MAXK)
REAL ZTOP(MAXK), TL(MAXK)
REAL ZF(0:MAXK), QL(MAXK), PL(MAXK)
C
C** local scalar variables
REAL UNS, VNS ! wind component
100