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

C
write(*,*) ' X3(1) in getmet = ', X3HTF (1,1,1), X3HTM(1,1,1)
c write(*,*) ' Computing contravariant wind components ..'
C** Compute necessary 3-D met parameters using 'X' arrays
C compute arrays such as density
! CALL MET3DSUP ( SDATE, STIME, ISTAT, MOUTFRQ, INHYD )
C
C** Set up contravariant velocity components
C
C U- & V- Component
C (Note: XUU and XVV correspond to actual wind components at DOT)
C
DO K = 1, METLAY
DO C = 1, NCOLS_X
DO R = 1, NROWS_X
C
XUHAT (C,R,K) = XMAPD (C,R) * XUU (C,R,K)
XVHAT (C,R,K) = XMAPD (C,R) * XVV (C,R,K)
XWHAT (C,R,K) = XWWIND(C,R,K) ! JACOBIAN = 1.0 here
END DO
END DO
END DO
XFLAGS (LXUHAT, 2) = 'YES'
XFLAGS (LXUHAT, 3) = PNAME
XFLAGS (LXVHAT, 2) = 'YES'
XFLAGS (LXVHAT, 3) = PNAME
WRITE(*,*) 'RETURNING TO MCIP from GETMET_CALMET'
RETURN
C
1010 CONTINUE
C--------------------------------------------------------------
WRITE(*,*) 'TERMINATION OF GETMET_CALMET'
STOP
END
78