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

C
XFLAGS (LXRNET, 2) = 'YES'
XFLAGS (LXRNET, 3) = PNAME
DO ILU = 1, LURADM
DO C = 1, NCOLS_X
DO R = 1, NROWS_X
Z0 (C,R,ILU) = Z00 (ILU,ISESN)
ENDDO
ENDDO
ENDDO
C
C Compute the logarithmic average surface roughness for each grid area
C Byun and Wesley's formula (in preparation)
DO 500 IZ=1, ITER ! ITER = 6
DO C = 1, NCOLS_X
DO R = 1, NROWS_X
XZRUF (C,R) = 0.0
Z0B1 (C,R) = 0.0
C XRADYN (C,R,0) = 1.0/SMALL ! initially, for avg
XRADYN (C,R,0) = 0.0 ! for inverse !
ENDDO
ENDDO
C
DO ILU = 1, LURADM
DO C = 1, NCOLS_X
DO R = 1, NROWS_X
Z0B1 (C,R) = Z0B1 (C,R) + XLUSE (C,R,ILU) *
& SQRT( LOG( ZREF / Z0 (C,R,ILU) ) )
ENDDO
ENDDO
ENDDO
C
DO C = 1, NCOLS_X
DO R = 1, NROWS_X
XZRUF (C,R) = ZREF * EXP( -Z0B1 (C,R)**2 )
ENDDO
ENDDO
C
C Now, XZRUF (C,R) contains regular roughness length in m
C Recompute sea-surface rougness using Charnock's relation
C Recompute cell-average surface fluxes
C --- this part is not land use dependent ---
C --- except for adjusting for water ---
C
C Rough estimate of Bowen Ratio for computation of XQFX
DO C=1, NCOLS_X
DO R=1,NROWS_X
BRATIO (C,R) = 0.10
IF(ILANDU(C,R).EQ.70) BRATIO(C,R) = 0.03
IF(ILANDU(C,R).GE.49.AND.ILANDU(C,R).LE.56) BRATIO(C,R)=0.15
ENDDO
106