CALPUFF and Postprocessors
CALPUFF and Postprocessors CALPUFF and Postprocessors
This process is simplified somewhat if a batch file is used to manage the filenames. One such batch file for DOS (RUNOPT.BAT) is included with the testcase. It requires three filenames as arguments: RUNOPT file1 file2 file3 where file1 OPTHILL.EXE executable program file file2 user.inp input file file3 user.out output file The batch file copies file2 to OPTHILL.INP, runs OPTHILL.EXE which creates OPTHILL.LST, then renames OPTHILL.LST to the name supplied as file3. For the example above, axis #1 would be processed by typing the command: RUNOPT OPTHILL.EXE AXIS1.INP AXIS1.LST With these results, hill information that is independent of the choice of coordinate system and the modeling grid for the wind model can be specified: xc,yc (m) thetah (deg) (depends on choice of coordinates) 69E zgrid (m) (depends on grid for wind model) relief (m) 300. expo (1) 1.91 expo (2) 1.24 scale (1) (m) 1523. scale (2) (m) 2896. axmax (1) (m) 2000. axmax (2) (m) 1500. Note that scale(2) is almost twice scale(1), even though axis 1 corresponds to the longer axis of the hill. This can occur because the "scale" parameter is a property of the entire inverse-polynomial function (Equation F-1), rather than just the portion of the function that is fit to the profile of the terrain. In Figure F-1, the shape of the terrain might best conform to the upper 10% of the polynomial function, in which case the "scale" parameter would exceed "axmax." In this example application of the OPTHILL program, we see that axmax(2) is substantially less than axmax(1), whereas scale(2) exceeds scale(1), indicating that a comparatively smaller portion of the polynomial function represents the terrain profile along the minor axis. MAR 2006 – CALPUFF F-283
G. POSTPROCESSORS The CALMET meteorological model generates a large unformatted meteorological data file which includes hourly gridded wind fields and temperature at multiple levels and hourly gridded surface meteorological fields such as PGT (Pasquill-Gifford-Turner) stability class, friction velocity, Monin- Obukhov length, mixing height, convective velocity scale, and precipitation rate. A postprocessor is designed to access this file: PRTMET is a postprocessor intended to aid in the analysis of the CALMET output data file by allowing the user to display selected portions of the meteorological data. After making one or more CALPUFF simulations, hourly concentrations and/or deposition fluxes for each species at each receptor exist in several unformatted data files. A single CALPUFF application can produce four such files: CONC.DAT (concentrations in g/m 3 ); WFLUX.DAT (wet deposition fluxes in g/m 2 /s); DFLUX.DAT (dry deposition fluxes in g/m 2 /s); and VISB.DAT (relative humidity for visibility analyses). When a period is simulated as a sequence of shorter-period CALPUFF applications, as when a year is simulated in chunks of about four weeks, for example, each of the shorter runs produces its own set of files. When certain groups of sources need to be characterized separately, individual CALPUFF applications are required for each group, again producing a set of data files for each application. Similarly, certain species may be modeled separately, producing more output data files. Data in all of these files must be processed to obtain results that can be used to characterize air quality impacts in terms of multiple-hour averages, increment consumption, threshold exceedences, visibility reduction, total deposition, and so forth. Postprocessors designed for this work include: APPEND CALSUM POSTUTIL CALPOST is a postprocessor which appends two or more sequential CALPUFF concentration, wet flux, dry flux or relative humidity (visibility) files in time. is a postprocessor which sums and scales concentrations or wet/dry fluxes from two or more source groups from different CALPUFF runs. is a postprocessor which operates on one or more CALPUFF concentration and wet/dry flux files to create new species as weighted combinations of modeled species; to sum wet and dry deposition fluxes; to merge species from different runs into a single output file; to sum and scale results from different runs; to repartition nitric acid/nitrate based on total available sulfate and ammonia; and to add time/space-varying background. is a postprocessor which operates on one CALPUFF concentration or wet/dry deposition flux file to perform visibility calculations; to average and summarize concentrations and deposition fluxes; to determine ranked concentration/ flux/light extinction values; and to create list files and plot files. MAR 2006 -- PostProc G-1
- Page 236 and 237: DFLX.DAT File - Data Records The DF
- Page 238 and 239: Table F-55 (Concluded) Unformatted
- Page 240 and 241: Character*132 COMMENT Character*80
- Page 242 and 243: Table F-56 (Continued) Unformatted
- Page 244 and 245: Table F-56 (Continued) Unformatted
- Page 246 and 247: Table F-56 (Continued) Unformatted
- Page 248 and 249: Table F-56 (Continued) Unformatted
- Page 250 and 251: Table F-56 (Concluded) Unformatted
- Page 252 and 253: Table F-57 Unformatted WFLX.DAT Fil
- Page 254 and 255: F.15.6 Relative Humidity File for V
- Page 256 and 257: Table F-58 Unformatted VISB.DAT fil
- Page 258 and 259: Table F-58 (Continued) Unformatted
- Page 260 and 261: Table F-58 (Continued) Unformatted
- Page 262 and 263: Table F-58 (Continued) Unformatted
- Page 264 and 265: Table F-58 (Continued) Unformatted
- Page 266 and 267: VISB.DAT File - Data Records The VI
- Page 268 and 269: Table F-59 (Concluded) Unformatted
- Page 270 and 271: Table F-60 Sample Debug Puff-Tracki
- Page 272 and 273: F.15.8 Mass Flux List File (MASSFLX
- Page 274 and 275: F.15.9 Mass Balance List File (MASS
- Page 276 and 277: F.16 OPTHILL When the subgrid scale
- Page 278 and 279: The following procedure is recommen
- Page 280 and 281: Table F-64 OPTHILL Input and Output
- Page 282 and 283: Figure F-2 Map of terrain surroundi
- Page 284 and 285: Table F-67 OPTHILL Output File for
- Page 288 and 289: G.1 PRTMET Meteorological Display P
- Page 290 and 291: Table G-1 PRTMET Input and Output F
- Page 292 and 293: Table G-2 (Continued) PRTMET Contro
- Page 294 and 295: Table G-2 (Concluded) PRTMET Contro
- Page 296 and 297: Table G-3 Sample PRTMET Control Fil
- Page 298 and 299: Table G-3 Sample PRTMET Control Fil
- Page 300 and 301: Table G-4 Sample PRTMET List File (
- Page 302 and 303: Table G-4 Sample PRTMET List File (
- Page 304 and 305: Table G-4 Sample PRTMET List File (
- Page 306 and 307: Table G-4 Sample PRTMET List File (
- Page 308 and 309: Table G-4 Sample PRTMET List File (
- Page 310 and 311: Table G-6 Sample vector plot file x
- Page 312 and 313: G.2 APPEND The APPEND program is de
- Page 314 and 315: Table G-9 APPEND Control File Struc
- Page 316 and 317: Table G-11 Sample APPEND List File
- Page 318 and 319: G.3 CALSUM The CALSUM program is de
- Page 320 and 321: Table G-13 CALSUM Control File Stru
- Page 322 and 323: Table G-15 Sample CALSUM List File
- Page 324 and 325: G.4 POSTUTIL POSTUTIL operates on o
- Page 326 and 327: The name and full path of each file
- Page 328 and 329: Table G-17 (Continued) POSTUTIL Con
- Page 330 and 331: Table G-18 Sample POSTUTIL Control
- Page 332 and 333: Table G-18 (Continued) Sample POSTU
- Page 334 and 335: INPUT GROUP: 2 -- Species Processin
This process is simplified somewhat if a batch file is used to manage the filenames. One such batch file<br />
for DOS (RUNOPT.BAT) is included with the testcase. It requires three filenames as arguments:<br />
RUNOPT file1 file2 file3<br />
where<br />
file1 OPTHILL.EXE executable program file<br />
file2 user.inp input file<br />
file3 user.out output file<br />
The batch file copies file2 to OPTHILL.INP, runs OPTHILL.EXE which creates OPTHILL.LST, then<br />
renames OPTHILL.LST to the name supplied as file3. For the example above, axis #1 would be<br />
processed by typing the comm<strong>and</strong>:<br />
RUNOPT OPTHILL.EXE AXIS1.INP AXIS1.LST<br />
With these results, hill information that is independent of the choice of coordinate system <strong>and</strong> the<br />
modeling grid for the wind model can be specified:<br />
xc,yc (m)<br />
thetah (deg)<br />
(depends on choice of coordinates)<br />
69E<br />
zgrid (m)<br />
(depends on grid for wind model)<br />
relief (m) 300.<br />
expo (1) 1.91<br />
expo (2) 1.24<br />
scale (1) (m) 1523.<br />
scale (2) (m) 2896.<br />
axmax (1) (m) 2000.<br />
axmax (2) (m) 1500.<br />
Note that scale(2) is almost twice scale(1), even though axis 1 corresponds to the longer axis of the hill.<br />
This can occur because the "scale" parameter is a property of the entire inverse-polynomial function<br />
(Equation F-1), rather than just the portion of the function that is fit to the profile of the terrain. In Figure<br />
F-1, the shape of the terrain might best conform to the upper 10% of the polynomial function, in which<br />
case the "scale" parameter would exceed "axmax." In this example application of the OPTHILL program,<br />
we see that axmax(2) is substantially less than axmax(1), whereas scale(2) exceeds scale(1), indicating<br />
that a comparatively smaller portion of the polynomial function represents the terrain profile along the<br />
minor axis.<br />
MAR 2006 – <strong>CALPUFF</strong><br />
F-283