A-1 APPENDIX A -- CALMET SCREENING CONFIGURATION

The tables below identify the CALMET configura tions used in developing the processed 6 km meteorological fields over the three CENRAP BART modeling domains shown in Figures 5-1 through 5-4. Also included in the tables belo w are the default CALMET options and parameter settings recommended in the IWAQM Phase 2 Report (EPA, 1998).

Table A-1. Input Groups in the CALMET Control File.

Input Group Description Applicable to CENRAP BART 0 Input and output file names Yes 1 General run control parameters Yes 2 Map Projection and Grid Control Parameters Yes 3 Output Options Yes 4 Meteorological Data Options Yes 5 Windfield Options and Parameters Yes 6 Mixing Height, Temperature and Precipitation Parameters Yes

7 Surface Meteorological Station Parameters Yes 8 Upper Air Meteorological Station Parameters Yes 9 Precipitation Station Parameters Yes

Table A-2. CALMET Model Input Group 0: Input and Output File Names.

Parameter Default CENRAP Comments Input GEO.DAT GEO.DAT Input SURF.DAT SURF.DAT Input CLOUD.DAT CLOUD.DAT Input PRECIP.DAT PRECIP.DAT Input MM4.DAT MM4.DAT Input WT.DAT WT.DAT Output CALMET.LST CALMET.LST Output CALMET.DAT CALMET.DAT Output PACOUT.DAT PACOUT.DAT NUSTA -- 0 Number of upper air stations NOWSTA -- 0 Number of over water met stations Input UP1.DAT UP1.DAT Input UP2.DAT UP2.DAT Input UP3.DAT UP3.DATA-2 Input SEA1.DAT SEA1.DAT Input DIAG.DAT DIAG.DAT Input PROG.DAT PROG.DAT Output TEST.PRT TEST.PRT Output TEST.OUT TEST.OUT Output TEST.KIN TEST.KIN Output TEST.FRD TEST.FRD Output TEST.SLP TEST.SLP

Table A-3. CALMET Model Input Group 1: General Run Control Parameters.

Parameter Default CENRAP Comments IBYR - 2001 Starting year IBMO - 1 Starting month IBDY - 1 Starting day IBHR - 1 Starting hour IBTZ - 6 Base time zone IRLG - 8736 Length of run IRTYPE 1 1 Run type (must = 1 to run CALPUFF) LCALGRD T F Compute CALGRID data fields ITEST 2 2 Stop run after SETUP to do input QA

Table A-4. CALMET Model Input Group 2: Map Projection and Grid Control Parameters.

Parameter Default CENRAP Comments PMAP UTM LCC Map Projection RLATO -- 40N Latitude (dec. degrees) of projection origin RLONO -- 97W Longitude (dec. degrees) of projection origin XLAT1 -- 33N Matching parallel(s) of latitude for projection XLAT2 -- 45N Matching parallel(s) of latitude for projection DATUM WGS-G WGS-G NX -- 300 Number of X grid cells in meteorological grid NY -- 192 Number of Y grid cells in meteorological grid DGRIDKM -- 6.0 Grid spacing, km XORIGKM -- -1008. Ref. Coordinate of SW corner of grid cell (1,1) YORIGKM -- 0.0 Ref. Coordinate of SW corner of grid cell (1,1) NZ -- 10 No. of vertical layers ZFACE -- 0, 20 40, Cell face heights in arbitrary vertical grid, mA-3 80, 160, 320, 640, 1200, 2000, 3000, 4000

Table A-5. CALMET Model Input Group 3: Output Options.

Parameter Default CENRAP Comments LSAVE T T Disk output option IFORMO 1 1 Type of unformatted output file LPRINT F F Print met fields IPRINF 1 1 Print intervals IUVOUT(NZ) NZ*0 NZ*0 Specify layers of u,v wind components to print IWOUT(NZ) NZ*0 NZ*0 Specify layers of w wind component to print ITOUT(NZ) NZ*0 NZ*0 Specify levels of 3-D temperature field to print LDB F F Print input met data and variables NN1 1 1 First time step for debug data to be printed NN2 1 1 Last time step for debug data to be printed IOUTD 0 0 Control variable for writing test/debug wind fields NZPRN2 1 0 Number of levels starting at surface to print IPR0 0 0 Print interpolated wind components IPR1 0 0 Print terrain adjusted surface wind components IPR2 0 0 Print initial divergence fields IPR3 0 0 Print final wind speed and direction IPR4 0 0 Print final divergence fields IPR5 0 0 Print winds after kinematic effects IPR6 0 0 Print winds after Froude number adjustment IPR7 0 0 Print winds after slope flows are addedA-4 IPR8 0 0 Print final wind field components

Table A-6. CALMET Model Input Group 4: Meteorological Data Options.

Parameter Default CENRAP Comments NOOBS 0 2 2 = No surface, overwater, or upper air observations; use MM5 for surface, overwater, and upper air data NSSTA -- 0 Number of meteorological surface stations NPSTA -- 0 Number of precipitation stations ICLOUD -- 3 Gridded cloud fields IFORMS 2 2 Formatted surface meteorological data file IFORMP 2 2 Formatted surface precipitation data file IFORMC 2 2 Formatted cloud data file

Table A-7. CALMET Model Input Group 5: Windfield Options and Parameters

.

Parameter Default CENRAP Comments IWFCOD 1 1 Model selection variable IFRADJ 1 1 Compute Froude number adjustment effects? IKINE 0 0 Compute kinematic effects? IOBR 0 0 Use O’Brien (1970) vertical velocity adjustment? ISLSOPE 1 1 Compute slope flow effects? IEXTRP -4 -1 Extrapolate surface wind obs to upper levels?

ICALM 0 0 Extrapolate surface winds even if calm?

BIAS NZ*0 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 Layer-dependent biases weighting aloft measurements

RMIN2 4. -1.0 Minimum vertical extrapolation distance IPROG 0 14 14 = Yes, use winds from MM5.DAT file as initial guess field [IWFCOD = 1 ISTEPPG 1 1 MM5 output timestep LVARY F T Use varying radius of influence RMAX1

-- 30. Maximum radius of influence over land in sfc layer RMAX2 -- 30. Maximum radius of influence over land aloftA-5 RMAX3 -- 50. Maximum radius of influence over water RMIN 0.1 0.1 Minimum radius of influence used anywhere TERRAD -- 12. Terrain features radius of influence R1 -- 1. Weighting of first guess surface field R2 -- 1. Weighting of first guess aloft field RPROG -- 0. MM5 windfield weighting parameter DIVLIM 5.E-6 5.E-6 Minimum divergence criterion NITER 50 50 Number of divergence minimization iterations NSMMTH 2, 4, 4, 4, 4, 4, 4 2, 4, 4, 4, 4, 4, 4 Number of passes through smoothing filter in each layer of CALMET

NITR2 99. 5, 5, 5, 5, 5, 5, 5, 5, 5, 5 Maximum number of stations used in each layer for the interpolation of data to a grid point

CRITFN 1.0 1.0 Critical Froude number ALPHA 0.1 0.1 Kinematic effects parameter FEXTR2 NZ*0.0 NZ*0.0 Scaling factor for extrapolating sfc winds aloft NBAR 0 0 Number of terrain barriers IDIOTP1 0 0 Surface temperature computation switch ISURFT -- 4 Number of sfc met stations to use for temp calcs IDIOPT2 0 0 Domain-averaged lapse rate switch IUPT 0 2 Upper air stations to use for lapse rate calculation ZUPT 200. 200. Depth through which lapse rate is calculated IDIOPT3 0 0 Domain-averaged wind component switch IUPWND -1 -1 Number of aloft stations to use for wind calc ZUPWND 1., 1000. 1., 1000. Bottom and top of layer through which the domain-scale winds are computed IDIOPT4 0 0 Observed surface wind component switch IDIOPT5 0 0 Observed aloft wind component switch LLBREZE F F Use Lake Breeze Module NBOX 0 0 Number of lake breeze regions NLB -- 0 Number of stations in the region METBXID(NLB) -- 0 Station ID’s in the region Table A-8. CALMET Model Input Group 6: Mixi ng Height, Temperature and Precipitation.

Parameter Default CENRAP CommentsA-6 CONSTB 1.41 1.41 Neutral stability mixing height coefficient CONSTE 0.15 0.15 Convective stability mixing height coefficient CONSTN 2400. 2400. Stable stability maxing height coefficient CONSTW 0.16 0.16 Overwater mixing height coefficient FCORIOL 1.E-4 1.E-4 Absolute value of Coriolis parameter IAVEZI 1 1 Conduct spatial averaging? Yes = 1 MNMDAV 1 10 Maximum search radius in averaging process HAFANG 30. 30. Half-angle of upwind looking cone for averaging ILEVZI 1 1 Layers of wind use in upwind averaging DPTMIN 0.001 0.001 Minimum potential temperature lapse rate in the stable layer above the current convective mixing ht DZZI 200. 200. Depth of layer above current conv. mixing height through which lapse rate is computed ZIMIN 50. 50. Minimum overland mixing height ZIMAX 3000. 3000. Maximum overland mixing height ZIMINW 50. 50. Minimum overwater mixing height ZIMAXW 3000. 3000. Maximum overwater mixing height ITPROG 0 2 3D temperature from observations or from MM5? IRAD 1 1 Type of interpolation; 1 = 1/R TRADKM 500. 36. Temperature interpolation radius of influence NUMTS 5 5 Max number of stations for temp interpolation IAVET 1 1 Spatially average temperatures? 1 = yes TGDEFB -.0098 -.0098 Temp gradient below mixing height over water TGDEFA -.0045 -.0045 Temp gradient above mixing height over water JWAT1 -- 55 Beginning land use categories over water JWAT2 -- 55 Ending land use categories for water NFLAGP 2 2 Precipitation interpolation flag; 2 = 1/R-squared

SIGMAP 100. 50. Radius of influence for precipitation interpolation CUTP 0.01 0.01 Minimum precipitation rate cutoff (mm/hr)