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

XFLAGS (LXDENSAM, 3) = PNAME
XFLAGS (LXDENSAF, 2) = 'YES'
XFLAGS (LXDENSAF, 3) = PNAME
XFLAGS (LX3JACOBM,2) = 'YES'
XFLAGS (LX3JACOBM,3) = PNAME
XFLAGS (LX3JACOBF,2) = 'YES'
XFLAGS (LX3JACOBF,3) = PNAME
C compute layer heights above ground - no change with CALMET vgrid
DO K = 1, METLAY
DO C = 1, NCOLS_X
DO R = 1, NROWS_X
X3HTF (C,R,K) = XTOPO(C,R) + XX3FACE (K)
X3HTM (C,R,K) = XTOPO(C,R) + XX3MIDL (K)
XDX3HTF(C,R,K) = X3HTF (C,R,K) - X3HTF (C,R,K-1)
ENDDO
ENDDO
ENDDO
write(*,*) ' X3(1) in getmet = ', X3HTF (1,1,1), X3HTM(1,1,1)
C
SECOND = .TRUE.
FIRST = .FALSE.
write(*,*) ' end of first time getmet_calmet call'
RETURN
END IF ! of first time
C
C-----------------------------------------------------------------------
C
IF( SECOND ) THEN
SECOND = .FALSE.
ELSE
CALL READCM
ENDIF
write(*,*) ' processing getmet_calmet ... wait '
C
C-----------------------------------------------------------------------
C
C *** Loop over read timestep starts here
C
C MAPPING F_ARRAYS INTO X-ARRYAS ON ALL COINCIDING GRID POINTS
C EXTRACTING INFO FROM THE MM DOMAIN FOR THE EXTENDED CTM domain
C
DO JW = 1,NCG_J
74