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

C
C
C
C
ZTEMP = 10.0 ! 10 meter
PTEMP = PL(1)+GRAV*XDENSAM(C,R,1)*(ZL(1)-ZTEMP)/100. !mb
CALL INTSFC ( UST, AMOL, THETA1,ZRUF, ZL(1),UL(1), THETA1
& ,ZTEMP, U2, THETA2)
XTEMP10(C,R) = THETA2*(PTEMP/1000.0)**.286
XSURF2 (C,R) = U2 ! for reference
ZTEMP = 1.5 ! 1.5 meter
PTEMP = PL(1)+GRAV*XDENSAM(C,R,1)*(ZL(1)-ZTEMP)/100. !mb
ZTEMP = 10.0 ! 10 meter
CALL INTSFC ( UST, AMOL, THETA1,ZRUF, ZL(1),UL(1), THETA1
& ,ZTEMP, U2, THETA2)
END DO
END DO
XSURF1 (C,R) = THETA2*(PTEMP/1000.0)**.286
C--- Compute W-Star
DO C = 1, NCOLS_X
DO R = 1,NROWS_X
XWSTAR (C,R) = XUSTAR (C,R,0) * ( XPBL (C,R)
& / ( VKAR * ABS( XMOL (C,R) ) ) )**0.3333333
c XWSTAR (C,R) = WSTR (C,R) ! directly from CALMET
XWSTAR (C,R) = CVMGP( 0.0, XWSTAR (C,R), XMOL (C,R) )
C
C--- Compute ALBEDO from MM5 values of incident and absorbed SW rad
IF(XRGRND(C,R).GT.0.) THEN
c XALBEDO(C,R) = 1. - XGSW(C,R) / XRGRND(C,R)
XALBEDO(C,R) = XALBEDO(C,R) ! from solar
ELSE
XALBEDO(C,R) = 0.
ENDIF
END DO
END DO
C
write(*,*) 'PBLSUP - RETURNING TO MCIP'
RETURN
C
1010 CONTINUE
STOP 'ABNORMAL TERMINATION IN PBLSUP'
END
110