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

C
C --- READ THE 2-D METEOROLOGICAL FIELDS FROM CALMET.DAT
C - PGT STABILITY CLASS,
C - FRICTION VELOCITY (M/S),
C - MIXING HEIGHT (M),
C - MONIN-OBUKHOV LENGTH (M),
C - CONVECTIVE VELOCITY SCALE (M/S),
C - PRECIPITATION RATE (MM/HR)
C --- (IF RUN TYPE NOT=0,I.E.,ONLY WINDS COMPUTED & STORED ON DISK)
IF(IRTYPE.EQ.1)THEN
CALL RDI2D(M1UNIT,IPGT,XBUF,MAXI,MAXJ,NX,NY,CLABEL,NDATHR)
CALL RDR2D(M1UNIT,USTR,XBUF,MAXI,MAXJ,NX,NY,CLABEL,NDATHR)
CALL RDR2D(M1UNIT,ZI,XBUF,MAXI,MAXJ,NX,NY,CLABEL,NDATHR)
CALL RDR2D(M1UNIT,EL,XBUF,MAXI,MAXJ,NX,NY,CLABEL,NDATHR)
CALL RDR2D(M1UNIT,WSTR,XBUF,MAXI,MAXJ,NX,NY,CLABEL,NDATHR)
IF(NPSTA.GT.0)CALL RDR2D(M1UNIT,RMM,XBUF,MAXI,MAXJ,NX,NY,CLABEL
& ,NDATHR)
C
C --- READ THE 1-D METEOROLOGICAL FIELDS FROM CALMET.DAT
C - AIR TEMPERATURE (DEG. K),
C - AIR DENSITY (KG/M**3),
C - SHORT-WAVE SOLAR RADIATION (W/M**2),
C - RELATIVE HUMIDITY (PERCENT),
C - PRECIPITATION CODE
CALL RDR1D(M1UNIT,TEMPK,NSSTA,CLABEL,NDATHR)
CALL RDR1D(M1UNIT,RHO,NSSTA,CLABEL,NDATHR)
CALL RDR1D(M1UNIT,QSW,NSSTA,CLABEL,NDATHR)
CALL RDI1D(M1UNIT,IRH,NSSTA,CLABEL,NDATHR)
IF(NPSTA.GT.0)CALL RDI1D(M1UNIT,IPCODE,NSSTA,CLABEL,NDATHR)
ENDIF
C Compute 2d-pressure, 2d-density, 3d-specific humidity
CALL INTCAL(nssta,xssta,yssta,rho,irh,tempk,
& elev,zface,ztemp,RHO_GRD,PSFC,PPRES,QQ)
C----------------------------------------------------------------------
C FILL UP THE ARRAYS FOR WHICH WE HAVE INFORMATION.
C NOTE: IN ORDER TO HAVE CALMET VARIABLES AGREE WITH THE WAY MCIP
C HANDLES I,J,K, WE NEED TO FLIP THEM !!!!!!!!! THIS IS DONE
C FOR ZTEMP,U,V,W,Z0,ELEV,ZI,EL,USTR,PSFC,RHO_GRD,PPRES
C ASSIGN 2D AND 3D PREDICTED FIELDS
DO J = 1, METCOL
DO I = 1, METROW
DO K = 1, METLAY
F_QQ (I,J,K) = QQ(J,I,K) ! from INTCAL subroutine
F_PP (I,J,K) = PPRES(J,I,K) ! from INTCAL subroutine
F_TT (I,J,K) = ZTEMP(J,I,K) ! from calmet.dat
ENDDO
85