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

C
C
VNS = 0.25 * ( XVV(C, R,K) + XVV(IP1, R,K)
& + XVV(C,JP1,K) + XVV(IP1,JP1,K) )
ENDDO
ENDDO
ENDDO
WSPEED (C, R, K) = MAX( 0.5, SQRT( UNS**2 + VNS**2 ) )
C
DO C = 1, NCOLS_X
DO R = 1,NROWS_X
C
C saturation vapor pressure (mb) over water
C Estimating water vapor at surface with TSURF = 0.5*(TEMPM(1)+TEMPG)
C assuming it has the same relative humidity as the layer 1 air,
C but at the reference temperature
C
C R.H = vapor_press / sat_vapor_press
C Use AERK formula
C
TSURF (C,R) = XTEMPM (C,R,1) ! test
VSAT = E_AERK( TSURF (C,R) - STDTEMP ) ! in Pascal
VPRESS = XDENSWM (C,R,1) * RWVAP * XTEMPM (C,R,1)
RHUMID = VPRESS / VSAT
PSURF = XPRESM(C,R,1) + GRAV*XDENSAM(C,R,1)*X3HTM(C,R,1)
VPRESS = VSAT * RHUMID
DENSW = VPRESS / ( RWVAP * XTEMPM (C,R,1) )
DENSD = XDENSAF(C,R,0)
QSURF (C,R) = DENSW / DENSD
C
C virtual temperature of ground and air above ground
VPTSURF = F_VTEMP( TSURF (C,R), QSURF (C,R) )
VPTSURF = F_PTEMP( VPTSURF, PSURF )
C
QRATIO(C,R) = XDENSWM (C,R,1) /
& ( XDENSAM (C,R,1) - XDENSWM (C,R,1) )
VPTEMP(C,R) = F_VTEMP( XTEMPM (C,R,1), QRATIO(C,R) ) !
VPTEMP(C,R) = F_PTEMP( VPTEMP (C,R), XPRESM (C,R,1) ) !
DTMPV (C,R) = VPTSURF - VPTEMP (C,R)
IF ( DTMPV(C,R) .LT. DTEMPMIN ) DTMPV (C,R) = DTEMPMIN
IF ( DTMPV(C,R) .GT. DTEMPMAX ) DTMPV (C,R) = DTEMPMAX
C ! limit for stable & unstable case differently
DTMPV (C,R) = CVMGZ( 1.E-10, DTMPV (C,R), DTMPV (C,R) )
ENDDO
ENDDO
104