A-1 

APPENDIX A -- CALMET SCREENING CONFIGURATION

The tables below identify the CALMET configurations 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 below 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.DAT 
 

A-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, m 
 

A-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 added 
 

A-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 aloft
 

A-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: Mixing Height, Temperature and Precipitation. 

Parameter Default CENRAP Comments 
 

A-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)