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

C NEW RECORD -- #3 - ADDITIONAL RUN CONTROL DATA -- 8 WORDS
READ(LEVMET(1:6),'(I6)') ILEVMET
IF(ILEVMET.GE.980304) THEN
C NEW HEADER RECORD FORMAT
READ(M1UNIT)XLAT0M,XLON0M,LLCONFM,CONECM,XLAT1M,XLAT2M,
& RLAT0M,RLON0M
ENDIF
C
C CHECK THAT ARRAY MAX DIMENSIONS HAVE BEEN SIZED PROPERLY
IF ( NX.GT.MAXI.OR.NY.GT.MAXJ.OR.NZ.GT.MAXK.OR.
& NSSTA.GT.MAXSS.OR.NUSTA.GT.MAXUS.OR.NPSTA.GT.MAXPS)THEN
WRITE(*,2354)NX,MAXI,NY,MAXJ,NZ,MAXK,NSSTA,MAXSS,
& NUSTA,MAXUS,NPSTA,MAXPS
2354 FORMAT(/1X,'ERROR -- ARRAY DIMENSIONS ARE TOO SMALL FOR ',
& 'THE DATA IN THE CALMET.DAT FILE'/
& 5X,'NX = ',I5,4X,'MAXI = ',I5/
& 5X,'NY = ',I5,4X,'MAXJ = ',I5/
& 5X,'NZ = ',I5,4X,'MAXK = ',I5/
& 5X,'NSSTA = ',I5,4X,'MAXSS = ',I5/
& 5X,'NUSTA = ',I5,4X,'MAXUS = ',I5/
& 5X,'NPSTA = ',I5,4X,'MAXPS = ',I5)
GOTO 99 ! STOP
ENDIF
C
C RECORD #4 - CELL FACE HEIGHTS (NZ + 1 WORDS)
NZP1=NZ+1
CALL RDR1D(M1UNIT,ZFACE,NZP1,CLABEL,IDUM)
C
C RECORDS #5 & 6 - X, Y COORDINATES OF SURFACE STATIONS
C RELATIVE TO SW CORNER OF DOMAIN IN METERS
C CONVERT TO ABSOLUTE COORDINATES:
C X/Y(ORIGINAL)=X/Y(RELATIVE)+XORIGR/YORIGR
C (NSSTA WORDS EACH RECORD)
IF(NSSTA.GT.0)THEN
CALL RDR1D(M1UNIT,XSSTA,NSSTA,CLABXS,IDUM)
CALL RDR1D(M1UNIT,YSSTA,NSSTA,CLABYS,IDUM)
DO M=1,NSSTA
XSSTA(M) = XSSTA(M) + XORIGR
YSSTA(M) = YSSTA(M) + YORIGR
ENDDO
ENDIF
C
C RECORDS #7 & 8 - X, Y COORDINATES OF UPPER AIR STATIONS
C (NUSTA WORDS EACH RECORD)
IF(NUSTA.GT.0)THEN
CALL RDR1D(M1UNIT,XUSTA,NUSTA,CLABXU,IDUM)
CALL RDR1D(M1UNIT,YUSTA,NUSTA,CLABYU,IDUM)
ENDIF
C
C RECORDS #9 & 10 - X, Y COORDINATES OF PRECIPITATION STATIONS
81