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

C
640 CONTINUE
C
C** compute cell average (using inverse directly)
DO C = 1, NCOLS_X
DO R = 1, NROWS_X
XRADYN (C,R,0) = XRADYN (C,R,0) +
& XLUSE (C,R,ILU) / XRADYN (C,R,ILU)
END DO
END DO
600 CONTINUE ! end of 'DO 600' landuse loop
C** now XRADYN has real dynamic resistance
C** ?? probably, we do not want to replace the original value ??
DO C = 1, NCOLS_X
DO R = 1, NROWS_X
XRADYN (C,R,0) = 1. / XRADYN (C,R,0)
END DO
END DO
XFLAGS (LXWSTAR,2) = 'YES'
XFLAGS (LXWSTAR,3) = PNAME
XFLAGS (LXRADYN,2) = 'YES'
XFLAGS (LXRADYN,3) = PNAME
C** Estimate temperatures at 1.5 m and 10 m
DO C = 1, NCOLS_X
DO R = 1, NROWS_X
ZF(0) = X3HTF(C,R,0)
C
DO K = 1, NLAYS
ZL (K) = X3HTM (C,R,K)
ZF (K) = X3HTF (C,R,K)
UL (K) = WSPEED (C,R,K)
TL (K) = XTEMPM (C,R,K)
QL (K) = XDENSWM (C,R,K) /
& ( XDENSAM (C,R,K) - XDENSWM (C,R,K) )
PL (K) = XPRESM (C,R,K) / 100.0 ! in hPascal = mb
END DO
C TSFC = 0.5 * ( XTEMPM (C,R,1) + XTEMPG (C,R) ) ! test
AMOL = XMOL (C,R)
ZRUF = XZRUF (C,R)
HPBL = XPBL (C,R)
C*bdx CPAIR = 1004.67 * (1.0 + 0.84 * QL(1)) ! J/(K KG)
CPAIR = CPD * (1.0 + 0.84 * QL(1)) ! J/(K KG)
UST = XUSTAR (C,R,0)
CBDX120197 PSFC = ( XPSTAR (C,R) + X3TOP ) / 100.0 !in hPascal=mb
PSURF=XPRESM(C,R,1)+GRAV*XDENSAM(C,R,1)*X3HTM(C,R,1)
C ! in Pascal
PSFC = PSURF / 100. ! in hPascal= mb
THETA1 = TL(1) *(1000.0/PL(1))**.286
109