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

C (NPSTA WORDS EACH RECORD)
IF(NPSTA.GT.0)THEN
CALL RDR1D(M1UNIT,XPSTA,NPSTA,CLABXP,IDUM)
CALL RDR1D(M1UNIT,YPSTA,NPSTA,CLABYP,IDUM)
ENDIF
C
C RECORD #11 - SURFACE ROUGHNESS LENGTHS (NX * NY WORDS)
CALL RDR2D(M1UNIT,Z0,XBUF,MAXI,MAXJ,NX,NY,CLABZ0,IDUM)
C
C RECORD #12 - LAND USE CATEGORIES (NX * NY WORDS)
CALL RDI2D(M1UNIT,ILANDU,XBUF,MAXI,MAXJ,NX,NY,CLABLU,IDUM)
C
C RECORD #13 - ELEVATIONS (NX * NY WORDS)
CALL RDR2D(M1UNIT,ELEV,XBUF,MAXI,MAXJ,NX,NY,CLABTE,IDUM)
C
C RECORD #14 - LEAF AREA INDEX (NX * NY WORDS)
CALL RDR2D(M1UNIT,XLAI,XBUF,MAXI,MAXJ,NX,NY,CLABLAI,IDUM)
C
C RECORD #15 - NEAREST SURFACE STATION TO EACH GRID POINT
C (NX * NY WORDS)
CALL RDI2D(M1UNIT,NEARS,XBUF,MAXI,MAXJ,NX,NY,CLABNSS,IDUM)
C
C ASSIGN HEADER INFO TO MCIP VARIABLES
IF(LLCONF) THEN
MAPPRO = 1
ELSE
MAPPRO = 3
ENDIF ! PROJECTION: 1=LAMBERT,3=MERCATOR
METCOL = NX ! CALMET GRID X DIMENSION
METROW = NY ! CALMET GRID Y DIMENSION
METLAY = NZ ! CALMET GRID Z DIMENSION
C ! START JULIAN DATE(YYYY,DDD,HH(GMT),MINS)
YYYY = IBYR
DDD = JULDAY(YYYY,IBMO,IBDY)
HH = IBHR + IBTZ
IF (HH.GT.24) THEN
HH = HH -24
DDD = DDD +1
ENDIF
IF (MOD(YYYY,4).EQ.0) THEN
MXDAYS = 366
ELSE
MXDAYS = 365
ENDIF
IF (DDD.GT.MXDAYS) THEN
DDD = DDD - MXDAYS
YYYY = YYYY + 1
ENDIF
MINS = 0
INAV = 0 ! 1=TURBULENT KINETIC ENERGY
82