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

ENDIF
CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XGLW)
CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XGSW)
CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XHFX)
CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XQFX)
CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XPBL)
CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XZRUF)
CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XCAPG)
CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XMOL)
CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XRGRND)
CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XRNET)
CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XFSOIL)
CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XCFRACH)
CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XCFRACM)
CALL BILIN2D (NDX,NCOLS_X,NROWS_X,XCFRACL)
C SECOND, BI-LINEAR INTERPOLATION OF 3-D FIELDS
CALL BILIN3D (NDX,NCOLS_X,NROWS_X,0,LURADM, XUSTAR)
CALL BILIN3D (NDX,NCOLS_X,NROWS_X,1,METLAY, XUU)
CALL BILIN3D (NDX,NCOLS_X,NROWS_X,1,METLAY, XVV)
CALL BILIN3D (NDX,NCOLS_X,NROWS_X,0,METLAY, XWWIND)
CALL BILIN3D (NDX,NCOLS_X,NROWS_X,1,METLAY, XPRESM)
CALL BILIN3D (NDX,NCOLS_X,NROWS_X,1,METLAY, XTEMPM)
CALL BILIN3D (NDX,NCOLS_X,NROWS_X,1,METLAY, XWVAPOR)
CALL BILIN3D (NDX,NCOLS_X,NROWS_X,1,METLAY, XCLDWTR)
CALL BILIN3D (NDX,NCOLS_X,NROWS_X,1,METLAY, XRANWTR)
ENDIF
C-----------------------------------------------------------------------
C
C CALCULATE DERIVED QUANTITIES SUCH AS LAYER HEIGHTS, DENSITY, etc.
c write(*,*) ' calling DERIVED '
c CALL DERIVED
DO R = 1, NROWS_X
DO C = 1, NCOLS_X
XDENSS(C, R) = XDENSAF(C, R, 0)
XDLUSE(C, R) = 1.
DO ILU = 2, LURADM
IF(XLUSE(C, R,ILU).GT.XLUSE(C, R,ILU-1)) XDLUSE(C, R) = ILU
ENDDO
ENDDO
ENDDO
C----------------------------------------------------------------------
77