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

! 0=NO TURBULENT KINETIC ENERGY
! FROM CALMET - INAV MUST = 0
YCENT_D= RLAT0M ! MET. DOMAIN CENTER LATITUDE
XCENT_D= -RLON0M ! MET. DOMAIN CENTER LONGITUDE
C ! (negative in Western Hemisphere)
P_BET_D= XLAT2M ! TRUE LATITUDE 1 (DEGREES)
P_ALP_D= XLAT1M ! TRUE LATITUDE 2 (DEGREES)
P_GAM_D= 90. ! POLE POSITION IN DEGREES LATITUDE???????
C
INHYD = 0 ! 1=NONHYDROSTATIC RUN, 0=HYDOSTATIC
C
RESOLN = DGRID ! DOMAIN GRID RESOLUTION IN METERS
I_11 = 1 ! I LOC. IN COARSE DOMAIN OF POINT (1,1)
J_11 = 1 ! J LOC. IN COARSE DOMAIN OF POINT (1,1)
C
C
F_ZTOP = ZFACE(NZP1) ! HEIGHT OF THE MODEL TOP (METERS)
TAPFRQ = 60.*60. ! TIME INTERVAL FOR DATA OUTPUT (MINS->S)
! CALMET.DAT OUTPUT FREQUENCY MUST BE 60min.
IF (METCOL.LT.(NROWS+2) .OR. METROW.LT.(NCOLS+2)
& .OR. METLAY.LT.NLAYS) THEN
CALL ERRMSG (PNAME,
& 'Met. data actual dimensions disagree with MCIPPARM')
WRITE(*,*)' Actual are ',METCOL,' * ',METROW,' * ',METLAY
WRITE(*,*)' MCIPPARM include file specifies NROWS,NCOLS,NLAYS',
& 'as', NROWS,' * ',NCOLS,' * ',NLAYS
WRITE(*,*)' To process met. data correctly, we should have '
WRITE(*,*)' METCOL >= NROWS+2, METROW >= NCOLS+2, and',
& ' METLAY >= NLAYS '
GOTO 99
ENDIF
NROW1 = NROWS
NCOL1 = NCOLS
IF(NDX.GT.1) NROW1 = NROW1 + 2
IF(NDX.GT.1) NCOL1 = NCOL1 + 2
IF (NROW1/NDX+2*NTHIK+1.GT.METCOL .OR.
& NCOL1/NDX+2*NTHIK+1.GT.METROW .OR.
& NLAYS.GT.METLAY) THEN
CALL ERRMSG (PNAME,
& 'Requested CTM dim. are larger than the actual Met. dim.')
WRITE(*,*)' Requested CTM dim.: ',NROWS,' x ',NCOLS,' x ',NLAYS
WRITE(*,*)' Met. dim.: ', METCOL,' x ',METROW,' x ',METLAY
WRITE(*,*)' Max. allowable dim.: ',METCOL-2*NTHIK-1,' x ',
& METROW-2*NTHIK-1,' x ',NLAYS
GOTO 99
ENDIF
FIRST = .FALSE.
83