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

RAUS( A1, B1, CKUST ) = PRO * LOG(
& ( ( 2.0 + A1 ) - 2.0 * SQRT( 1.0 + A1 ) ) *
* ( ( 2.0 + B1 ) + 2.0 * SQRT( 1.0 + B1 ) ) /
& ( A1 * B1 ) ) / CKUST
C
C aerodynamic resistance for a stable pbl
RASP( X1, X2, BETAB, CKUST ) = PRO * ( 2.0 * ( BETAB + 1.0 ) *
& ( 1.0 / SQRT( 1.0 - X2 ) - 1.0 / SQRT( 1.0 - X1 ) ) +
& LOG( ABS( ( -1.0 + SQRT( 1.0 - X2 ) )
& * ( 1.0 + SQRT( 1.0 - X1 ) ) /
& ( ( 1.0 + SQRT( 1.0 - X2 ) )
& * ( -1.0 + SQRT( 1.0 - X1 ) ) ) ) ) ) / CKUST
C
C aerodynamic resistance for a strongly stable PBL
C?? do we need a new formula??? - may be not??
C aerodynamic resistance for a mixed layer
RAUM( X1, X2, CKWST ) = LOG( X2 * ( 1.0 - X1 ) /
& ( X1 * ( 1.0 - X2 ) ) ) / CKWST
C ------------------------------------------------------------------
C saturation vapor pressure in pascal ! VP0 = 611.29
C TEMPC in C
E_AERK( TEMPC ) = VP0 *
& EXP( 17.625 * TEMPC / ( 243.04 + TEMPC ) )
C virtual temperature
C TEMP in K, QMIX is specific humidity
F_VTEMP( TEMP, QMIX ) = TEMP * ( 1.0 + 0.6077 * QMIX )
C ! where do we get 0.6077?
C potential temperature
C TEMP in K, PRESS is in Pascal
F_PTEMP( TEMP, PRESS ) = TEMP * ( 100000.0 / PRESS )**(2.0/7.0)
C
C***********************************************************************
C EXECUTABLE CODE STARTS HERE
C***********************************************************************
C ------------------------------------------------------------------
C
C INITIALIZE SOME VARIABLES
C
C ------------------------------------------------------------------
C
C** Compute wind speed profile from true horizontal wind data
DO K = 1, NLAYS
DO C = 1, NCOLS_X
IP1 = MIN( C+1, NCOLS_X )
DO R = 1,NROWS_X
JP1 = MIN( R+1, NROWS_X )
UNS = 0.25 * ( XUU(C, R,K) + XUU(IP1, R,K)
& + XUU(C,JP1,K) + XUU(IP1,JP1,K) )
103