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

REAL ELEV_STA(MAXSS), ZMID(MAXK)
REAL DIST2
REAL X0, Y0
REAL SUM_DIST2,D2
REAL DELTA(MAXI,MAXJ), DELTAMIN
REAL DELTAX2, DELTAY2
REAL TEMP_0(MAXI,MAXJ)
REAL ESAT(MAXI,MAXJ,MAXK), QQSAT(MAXI,MAXJ,MAXK)
REAL T0(MAXSS)
REAL PWR, GAMMA_A, GAMMA, RDGASS
REAL PWR1, AA,BB,CC,TS0,P00,LNP
DATA GAMMA_A /0.0342/ !deg./m - adiabatic lapse rate
DATA GAMMA /0.0065/ !deg./m - lapse rate
C ------------------------------------------------------------
C SINCE WE DON'T HAVE THE ELEVATIONS OF SURFACE STATIONS,
C WE NEED TO FIND THE CLOSEST GRID POINT TO THEIR COORDINATES
C AND ASSIGN ITS ELEVATION AS AN ELEVATION OF THE SFC. STATION
C ------------------------------------------------------------
X0 = XORIGR - 0.5*RESOLN !The point to the west of i=1
Y0 = YORIGR - 0.5*RESOLN !The point to the south of j=1
C Initialize PSFC(i,j), RHO_GRD(i,j) and RH_GRD(i,j)
DO J=1,NY
DO I=1,NX
PSFC(I,J) = 0.0
RHO_GRD(I,J) = 0.0
RH_GRD(I,J) = 0.0
ENDDO
ENDDO
C Initialize qq(i,j,k)
DO J=1,NY
DO I=1,NX
DO K=1,NZ
QQ(I,J,K) = 0.0
ENDDO
ENDDO
ENDDO
C Find grid coordinates closest to station coordinates
DO M = 1,NSSTA
C Set 'deltamin' to maximum distance in domain
DELTAMIN = (FLOAT(MAXI)*RESOLN)**2 +
& (FLOAT(MAXJ)*RESOLN)**2
IMIN = -999
JMIN = -999
DO J = 1,NY
DO I = 1,NX
92