ADDENDUM USER'S GUIDE FOR THE AMS/EPA REGULATORY MODEL - AERMOD (EPA-454/B-03-001, September 2004) U.S. ENVIRONMENTAL PROTECTION AGENCY Office of Air Quality Planning and Standards Air Quality Assessment Division Research Triangle Park, North Carolina 27711 September 2009 ------- PREFACE This document provides updated user instructions for the AERMOD dispersion model, including modifications introduced with version 06341 and later. This addendum supplements and updates the information contained in the current AERMOD User's Guide (EPA, 2004a). 11 ------- ACKNOWLEDGEMENTS Portions of this Addendum to the "User's Guide for the AMS/EPA Regulatory Model - AERMOD" related to the PSDCREDIT option for the Plume Volume Molar Ratio Method (PVMRM) were prepared by MACTEC Federal Programs, Inc., Research Triangle Park, North Carolina. This effort was funded by the Environmental Protection Agency, Region 10, under Contract No. EP-D-05-096, with Herman Wong as the Work Assignment Manager (WAM). in ------- CONTENTS PREFACE ii ACKNOWLEDGEMENTS iii TABLES vi 1.0 INTRODUCTION 1 1.1 OVERVIEW OF AERMOD REVISIONS 1 1.2 BACKGROUND ON DEPOSITION ALGORITHMS 3 2.0 USER INSTRUCTIONS 6 2.1 DISPERSION MODELING OPTIONS 6 2.1.1MODELOPT Keyword 6 2.1.2 BETA Test Options 10 2.1.3 Processing for Paniculate Matter (PM) NAAQS 11 2.1.4 Specifying Multiple Urban Areas 14 2.1.5 BETA Options for Capped and Horizontal Stack Releases 15 2.2 DEPOSITION ALGORITHM INPUTS AND OPTIONS 17 2.2.1 Definition of Seasons for Gas Dry Deposition 17 2.2.2 Definition of Land Use Categories for Gas Dry Deposition 18 2.2.3 Option for Overriding Default Parameters for Gas Dry Deposition 19 2.2.4 Specifying Source Parameters for Gas Deposition (Dry and/or Wet) 20 2.2.5 Option for Specifying the Deposition Velocity for Gas Dry Deposition 20 2.2.6 Specifying Source Parameters for Particle Deposition 21 2.2.7 Specifying Emission and Output Units 23 2.2.8 Deposition Velocity and Resistance Outputs 24 2.2.9 Meteorological Data for Deposition Algorithms 25 2.3 OPEN PIT SOURCE OPTION 26 2.4 PVMRM AND OLM OPTIONS FOR MODELING NO2 28 2.4.1 Specifying Ozone Concentrations for PVMRM and OLM Options 29 2.4.2 Specifying the Ambient Equilibrium NO2/NOx Ratio for PVMRM 30 2.4.3 Specifying the Default In-stack NO2/NOx Ratio for PVMRM and OLM 31 2.4.4 Specifying In-stack NO2/NOx Ratios by Source for PVMRM and OLM 31 2.4.5 Specifying Combined Plumes for OLM 32 2.4.6 Modeling NO2 Increment Credits with PVMRM 34 2.5 VARIABLE EMISSION RATES 39 2.5.1 Specifying Variable Emission Factors (EMISFACT) 39 2.5.2 Specifying an Hourly Emission Rate File (HOUREMIS) 43 2.6 OUTPUT FILE OPTIONS 45 2.7 MODEL STORAGE LIMITS 46 3.0 REFERENCES 48 APPENDIX A. ALPHABETICAL KEYWORD REFERENCE A-l IV ------- APPENDIX B. FUNCTIONAL KEYWORD/PARAMETER REFERENCE B-l APPENDIX C. LIST OF ERROR/WARNING MESS AGES C-l APPENDIX D. EPA Model Clearinghouse Memorandum Dated July 9, 1993 D-1 ------- TABLES Table Page 1-1. Summary of Deposition Options in AERMOD 5 B-l. DESCRIPTION OF CONTROL PATHWAY KEYWORDS B-3 B-2. DESCRIPTION OF CONTROL PATHWAY KEYWORDS AND PARAMETERS... B-5 B-3. DESCRIPTION OF SOURCE PATHWAY KEYWORDS B-9 B-4. DESCRIPTION OF SOURCE PATHWAY KEYWORDS AND PARAMETERS .... B-11 B-5. DESCRIPTION OF RECEPTOR PATHWAY KEYWORDS B-16 B-6. DESCRIPTION OF RECEPTOR PATHWAY KEYWORDS AND PARAMETERS B-17 B-7. DESCRIPTION OF METEOROLOGY PATHWAY KEYWORDS B-20 B-8. DESCRIPTION OF METEOROLOGY PATHWAY KEYWORDS AND PARAMETERS B-21 B-9. DESCRIPTION OF EVENT PATHWAY KEYWORDS B-23 B-10. DESCRIPTION OF EVENT PATHWAY KEYWORDS AND PARAMETERS B-24 B-ll. DESCRIPTION OF OUTPUT PATHWAY KEYWORDS B-25 B-12. DESCRIPTION OF OUTPUT PATHWAY KEYWORDS AND PARAMETERS .... B-26 VI ------- 1.0 INTRODUCTION 1.1 OVERVIEW OF AERMOD REVISIONS This document provides user instructions for revisions to the AERMOD dispersion model. The discussion provided here supplements and updates the information contained in the current AERMOD User's Guide (EPA, 2004a), and it is assumed that the reader is already familiar with the contents of that document. Note that portions of the current AERMOD User's Guide (EPA, 2004a) are no longer valid or applicable. Four sets of AERMOD revisions are included in this Addendum: 1. The first set of revisions, first introduced with version 03273 of AERMOD, includes dry and wet deposition algorithms for both particulate and gaseous emissions (see Sections 1.2 and 2.2 for more details), and the OPENPIT source option, originally incorporated in the ISCST3 model (EPA, 1995a), for modeling particulate emissions from open pit (below grade) sources, such as surface coal mines and rock quarries; 2. The second set of revisions, first introduced with version 04300 of AERMOD, includes two non-DFAULT options for modeling conversion of NOX to NO2: 1) the Plume Volume Molar Ratio Method (PVMRM) (Hanrahan, 1999a and 1999b); and 2) the Ozone Limiting Method (OLM); 3. The third set of revisions, first introduced with version 06341 of AERMOD, includes the following: a. A new "BETA" option on the CO MODELOPT card to allow for new features to be added to AERMOD that are still in BETA-test status; b. A BETA option for incorporating NO/NO2 chemistry for NO2 increment consumption calculations with PSD credits using the PVMRM option; c. BETA options for capped and horizontal stack releases; d. An option to specify an initial default in-stack NO2/NOX ratio for the PVMRM and OLM options; e. New options for varying emissions by month, hour-of-day, and day-of-week (MHRDOW and MHRDOW7); f. An option to allow multiple urban areas to be defined in a single model run; ------- g. Updated processing to support modeling demonstrations for the National Ambient Air Quality Standards (NAAQS) for PM, including the 24-hour average design value for PM-2.5 impacts; and h. Use of dynamic array allocation for AREAPOLY sources to allocate array limits for the maximum number of vertices at model runtime, replacing the previous fixed array limit of 20 vertices. 4. The fourth set of revisions, first introduced with version 09260 of AERMOD, includes the following (additional information is provided in Model Change Bulletin (MCB) #3 provided on the SCRAM AERMOD webpage): a. New options for varying emissions by hour-of-day and day-of-week (HRDOW and HRDOW7); b. Modification of the regulatory default option (DFAULT) on the CO MODELOPT card to impose a restriction on the urban roughness length parameter to be 1 meter for regulatory default applications. Any value other than 1 meter for the urban roughness length option on the CO URBANOPT card will be treated as a non-DFAULT option; c. Removal of the TOXICS option from the MODELOPT key word. Options formerly associated with the TOXICS option are still considered non-DFAULT options within AERMOD. The area source optimizations formerly associated with the TOXICS option are now selected using the new non-DFAULT FASTAREA option on the MODELOPT keyword; d. A new non-DFAULT option for optimizing runtime for POINT and VOLUME sources based on an alternative implementation of the horizontal meander algorithm has been incorporated through the FASTALL option on the MODELOPT keyword. The FASTALL option also includes the FASTAREA optimizations if area sources are included in the model inputs; e. The option for specifying hourly emissions from a separate file through the HOUREMIS keyword has been enhanced to allow the use of hourly varying release heights and initial dispersion coefficients for VOLUME and AREA/AREAPOLY/AREACIRC sources; f The OPENPIT source option has been modified to allow for use of the OPENPIT source for gaseous (non-particulate) emissions; g. A new SUMMFILE option has been included on the OU pathway to output the summary of high ranked values to a separate file; h. The non-DFAULT option of FLAT terrain can now be specified on a source-by- source basis, allowing both FLAT and ELEV terrain treatments within the same model run (see Section 4.1 of'the AERMOD Implementation Guide regarding modeling of sources with terrain-following plumes in sloped terrain); i. A non-DFAULT option for a user-specified dry deposition velocity for gaseous emissions has been added under the GASDEPVD keyword on the CO pathway; ------- j. An option to use scientific notation for output result files has been added through the FILEFORM keyword on the OU pathway. The FILEFORM option is applicable to PLOTFILEs, plot-formatted POSTFILEs, MAXIFILEs, RANKFILEs, and SEASONHR files; k. An option (WARNCHKD) has been added to the MODELOPT keyword to allow issuing of warning messages rather than fatal errors for non-sequential meteorological data files, in order to allow use of multi-year meteorological data files that may contain gaps between years of data under the DF AULT option; and 1. The maximum length of filenames specified in the 'aermod.inp' file has been increased to 200 (controlled by the ILEN_FLD parameter in modules.f), and the maximum input string length to 512 (controlled by the ISTRG parameter in modules.f). Double quotes (") are also allowed as field delimiters in the 'aermod.inp' file to support filenames with embedded spaces. 1.2 BACKGROUND ON DEPOSITION ALGORITHMS The deposition algorithms incorporated into AERMOD are based on the draft Argonne National Laboratory (ANL) report (Wesely et al., 2002), with modifications based on peer review. Treatment of wet deposition was revised from Wesely et al. (2002) based on recommendations by peer review panel members (Walcek et al., 2001). A full technical description of the deposition algorithms implemented in AERMOD is provided in an EPA report (EPA, 2003). The deposition algorithms based on the ANL report were initially implemented in the AERMOD model under the non-DFAULT TOXICS option, which was selected by including the TOXICS keyword on the CO MODELOPT card. Beginning with version 09260 of AERMOD, the TOXICS option has been removed from the model. Those options in AERMOD that were formerly associated with the TOXICS option are still considered non-DFAULT options but no longer require the specification of the TOXICS option to allow their use. The other changes to the AERMOD model inputs associated with the deposition algorithms are limited to the CO (control) pathway and the SO (source) pathway. For gaseous dry deposition based on the ANL algorithms, the user must define the seasonal categories based on the ANL report for each of the calendar months, and must also define the land use category and three pollutant-specific physical parameters that are provided in the appendices of the ANL report. An optional keyword is also ------- provided to override default values for three parameters used in the gas deposition algorithm. The input requirements for "Method 1" particle deposition in AERMOD are the same as for the particle deposition algorithm in the ISCST3 model and are described below in Section 2.2. For "Method 2" particle deposition, the user must define the fraction of the particle mass in the fine particle category (less than 2.5 microns) and a representative mass mean diameter for the particles, which are also provided for selected pollutants in Appendix B of the ANL report. The keywords used to define these inputs and the meteorological data requirements for deposition are described in Section 2.2. Consistent with Section 7.2.7(b) of the Guideline on Air Quality Models (40 CFR Part 51, Appendix W), the use of "Method 1" for particle deposition in AERMOD is allowed under the regulatory DFAULT option. However, use of the "Method 1" particle deposition algorithms require quantification of the particle size distribution and should be done in coordination with the appropriate reviewing authority. Use of the "Method 2" option for particle deposition and the gas deposition algorithms (both dry and wet) are considered non-DFAULT options in AERMOD. Table 1-1 summarizes the required keywords for the various deposition options within AERMOD and whether they are allowed under the DFAULT option. ------- Table 1-1. Summary of Deposition Options in AERMOD Pollutant Type Gaseous Gaseous Gaseous Particulate ("Method 1") Particulate ("Method 2") Model Output Type CONC w/dry depletion DDEP CONC w/wet depletion WDEP CONC w/dry & wet depletion DEPOS CONC w/dry and/or wet depletion DEPOS DDEP WDEP CONC w/dry and/or wet depletion DEPOS DDEP WDEP Required Keywords CO GASDEPVD or CO GDSEASON, CO GDLANUSE, and SO GASDEPOS SO GASDEPOS CO GDSEASON, CO GDLANUSE, and SO GASDEPOS SO PARTDIAM, SO PARTDENS, and SO MASSFRAX SO METHOD 2 Allowed under DFAULT? No No No Yes No ------- 2.0 USER INSTRUCTIONS 2.1 DISPERSION MODELING OPTIONS 2.1.1 MODELOPT Keyword The dispersion options are controlled by the MODELOPT keyword on the CO pathway. The syntax, type, and order of the MODELOPT keyword are summarized below: Syntax: CO MODELOPT Type: Order: DFAULT BETA CONG AREADPLT FLAT NOSTD NOCHKD NOWARN SCREEN DEPOS and/or or DDEP ELEV WARNCHKD and/or WDEP DRYDPLT WETDPLT or or NODRYDPLT NOWETDPLT Mandatory, Non-repeatable Must precede POLLUTID, HALFLIFE and DCAYCOEF SCIM PVMRM PSDCREDIT FASTALL or or OLM FAS TARE A where: DFAULT- BETA- CONC- DEPOS- DDEP- WDEP- Specifies that the regulatory default options will be used; note that specification of the DFAULT option will override some non-DFAULT options that may be specified in the input file, while other non- DFAULT options will cause fatal errors when DFAULT is specified (see below for details); Non-DFAULT option that allows for draft, "Beta" test options to be used; currently includes the PSDCPxEDIT option and options for capped and horizontal stack releases (see Section 2.1.2 for more details); Specifies that concentration values will be calculated; Specifies that total deposition flux values (both dry and wet) will be calculated; Specifies that dry deposition flux values will be calculated; Specifies that wet deposition flux values will be calculated; ------- AREADPLT - FLAT- ELEV- NOSTD - NOCHKD- WARNCHKD- NOWARN- SCREEN - SCIM- PVMRM- OLM- PSDCREDIT - Specifies that a non-DFAULT method for optimized plume depletion due to dry removal mechanisms will be included in calculations for area sources; Specifies that the non-DFAULT option of assuming flat terrain will be used; Note that FLAT and ELEV may be specified in the same model run to allow specifying the non-DFAULT FLAT terrain option on a source-by-source basis; Specifies that the default option of assuming elevated terrain will be used; Note that FLAT and ELEV may be specified in the same model run to allow specifying the non-DFAULT FLAT terrain option on a source-by-source basis; Specifies that the non-DFAULT option of no stack-tip downwash will be used; Specifies that the non-DFAULT option of suspending date checking will be used for non-sequential meteorological data files; Specifies that the option of issuing warning messages rather than fatal errors will be used for non-sequential meteorological data files; Specifies that the option of suppressing the detailed listing of warning messages in the main output file will be used (the number of warning messages is still reported, and warning messages are still included in the error file controlled by the CO ERRORFIL keyword); Specifies that the non-DFAULT option for running AERMOD in a screening mode will be used; Sampled Chronological Input Model - used only with the ANNUAL average option to reduce runtime by sampling meteorology at a user- specified regular interval; SCEVI sampling parameters must be specified on the ME pathway; Specifies that the non-DFAULT Plume Volume Molar Ratio Method (PVMRM) for NC>2 conversion will be used; Specifies that the non-DFAULT Ozone Limiting Method (OLM) for NO2 conversion will be used; Specifies that the non-DFAULT BETA test option will be used to calculate the increment consumption with PSD credits using the PVMRM option; 7 ------- FASTALL - Non-DFAULT option to optimize model runtime through use of alternative implementation of horizontal meander for POINT and VOLUME sources; also optimizes model runtime for AREA/ AREAPOLY/AREACIRC and OPENPIT sources through hybrid approach (formerly associated with TOXICS option, now controlled by FASTAREA option); FASTAREA - Non-DFAULT option to optimize model runtime through hybrid approach for AREA/ AREAPOLY/AREACIRC and OPENPIT sources (formerly associated with TOXICS option); DRYDPLT - Option to incorporate dry depletion (removal) processes associated with dry deposition algorithms; this requires specification of dry deposition source parameters and additional meteorological variables; dry depletion will be used by default if dry deposition algorithms are invoked; NODRYDPLT - Option to disable dry depletion (removal) processes associated with dry deposition algorithms; WETDPLT - Option to incorporate wet depletion (removal) processes associated with wet deposition algorithms; this requires specification of wet deposition source parameters and additional meteorological variables; wet depletion will be used by default if wet deposition algorithms are invoked; and NOWETDPLT - Option to disable wet depletion (removal) processes associated with wet deposition algorithms. The DFAULT option has been modified beginning with version 09260 to impose a restriction on the optional urban roughness length parameter to be 1 meter for regulatory default applications. The user may select any or all of the output types (CONC, DEPOS, DDEP and/or WDEP) to be generated in a single model run. The order of these secondary keywords on the MODELOPT card has no effect on the order of results in the output files - the outputs will always be listed in the order of CONC, DEPOS, DDEP, and WDEP. Appropriate deposition parameters must be specified in order to output deposition fluxes using the DEPOS, DDEP, and/or WDEP keywords (see Sections 1.2 and 2.2 for more details). ------- Beginning with version 04300, the dry and/or wet removal (depletion) mechanisms (the DRYDPLT and WETDPLT options in earlier versions of AERMOD) will automatically be included in the calculated concentrations or deposition flux values if the dry and/or wet deposition processes are considered, unless the user specifies the NODRYDPLT and/or NOWETDPLT options. Note that dry and wet removal effects on calculated concentration values can be included even if deposition flux values are not being calculated. However, the additional data requirements for dry and wet deposition, described in Section 2.2, must be met in order for dry and wet removal to be included in the concentration calculations. The use of the NODRYDPLT and/or NOWETDPLT options will result in a more conservative estimate of concentrations and/or deposition fluxes for applications involving deposition processes, but the degree of additional conservatism will vary depending on the source characteristics, meteorological conditions, receptor locations and terrain influences. However, the inclusion of particle deposition effects may increase ground-level concentrations for some sources compared to the same source modeled as a gaseous emission, due to the effect of gravitational settling on the particulate plume. The magnitude of this effect will depend on the source characteristics (elevated or low-level) and particle size distribution. The PVMRM and OLM options for modeling NO2 conversion are non-DFAULT options, and only one of these options can be specified for a given model run. Both options require that the pollutant ID be specified as 'NO2' on the CO POLLUTE) card (see Section 3.2.5 of the AERMOD User's Guide). These options have additional input requirements as described in Section 2.4. Beginning with version 09260, the TOXICS option in no longer used in AERMOD and the FASTAREA option on the MODELOPT is now used to select the non-DFAULT option to optimize model runtime for AREA sources (including AREA, AREAPOLY, AREACIRC and OPENPIT source types). When the FASTAREA option is specified, the area source integration routine is optimized to reduce model runtime by incorporation of a three-tiered approach using the Romberg numerical integration, a 2-point Gaussian Quadrature routine for numerical integration, or a point source approximation, depending on the location of the receptor relative to the source. In the regulatory default mode the Romberg numerical integration is utilized for all ------- receptors. Also beginning with version 09260, a non-DFAULT option to optimize model runtime for POINT and VOLUME sources has been included, which is selected by the FASTALL option on the MODELOPT keyword. Specification of the FASTALL option also activates the FASTAREA option if AREA sources are including in the model inputs. Both FASTALL and FASTAREA skip receptors that are more than 80 kilometers from the source. The FASTALL option for POINT and VOLUME sources uses an alternative implementation of the horizontal meander algorithm based on an effective horizontal dispersion coefficient (oyeff) that replicates the centerline concentration based on the full meander approach. Use of the effective oy allows the model runtime to be optimized by skipping receptors that are more than 4ayeff off the plume centerline. Based on tests conducted to date, comparisons of concentrations based on the FASTALL option for POINT and VOLUME sources with concentrations based on the DFAULT option are similar to comparisons of concentrations for AREA sources using the FASTAREA option. The average ratio of FASTALL concentrations to DFAULT values is about 1.02 for high ranked values, showing a slight bias toward overprediction for the FASTALL option. However, the range of ratios for high ranked values shows both overpredictions and underpredictions relative the DFAULT option, and differences at specific receptors may be much larger. 2.1.2 BETA Test Options A 'BETA' test option switch is included on the CO MODELOPT keyword to identify and allow for new features to be added to the model that are still in a draft BETA-test status. The BETA option is a non-DFAULT option, and will result in a fatal error if the DFAULT option is also specified. Two draft enhancements are included in AERMOD under the BETA option, beginning with version 06341: 1) Options for capped stacks (source type = POINTCAP) and for horizontal releases (source type = POINTHOR); and 2) The PSDCREDIT option for PVMRM to account for NO/NO2 chemistry of combined plumes in the computation of increment consumption with PSD credits. 10 ------- Inclusion of these draft BETA-test options does not imply any endorsement of their use for regulatory or non-regulatory applications of the model. In addition, the designation of BETA-test to these draft enhancements does not imply that these options have completed rigorous internal ("Alpha") testing prior to being included in a public release of the model. 2.1.3 Processing for Particulate Matter (PM) NAAQS 2.1.3.1 Processing for Fine Particulate Matter (PM-2.5) A NAAQS for fine particulate matter, with aerodynamic particle diameters of 2.5 microns or less (PM-2.5), was promulgated in 1997, and the 24-hour standard was revised in December 2006. For attainment demonstrations, the PM-2.5 standard is based on a 3-year average of the 98th percentile 24-hour average and a 3-year average of the annual mean at each ambient monitor. For purposes of modeling demonstrations of compliance with the NAAQS, the eighth-highest value is an unbiased surrogate for the 98th percentile 24-hour average concentration at a particular receptor over a one-year period. Since the Guideline on Air Quality Models (40 CFR Part 51, Appendix W) prescribes the use of a 5-year data set for off-site National Weather Service (NWS) meteorological data or at least one year of site-specific data, a policy was established to utilize all available meteorological data (both single and multiple years of data) as an unbiased estimate of the 3-year averages for purposes of modeling demonstrations of compliance with the NAAQS. Based on this policy, the 24-hour design value for purposes of modeling demonstrations of compliance with the PM-2.5 NAAQS is based on the highest of the eighth-highest (H8H) concentrations at each receptor, if one year of site-specific meteorological data is input to the model, or the highest of the multi-year average of the eighth-highest concentrations at each receptor, if more than one year of meteorological data is input to the model. In other words, the model calculates the eighth-highest concentration at each receptor for each year modeled, averages those eighth-highest concentrations at each receptor across the number of years of 11 ------- meteorological data, and then selects the highest, across all receptors, of the N-year averaged eighth-highest values. Similar to the 24-hour averages, an unbiased estimate of the 3-year average annual mean is simply the annual mean, if only one year of site-specific meteorological data is input to the model, or the multi-year average of the annual means if multiple years of meteorological data are used. The annual design value for PM-2.5 is then based on the highest annual average across the receptor domain for single-year meteorological data input, or the highest of the multi-year averaged annual means across the receptor domain for multi-year meteorological data input. The special processing of the 24-hour and annual averages for the PM-2.5 NAAQS is triggered by specifying a pollutant ID of 'PM25' or 'PM-2.5' on the CO POLLUTE) card. In this case, the model will compute the 24-hour and annual average design values as described in the previous paragraphs. In order for the PM-2.5 processing to work correctly for multiple year periods, the yearly meteorological data files must be concatenated into a single multi-year file for input into the model. Multi-year meteorological data files can also be generated by processing multi-year inputs in AERMET, the meteorological processor for AERMOD. There is no requirement to remove the header records between concatenated surface meteorological data files prior to running the model. Processing the average of the individual annual mean values across multiple years for PM-2.5 also requires use of the ANNUAL average option on the AVERTEVIE keyword, rather than PERIOD average. The PERIOD option computes a single multi-year average concentration for each receptor, which may give slightly different results than the multi-year average of individual ANNUAL mean concentrations due to differences in the number of calms and/or missing data from year to year. In order to comply with these processing requirements, the following restrictions which are applied to the PM-2.5 NAAQS processing whenever a pollutant ID of 'PM25' or 'PM-2.5' is specified: 1. The averaging periods on the AVERTIME keyword are limited to the 24-hour and ANNUAL averages. Use of the PERIOD average or use of a short-term average other than 24-hour will result in a fatal error message being generated. 12 ------- 2. Only the EIGHTH (or 8TH) highest value may be requested on the RECTABLE keyword for 24-hour averages. Specifying another high value on the RECTABLE card will result in a fatal error message being generated. 3. The model will only process meteorological data for periods of record that span complete years, although the meteorological data period does not need to follow calendar years (i.e., the data period does not need to start on January 1, hour 1). If the period of record spans less than one complete year of data, a fatal error message will be generated and the model run will be unsuccessful. If additional meteorological data remains after the end of the last complete year of data, the remaining data will be ignored, and a non-fatal warning message will be generated specifying the number of hours ignored. 4. The MULTYEAR card on the CO pathway cannot be used to calculate multi-year averages for the PM-2.5 NAAQS. Multiple year analyses must be accomplished by including the multiple years of meteorology in a single data file. 5. Since the 24-hour average design values for PM-2.5 analyses, based on the H8H averaged over N years, may consist of averages over a multi-year period, they are incompatible with the EVENT processor. If the MAXIFILE option is used to output 24- hour average threshold violations, these may be used with the EVENT processor. Therefore, if the EVENTFIL option is used without the MAXIFILE option for PM-2.5 analyses, a non-fatal warning message will be generated, and the EVENTFIL option will be ignored. 2.1.3.2 Processing for Particulate Matter of 10 Microns or Less (PM-10) The 24-hour NAAQS for parti culate matter with aerodynamic particle diameters of 10 microns or less (PM-10) is in the form of an expected exceedance value, which cannot be exceeded more than once per year on average over a three year period for purposes of attainment demonstrations. Modeling demonstrations of compliance with the PM-10 NAAQS are based on the High-N+1-High value over N years, or in the case of five years of NWS meteorological data, the High-6th-High (H6H) over five years. In the AERMOD model, the H6H 24-hour average over five years can be modeled in one of two ways: 1) running five individual years and combining the results using the CO MULTYEAR option, as described in Section 3.2.12 of the AERMOD User's Guide; or 2) using a single five-year meteorological data file and specifying the SIXTH (or 6TH) highest value on the OU RECTABLE card. If applied properly, the results of these two approaches will be equivalent. The special processing consisting of the 99th percentile 24-hour value averaged over N years for PM-10 in previous 13 ------- versions of AERMOD, referred to as the "Post-1997" PM-10 option, has been removed since that standard was vacated. 2.1.4 Specifying Multiple Urban Areas The AERMOD model (beginning with the version dated 06341) includes the option to specify multiple urban areas within the same model run. This option may be applicable for large domains that encompass more than one identifiable urban area where the separation is large enough to warrant separate treatment of the urban boundary layer effects. Use of the option for multiple urban areas eliminates the need for post-processing for such applications. The multiple urban areas are defined using multiple CO URBANOPT cards. The syntax of the modified URBANOPT keyword is as follows: Syntax: For Multiple Urban Areas: CO URBANOPT UrbanID UrbPop (UrbName) (UrbRoughness) For Single Urban Areas: CO URBANOPT UrbPop (UrbName) (UrbRoughness) Type: Optional, Repeatable for multiple urban areas where the UrbanID parameter is the alphanumeric urban ID defined by the user (up to eight characters) when multiple urban areas are defined, the UrbPop parameter specifies the population of the urban area, the optional UrbName parameter may be used to identify the name of the urban area, and the optional UrbRoughness parameter may be used to specify the urban surface roughness length. Note the UrbName must be specified if the user wants to specify the urban roughness length. A default value of 1.0 meter will be used for the urban roughness length if the UrbRoughness parameter is omitted. Beginning with version 09260, any value for the urban roughness length other than 1.0 meter will be treated as a non-DFAULT option. Caution should be used when specifying a non-default urban roughness length, and use of a non-default value should be clearly documented and justified. Note that the syntax of the URBANOPT keyword for single urban areas has not changed from the previous version of AERMOD, so that existing input files will not require modification. 14 ------- The syntax of the URBANSRC keyword on the SO pathway has also been modified to allow for the option of specifying multiple urban areas. The syntax of the modified URBANSRC keyword is as follows: Syntax: For Multiple Urban Areas: SO URBANSRC UrbanID SrcID's and/or SrcRng's For Single Urban Areas: SO URBANSRC SrcID's and/or SrcRng's Type: Optional, Repeatable Order: Must follow the LOCATION card for each source input where the UrbanID parameter is the alphanumeric urban ID (up to eight characters) defined by the user on the CO URBANOPT keyword when multiple urban areas are defined, and the SrcID's and SrcRng's are the individual source IDs and/or source ID ranges that are to be modeled with urban effects. Source ranges are described in more detail in Section 3.3.3 of the AERMOD User's Guide (EPA, 2004a). As with the URB ANOPT keyword, the syntax of the URBANSRC keyword for applications with single urban areas has not changed from the previous version of AERMOD, so that existing input files will not require modification. 2.1.5 BETA Options for Capped and Horizontal Stack Releases Draft BETA test options are included in AERMOD (beginning with the version dated 06341) for modeling releases from capped and horizontal stacks. For sources that are not subject to building downwash influences, the plume rise for these capped and horizontal stacks is simulated based on an EPA Model Clearinghouse Memorandum, dated July 9, 1993, included in Appendix A of this Addendum. The Model Clearinghouse procedure for these sources entails setting the exit velocity very low (0.001 m/s) to account for suppression of vertical momentum for the plume, and using an effective stack diameter that maintains the actual flow rate of the plume. Maintaining the flow rate will also serve to maintain the buoyancy of the plume in order 15 ------- to provide a more realistic estimate of plume rise. The Model Clearinghouse procedure also addresses the issue of stack-tip downwash for these cases. The Model Clearinghouse procedure is not considered appropriate for sources subject to building downwash influences with the PRIME downwash algorithm for the following reason. The PRIME algorithm uses the specified stack diameter to define the initial radius of the plume for the numerical plume rise calculation; use of an effective diameter adjusted to maintain flow rate is not appropriate and could produce unrealistic results. For PRIME downwash sources modeled using the BETA options for capped and horizontal releases, the basic premise of the Model Clearinghouse procedure, i.e. that the vertical momentum is suppressed while the buoyancy of the plume is conserved, have been adapted for the PRIME numerical plume rise formulation. However, this adaptation of the Model Clearinghouse procedure to PRIME downwash sources has not been validated by field tracer or wind tunnel data. The user selects the BETA options for capped and/or horizontal releases by specifying one of the new source types on the SO LOCATION card: POINTCAP for capped stacks, and POINTHOR for horizontal releases. For each of these options, the user specifies the actual stack parameters [release height (m), exit temperature (K), exit velocity (m/s), and stack diameter (m)] using the SO SRCPARAM card as if the release were a non-capped vertical point source. The syntax of the SO LOCATION and SRCPARAM keywords is described in Sections 3.3.1 and 3.3.2, respectively, of the AERMOD User's Guide (EPA, 2004a). The AERMOD model performs the necessary adjustments internally to account for plume rise and stack-tip downwash. For horizontal releases, the model currently assumes that the release is oriented with the wind direction, and the model does not account for directional effects that may occur with horizontal releases. The model also does not account for stacks oriented at a non-horizontal angle relative to vertical. For PREVIE-downwashed sources, the user-specified exit velocity for horizontal releases is treated initially as horizontal momentum in the downwind direction. More details regarding the BETA options for capped and horizontal releases will be provided later, as warranted based on further testing and evaluation. 16 ------- 2.2 DEPOSITION ALGORITHM INPUTS AND OPTIONS The AERMOD model includes algorithms for both dry and wet deposition of both particulate and gaseous emissions. Based on the guidance provided for application of the AERMOD model in Appendix W, and based on the history of deposition algorithms in the AERMOD and ISC models, the particle deposition algorithms with a user-specified particle size distribution (referred to below as "Method 1") can be applied under the regulatory default option. This option is comparable to the particle deposition algorithm in the ISCST3 model (EPA, 1995a). The gas deposition algorithms and the "Method 2" option for particle deposition based on the ANL draft report (Wesely, et al, 2002) are considered to be non-DFAULT options in AERMOD, and the model will issue a fatal error message and abort processing if the DFAULT option is specified with the gas deposition or Method 2 particle deposition options. As discussed above in Section 1.2, the TOXICS option formerly associated with the options for gas deposition and Method 2 for particle deposition is no longer used in the AERMOD model. No additional option switches are required to allow use of these non-DFAULT options. Table 1-1 in Section 1.2 summarizes the required keywords and regulatory status of various deposition options within AERMOD. The remainder of this section provides a detailed description of the model input parameters associated with the deposition algorithms. 2.2.1 Definition of Seasons for Gas Dry Deposition The gas deposition algorithms in AERMOD include land use characteristics and some gas deposition resistance terms based on five seasonal categories, defined in Table 2 of the ANL report as: Seasonal Category 1: Midsummer with lush vegetation Seasonal Category 2: Autumn with unharvested cropland Seasonal Category 3: Late autumn after frost and harvest, or winter with no snow Seasonal Category 4: Winter with snow on ground (with generally continuous snow cover) Seasonal Category 5: Transitional spring with partial green coverage or short annuals 17 ------- The user correlates these seasonal definitions to calendar months through the GDSEASON keyword on the CO pathway. The syntax and type of the GDSEASON keyword are: Syntax: CO GDSEASON Jan Feb Mar ... Dec Type: Optional, Non-repeatable where a numeric value from 1 to 5 is entered for each of the twelve calendar months to associate it with the seasonal definitions given above. This keyword is optional for the model, but mandatory when applying the gas deposition algorithms, unless the GASDEPVD option for user- specified dry deposition velocity on the CO pathway is used, described below in Section 2.2.5. Note that some of the seasonal categories defined above may not apply for certain regions, such as Category 4, winter with continuous snow cover, for moderate climates. 2.2.2 Definition of Land Use Categories for Gas Dry Deposition The gas deposition algorithms include some gas deposition resistance terms based on five seasonal categories, defined above, and nine land use categories as follows (from Table 1 of the ANL report): Land Use Category Description 1 Urban land, no vegetation 2 Agricultural land 3 Rangeland 4 Forest 5 Suburban areas, grassy 6 Suburban areas, forested 7 Bodies of water 8 Barren land, mostly desert 9 Non-forested wetlands The user defines the land use categories by direction sector through the GDLANUSE keyword on the CO pathway. The syntax and type of the GDLANUSE keyword are: 18 ------- Syntax: CO GDLANUSE Seel Sec2 Sec3 ... Sec36 Type: Optional, Non-repeatable where a numeric value from 1 to 9 is entered for each of the 36 direction sectors (every 10 degrees) to associate it with the land use definitions given above. This keyword is optional for the model, but mandatory when applying the gas deposition algorithms, unless the GASDEPVD option for user-specified deposition velocity is used. The first value, Seel, corresponds with the land use category, downwind of the application site, for winds blowing toward 10 degrees, plus or minus 5 degrees. The downwind sectors are defined in clockwise order, with Sec36 corresponding to winds blowing toward 360 degrees (North), and should generally reflect conditions downwind relative to the source location. The user can specify "repeat values" by entering a field such as "36*3" as a parameter for the GDLANUSE keyword. The model will interpret this as "36 separate entries, each with a value of 3." Since the model must identify this as a single parameter field, there must not be any spaces between the repeat-value and the value to be repeated. 2.2.3 Option for Overriding Default Parameters for Gas Dry Deposition An optional keyword is available on the Control (CO) pathway to allow the user to override the default values of the reactivity factor (f0\ and the fraction (F) of maximum green leaf area index (LAI) for seasonal categories 2 (autumn/unharvested cropland) and 5 (transitional spring), for use with the gas dry deposition algorithms. The syntax and type of the GASDEPDF keyword are summarized below: Syntax: CO GASDEPDF React F_Seas2 F_Seas5 (Refpoll) Type: Optional, Non-repeatable where the parameter React is the value for pollutant reactivity factor (f0), and F_Seas2 and F_Seas5 are the fractions (F) of maximum green LAI for seasonal categories 2 and 5, respectively. The parameter Refpoll is the optional name of the pollutant. If the optional GASDEPDF keyword is omitted, then the default value of 0 is used for React, and default values 19 ------- of 0.5 and 0.25 are used for F_Seas2 and F_Seas5, respectively. A value of F=1.0 is used for seasonal categories 1, 3, and 4. A reactivity factor value of 1 should be input for ozone (O3), titanium tetrachloride (TiCU), and divalent mercury (Hg2+), and a value of 0.1 should be input for nitrogen dioxide (NC^). 2.2.4 Specifying Source Parameters for Gas Deposition (Dry and/or Wet) The input of source parameters for dry and wet deposition of gaseous pollutants is controlled by the GASDEPOS keyword on the SO pathway. The gas deposition variables may be input for a single source, or may be applied to a range of sources. The syntax, type, and order for the GASDEPOS keyword are summarized below: Syntax: SO GASDEPOS Srcid (or Srcrng) Da Dw rcl Henry Type: Optional, Repeatable Order: Must follow the LOCATION card for each source input where the Srcid or Srcrng identify the source or sources for which the inputs apply, the parameter Da is the diffusivity in air for the pollutant being modeled (cm2/s), Dw is the diffusivity in water for the pollutant being modeled (cm2/s), rcl is the cuticular resistance to uptake by lipids for individual leaves (s/cm), and Henry is the Henry's Law constant (Pa m3/mol). Values of the physical parameters for several common pollutants may be found in the appendices to the ANL report (Wesely, et. al, 2002). 2.2.5 Option for Specifying the Deposition Velocity for Gas Dry Deposition An optional keyword is available on the Control (CO) pathway to allow the user to specify the dry deposition velocity for gaseous emissions. A single dry deposition velocity can be input for a given model run, and is used for all sources of gaseous pollutants. Selection of this option will by-pass the algorithm for computing deposition velocities for gaseous pollutants, and should only be used when sufficient data to run the algorithm are not available. Results of the AERMOD model based on a user-specified deposition velocity should be used with extra 20 ------- caution. The syntax and type of the GASDEPVD keyword are summarized below: Syntax: CO GASDEPVD Uservd Type: Optional, Non-repeatable where the parameter Uservd is the gaseous dry deposition velocity (m/s). A non-fatal warning message is generated by the model if a value of Uservd greater than 0.05 m/s (5 cm/s) is input by the user. When the GASDEPVD keyword is used, the GDSEASON, GDLANUSE, and GASDEPRF keywords for the CO pathway, and the GASDEPOS keyword for the SO pathway, are no longer applicable and cannot be used in the same model run. As a result, gas wet deposition processes (DEPOS, WDEP, and WETDPLT) cannot be simulated with the GASDEPVD option is used. 2.2.6 Specifying Source Parameters for Particle Deposition The AERMOD model includes two methods for handling dry and/or wet deposition of particulate emissions. Method 1 is used when a significant fraction (greater than about 10 percent) of the total paniculate mass has a diameter of 10 um or larger, or when the particle size distribution is known. The particle size distribution must be known reasonably well in order to use Method 1. Method 2 may be used when the particle size distribution is not well known and when a small fraction (less than 10 percent of the mass) is in particles with a diameter of 10 um or larger. The deposition velocity for Method 2 is calculated as the weighted average of the deposition velocity for particles in the fine mode (i.e., less than 2.5 um in diameter) and the deposition velocity for the coarse mode (i.e., greater than 2.5 um but less than 10 um in diameter). As described in Sections 1.2 and 2.2, use of the Method 2 option is considered non- DFAULT. 21 ------- 2.2.6.1 Specifying Particle Inputs for Method 1 The input of source variables for particle deposition using Method 1 is controlled by three keywords on the SO pathway, PARTDIAM, MASSFRAX, and PARTDENS. These inputs are comparable to the particulate inputs used in the ISCST3 model (EPA,1995a). The particle variables may be input for a single source, or may be applied to a range of sources. The syntax, type and order for these three keywords are summarized below: Syntax: SO PARTDIAM Srcid (or Srcrng) Pdiam(i), i=l,Npd SO MASSFRAX Srcid (or Srcrng) Phi(i), i=l,Npd SO PARTDENS Srcid (or Srcrng) Pdens(i), i=l,Npd Type: Optional, Repeatable Order: Must follow the LOCATION card for each source input where the Srcid or Srcrng identify the source or sources for which the inputs apply, and where the Pdiam array consists of the mass-mean aerodynamic particle diameter (microns) for each of the particle size categories, the Phi array is the corresponding mass fractions (between 0 and 1) for each of the categories, and the Pdens array is the corresponding particle density (g/cm3) for each of the categories. The number of particle size categories for a particular source is Npd. The user does not explicitly tell the model the number of categories being input, but if continuation cards are used to specify particle size variables, all inputs of a keyword for a particular source or source range must be contiguous, and the number of categories must agree for each of the three keywords input for a particular source. As many continuation cards as needed may be used to define the inputs for a particular keyword. The model checks the inputs to ensure that the mass fractions sum to 1.0 (within 2 percent) for each source input, and issues a warning message if that range is exceeded. The model also ensures that mass fractions for each particle size category are within the proper range (between 0 and 1), and issues fatal error messages for any value exceeded that range. 22 ------- 2.2.6.2 Specifying Particle Inputs for Method 2 The Method 2 particle information is input through the METHOD_2 keyword on the SO pathway. The syntax, type, and order for the METHOD_2 keyword are summarized below: Syntax: SO METHOD_2 Srcid (or Srcrng) FineMassFraction Dmm Type: Optional, Repeatable Order: Must follow the LOCATION card for each source input where the Srcid or Srcrng identify the source or sources for which the inputs apply, the parameter FineMassFraction is the fraction of particle mass emitted in the fine mode, less than 2.5 microns, and Dmm is the representative mass-mean aerodynamic particle diameter in microns. Estimated values of fine particle fractions and mass mean diameters for various pollutants are provided in Appendix B of the ANL report (Wesely, et al, 2002). 2.2.7 Specifying Emission and Output Units Since the AERMOD model allows for both concentration and deposition to be output in the same model run, the EMISUNIT keyword (see Section 3.3.6 of the AERMOD User's Guide (EPA, 2004a)) cannot be used to specify emission unit factors if more than one output type is being generated. The AERMOD model therefore allows for concentration and deposition units to be specified separately through the CONCUNIT and DEPOUNIT keywords, respectively. The syntax and type of the CONCUNIT keyword are summarized below: Syntax: SO CONCUNIT Emifac Emilbl Conlbl Type: Optional, Non-repeatable where the parameter Emifac is the factor to convert emission rate input units to the desired output units, Emilbl is the label for the emission input units (up to 40 characters), and Conlbl is the output unit label (up to 40 characters) for concentration calculations. The syntax and type of the DEPOUNIT keyword are summarized below: 23 ------- Syntax: SO DEPOUNIT Emifac Emilbl Deplbl Type: Optional, Non-repeatable where the parameter Emifac is the factor to convert emission rate input units to the desired output units, Emilbl is the label for the emission input units (up to 40 characters), and Deplbl is the output unit label (up to 40 characters) for deposition calculations. 2.2.8 Deposition Velocity and Resistance Outputs In order to facilitate review and testing of the deposition algorithms in the AERMOD model, the model includes an option to output the main resistance terms and deposition velocities for gaseous and particle sources. These optional outputs are generated if the user specifies the 'CO DEBUGOPT MODEL' option described in Section 3.2.12 of the AERMOD User's Guide (EPA, 2004a). The gas deposition data are written to a file called GDEP.DAT, which includes the values of Ra, Rb, Rc, and Vdg (see Wesely, et al, 2002, for definitions) for each source and for each hour modeled. A header record is included to identify the columns. The particle deposition data are written to a file called PDEP.DAT, which includes the values of Ra, Rp, Vg, and Vd for each source and for each hour modeled. The particle outputs are labeled as being based on either Method 1 or Method 2. For Method 1, results are output for each particle size category. The filename and file units for these data files are hardcoded in the model, and the files are overwritten each time the model is executed. Since these files include data for each source for each hour, file sizes may become large. 24 ------- 2.2.9 Meteorological Data for Deposition Algorithms The AERMET meteorological processor was modified (beginning with the version dated 04300) to output additional meteorological parameters needed for the deposition algorithms in AERMOD. The additional variables include the precipitation code, precipitation rate, relative humidity, surface pressure, and cloud cover. These additional variables are automatically included after the standard variables for each hour, and do not require any additional user input. The precipitation data needed for wet deposition calculations in AERMOD can be obtained from the SAMSON, HUSWO or ISHD (TD-3505) formats currently supported by AERMET (EPA, 2004b). The input meteorological data file consists of a header record that includes the latitude and longitude, surface station ID (e.g., WBAN number), upper air station ID, the on-site station ID, and the AERMET version date (see Section D.I of the AERMOD model user's guide (EPA, 2004a)). The meteorological data file for the deposition algorithms is read as a free format file, i.e., each field on a record is separated from adjacent fields by a comma or by one or more spaces. The subsequent data records contain the following variables in the order listed: Variable Description Format Year Integer Month Integer Day of Month Integer Julian Day (Day of Year) Integer Hour of Day Integer Heat Flux (W/m2) Real Surface Friction Velocity, u* (m/s) Real Convective Velocity Scale, w* (m/s) Real Lapse Rate Above Mixing Height (K/m) Real Convective Mixing Height (m) Real Mechanical Mixing Height (m) Real Monin-Obukhov Length, L (m) Real Surface Roughness Length, z0 (m) Real Bowen Ratio Real Albedo Real Reference Wind Speed (m/s) Real Reference Wind Direction (degrees) Real Reference Height for Wind (m) Real 25 ------- Ambient Temperature (K) Real Reference Height for Temperature (m) Real Precipitation Code (0-45) Integer Precipitation Amount (mm) Real Relative Humidity (%) Real Surface Pressure (mb) Real Cloud Cover (tenths) Integer 2.3 OPEN PIT SOURCE OPTION The open pit source option is invoked by specifying a source type of OPENPIT on the source location (SO LOCATION) card. The OPENPIT source algorithm can be used to model particulate or gaseous emissions from open pits, such as surface coal mines and rock quarries. If particulate emissions modeled, the user must specify the particle size information through the appropriate keywords described in Section 2.2. The OPENPIT algorithm uses an effective area for modeling pit emissions, based on meteorological conditions, and then utilizes the numerical integration area source algorithm to model the impact of emissions from the effective area sources. A complete technical description of the OPENPIT source algorithm is provided in the ISC3 Model User's Guide - Volume II (EPA, 1995b). The AERMOD model accepts rectangular pits with an optional rotation angle specified relative to a north-south orientation. The rotation angle is specified relative to the vertex used to define the source location on the SO LOCATION card (e.g., the southwest corner). The syntax, type and order for the SRCPARAM card for OPENPIT sources are summarized below: Syntax: SO SRCPARAM Srcid Opemis Relhgt Xinit Yinit Pitvol (Angle) Type: Optional, Repeatable Order: Must follow the LOCATION card for each source input where the Srcid parameter is the same source ID that was entered on the LOCATION card for a particular source, and the other parameters are as follows: Opemis - open pit emission rate in g/(s-m2), 26 ------- Relhgt - average release height above the base of the pit in meters, Xinit - length of X side of the open pit (in the east-west direction if Angle is 0 degrees) in meters, Yinit - length of Y side of the open pit (in the north-south direction if Angle is 0 degrees) in meters, Pitvol - volume of open pit in cubic meters, and Angle - orientation angle for the rectangular open pit in degrees from North, measured positive in the clockwise direction (optional). The same emission rate is used for both concentration and deposition calculations in the AERMOD model. It should also be noted that the emission rate for the open pit source is an emission rate per unit area, which is different from the point and volume source emission rates, which are total emissions for the source. The Relhgt parameter cannot exceed the effective depth of the pit, which is calculated by the model based on the length, width and volume of the pit. A Relhgt of 0.0 indicates emissions that are released from the base of the pit. If the optional Angle parameter is input, and the value does not equal 0.0, then the model will rotate the open pit clockwise around the vertex defined on the SO LOCATION card for this source. The relationship between the Xinit, Yinit, and Angle parameters and the source location, (Xs,Ys), for a rotated pit is the same as for rectangular area sources. The Xinit dimension is measured from the side of the area that is counterclockwise along the perimeter from the vertex defined by (Xs,Ys), while the Yinit dimension is measured from the side of the open pit that is clockwise along the perimeter from (Xs,Ys). Unlike the area source inputs, the Yinit parameter is not optional for open pit sources. The Angle parameter is measured as the orientation relative to North of the side that is clockwise from (Xs,Ys), i.e. the side with length Yinit. The Angle parameter may be positive (for clockwise rotation) or negative (for counterclockwise rotation), and a warning message is generated if the absolute value of Angle is greater than 180 degrees. The selection of the vertex to use for the source location is not critical, as long as the relationship described above for the Xinit, Yinit, and Angle parameters is maintained. 27 ------- The aspect ratio (i.e., length/width) of open pit sources should be less than 10 to 1. However, since the pit algorithm generates an effective area for modeling emissions from the pit, and the size, shape and location of the effective area is a function of wind direction, an open pit cannot be subdivided into a series of smaller sources. Aspect ratios of greater than 10 to 1 will be flagged by a warning message in the output file, and processing will continue. Since open pit sources cannot be subdivided, the user should characterize irregularly-shaped pit areas by a rectangular shape of equal area. Receptors should not be located within the boundaries of the pit; concentration and/or deposition at such receptors will be set to zero. Such receptors will be identified during model setup and will be flagged in the summary of inputs. An example of a valid SRCPARAM input card for an open pit source is given below: SOSRCPARAM NORTHPIT 1.15E-4 0.0 150.0 500.0 3.75E+6 30.0 where the source ID is NORTHPIT, the emission rate is 1.15E-4 g/(s-m2), the release height is 0.0 m, the X-dimension is 150.0 m, the Y-dimension is 500.0 m, the pit volume is 3.75E+6 cubic meters (corresponding to an effective pit depth of about 50 meters) and the orientation angle is 30.0 degrees clockwise from North. 2.4 PVMRM AND OLM OPTIONS FOR MODELING NO2 This section provides a description of the inputs related to the non-DFAULT PVMRM and OLM options for modeling the conversion of NOX to NO2. A technical description of the PVMRM algorithm is provided in an Addendum to the AERMOD Model Formulation Document (Cimorelli, et a/., 2004). Background on the original development of the PVMRM option is provided by Hanrahan (1999a and 1999b). The PVMRM and OLM algorithms have been implemented as non-DFAULT options, which means that the PVMRM and OLM options cannot be used if the DFAULT keyword is included on the CO MODELOPT card. As described in Section 2.1.1, a BETA-test draft model option, PSDCREDIT, has been added for use when an application is for increment consumption 28 ------- with PSD credits using PVMRM. The special source grouping required for the PSDCREDIT option is described below in Section 2.4.6.1. 2.4.1 Specifying Ozone Concentrations for PVMRM and OLM Options The background ozone concentrations can be input as a single value through the OZONEVAL keyword on the CO pathway, or may be input as hourly values through a separate data file specified through the OZONEFIL keyword on the CO pathway. The user must specify either the OZONEVAL or OZONEFIL keywords, or both, in order to use the PVMRM or OLM options. If both keywords are entered, then the value entered on the OZONEVAL keyword will be used to substitute for hours with missing ozone data in the ozone data file. The syntax of the OZONEVAL keyword is as follows: Syntax: CO OZONEVAL OSValue (OSUnits) Type: Optional, Non-repeatable where the OSValue parameter is the background ozone concentration in the units specified by the optional OSUnits parameter (PPM, PPB, or UG/M3). If the optional OSUnits parameter is missing, then the model will assume units of micrograms/cubic-meter (UG/M3) for the background ozone values. If units of PPM or PPB are used, then the model will convert the concentrations to micrograms/cubic-meter based on the reference ambient temperature and the base elevation specified on the ME PROFBASE card. The OZONEVAL keyword word is optional and non-repeatable. The syntax of the OZONEFIL keyword is as follows: Syntax: CO OZONEFIL OSFileName (OSUnits) (Format) Type: Optional, Non-repeatable where the OSFileName parameter is the filename for the hourly ozone concentration file, the optional OSUnits parameter specifies the units of the ozone data (PPM, PPB, or UG/M3, with 29 ------- UG/M3 as the default), and the optional Format parameter specifies the Fortran FORMAT to read the ozone data. The OSFileName can be up to 200 characters in length based on the default parameters in AERMOD. Double quotes (") at the beginning and end of the filename can also be used as field delimiters to allow filenames with embedded spaces. If the optional Format parameter is missing, then the model will read the ozone data using a Fortran free format, i.e., assuming that commas or spaces separate the data fields. The contents of the ozone data file must include the year (2-digits), month, day, hour and ozone value in that order (unless specified differently through the Format parameter). The date sequence in the ozone data file must match the date sequence in the hourly meteorological data files. As with the OZONEVAL keyword, if units of PPM or PPB are used, then the model will convert the concentrations to micrograms/cubic-meter based on the reference ambient temperature and the base elevation specified on the ME PROFBASE card. Values of ozone concentrations in the ozone data file that are less than zero or greater than or equal to 900.0 will be regarded as missing. If a background ozone value has been specified using the OZONEVAL keyword, then that value will be used to substitute for missing ozone data from the ozone file. If no OZONEVAL keyword is used, then the model will assume full conversion for hours with missing ozone data. 2.4.2 Specifying the Ambient Equilibrium NO2/NOx Ratio for PVMRM The PVMRM option for modeling conversion of NO to NO2 incorporates a default NO2/NOX ambient equilibrium ratio of 0.90. A NO2/NOX equilibrium ratio other than 0.90 can be specified through the optional NO2EQUIL keyword on the CO pathway. The syntax of the NO2EQUIL keyword is as follows: Syntax: CO NO2EQUIL NO2Equil Type: Optional, Non-repeatable where the NO2Equil parameter is the NO2/NOX equilibrium ratio and must be between 0.10 and 1.0, inclusive. 30 ------- 2.4.3 Specifying the Default In-stack NO2/NOx Ratio for PVMRM and OLM The PVMRM and OLM options for modeling conversion of NO to NO2 require that an in-stack NO2/NOX ratio be specified. A default value of 0.10 for the in-stack NO2/NOX ratio will be used for OLM for all sources unless a user-specified value is provided. However, the PVMRM option requires that the user specify the in-stack NO2/NOX ratio for each source. The in-stack NO2/NOX ratio can be specified for the PVMRM or OLM options by using either the CO NO2STACK card to specify a default value to be used for all sources, or by using the SO NO2RATIO card to specify a value on a source-by-source basis. The SO NO2RATIO card can also be used to override the default value if the CO NO2STACK card has been specified. The syntax of the NO2STACK keyword is as follows: Syntax: CO NO2STACK NO2Ratio Type: Optional, Non-repeatable where the NO2Ratio parameter is the default in-stack NO2/NOX ratio that will be used, unless overridden on a source-by-source basis by the SO NO2RATIO card (described below). The value of NO2Ratio must be between 0.0 and 1.0, inclusive. Users should note that while CO NO2STACK is an optional keyword, the PVMRM option requires the user to specify an in-stack NO2/NOX ratio for each source, using either the CO NO2STACK or SO NO2RATIO cards (described in Section 2.4.4), or both. 2.4.4 Specifying In-stack NO2/NOx Ratios by Source for PVMRM and OLM The PVMRM and OLM options for modeling NO2 conversion assume a default in-stack NO2/NOX ratio of 0.10 (or 10 percent). The user can specify in-stack NO2/NOX ratios through the optional NO2RATIO keyword on the SO pathway. The syntax of the NO2RATIO keyword is as follows: ------- Syntax: SO NO2RATIO SrcID or SrcRange NO2Ratio Type: Optional, Repeatable Order: Must follow the LOCATION card for each source input where the SrcID or SrcRange identify the source or sources for which the inputs apply, and where the NO2Ratio parameter specifies the in-stack ratio. In this way, the user can specify a single in-stack NO2/NOX ratio for a group of stacks. For example, the following input: SONO2RATIO STACK1-STACK10 0.15 will apply the in-stack ratio of 0.15 to sources with IDs falling within the range STACK1 to STACK10. Any value specified on the SO NO2RATIO card will override the default ratio, if any, specified on the CO NO2STACK card. Users should note that while SO NO2RATIO is an optional keyword, the PVMRM option requires the user to specify an in-stack NO2/NOX ratio for each source, using either the CO NO2STACK (described in Section 2.4.3) or SO NO2RATIO cards, or both. 2.4.5 Specifying Combined Plumes for OLM The OLM option for modeling NO2 conversion includes an option for specifying which sources are to be modeled as combined plumes. Sources which are not specified for modeling as combined plumes will be modeled as individual plumes. The selection of individual or combined plume option for OLM is specified through the OLMGROUP keyword on the SO pathway. The syntax of the OLMGROUP card is as follows: Syntax: Type: Order: SO OLMGROUP OLMGrpID SrcID's and/or SrcRange's Optional, Repeatable Must follow the LOCATION card for each source input where OLMGrpID identifies a group to be treated as a combined plume with OLM, and the SrcID's and/or SrcRange's identify the sources to be included in the OLM group. As with the SO 32 ------- SRCGROUP card, individual source IDs and source ranges may be used on the same card, and if more than one input card is needed to define the sources for a particular OLM group, then additional cards may be input by repeating the pathway, keyword and OLM group ID. A user can also specify an OLMGrpID of ALL, which means that OLM will be applied on a combined plume basis to all sources. However, unlike the SO SRCGROUP card, the results will not be output for an OLM group unless the same group of sources is also identified on a SRCGROUP card. Another constraint is that a source cannot be included in more than one OLM group. If a source is not selected for an OLMGROUP card, then OLM will be applied to that source as an individual plume. Other than the similarity in syntax, there is no connection in the model between the groups defined on the OLMGROUP card and groups defined on the SRCGROUP card. The OLMGROUP card relates to how the results are processed within the model for the OLM model, and the SRCGROUP card simply controls how source impacts are grouped in the model outputs. If the user identifies one or more groups of sources to apply OLM on a combined plume basis using the OLMGROUP card, the model will still need to calculate the concentration for individual plumes within the OLM group in order for the model to sum the results for the sources listed on the SRCGROUP card(s). The individual source concentrations are calculated by applying the ratio of the combined concentration for the OLM group with and without OLM to each source within the OLM group. 33 ------- 2.4.6 Modeling NO2 Increment Credits with PVMRM NO Due to the ozone-limiting effects of the PVMRM option, the predicted concentrations of are not linearly proportional to the emission rate. Therefore, the approach of modeling NC>2 increment consumption with PSD credits through the use of a negative emission rate for credit sources cannot be used with the PVMRM option. However, the draft PSDCREDIT option allows modeling PSD increment credits for NC>2 when the PVMRM option is specified. The PSDCREDIT option is currently implemented as a BETA-test option, and requires that the PVMRM and BETA options be specified. The PSDCREDIT option utilizes a new PSDGROUP keyword, described below, to identify which sources consume or expand increment. This option is not valid if the OLM option is specified, and no comparable option is available for modeling increment credits with the OLM option. The user should check with the appropriate reviewing authority for further guidance on modeling increment credits for NO2. A general discussion of concepts related to modeling increment consumption is provided below, followed by a description of inputs required to use the BETA-test PSDCREDIT option for PVMRM. 2.4.6.1 Increment Consuming and Baseline Sources Increment is the maximum allowable increase in concentration of a pollutant above a baseline concentration for an area defined under the Prevention of Significant Deterioration (PSD) regulations. The PSD baseline area can be an entire State or a subregion of a State such as a county or group of counties. Increment standards exist for three pollutants: SO2 (3-hr, 24-hr, and annual averages), NO2 (annual average), and PM-10 (24-hr and annual average). Increment consumption is the additional air quality impact above a baseline concentration. The baseline concentration is the ambient concentration of the pollutant that existed in the area at the time of the submittal of the first complete permit application by any source in that area subject to PSD regulations. A baseline source is any source that existed prior to that first application and the baseline date is the date of the PSD application. This baseline date is referred 34 ------- to as the minor source baseline date in PSD regulations. By definition, baseline sources do not consume increment. However, any baseline source that retires from service after the baseline date expands the increment available to new sources. Therefore, a PSD modeling analysis performed for a new source may need to account for this increment expansion. Such an analysis may therefore involve identification of three groups of sources: 1) increment-consuming sources; 2) retired (increment-expanding) baseline sources; and 3) existing, non-retired, baseline sources. 2.4.6.2 Calculating Increment Consumption under the PSDCREDIT Option Calculating increment consumption under the PSDCREDIT option in AERMOD is not a simple arithmetic exercise involving the three groups of sources defined above. Since the amount of ozone available in the atmosphere limits the conversion of NO to NO2, interactions of plumes from the existing and retired baseline sources with those from the increment consuming sources must be considered as part of the calculation of net increment consumption. Without the PSDCREDIT option, properly accounting for the potential interaction of plumes among the different source categories would require post-processing of results from multiple model runs. Internal "posf'-processing algorithms have been incorporated in AERMOD under the PSDCREDIT option to account for the apportioning of the three groups of sources to properly calculate increment consumption from a single model run. Define the following three source groupings for the discussion that follows: A = increment-consuming sources; B = non-retired baseline sources; and C = retired baseline, increment-expanding sources. The calculation of the amount of increment consumption by the A sources cannot simply be estimated by modeling the A sources alone because of the possible interaction of those plumes with the plumes from B sources. The PVMRM algorithm is designed to account for such plume interactions and calculate the total NO to NO2 conversion in the combined plumes based on the amount of ozone available. Therefore, the total increment consumption by the A sources is 35 ------- given by the difference between (1) the total future impact of increment consuming sources and non-retired baseline sources (A+B) and (2) the total current impact (B), which can be expressed as (A+B) - (B). Here (A+B) represents the value that would be compared against the National Ambient Air Quality Standard (NAAQS) for NC>2 during PSD review of the A sources. In a case where some of the baseline sources have been retired from service (C sources), the PSD regulations allow the consideration of increment expansion when assessing compliance with the PSD increment. However, the amount of increment expansion cannot be estimated by simply modeling the C sources alone because of the possible interaction of those plumes with the plumes from B sources. Therefore, the total increment expansion, i.e., PSD credit, is calculated as the difference between (1) the total impact prior to the retirement of C sources, i.e. (B+C), and (2) the total impact from existing (non-retired) baseline sources (B), which can be expressed as (B+C) - (B). Finally, the net increment consumption is given by the difference between total increment consumption and the total increment expansion, or [(A+B)-(B)]-[(B+C)-(B)] (1) Note that in the absence of any increment expansion, the net increment consumption is equal to the total increment consumption [(A+B) - (B)], as described above. These expressions of net increment consumption and expansion cannot be interpreted as algebraic equations. Instead, the terms within parentheses represent the results of separate model runs that account for the combined effects of NOX conversion chemistry on specific groups of sources. The expression shown in Equation 1 above represents four model simulations: (A+B), (B), (B+C), and (B) again. In this case, the two (B) terms do cancel each other and we are left with: [(A+B)] - [(B+C)] (2) 36 ------- The expression presented in Equation 2 summarizes how the net increment consumption calculation is performed under the PSDCREDIT option. Under this option, AERMOD first models the A and B groups together, then models the B and C groups together, and finally computes the difference to obtain the desired result, i.e., the value to compare to the PSD increment standard. In order for AERMOD to perform the special processing associated with this option, the user must define which sources belong to each of the groupings defined above. The next section describes how this is accomplished. 2.4.6.3 Specifying Source Groups under the PSDCREDIT Option The PSDCREDIT option introduces limitations on grouping sources in order to calculate increment consumption as described in the previous section. A new keyword, PSDGROUP, is used to group the sources to correctly calculate the increment consumption. The syntax, type, and order are similar to the regular SRCGROUP keyword and are summarized below: Syntax: SO PSDGROUP Grpid Srcid's and/or Srcrng's Type: Mandatory for PSDCREDIT option, Repeatable Order: Must follow the last keyword in the SO pathway before FINISHED If the PSDCREDIT model option is specified, the PSDGROUP keyword must be used. The SRCGROUP keyword cannot be used under the PSDCREDIT option since results from other groupings beyond these three do not have any meaning when the PSDCREDIT option is invoked and sources are allocated to the calculation of increment consumption. Special source groups for outputting model results are defined within AERMOD for the PSDCREDIT option, as described in the next section. Only the following special PSD group ID's can be used. Failure to use these group ID's will result in a fatal error message during setup processing by AERMOD. The group ID's are: INCRCONS - increment-consuming sources (group A above); these can be new sources or modifications to existing sources; NONRBASE - existing, non-retired baseline sources (group B above); and 37 ------- RETRBASE - retired (increment-expanding or PSD credit) baseline sources (group C above). It is important to note that the source emission inputs for sources included in the RETRBASE PSD group must be entered as positive numbers, unlike other types of PSD credit modeling where negative emissions are input to simulate the impact of the credit sources on the increment calculation. The increment-expanding contribution from RETRBASE sources is accounted for within the AERMOD model under the PSDCREDIT option. The group ID's can appear in any order, but these are the only three that can be specified. If there are no retired baseline sources (i.e., no baseline sources are retired), the keyword RETRBASE can be omitted. Likewise, if there are no non-retired baseline sources (i.e., all baseline sources have been retired), the NONRBASE keyword can be omitted. The special group ID 'ALL' that can be use with the SRCGROUP keyword cannot be used with the PSDGROUP keyword. As with the SRCGROUP keyword for non-PSDCREDIT applications, the group ID's are repeatable and they must be the last keyword before FINISHED on the SO pathway when the PSDCREDIT option is specified. Source ranges, which are described in more detail in Section 3.3.3 of the AERMOD User's Guide (EPA, 2004a), are input as two source IDs separated by a dash, e.g., STACK1- STACK10. Individual source IDs and source ranges may be used on the same card. If more than one input card is needed to define the sources for a particular group, then additional cards may be input, repeating the pathway, keyword and group ID. A source can appear in only one of these source groups, and must be assigned to one of the groups. The requirements for specifying sources and source groups under the PSDCREDIT option are summarized below: • The SRCGROUP keyword cannot be used with the PSDCREDIT option; • Special PSD group ID's must be used with the PSDGROUP keyword; • The group ID ALL is not allowed when the PSDCREDIT option is specified; • A source must appear in one, and only one, of the PSDGROUPs; and 38 ------- • Emission rates for increment-expanding (RETRBASE) sources must be entered as positive values. 2.4.6.4 Model Outputs under the PSDCREDIT Option Unlike the regular SRCGROUP keyword, the PSDGROUP keyword does not define how the source impacts are grouped for model output. As described in the previous sections, the PSDGROUP keyword defines the different categories of sources needed in order to properly account for NOX conversion chemistry under the PVMRM option. The model outputs under the PSDCREDIT option in AERMOD are based on demonstrating compliance with the air quality standards, i.e., the NAAQS and PSD increment for NO2. As a result, AERMOD uses hardcoded "SRCGROUP" names of 'NAAQS' and 'PSDINC' to label these two types of outputs. The results output under the 'NAAQS' source group label are based on the calculation of (A+B) as described above in Section 2.4.6.2. The results reported under the 'PSDINC' source group label are based on the expression presented above in Equation 2. 2.5 VARIABLE EMISSION RATES 2.5.1 Specifying Variable Emission Factors (EMISFACT) The AERMOD model provides the option of specifying variable emission rate factors for individual sources or for groups of sources. The syntax, type and order of the EMISFACT keyword are summarized below: Syntax: SO EMISFACT SrcID or SrcRange Qflag Qfact(i), i=l,n Type: Optional, Repeatable Order: Must follow the LOCATION card for each source input where the SrcID parameter is the same source ID that was entered on the LOCATION card for a particular source. The user also has the option of using the SrcRange parameter for specifying a 39 ------- range of sources for which the emission rate factors apply, instead of identifying a single source. This is accomplished by two source ID character strings separated by a dash, e.g., STACK1- STACK10. The use of the SrcRange parameter is explained in more detail in the description of the BUILDHGT keyword (see Section 3.3.3 of the AERMOD User's Guide). The parameter Qflag is the variable emission rate flag, and must be specified as one of the following secondary keywords: SEASON - emission rates vary seasonally (n=4), MONTH - emission rates vary monthly (n=12), HROFDY - emission rates vary by hour-of-day (n=24), WSPEED - emission rates vary by wind speed (n=6), SEASHR - emission rates vary by season and hour-of-day (n=96), HRDOW - emission rates vary by hour-of-day, and day-of-week [M-F, Sat, Sun] (n=72), FIRDOW7 - emission rates vary by hour-of-day, and the seven days of the week [M, Tu, W, Th, F, Sat, Sun] (n=168), SHRDOW - emission rates vary by season, hour-of-day, and day-of-week [M-F, Sat, Sun] (n=288), SHRDOW7 - emission rates vary by season, hour-of-day, and the seven days of the week [M, Tu, W, Th, F, Sat, Sun] (n=672), MHRDOW - emission rates vary by month, hour-of-day, and day-of-week [M-F, Sat, Sun] (n=864), and MHRDOW7 - emission rates vary by month, hour-of-day, and the seven days of the week [M, Tu, W, Th, F, Sat, Sun] (n=2,016). The Qfact array is the array of factors, where the number of factors is shown above for each Qflag option. The seasons are defined in the following order: Winter (Dec., Jan., Feb.), Spring (Mar., Apr., May), Summer (Jun., Jul., Aug.), and Fall (Sep., Oct., Nov.). The wind speed categories used with the WSPEED option may be defined using the ME WINDCATS keyword. 40 ------- If the WINDCATS keyword is not used, the default wind speed categories are defined by the upper bound of the first five categories as follows (the sixth category is assumed to have no upper bound): 1.54, 3.09, 5.14, 8.23, and 10.8 m/s. The EMISFACT card may be repeated as many times as necessary to input all of the factors, and repeat values may be used for the numerical inputs. Examples for the more recent HRDOW and MHRDOW options are presented below, with column headers to indicate the order in which values are to be to input: 41 ------- ** Weekdays: Mrs: 1-5 6 7-17 18 19-24 SO EMISFACT STK1 HRDOW 5*0.3 0.5 11*1.0 0.5 6*0.3 ** Saturdays: Mrs: 1-5 6 7-17 18 19-24 SO EMISFACT STK1 HRDOW 5*0.3 0.5 11*1.0 0.5 6*0.3 ** Sundays: Mrs: 1-5 6 7-17 18 19-24 SO EMISFACT STK1 HRDOW I SO EMISFACT STK1 HRDOW7 enter 24 hourly scalars for each of the "days", first for Mondays, then for Tuesdays, ..., then for Saturdays, and finally for Sundays, e.g., ** Mondays: Mrs: 1-5 6 7-17 18 19-24 SO EMISFACT STK1 HRDOW7 5*0.3 0.5 11*1.0 0.5 6*0.3 ** Tuesdays: Mrs: 1-5 6 7-17 18 19-24 SO EMISFACT STK1 HRDOW7 5*0.3 0.5 11*1.0 0.5 6*0.3 ** Saturdays: Mrs: 1-5 6 7-17 18 19-24 SO EMISFACT STK1 HRDOW7 5*0.3 0.5 11*1.0 0.5 6*0.3 ** Sundays: Mrs: 1-5 6 7-17 18 19-24 SO EMISFACT STK1 HRDOW7 5*0.3 0.5 11*1.0 0.5 6*0.3 SO EMISFACT STK1 MHRDOW enter 24 hourly scalars for each of the twelve months, first for Weekdays (Monday-Friday), then for Saturdays, and finally for Sundays, e.g. ** Weekdays JAN FEE MAR APR MAY JUN . . . NOV DEC SO EMISFACT STK1 MHRDOW 24*1.0 24*0.8 24*0.6 24*0.8 24*1.0 24*0.8 24*0.6 24*0.8 ** Saturdays: SO EMISFACT STK1 MHRDOW 24*1.0 24*0.8 24*0.6 24*0.8 24*1.0 24*0.8 24*0.6 24*0.8 * * Sundays: SO EMISFACT STK1 MHRDOW 24*1.0 24*0.8 24*0.6 24*0.8 24*1.0 24*0.8 24*0.6 24*0.8 SO EMISFACT STK1 MHRDOW7 enter 24 hourly scalars for each of the twelve months, first for Mondays, then for Tuesdays, ..., then for Saturdays, and finally for Sundays, e.g., ** Mondays JAN FEE MAR APR MAY JUN . . . NOV DEC SO EMISFACT STK1 MHRDOW7 24*1.0 24*0.8 24*0.6 24*0.8 24*1.0 24*0.8 24*0.6 24*0.8 ** Tuesdays JAN FEE MAR APR MAY JUN . . . NOV DEC SO EMISFACT STK1 MHRDOW7 24*1.0 24*0.8 24*0.6 24*0.8 24*1.0 24*0.8 24*0.6 24*0.8 ** Saturdays: SO EMISFACT STK1 MHRDOW7 24*1.0 24*0.8 24*0.6 24*0.8 24*1.0 24*0.8 24*0.6 24*0.8 ** Sundays: SO EMISFACT STK1 MHRDOW7 24*1.0 24*0.8 24*0.6 24*0.8 24*1.0 24*0.8 24*0.6 24*0.8 42 ------- 2.5.2 Specifying an Hourly Emission Rate File (HOUREMIS) The source (SO) pathway includes an option for inputting hourly emission rates for the AERMOD model, controlled by the HOUREMIS keyword. AERMOD currently allows for a single hourly emission file to be used with each model run. The syntax, type and order for this keyword are summarized below: Syntax: SO HOUREMIS Emifil Srcid's (and/or Srcrng's) Type: Optional, Repeatable Order: Must follow the LOCATION card for each source input where the Emifil parameter specifies the filename for the hourly emission file, and Srcid or Srcrng identify the source or sources for which hourly emission rates are included. The Emifil filename can be up to 200 characters in length based on the default parameters in AERMOD. Double quotes (") at the beginning and end of the filename can also be used as field delimiters to allow filenames with embedded spaces. Source ranges, which are described in more detail in Section 3.3.3 of the AERMOD User's Guide (EPA, 2004a), are input as two source IDs separated by a dash, e.g., STACK1-STACK10. The user may include more than one HOUREMIS card in a runstream file, if needed to specify additional sources, but there can be only one hourly emissions file, and therefore the filename must be the same on all HOUREMIS cards. The format of each record of the hourly emissions file includes a pathway and keyword (SO HOUREMIS), followed by the Year, Month, Day, Hour, Source ID, and emission rate (in the appropriate units). For POINT sources, the stack gas exit temperature (K), and stack gas exit velocity (m/s) are also specified. Beginning with version 09260, the release heights and initial dispersion coefficients can also be varied on an hourly basis for AREA, AREAPOLY, AREACIRC, and VOLUME sources using the HOUREMIS option. The user selects this enhanced option by including the additional source parameters in the hourly emissions file. AERMOD determines whether hourly release heights and initial dispersion coefficients are being used based on the first HOUREMIS record for each source, and these additional 43 ------- parameters must be included on all HOUREMIS records unless the emissions are missing, which is indicated but leaving the emission rate and all fields beyond the source ID blank. The hourly emissions file is processed using the same routines used to process the runstream input file, therefore each of the parameters must be separated by at least one space, but otherwise the format is flexible. It is also not necessary to include the SO HOUREMIS on each line, as long as the parameters (Year, Month, etc.) do not begin before column 13. The data in the hourly emission file must also include the exact same dates as are included in the meteorological input files, and the source IDs must correspond to the source IDs defined on the SO LOCATION cards and be in the same order as defined in the 'aermod.inp' file. The model will check for a date mismatch between the hourly emissions file and the meteorological data, and also for a source ID mismatch. However, it is not necessary to process the entire hourly emissions file on each model run, i.e., the correct emissions data will be read if the ME DAYRANGE or the ME STARTEND cards (see Section 3.5.4 of the AERMOD User's Guide) are used, as long as all the dates (including those that are processed and those that are skipped) match the meteorological data files. An example of several lines from an hourly emissions file for two point sources is provided below: SO HOUREMIS 88 8 16 1 STACK1 52.467 382.604 12.27 SO HOUREMIS 88 8 16 1 STACK2 44.327 432.326 22.17 SO HOUREMIS 88 8 16 2 STACK1 22.321 377.882 9.27 SO HOUREMIS 88 8 16 2 STACK2 42.166 437.682 19.67 SO HOUREMIS 88 8 16 3 STACK1 51.499 373.716 11.87 SO HOUREMIS 88 8 16 3 STACK2 41.349 437.276 18.77 SO HOUREMIS 88 8 16 4 STACK1 36.020 374.827 9.63 SO HOUREMIS 88 8 16 4 STACK2 43.672 437.682 18.23 The use of hourly varying release heights and initial dispersion coefficients for VOLUME and AREA sources is illustrated in the following example: 44 ------- SO HOUREMIS 88 3 1 1 VOL1 500.0 2.0 2.0 2.0 SO HOUREMIS 88 3 1 1 AREA1 5.000 2.0 2.0 SO HOUREMIS 88 3 1 2 VOL1 500.0 2.0 2.0 3.0 SO HOUREMIS 88 3 1 2 AREA1 5.000 2.0 3.0 SO HOUREMIS 88 3 1 3 VOL1 500.0 2.0 2.0 4.0 SO HOUREMIS 88 3 1 3 AREA1 5.000 2.0 4.0 For POINT sources, the model will use the stack release height and stack inside diameter defined on the SO SRCPARAM card, but will use the emission rate, exit temperature and exit velocity from the hourly emission file. As noted above regarding VOLUME and AREA sources, if the emission rate, exit temperature and exit velocity are not included for a particular hour, i.e, any or all of those fields are blank, the model will interpret emissions data for that hour as missing and will set the parameters to zero. Since the emission rate will be zero, there will be no calculations made for that hour and that source. 2.6 OUTPUT FILE OPTIONS Two new options have been incorporated on the OU pathway, beginning with version 09260 of AERMOD. The optional SUMMFILE keyword can be used to generate a separate formatted output file containing the summary of high ranked values included at the end of the standard 'aermod.ouf file. The optional FILEFORM keyword can be used to specify the use of exponential notation, rather than fixed format as currently used, for results that are output to separate result files. The syntax, type and order of the optional SUMMFILE keyword are summarized below: Syntax: OU SUMMFILE SummFileName Type: Optional, Non-repeatable where the SummFileName is the name of the external file containing the summary of high ranked values. The SUMMFILE filename can be up to 200 characters in length based on the default parameters in AERMOD. Double quotes (") at the beginning and end of the filename can also be used as field delimiters to allow filenames with embedded spaces. In addition to the 45 ------- summary of high ranked values, the SUMMFILE also includes the "MODEL SETUP OPTIONS SUMMARY" page from the main 'aermod.out' file. The syntax, type and order of the optional FILEFORM keyword are summarized below: Syntax: OU FILEFORM EXP or FIX Type: Optional, Non-repeatable where the EXP parameter specifies that output results files will use exponential-formatted values, and the FIX parameter specifies that the output results files will use fixed-formatted values. The default option is to use fixed-formatted results, so use of FILEFORM = 'FIX' is extraneous. Note that AERMOD only examines the first three characters of the input field, so that the full terms of 'EXPONENTIAL' or 'FIXED' can also be used. The format specified on this optional keyword is applicable to PLOTFILEs, plot-formatted POSTFILEs, MAXIFILEs, RANKFILEs, and SEASONHR files, but will not affect the format of results in the standard 'aermod.out' file or the optional SUMMFILE. The FILEFORM optional may be useful to preserve precision in applications with relatively small impacts, especially for the purpose of post-processing hourly concentrations using the POSTFILE option. The option may also be useful for applications with relatively large impacts that may overflow the Fortran format specifier of F13.5 used for fixed-formatted outputs. AERMOD will issue a warning message if values that exceed the range allowed for fixed-format are detected unless the FILEFORM EXP option has been selected. 2.7 MODEL STORAGE LIMITS The AERMOD model has been designed using a dynamic storage allocation approach, where the model allocates data storage as needed based on the number of sources, receptors, source groups, and input requirements, up to the maximum amount of memory available on the computer being used. The AERMOD model uses dynamic arrays to allocate data storage at model runtime rather than at compile time. The AERMOD model preprocesses the model runstream input file to determine the data storage requirements for a particular model run, and 46 ------- then allocates the input data arrays before processing the setup data. Once the setup processing is completed, the model allocates storage for the result arrays. When allocating data storage, the AERMOD model traps for errors, e.g., not enough memory available to allocate. If the allocation is unsuccessful, then an error message is generated by the model and further processing is prevented. If the CO RUNORNOT NOT option is selected, the model will still go through all array allocations so that the user can determine if sufficient memory is available to complete the run. Also, a rough estimate of the total amount of memory needed for a particular run is printed out as part of the first page of printed output. The storage parameters that are established at model runtime are as follows: NSRC = Number of Sources NREC = Number of Receptors NGRP = Number of Source Groups NOLM = Number of OLM Groups (OLMGROUP Keyword) NAVE = Number of Short Term Averaging Periods NVAL = Number of High Values by Receptor (RECTABLE Keyword) NTYP = Number of Output Types (CONC currently is the only output type) NMAX = Number of Overall Maximum Values (MAXTABLE Keyword) NQF = Number of Variable Emission Rate Factors per Source NPDMAX = Number of Particle Diameter Categories per Source NVMAX = Number of Vertices for Area Sources (including AREA, AREACIRC, and AREAPOLY source types) and/or OPENPIT Sources NSEC = Number of Sectors for Building Downwash Parameters (set to 36 if downwash sources are included) NURB = Number of Urban Areas (URB ANOPT Keyword) NNET = Number of Cartesian and/or Polar Receptor Networks IXM = Number of X-coord (Distance) Values per Receptor Network IYM = Number of Y-coord (Direction) Values per Receptor Network NARC = Number of Receptor Arcs Used with EVALCART Keyword NEVE = Number of Events for EVENT processing 47 ------- 3.0 REFERENCES Cimorelli, A. 1, S. G. Perry, A. Venkatram, J. C. Weil, R. J. Paine, R. B. Wilson, R. F. Lee, W. D. Peters, R. W. Erode, and J. O. Paumier, 2004: AERMOD: Description of Model Formulation. EPA 454/R-03-004. U. S. Environmental Protection Agency, Research Triangle Park, North Carolina 27711. EPA, 1995a: User's Guide for the Industrial Source Complex (ISC3) Dispersion Models, Volume I - User Instructions. EPA-454/B-95-003a. U.S. Environmental Protection Agency, Research Triangle Park, North Carolina 27711. EPA, 1995b: User's Guide for the Industrial Source Complex (ISC3) Dispersion Models, Volume II - Description of Model Algorithms. EPA-454/B-95-003b. U.S. Environmental Protection Agency, Research Triangle Park, North Carolina 27711. EPA, 2003: AERMOD Deposition Algorithms - Science Document (Revised Draft). U.S. Environmental Protection Agency, Research Triangle Park, North Carolina 27711. EPA, 2004a: User's Guide for the AMS/EPA Regulatory Model - AERMOD. EPA- 454/B-03-001. U.S. Environmental Protection Agency, Research Triangle Park, North Carolina 27711. EPA, 2004b: User's Guide for the AERMOD Meteorological Processor (AERMET). EPA- 454/B-03-002. U.S. Environmental Protection Agency, Research Triangle Park, North Carolina 27711. Hanrahan, P.L., 1999a. "The plume volume molar ratio method for determining NO2/NOx ratios in modeling. Part I: Methodology," J. Air & Waste Manage. Assoc., 49, 1324-1331. Hanrahan, P.L., 1999b. "The plume volume molar ratio method for determining NO2/NOx ratios in modeling. Part II: Evaluation Studies," J. Air & Waste Manage. Assoc., 49, 1332-1338. Walcek, C., G. Stensland, L. Zhang, H. Huang, J. Hales, C. Sweet, W. Massman, A. Williams, J, Dicke, 2001: Scientific Peer-Review of the Report "Deposition Parameterization for the Industrial Source Complex (ISC3) Model." The KEVRIC Company, Durham, North Carolina. Wesely, M.L, P. V. Doskey, and J.D. Shannon, 2002: Deposition Parameterizations for the Industrial Source Complex (ISC3) Model. Draft ANL report ANL/ER/TR-01/003, DOE/xx- nnnn, Argonne National Laboratory, Argonne, Illinois 60439. Note: All references listed here, with the exception of Hanrahan (1999a and 1999b), can be found on the U.S. EPA SCRAM website at the following url: http://www.epa.gov/scram001/ 48 ------- APPENDIX A. ALPHABETICAL KEYWORD REFERENCE This appendix provides an alphabetical listing of all of the keywords used by the AERMOD model. Each keyword is identified as to the pathway for which it applies, the keyword type (either mandatory or optional, and either repeatable or non-repeatable), and with a brief description of the function of the keyword. For a more complete description of the keywords, including a list of associated parameters, refer to Section 2.0 of this Addendum, the Functional Keyword/Parameter Reference in Appendix B of this Addendum, and/or the Detailed Keyword Reference in Section 3 of the AERMOD User's Guide (EPA, 2004a). A-l ------- Keyword AREA VERT AVERTIME BUILDHGT BUILDLEN BUILDWID CONCUNIT DAYRANGE DAYTABLE DCAYCOEF DEBUGOPT DEPOUNIT DISCCART DISCPOLR ELEVUNIT EMISFACT EMISUNIT ERRORFIL EVALCART EVALFILE EVENTFIL EVENTLOC EVENTOUT EVENTPER FILEFORM Path SO CO so so so CO ME OU CO CO SO RE RE SO RE SO SO CO RE OU CO EV OU EV OU Type M-R M-N O-R O-R O-R O-N O-R O-N O-N O-N O-N O-R O-R O-N O-N O-R O-N O-N O-R O-R O-N M-R M-N M-R O-N Keyword Description Specifies location of vertices for an AREAPOLY source type (mandatory if AREAPOLY source is used) Averaging time(s) to process Building height values for each wind sector Building projected length values for each wind sector Building projected width values for each wind sector Optional conversion factors for emission input units and concentration output units Specifies days or ranges of days to process (default is to process all data) Option to provide summaries for each averaging period for each day processed. Optional decay coefficient for exponential decay Option to generate detailed result and meteorology files for debugging purposes Optional conversion factors for emission input units and deposition output units Defines discretely placed receptors referenced to a Cartesian system Defines discretely placed receptors referenced to a polar system Defines input units for receptor elevations (RE path), or source elevations (SO path) (defaults to meters) Optional input for variable emission rate factors Optional conversion factors for emission units and concentration units Option to generate detailed error listing file (error file is mandatory for CO RUNORNOT NOT case) Defines discretely placed receptor locations referenced to a Cartesian system, grouped by arc for use with the EVALFILE output option Option to output file of normalized arc maxima for model evaluation studies Specifies whether to generate an input file for EVENT model Describes receptor location for an event Specifies level of output information provided by the EVENT model Describes data and averaging period for an event Specify fixed or exponential format for output results files A-2 ------- Keyword FINISHED FLAGPOLE GASDEPDF GASDEPOS GASDEPVD GDLANUSE GDSEASON GRIDCART GRIDPOLR HALFLIFE HOUREMIS INCLUDED INITFILE LOCATION MASSFRAX MAXIFILE MAXTABLE METHOD 2 MODELOPT MULTYEAR NO2EQUIL NO2RATIO NO2STACK OLMGROUP OZONEFIL Path ALL CO CO SO CO CO CO RE RE CO SO SO, RE, EV CO SO SO OU OU SO CO CO CO SO CO so CO Type M-N O-N O-N O-R O-N O-N O-N O-R O-R O-N O-R O-R O-N M-R O-R O-R O-R O-R M-N O-N O-N O-R O-N O-R O-N Keyword Description Identifies the end of inputs for a particular pathway Specifies whether to accept receptor heights above local terrain (m) for use with flagpole receptors, and allows for a default flagpole height to be specified Option to override default parameters for gas dry deposition Specify source parameters for gas deposition algorithms Option to specify deposition velocity for gas dry deposition Specify land use categories by sector for gas dry deposition Specify seasonal definitions for gas dry deposition Defines a Cartesian grid receptor network Defines a polar receptor network Optional half-life for exponential decay Option for specifying hourly emission rates in a separate file Option to include input data from a separate file in the runstream for the SO and/or RE pathways, or for the EV pathway for EVENTs Option to initialize model from file of intermediate results generated by SAVEFILE option Identifies coordinates for particular source location Optional input of mass fraction for each particle size category Option to list events exceeding a threshold value to file (if CO EVENTFIL option is used, these events are included in the input file generated for the EVENT model) Option to summarize the overall maximum values Specify optional source parameters for METHOD 2 option for particle deposition Job control and dispersion options Specifies that run is part of a multi-year run, e.g., for PM-10 H6H in five years Option to override default NO2/NOx equilibrium ratio for PVMRM Option to specify in-stack NO2/NOx equilibrium ratio for OLM and PVMRM options by source Option to specify default in-stack NO2/NOx equilibrium ratio for OLM and PVMRM options; may be overridden by NO2RATIO Specifies sources to combine for OLM option for merging plumes Specifies hourly ozone file for OLM and PVMRM options A-3 ------- Keyword OZONEVAL PARTDENS PARTDIAM PLOTFILE POLLUTID POSTFILE PROFBASE PROFFILE PSDGROUP RANKFILE RECTABLE RUNORNOT SAVEFILE SCIMBYHR SEASONHR SITEDATA SRCGROUP SRCPARAM STARTEND STARTING SUMMFILE SURFDATA SURFFILE TITLEONE TITLETWO TOXXFILE UAIRDATA URBANOPT URBANSRC WDROTATE WINDCATS Path CO SO SO OU CO OU ME ME SO OU OU CO CO ME OU ME SO SO ME ALL OU ME ME CO CO OU ME CO SO ME ME Type O-N O-R O-R O-R M-N O-R M-N M-N O-R O-R O-R M-N O-N O-N O-R O-N M-R M-R O-N M-N O-N M-N M-N M-N O-N O-R M-N O-R O-R O-N O-N Keyword Description Specifies background value of ozone for OLM and PVMRM options Specifies particle density by size category for particle deposition Specifies particle diameters by size category for particle deposition Option to write certain results to a storage file suitable for input to plotting routines Identifies pollutant being modeled Option to write results to a mass storage file for postprocessing Specifies the base elevation for the potential temperature profile Describes input profile meteorological data file Specifies source groups for PSDCREDIT option with PVMRM Option to produce output file of ranked values for Q-Q plots Option to output high ranked value (s) by receptor Identifies whether to run model or process setup information only Option to store intermediate results for later restart of the model after user or system interrupt Specifies sampling parameters for the SCIM option Option to output values by season and hour-of-day Describes on-site meteorological station Identification of source groups Identifies source parameters for a particular source Specifies start and end dates to be read from input meteorological data file (default is to read entire file) Identifies the start of inputs for a particular pathway Option to output summary of high ranked values to separate file Surface meteorological station Describes input surface meteorological data file First line of title for output Optional second line of output title Creates output file for use with TOXX model component of TOXST Upper air meteorological station Option to specify population for urban option Option to specify use of urban option by source Wind direction rotation adjustment Upper bound of wind speed categories A-4 ------- Keyword XBADJ YBADJ Path SO SO Type O-R O-R Keyword Description Along -flow distances from the stack to the center of the upwind face of the projected building Across-flow distances from the stack to the center of the upwind face of the projected building A-5 ------- APPENDIX B. FUNCTIONAL KEYWORD/PARAMETER REFERENCE This appendix provides a functional reference for the keywords and parameters used by the input runstream files for the AERMOD model. The keywords are organized by functional pathway, and within each pathway the order of the keywords is based on the function of the keyword within the preprocessor. The pathways used by the preprocessor are as follows, in the order in which they appear in the runstream file and in the tables that follow: CO - for specifying overall job COntrol options; and SO - for specifying SOurce location information (optional); D1? - for specifying REceptor information; and - for specifying MEteorology information and options; - for specifying EVent information and options; - for specifying OUtput file information. RE ME EV- OU The pathways and keywords are presented in the same order as in the Detailed Keyword Reference in Section 3 of the AERMOD User's Guide (EPA, 2004a), and in the Quick Reference at the end of the manual, with the exception of new keywords that are not reflected in the original user's guide. Two types of tables are provided for each pathway. The first table lists all of the keywords for that pathway, identifies each keyword as to its type (either mandatory or optional and either repeatable or non-repeatable), and provides a brief description of the function of the keyword. The second type of table presents the parameters for each keyword, in the order in which they should appear in the runstream file where order is important, and describes each parameter in detail. The following convention is used for identifying the different types of input parameters. Parameters corresponding to secondary keywords which should be input "as is" are listed on the tables with all capital letters and are underlined, although none of the inputs to AERMAP are treated as case-sensitive. Other parameter names are given with an initial capital letter and are not input "as is." In all cases, the parameter names are intended to be descriptive of the input variable being represented, and they often correspond to the Fortran variable names used in the B-l ------- preprocessor code. Parentheses around a parameter indicate that the parameter is optional for that keyword. The default that is taken when an optional parameter is left blank is explained in the discussion for that parameter. B-2 ------- TABLE B-l DESCRIPTION OF CONTROL PATHWAY KEYWORDS CO Keywords STARTING TITLEONE TITLETWO MODELOPT AVERTIME URBANOPT POLLUTID HALFLIFE DCAYCOEF GASDEPDF GASDEPVD GDLANUSE GDSEASON NO2EQUIL NO2STACK OZONEFIL OZONEVAL FLAGPOLE RUNORNOT EVENTFIL2 SAVEFILE3 INITFILE3 MULTYEAR3 DEBUGOPT ERRORFIL FINISHED Type M-N M-N O-N M-N M-N O-R M-N O-N1 O-N1 O-N O-N O-N O-N O-N O-N O-N O-N O-N M-N O-N O-N O-N O-N O-N O-N M-N Keyword Description Identifies the start of CONTROL pathway inputs First line of title for output Optional second line of title for output lob control and dispersion options Averaging time(s) to process Specifies parameters for urban dispersion option Identifies type of pollutant being modeled Optional half life used for exponential decay Optional decay coefficient Option to override default parameters for gas dry deposition Option to specify deposition velocity for gas dry deposition Specify land use categories by sector for gas dry deposition Specify seasonal definitions for gas dry deposition Option to override default NO2/NOx equilibrium ratio for PVMRM Option to specify default in-stack NO2/NOx equilibrium ratio for OLM and PVMRM options; may be overridden by NO2RATIO option on SO pathway Specifies filename for hourly ozone file for use with OLM and PVMRM options Specifies background value of ozone for use with OLM and PVMRM options Specifies whether to accept receptor heights above local terrain (m) for use with flagpole receptors, and allows for default flagpole height to be specified Identifies whether to run model or process setup information only Specifies whether to generate an input file for EVENT model Option to store intermediate results for restart of model after user or system interrupt Option to initialize model from intermediate results generated by SAVEFILE option Option to process multiple years of meteorological data (one year per run) and accumulate high short term values across years Option to generate detailed result and meteorology files for debugging purposes Option to generate detailed error listing file Identifies the end of CONTROL pathway inputs ------- Type: M - Mandatory N - Non-Repeatable O - Optional R - Repeatable 1) Either HALFLIFE or DCAYCOEF may be specified. If both cards appear a warning message will be issued and the first value entered will be used in calculations. The DFAULT option assumes a half life of 4 hours for S02 modeled in urban mode. 2) The EVENTFIL keyword controls whether or not to generate an input file for EVENT processing. The primary difference between AERMOD regular processing and EVENT processing by AERMOD is in the treatment of source group contributions. The AERMOD model treats the source groups independently, whereas EVENT processing determines individual source contributions to particular events, such as the design concentrations determined from AERMOD, or user-specified events. By specifying the EVENTFIL keyword, an input runstream file will be generated that can be used directly for EVENT processing. The events included in the generated EVENT processing input file are defined by the RECTABLE and MAXIFILE keywords on the OU pathway, and are placed in the EVent pathway. 3) The SAVEFILE and INITFILE keywords work together to implement the model's re-start capabilities. Since the MULTYEAR option utilizes the re-start features in a special way to accumulate high short term values from year to year, it cannot be used together with the SAVEFILE or INITFILE keyword in the same model run. B-4 ------- TABLE B-2 DESCRIPTION OF CONTROL PATHWAY KEYWORDS AND PARAMETERS Keyword Parameters TITLEONE Title 1 where: Title 1 First line of title for output, character string of up to 68 characters TITLETWO Title2 where: Title2 Optional second line of title for output, character string of up to 68 characters MODELOPT DFAULT BETA CONG AREADPLT FLAT NOSTD NOCHKD NOWARN SCREEN SCIM PVMRM PSDCREDIT DDEP and/or WDEP DRYDPLT or FASTAREA and/or ELEV WETDPLT or NOWETDPLT where: DFAULT BETA CONC DEPPS DDEP WDEP AREADPLT FLAT ELEV NOSTD NOCHKD WARNCHKD NOWARN SCREEN Specifies that regulatory default options will be used; specification of DFAULT option will override non-DFAULT options that may be specified; Non-DFAULT option that allows for draft, "Beta" test options to be used; includes PSDCREDIT option and capped and horizontal stack releases; Specifies that concentration values will be calculated; Specifies that total deposition flux values will be calculated; Specifies that dry deposition flux values will be calculated; Specifies that wet deposition flux values will be calculated; Specifies use of non-DFAULT method for optimized plume depletion due to dry removal mechanisms for area sources; Non-DFAULT option of assuming flat terrain will be used; Default option of assuming elevated terrain will be used; Note that FLAT and ELEV may be specified in the same model run to allow specifying the non-DFAULT FLAT terrain option on a source-by-source basis; Non-DFAULT option of no stack-tip downwash will be used; Non-DFAULT option of suspending date checking will be used for non-sequential meteorological data files, also implemented when SCREEN option is specified; Specifies option for issuing warning messages rather than fatal errors for non-sequential meteorological data files Option to suppress detailed listing of warning messages in the main output file will be used; Non-DFAULT option for running AERMOD in a screening mode for AERSCREEN will be used; B-5 ------- Keyword AVERTIME where: URBANOPT where: POLLUTID where: HALFLIFE where: DCAYCOEF where: Parameters SCIM PVMRM OLM PSDCREDIT FASTALL FASTAREA DRYDPLT NODRYDPLT WETDPLT NOWETDPLT Sampled Chronological Input Model - ANNUAL averages only; SCIM sampling parameters must be specified on the ME pathway; Non-DFAULT Plume Volume Molar Ratio Method (PVMRM) for NO2 conversion will be used; Non-DFAULT Ozone Limiting Method (OLM) for NO2 conversion will be used; Specifies use of BETA test option to calculate the increment consumption with PSD credits using the PVMRM option; Non-DFAULT option to optimize model runtime for POINT, VOLUME and AREA sources (AREA optimizations formerly associated with TOXICS option); Non-DFAULT option to optimize model runtime for AREA sources (formerly associated with TOXICS option); Option to incorporate dry depletion (removal) processes associated with dry deposition algorithms; dry depletion will be used by default if dry deposition algorithms are invoked; Option to disable dry depletion (removal) processes; Option to incorporate wet depletion (removal) processes associated with wet deposition algorithms; wet depletion will be used by default if wet deposition algorithms are invoked; and Option to disable wet depletion (removal) processes. Timel Time2 . . . TimeN MONTH PERIOD or ANNUAL TimeN MONTH PERIOD ANNUAL Nth optional averaging time (123468 12 or 24-hr) Option to calculate MONTHly averages (treated as Option to calculate averages for the entire data PERIOD Option to calculate ANNUAL averages (assumes complete years) Urbpop (Urbname) (UrbRoughness) UrbPop UrbName UrbRoughness Specifies the population of the urban area Specifies the name of the urban area (optional) Specifies the urban surface roughness length, meters (optional, defaults to 1.0m; value other than 1.0m treated as non-DFAULT) Pollut Pollut Identifies type of pollutant being modeled. Any name of up to eight characters mav be used. e.e.. SO2. NOX. CO. PM25. PM-2.5. PM10. PM-10. TSP or OTHER. Use of PM10. PM-10 or OTHER allows for the use of the MULTYEAR option. Haflif Haflif Half life used for exponential decay (s) Decay Decay Decay coefficient for exponential decay (s"1) = 0.693/HAFLIF B-6 ------- Keyword GASDEPDF where: GASDEPVD where: GDLANUSE where: GDSEASON where: NO2EQUIL where: NO2STACK where: OZONEFIL where: OZONEVAL where: FLAGPOLE where: Parameters React F_Seas2 F_Seas5 (Refpoll) React F Seas2 F SeasS (Refpoll) Value for pollutant reactivity factor (Q Fraction (F) of maximum green LAI for seasonal category 2 Fraction (F) of maximum green LAI for seasonal category 5 Optional name of reference pollutant Uservd Uservd User-specified dry deposition velocity (m/s) for gaseous pollutants Seel Sec2 ... Sec36 Seel Sec2 Sec36 Land use category for winds blowing toward sector 1(10 degrees) Land use category for winds blowing toward sector 2 (20 degrees) Land use category for winds blowing toward sector 36 (360 degrees) Jan Feb ... Dec Jan Dec Seasonal category for January: 1 = Midsummer/Lush vegetation; 2 = Autumn/Unharvested cropland; 3 = Late autumn after harvest or Winter with no snow; 4 = Winter with continuous snow cover; or 5 = Transitional spring/partial green coverage/short annuals) Seasonal category for December NO2Equil NO2Equil Equilibrium ratio for PVMRM option (default is 0.9) NO2Ratio NO2Ratio Default in-stack ratio of NO2/NOx for OLM and PVMRM options; default is 0. 1 for OLM, which may be overridden by NO2RATIO keyword on SO pathway OSFileName (OSUnits) (Format) OSFileName (OSUnits) (Format) Filename for hourly ozone data file (YR, MN, DY, HR, O3Value) Units of ozone data (PPM, PPB, or UG/M3); default is UG/M3 Fortran format statement to read ozone file; default is free-format, i.e., comma or space-delimited data fields OSValue (OSUnits) OSValue (OSUnits) Background ozone concentration; also used to substitute for missing data in OZONEFIL Units of ozone value (PPM, PPB, or UG/M3); default is UG/M3 (Flagdf) Flagdf Default value for height of (flagpole) receptors above local ground, a B-7 ------- Keyword RUNORNOT where: EVENTFIL where: SAVEFILE where: INITFILE where: MULTYEAR where: DEBUGOPT where: ERRORFIL where: Parameters RUN or NOT RUN NOT default value of 0.0 m is used if this optional parameter is omitted Indicates to run full model calculations Indicates to process setup data and report errors, but to not run full model calculations (Evfile) (Evopt) Evfile Evopt Identifies the filename to be used to generate a file for input to EVENT model (Default=EVENTFIL.INP) Optional parameter to specify the level of output detail selected for the EVENT model: either SOCONT or DETAIL (default is DETAIL if this parameter is omitted) (Savfil) (Dayinc) (Savfl2) Savfil Dayinc Savfl2 Specifies name of disk file to be used for storing intermediate results (default = TMP.FIL) file is overwritten after each dump) Number of days between dumps (optional: default is 1) Optional second disk filename to be used on alternate dumps - eliminates risk of system crash during the dump. If blank, file is overwritten each time. (Inifil) Inifil H6H Savfil dm H6H Savfil Inifil MODEL (Dbafi MODEL (Dbgfil) METEOR (Dbmfil) Specifies name of disk file of intermediate results to be used for initializing run (default = TMP.FIL) fil) Specifies the High- Sixth-High is being calculated for use in PM10 processing Specifies name of file to be used for storing results at end of the year Optional name of file used for initializing the results arrays from previous year(s). The Inifil parameter is not used for the first year in the multi-year run. 1) and/or METEOR (Dbmfil) Specifies that MODEL debugging output will be generated Optional filename for the model calculation debug file Specifies that METEORological profile data file will be generated Optional filename for the meteorological profile data file (Errfil) Errfil Specifies name of detailed error listing file (default = ERRORS. LST) ------- TABLE B-3 DESCRIPTION OF SOURCE PATHWAY KEYWORDS SO Keywords STARTING ELEVUNIT LOCATION SRCPARAM BUILDHGT BUILDLEN BUILDWID XBADI YBADI AREA VERT URBANSRC EMISFACT EMISUNIT CONCUNIT DEPOUNIT PARTDIAM MASSFRAX PARTDENS METHOD 2 GASDEPOS NO2RATIO HOUREMIS INCLUDED OLMGROUP Type M-N O-N M-R M-R O-R O-R O-R O-R O-R M-R O-R O-R O-N O-N O-N O-R O-R O-R O-R O-R O-R O-R O-R O-R Keyword Description Identifies the start of SOURCE pathway inputs Defines input units for source elevations (defaults to meters), must be first keyword after SO STARTING if used. Identifies coordinates for particular source Identifies source parameters for a particular source Building height values for each wind sector Building projected length values for each wind sector Building projected width values for each wind sector Along -flow distances from the stack to the center of the upwind face of the projected building Across-flow distances from the stack to the center of the upwind face of the projected building Specifies location of vertices for an AREAPOLY source type (mandatory if AREAPOLY source is used) Identifies which sources to model with urban effects Optional input for variable emission rate factors Optional unit conversion factors for emissions, concentrations Optional conversion factors for emissions and concentrations Optional conversion factors for emissions and depositions Input variables for optional input of particle size (microns) Optional input of mass fraction for each particle size category Optional input of particle density (g/cm3) for each size category Optional input of parameters for METHOD_2 particle deposition Optional input of gas deposition parameters Option to specify in-stack NO2/NOx equilibrium ratio for OLM and PVMRM options by source Option for specifying hourly emission rates in a separate file Option to include data from a separate file in the runstream Specifies sources to combine for OLM option to account for merging plumes B-9 ------- SO Keywords PSDGROUP1 SRCGROUP1 FINISHED Type O-R M-R M-N Keyword Description Specifies source groups for PSDCREDIT option with PVMRM Identification of source groups Identifies the end of SOURCE pathway inputs 1) The PSDGROUP or SRCGROUP keywords must be the last keyword within the SO pathway before the FINISHED keyword. The SRCGROUP keyword is mandatory, unless the PSDCREDIT option is used, which requires the PSDGROUP option instead. B-10 ------- TABLE B-4 DESCRIPTION OF SOURCE PATHWAY KEYWORDS AND PARAMETERS Keyword Parameters ELEVUNIT METERS or FEET where: METERS FEET Specifies input units for source elevations of meters (default if ELEVUNIT is omitted) Specifies input units for source elevations of feet Note: This keyword applies to source elevations only. LOCATION Srcid Srctyp Xs Ys (Zs) where: Srcid Srctyp Xs Ys Zs Source identification code (alphanumeric string of up to eight characters) Source type: POINT. VOLUME. AREA. AREAPOLY. AREACIRC. OPENPIT x-coord of source location, corner for AREA. AREAPOLY. and OPENPIT. center for AREACIRC (m) y-coord of source location, corner for AREA. AREAPOLY. and OPENPIT. center for AREACIRC (m) Optional z-coord of source location (elevation above mean sea level, defaults to 0.0 if omitted) SRCPARAM Srcid Ptemis Stkhgt Stktmp Stkvel Stkdia Vlemis Relhgt Syinit Szinit Aremis Relhgt Xinit (Yinit) (Angle) (Szinit) Aremis Relhgt Nverts (Szinit) Aremis Relhgt Radius (Nverts) (Szinit) Opemis Relhgt Xinit Yinit Pitvol (Angle) (POINT source) (VOLUME source) (AREA source) (AREAPOLY source) (AREACIRC source) (OPENPIT source) where: Srcid Emis Hgt Stktmp Stkvel Stkdia Syinit Szinit Xinit Yinit Angle Nverts Source identification code Source emission rate: in g/s for Ptemis or Vlemis, g/(s-m2) for Aremis for concentration Source physical release height above ground (center of height for VOLUME, height above base of pit for OPENPIT) Stack gas exit temperature (K) Stack gas exit velocity (m/s) Stack inside diameter (m) Initial lateral dimension of VOLUME source (m) Initial vertical dimension of VOLUME or AREA source (m) Length of side of AREA source in X-direction (m) Length of side of AREA source in Y-direction (m) (optional parameter, assumed to be equal to Xinit if omitted) Orientation angle (deg) of AREA or OPENPIT source relative to North measured positive clockwise, rotated around the source location, (XsYs) (optional parameter, assumed to be 0.0 if omitted) Number of vertices used for AREAPOLY or AREACIRC source (optional for AREACIRC sources) B-ll ------- Keyword BUILDHGT where: BUILDWID where: XBADJ where: YBADJ where: AREA VERT where: URBANSRC where: EMISFACT Parameters Radius Pitvol Radius of circular area for AREACIRC source (m) Volume of OPENPIT source (m3) Srcid (or Srcrng) Dsbh(i), i=l,36 Srcid Srcrng Dsbh Source identification code Range of sources (inclusive) for which building dimensions apply, entered as two alphanumeric strings separated by a '-' Array of direction-specific building heights (m) beginning with 10 degree flow vector and increment- ing by 10 degrees clockwise Srcid (or Srcrng) Dsbw(i), i=l,36 Srcid Srcrng Dsbw Source identification code Range of sources (inclusive) for which building dimensions apply Array of direction-specific building widths (m) beginning with 10 degree flow vector and increment- ing by 10 degrees clockwise Srcid (or Srcrng) Xbadj(i), i=l,36 Srcid Srcrng Xbadj(i) Source identification code Range of sources (inclusive) for which XBADJ distances apply Array of direction-specific along -wind distances beginning with 10 degree flow vector and incrementing by 10 degrees clockwise Srcid (or Srcrng) Ybadj(i), i=l,36 Srcid Srcrng Ybadj(i) Source identification code Range of sources (inclusive) for which YBADJ distances apply Array of direction-specific across-wind distances beginning with 10 degree flow vector and incrementing by 10 degrees clockwise Srcid Xv(l) Yv(l) Xv(2) Yv(2) ... Xv(i) Yv(i) Srcid Xv(l) Yv(l) Xv(i) Yv(i) Source identification code X-coordinate of the first vertex of an AREAPOLY source (must be the same as the value of Xs for that source defined on the SO LOCATION card) Y-coordinate of the first vertex of an AREAPOLY source (must be the same as the value of Ys for that source defined on the SO LOCATION card) X-coordinate for the ith vertex of an AREAPOLY source Y-coordinate for the ith vertex of an AREAPOLY source Srcid's and/or Srcrng's Srcid Srcrng Specifies which source(s) will be modeled with urban effects Specifies a range of sources that will be modeled with urban effects Srcid (or Srcrng) Qflag Qfact(i), i=l,n B-12 ------- Keyword where: EMISUNIT where: CONCUNIT where: DEPOUNIT where: PARTDIAM where: MASSFRAX where: PARTDENS Parameters Srcid Srcrng Qflag Qfact Source identification code Range of sources (inclusive) for which emission rate factors apply Variable emission rate flag: SEASON for seasona; MONTH for monthly; HROFDY for hour-of-day; WSPEED for wind speed category; SEASHR for season-by-hour; HRDOW for emission rates vary by hour-of-day. and day-of-week [M-F, Sat, Sun]; HRDOW7 for emission rates vary by hour-of-day, and the seven days of the week [M, Tu, W, Th, F, Sat, Sun]; SHRDOW for season by hour-of-day by day-of- week (M-F,Sat,Sun); SHRDOW7 for season by hour-of-day by day-of-week (M,Tu,W,Th,F,Sat,Sun); MHRDOW for month by hour-of-dav by dav-of-week (M-F.Sat.Sun); MHRDOW7 for month by hour-of-day by day-of-week (M,Tu,W,Th,F,Sat,Sun) Array of scalar emission rate factors, for: SEASON. n=4; MONTH. n=12; HROFDY, n=24; WSPEED. n=6; SEASHR. n=96; HRDOW. n=72; HRDQW7,n=168; SHRDOW, n=28 8; SHRDQW7, n=672; MHRDOW. n=864; MHRDOW7. n=2016 Emifac Emilbl Conlbl Emifac Emilbl Conlbl Emission rate factor used to adjust units of output (default value is 1 .OE06 for CONC for grams to micrograms) Label to use for emission units (default is grams/sec) Label to use for concentrations (default is micrograms/m3) Emifac Emilbl Conlbl Emifac Emilbl Conlbl Emission rate factor used to adjust units of output (default value is 1 .OE06 for concentration for grams to micrograms) Label to use for emission units (default is grams/sec) Label to use for concentrations (default is micrograms/m3) Emifac Emilbl Deplbl Emifac Emilbl Deblbl Emission rate factor used to adjust units of output for deposition (default value is 3600 for grams/sec to grams/m2/hr) Label to use for emission units (default is grams/sec) Label to use for deposition (default is grams/m2) Srcid (or Srcrng) Pdiam(i), i=l,Npd Srcid Srcrng Pdiam Source identification code Range of sources (inclusive) for which size categories apply Array of particle diameters (microns) Srcid (or Srcrng) Phi(i), i=l,Npd Srcid Srcrng Phi Source identification code Range of sources (inclusive) for which mass fractions apply Array of mass fractions for each particle size category Srcid (or Srcrng) Pdens(i), i=l,Npd B-13 ------- Keyword where: METHOD 2 where: GASDEPOS where: NO2RATIO where: HOUREMIS where: INCLUDED where: OLMGROUP where: PSDGROUP where: Parameters Srcid Srcrng Pdens Source identification code Range of sources (inclusive) for which particle densities apply Array of particle densities (g/cm3) for each size category Srcid (or Srcrng) FineMassFraction Dmm Srcid FineMassFractio n Dmm Source identification code Fraction of particle mass emitted in fine mode, less than 2.5 microns Representative mass mean particle diameter in microns SrcID (or SrcRange) Da Dw rcl Henry SrcID Da Dw rcl Henry Source identification code Diffusivity in air for the pollutant being modeled (cm2/s) Diffusivity in water for the pollutant being modeled (cm2/s) Cuticular resistance to uptake by lipids for individual leaves (s/cm) Henry's Law constant (Pa m3/mol) SrcID (or SrcRange) NO2Ratio SrcID SrcRange NO2Ratio Source identification code Source ID range for specified ratio In-stack ratio of NO2/NOx Emifil Srcid's Srcrng's Emifil Srcid's Srcrng's Specifies name of the hourly emission rate file Discrete source IDs that are included in the hourly emission file Source ID ranges that are included in the hourly emission file Incfil SrcIncFile Filename for the included source file, up to 200 characters in length; double quotes (") may be used as delimiters for the filename to allow for embedded spaces; and quotes don't count toward the limit of 200 OLMGrpID SrcID's SrcRange's OLMGrpID SrcID's SrcRange's Group ID (Grpid = ALL specifies group including all sources) Discrete source IDs to be included in group Source ID ranges to be included in group Note: Card may be repeated with same Grpid if more space is needed to specify sources PSDGrpID SrcID's SrcRange's PSDGrpID SrcID's SrcRange's PSD GroupID for PSDCREDIT option, must be one of the following: INCRCONS - increment-consuming sources, NONRBASE - non-retired baseline sources, or RETRBASE - retired (increment-expanding) baseline sources. Discrete source IDs to be included in group Source ID ranges to be included in group B-14 ------- Keyword SRCGROUP where: Parameters Note: Card may be repeated with same PSDGrpID if more space is needed to specify sources SrcGrpID SrcID's SrcRange's SrcGrpID SrcID's SrcRange's Group ID (Grpid = ALL specifies group including all sources) Discrete source IDs to be included in group Source ID ranges to be included in group Note: Card may be repeated with same Grpid if more space is needed to specify sources B-15 ------- TABLE B-5 DESCRIPTION OF RECEPTOR PATHWAY KEYWORDS RE Keywords STARTING ELEVUNIT GRIDCART GRIDPOLR DISCCART DISCPOLR EVALCART INCLUDED FINISHED Type M-N O-N O'-R O'-R O'-R O'-R O'-R O-R M-N Keyword Description Identifies the start of RECEPTOR pathway inputs Defines input units for receptor elevations (defaults to meters), keyword after RE STARTING if used. must be first Defines a Cartesian grid receptor network Defines a polar receptor network Defines the discretely placed receptor locations referenced to a system Defines the discretely placed receptor locations referenced to a Cartesian polar system Defines discrete Cartesian receptor locations for use with EVALFILE output option Identifies an external file containing receptor locations to be included in the inputs Identifies the end of RECEPTOR pathway inputs 1) At least one of the following must be present: GRIDCART, GRIDPOLR DISCCART, DISCPOLR or EVALCART, unless the INCLUDED keyword is used to include receptor inputs from an external file. Multiple receptor networks can be specified in a single run, including both Cartesian and polar. B-16 ------- TABLE B-6 DESCRIPTION OF RECEPTOR PATHWAY KEYWORDS AND PARAMETERS Keyword Parameters ELEVUNIT METERS or FEET where: METERS FEET Specifies input units for receptor elevations of meters Specifies input units for receptor elevations of feet Note: This keyword applies to receptor elevations only. GRIDCART Netid STA XYINC Xinit Xnum Xdelta Yinit Ynum Ydelta or XPNTS Gridxl Gridx2 GridxS .... GridxN, and YPNTS Gridyl Gridy2 GridyS .... GridyN ELEV Row Zelevl Zelev2 ZelevS ...ZelevN FLAG Row Zflagl Zflag2 ZflagS ...ZflagN END where: Netid STA XYINC Xinit Xnum Xdelta Yinit Ynum Ydelta XPNTS Gridxl GridxN YPNTS Gridyl GridyN ELEV Row Zelev FLAG Row Zflag END Receptor network identification code (up to eight alphanumeric characters) Indicates STArt of GRIDCART subpathway, repeat for each new Netid Keyword identifying grid network generated from x and y increments Starting local x-axis grid location in meters Number of x-axis receptors Spacing in meters between x-axis receptors Starting local y-axis grid location in meters Number of y-axis receptors Spacing in meters between y-axis receptors Keyword identifying grid network defined by a series of x and y coordinates Value of first x-coordinate for Cartesian grid Value of 'nth' x-coordinate for Cartesian grid Keyword identifying grid network defined by a series of x and y coordinates Value of first y-coordinate for Cartesian grid Value of 'nth' y-coordinate for Cartesian grid Keyword to specify that receptor elevations follow Indicates which row (y-coordinate fixed) is being input An array of receptor terrain elevations for a particular Row Keyword to specify that flagpole receptor heights follow Indicates which row (y-coordinate fixed) is being input An array of receptor heights above local terrain elevation for a particular Row (flagpole receptors) Indicates END of GRIDCART subpathway, repeat for each new Netid GRIDPOLR Netid STA ORIG Xinit Srcid Ringl or ORIG DIST Yinit, Ring2 RingS ... RingN B-17 ------- Keyword where: DISCCART where: DISCPOLR where: Parameters DDIR Dirl Dir2 Dir3 ... DirN. or GDIR Dirnum Dirini Dirinc ELEV Dir Zelevl Zelev2 ZelevS ... ZelevN FLAG Dir Zflael Zflae2 ZflaaS ... ZflaeN END Netid STA ORIG Xinit Yinit Srcid DIST Ringl RingN DDIR Dirl DirN GDIR Dirnum Dirini Dirinc ELEV Dir Zelev FLAG Dir Zflag END Receptor network identification code (up to eight alphanumeric characters) Indicates STArt of GRIDPOLR subpathwav, repeat for each new Netid Optional keyword to specify the origin of the polar network (assumed to be at x=0, y=0 if omitted) local x-coordinate for origin of polar network (m) local y-coordinate for origin of polar network (m) Source ID of source used as origin of polar network Keyword to specify distances for the polar network Distance to the first ring of polar coordinates (m) Distance to the 'nth' ring of polar coordinates (m) Keyword to specify discrete direction radials for the polar network First direction radial in degrees (1 to 360) The 'nth' direction radial in degrees (1 to 360) Keyword to specify generated direction radials for the polar network Number of directions used to define the polar system Starting direction of the polar system Increment (in degrees) for defining directions Keyword to specify that receptor elevations follow Indicates which direction is being input An array of receptor terrain elevations for a particular direction radial Keyword to specify that flagpole receptor heights follow Indicates which direction is being input An array of receptor heights above local terrain elevation for a particular direction (flagpole receptors) Indicates END of GRIDPOLR subpathwav, repeat for each new Netid Xcoord Ycoord (Zelev) (Zflag) Xcoord Ycoord Zelev Zflag local x-coordinate for discrete receptor location (m) local y-coordinate for discrete receptor location (m) Elevation above sea level for discrete receptor location (optional), used only for ELEV terrain Receptor height (flagpole) above local terrain (optional), used only with FLAGPOLE keyword Srcid Dist Direct (Zelev) (Zflag) Srcid Dist Direct Zelev Zflag Specifies source identification for which discrete polar receptor locations apply (used to define the origin for the discrete polar receptor) Downwind distance to receptor location (m) Direction to receptor location, in degrees clockwise from North Elevation above sea level for receptor location (optional), used onlv for ELEV terrain Receptor height (flagpole) above local terrain (optional), B-18 ------- Keyword EVALCART where: INCLUDED where: Parameters used only with FLAGPOLE keyword Xcoord Ycoord Zelev Zflag Arcid (Name) Xcoord Ycoord Zelev Zflag Arcid Name Local x-coordinate for discrete receptor location (m) Local y-coordinate for discrete receptor location (m) Elevation above sea level for discrete receptor location (optional), used only for ELEV terrain Receptor height (flagpole) above local terrain (optional), used only with FLAGPOLE keyword Receptor arc ID used to group receptors along an arc or other grouping (up to eight characters) Optional name for receptor (up to eight characters) RecIncFile RecIncFile Identifies the filename for the included receptor file, up to 200 characters in length; double quotes (") may be used as delimiters for the filename to allow for embedded spaces; and quotes don't count toward the limit of 200 B-19 ------- TABLE B-7 DESCRIPTION OF METEOROLOGY PATHWAY KEYWORDS ME Keywords STARTING SURFFILE PROFFILE SURFDATA UAIRDATA SITEDATA PROFBASE STARTEND DAYRANGE SCIMBYHR WDROTATE WINDCATS FINISHED Type M-N M-N M-N M-N M-N O-N M-N O-N O-R O-N O-N O-N M-N Keyword Description Identifies the start of METEOROLOGY pathway inputs Describes input meteorological surface data file Describes input meteorological profile data file Describes surface meteorological station Describes upper air meteorological station Describes on-site meteorological station Specifies the base elevation for the potential temperature profile Specifies start and end dates to be read from input meteorological data file (default is to read entire file) Specifies days or ranges of days to process (default is to process all data) Specifies the parameters for the SCIM (Sampled Chronological Input Model) option (see CO MODELOPT) May be used to correct for alignment problems of wind direction measurements, or to convert wind direction from to flow vector Input upper bounds of wind speed categories, five values input - sixth category is assumed to have no upper bound (used for WSPEED option on the EMISFACT keyword) Identifies the end of METEOROLOGY pathway inputs B-20 ------- TABLE B-8 DESCRIPTION OF METEOROLOGY PATHWAY KEYWORDS AND PARAMETERS Keyword SURFFILE where: PROFFILE where: SURFDATA where: UAIRDATA where: SITEDATA where: PROFBASE where: STARTEND where: Parameters Sfcfil (Format) Sfcfil Format Specify filename for surface meteorological input file Specify format for input file (Note: FREE format is used for SURFFILE reads beginning with version 09260) all Profil (Format) Profil Format Specify filename for profile meteorological input file Specify format for input file (Note: FREE format is used for all SURFFILE reads beginning with version 09260) Stanum Year (Name) (Xcoord Ycoord) Stanum Year Name Xcoord Ycoord Station number, e.g. 5 -digit WBAN number for NWS station Year of data being processed (four digits) Station name (optional) x-coordinate of station location (m) (optional) y-coordinate of station location (m) (optional) Stanum Year (Name) (Xcoord Ycoord) Stanum Year Name Xcoord Ycoord Station number, e.g. 5 -digit WBAN number for NWS station Year of data being processed (four digits) Station name (optional) x-coordinate of station location (m) (optional) y-coordinate of station location (m) (optional) Stanum Year (Name) (Xcoord Ycoord) Stanum Year Name Xcoord Ycoord Station number for on-site meteorological data station Year of data being processed (four digits) Station name (optional) x-coordinate of station location (m) (optional) y-coordinate of station location (m) (optional) BaseElev (Units) BaseElev Units Base elevation (above MSL) for the potential temperature prc Units of BaseElev: METERS or FEET (default is METERS) .file Strtyr Strtmn Strtdy (Strthr) Endyr Endmn Enddy (Endhr) Strtyr Strtmn Strtdy Strthr Endyr Endmn Year of first record to be read Month of first record to be read Day of first record to be read Hour of first record to be read (optional) Year of last record to be read Month of last record to be read B-21 ------- Keyword DAYRANGE where: SCIMBYHR where: WDROTATE where: WINDCATS where: Parameters Enddy Endhr Day of last record to be read Hour of last record to be read (optional) Note: File read begins with hour 1 of the start date and ends with hour 24 of the end date if Stahr and Endhr are omitted. Range 1 Range2 RangeS ... RangeN Range 1 RangeN First range of days to process, either as individual day (XXX) or as range (XXX-YYY); days may be input as Julian dates (XXX) or as month and day (XX/YY) The 'nth' range of days to process NRegStart NReglnt (SfcFilnam PflFilnam) NRegStart NReglnt SfcFilnam PflFilnam Specifies the first hour to be sampled with the SCIM option; required to have a value from 1 to 24 Specifies the sampling interval, in hours Optional output file name to list the surface meteorological data for the sampled hours Optional output file name to list the profile meteorological data for the sampled hours Rotang Rotang Specifies angle (in degrees) to rotate wind direction measurements to correct for alignment problems; value of Rotang is subtracted from WD measurements, i.e., rotation is counterclockwise Wsl Ws2 Ws3 Ws4 Ws5 Wsl Ws2 Ws3 Ws4 Ws5 Upper bound of first wind speed category (m/s) Upper bound of second wind speed category (m/s) Upper bound of third wind speed category (m/s) Upper bound of fourth wind speed category (m/s) Upper bound of fifth wind speed category (m/s) (sixth category is assumed to have no upper bound) B-22 ------- TABLE B-9 DESCRIPTION OF EVENT PATHWAY KEYWORDS EV Keywords STARTING EVENTPER EVENTLOC INCLUDED FINISHED Type M-N M-R M-R O-R M-N Keyword Description Identifies the start of EVENT pathway inputs Describes data and averaging period for an event Describes receptor location for an event Identifies an external file containing EVENT data to be included in the inputs Identifies the end of EVENT pathway inputs B-23 ------- TABLE B-10 DESCRIPTION OF EVENT PATHWAY KEYWORDS AND PARAMETERS Keyword EVENTPER where: EVENTLOC where: INCLUDED where: Parameters Evname Aveper Grpid Date Name Grpid Aveper Date Specify name of event to be processed (e.g. H2H24ALL), (up to eight alphanumeric characters) Specify source group ID for event Specify averaging period for event Specify data period for event (ending YYMMDDHH for averaging period) Evname XR=Xr YR=Yr (Zelev) (Zfiaz) or RNG= Rna DIR= Dir (Zelev) (ZQaz) Evname XR= YR= RNG= DIR= Zelev Zflag Specify name of event to be processed (e.g. H2H24ALL), (up to eight alphanumeric characters) X-coordinate for event (discrete Cartesian receptor) Y-coordinate for event (discrete Cartesian receptor) Distance range for event (discrete polar receptor) Radial direction for event (discrete polar receptor) Terrain elevation for event (optional) Receptor height above ground for event (optional) EventlncFile EventlncFile Identifies the filename for the included EVENT file, up to 200 characters in length; double quotes (") may be used as delimiters for the filename to allow for embedded spaces; and quotes don't count toward the limit of 200 EVENT locations can be input as either discrete Cartesian receptors (XR=, YR=) or as discrete polar receptors (RNG=, DIR=). Events that are specified in the file generated by the AERMOD model (CO EVENTFIL card) are always given as discrete Cartesian coordinates. Discrete polar receptors are assumed to be relative to an origin of (0,0) . B-24 ------- TABLE B-11 DESCRIPTION OF OUTPUT PATHWAY KEYWORDS OU Keywords STARTING RECTABLE MAXTABLE DAYTABLE MAXIFILE POSTFILE1 PLOTFILE1 TOXXFILE RANKFILE EVALFILE SEASONHR SUMMFILE FILEFORM EVENTOUT FINISHED Type M-N O-R O-R O-N O-R O-R O-R O-R O-R O-R O-R O-N O-N M-N M-N Keyword Description Identifies the start of OUTPUT pathway inputs Option to specify value(s) by receptor for output Option to summarize the overall maximum values Option to print summaries for each averaging period for each day processed. Option to list events exceeding a threshold value to file (if CO EVENTFIL option is used, these events are included in the input file generated for the EVENT model). Option to write results to a mass storage file for postprocessing. Option to write certain results to a storage file suitable for input to plotting routines Option to write results to a storage file suitable for input to the TOXX model component of TOXST or the RISK Option to output file of ranked values for Q-Q plots (must be used with the MAXTABLE keyword) Option to output file of normalized arc maxima from EVALCART receptors for model evaluation studies Option to output results by season and hour-of-day Option to output summary of high ranked values to separate file Specify fixed or exponential format for output results files Specifies the level of output information provided for EVENT Processing [EVENT Only] Identifies the end of OUTPUT pathway inputs 1) POSTFILE is used to output concurrent concentration values for particular source groups and averaging times across the receptor network suitable for postprocessing. PLOTFILE is used to output specific design values, such as second high concentrations, across the receptor network, suitable for plotting concentration contours. B-25 ------- TABLE B-12 DESCRIPTION OF OUTPUT PATHWAY KEYWORDS AND PARAMETERS Keyword RECTABLE where: MAXTABLE where: DAYTABLE where: MAXIFILE where: POSTFILE where: Parameters Aveper FIRST S Aveper 1ST 2> ECOND . . . SIXTH . . . TENTH or ID ... 6TH . . . 10TH Aveper FIRST SECOND SIXTH 1ST 2ND 6TH Averaging period to summarize with high values (kevword ALLAVE specifies all averaging periods) Select summaries of FIRST highest values by receptor Select summaries of SECOND highest values by receptor Select summaries of SIXTH highest values by receptor Select summaries of 1ST highest values by receptor Select summaries of 2ND highest values by receptor Select summaries of 6TH highest values by receptor Note: If two keywords are input separated by a dash (e.g. FIRST-THIRD) . then summaries of all high values in that range are provided. Also, if the CO EVENTFIL keyword is exercised, then the events generated by the RECTABLE keyword are included in the input file for EVENT model. Aveper Maxnum Aveper Maxnum Averaging period to summarize with overall maximum values (keyword ALLAVE specifies all averaging periods) Specifies number of overall maximum values to summarize Avperl Avper2 AvperS . . . Avperl Averaging period, e.g., 24 for 24-hr averages, to summarize with values by receptor for each day of data processed (keyword ALLAVE for first parameter specifies all averaging periods) Aveper Grpid Thresh Filnam (Funit) Aveper Grpid Thresh Filnam Funit Specifies averaging period for list of values equal to or exceeding a threshold value Specifies source group to be output to file Threshold value (e.g. NAAQS) for list of exceedances Name of disk file to store maximum values Optional parameter to specify the file unit Note: If the CO EVENTFIL keyword is exercised, then the events generated by the MAXIFILE keyword are included in the input file for EVENT processing. Aveper Grpid Format Filnam (Funit) Aveper Grpid Format Specifies averaging period to be output to file, e.g.. 24 for 24-hr averages, PERIOD for period averages Specifies source group to be output to file Specifies format of file, either UNFORM for B-26 ------- Keyword PLOTFILE where: TOXXFILE where: RANKFILE where: EVALFILE where: SEASONHR SUMMFILE FILEFORM Parameters Filnam Funit unformatted files or PLOT for formatted files for plotting Specifies filename for output file Optional parameter to specify the file unit Aveper Grpid Hivalu Filnam (Funit) (Short Term values) Aveper Grpid Filnam (Funit) (PERIOD or ANNUAL averages) Aveper Grpid Hivalu Filnam Funit Specifies averaging period to be output to file, e.g., 24 for 24-hr averages, PERIOD for period averages, etc. Specifies source group to be output to file Specifies hiah value summarv (e.e. FIRST. SECOND. 1ST. 2ND. etc.) to be output to file (must be selected on RECTABLE card) Specifies filename for output file Optional parameter to specify the file unit Aveper Cutoff Filnam (Funit) Aveper Cutoff Filnam Funit Specifies averaging period to be output to file, e.g., 1 for 1-hr averages. Specifies cutoff (threshold) value in g/m3 for outputting results for AERMOD model Specifies filename for output file Optional parameter to specify the file unit Aveper Hinum Filnam (Funit) Aveper Hinum Filnam Funit Specifies averaging period to be output to file, e.g., 24 for 24-hr averages Specifies the number of high values to be ranked Specifies filename for output file Optional parameter to specify the file unit Srcid Filnam (Funit) Srcid Filnam Funit Specifies the source ID to be output to file Specifies filename for output file Optional parameter to specify the file unit GroupID FileName (FileUnit) GroupID FileName (FileUnit) Specifies the source group ID to be output to file Specifies filename for output file Optional parameter to specify file unit SummFileName SummFileName EXP or FIX EXP FIX Specifies filename of output file Specifies that the output results files will use exponential-formatted values Specifies that the output results files will use fixed-formatted values B-27 ------- Keyword EVENTOUT where: Parameters SOCONT or DETAIL [EVENT Onlvl SOCONT DETAIL Provide source contribution information only in the event output Include hourly concentrations for each source and hourly meteorological data in the event output B-28 ------- APPENDIX C. LIST OF ERROR/WARNING MESSAGES This appendix provides a list of error, warning and informational messages used in the current version of the AERMOD model. Three types of messages can be produced by the model during processing of model inputs and during model calculations. These are described below: • Fatal Errors that will halt any further processing, except to identify additional error conditions (type E); • Warnings that do not halt processing but indicate possible errors or suspect conditions (type W); and • Informational messages that may be of interest to the user but have no direct bearing on the validity of the results (type I). A more detailed explanation of the error handling and reporting procedures used in AERMOD is provided Appendix C of the AERMOD User's Guide (EPA, 2004a). The three message types are identified with the letters E (errors), W (warnings), and I (informational messages). The 3- digit message numbers are generally grouped into categories corresponding to the different stages of the processing, although these distinctions may not always be obvious. Theses categories are: 100-199 Input Runstream Image Structure Processing 200 - 299 Parameter Setup Processing 300-399 Data and Quality Assurance Processing 400 - 499 Run Time Message Processing 500 - 599 Input/Output Message Processing The list provided below includes the message number, the main message text, and the "hint" field that may include additional details regarding the message. Note that the listing of messages in the AERMOD User's Guide (EPA, 2004a) is not up-to-date. Several additional messages have been added to AERMOD since that time, and the numbering of some messages may have changed. C-l ------- Message Number 100 105 110 115 120 125 130 135 140 141 142 143 144 145 146 147 149 150 152 154 155 156 157 158 159 160 170 175 180 185 189 190 191 192 Error/Warning Message Invalid Pathway Specified. The Troubled Pathway is Invalid Keyword Specified. The Troubled Keyword is Keyword is Not Valid for This Pathway. Keyword is STARTING or FINISHED Out of Sequence: Pathway = Pathway is Out of Sequence: Pathway = Missing FINISHED-Runstream File Incomplete: ISTAT= Missing Mandatory Keyword. The Missing Keyword is Nonrepeatable Keyword or Recursed INCLUDED: Keywrd Invalid Order of Keyword. The Troubled Keyword is Conflicting Options: PVMRM and OLM both specified Following Keyword Invalid Without PVMRM or OLM: Following Keyword Invalid Without PVMRM Option: Following Keyword Invalid Without OLM Option: Conflicting Options: MULTYEAR and Re-Start Option PSDGROUP Keyword Specified without PSDCREDIT Opt. Following Option is Invalid with PSDCREDIT Option: Conflicting options specified on MODELOPT keyword: Conflicting Options: MULTYEAR for Wrong Pollutant ELEVUNIT card must be first for this Pathway: Conflicting options: SCIM cannot be used with Conflicting Decay Keyword. Inputs Ignored for Option ignored - not valid with SCIM. Option = Wet SCIM Not Supported - Wet SCIM Inputs Ignored EMISUNIT Keyword Used With More Than 1 Output Type EMISUNIT Keyword Used With the Following Keyword: Duplicate ORIG Secondary Keyword for GRIDPOLR: Invalid Secondary Keyword for Receptor Grid: Missing Secondary Keyword END for Receptor Grid: Conflicting Secondary Keyword for Receptor Grid: Missing Receptor Keywords. No Receptors Specified. No Keywords for OU Path and No PERIOD/ANNUAL Aves. Incompatible Option Used With SAVEFILE or INITFILE PM-2.5 without MAXIFILE is incompatible with FASTALL option also implies use of FASTAREA option Hint Path Keyword Path Path Path Keyword Keyword Keyword Keyword Keyword Keyword Keyword Keyword Option Option Keyword Path ST AVEs Keyword DEPOS Keyword Keyword NETID NETID NETID NETID Keyword EVENTFIL C-2 ------- Message Number 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 Error/Warning Message Incompatible Keyword used with GASDEPVD option Gas deposition algorithms are non-DFAULT options METHOD_2 for participates is a non-DFAULT option TOXICS Option obsolete; see Users Guide Addendum Non-DFAULT BETA Option Required for Missing Parameter(s). No Options Specified For Not Enough Parameters Specified For the Keyword of Too Many Parameters Specified For the Keyword of Invalid Parameter Specified. Troubled Parameter: Regulatory DFAULT Conflicts with Non-DFAULT Option No Option Parameter Setting. Forced by Default to Regulatory DFAULT Overrides Non-DFAULT Option For No Parameters Specified. Default Values Will Used. Illegal Numerical Field Encountered in Negative Value Appears For Non-negative Variable. Number of Short Term Averages Exceeds Max: NAVE = Duplicate Averaging Period Specified for Keyword END Encountered Without (X,Y) Points Properly Set ELEV Input Inconsistent With Option: Input Ignored ELEV Input Inconsistent With Option: Defaults Used FLAG Input Inconsistent With Option: Input Ignored FLAG Input Inconsistent With Option: Defaults Used More Than One Delimiter In A Field for Keyword Number of (X,Y) Points Not Match With Number Of Number Of Receptors Specified Exceeds Max: NREC = Missing Origin (Use Default = 0,0) In GRIDPOLR Missing Distance Setting In Polar Network Missing Degree Or Dist Setting In Polar Network Missing Distance or Degree Field in Number of Receptor Networks Exceeds Max: NNET = Number of X-Coords Specified Exceeds Max: IXM = Number of Y-Coords Specified Exceeds Max: IYM = No Receptors Were Defined on the RE Pathway. Default(s) Used for Missing Parameters on Keyword Too Many Parameters - Inputs Ignored on Keyword Hint Keyword Keyword Option Keyword Keyword Keyword Parameter Option Option Keyword Variable Keyword NETID NETID NETID NETID NETID Keyword NETID NETID NETID Keyword Keyword C-2 ------- Message Number 230 231 232 233 234 235 236 237 239 240 241 242 243 244 245 246 247 248 249 250 252 254 256 259 260 262 264 265 266 270 275 280 281 282 Error/Warning Message Source ID field is too long (>8); first 8 chars: Too Many Numerical Values Specified for Number Of Specified Sources Exceeds Maximum: NSRC= Building Dimensions Specified for Non-POINT Source Too Many Sectors Input for Number of Source Groups Exceeds Maximum: NGRP = Not Enough BUILDHGTs Specified for SourcelD Not Enough BUILDWIDs Specified for SourcelD Not Enough QFACTs Specified for SourcelD Inconsistent Number of Particle Categories for Not Enough BUILDLENs Specified for SourcelD No Particle Cat. or Gas Depos. Specified for SRCID Wet depos (DEPOS, WDEP, WETDPLT) incompatible with Source parameters are missing or incomplete for No. of Particle Categories Exceeds Max: NPDMAX = Not Enough XBADJs Specified for SourcelD Not Enough YBADJs Specified for SourcelD No Sources Were Defined on the SO Pathway. Source elevation is missing (-9999.0); SRCID = Duplicate XPNT/DIST or YPNT/DIR Specified for GRID Duplicate Receptor Network ID Specified. NETID = Number of Receptor ARCs Exceeds Max: NARC = EVALFILE Option Used Without EVALCART Receptors Receptor elevation is missing (-9999.0); IREC = Number of Emission Factors Exceeds Max: NQF = First Vertex Does Not Match LOCATION for AREAPOLY Too Many Vertices Specified for AREAPOLY Source Not Enough Vertices Specified for AREAPOLY Source Invalid shape defined (area=0) for AREAPOLY source Number of High Values Specified Exceeds Max: NVAL= Number of Max Values Specified Exceeds Max: NMAX= Number of Output Types Specified Exceeds Max:NTYP= Number of OLMGROUPs Specified Exceeds Max: NOLM = Following SRCID Included in Multiple OLMGROUPs: Hint SRCID(1:8) EMISFACT SRCID Keyword SRCID SRCID SRCID SRCID SRCID SRCID GASDEPVD DRYDEP or WETDEP SRCID SRCID SRCID NETID NETID NUMARC=0 Rec Number SRCID SRCID SRCID SRCID SRCID C-4 ------- Message Number 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 310 313 315 317 318 319 320 322 324 325 330 332 334 335 Error/Warning Message Either OZONEVAL or OZONEFIL Card Needed for Option Invalid POLLUTID Specified for PVMRM/OLM; Must Use Number of Urban Areas Exceeds Maximum. NURB = Following SRCID Included in Multiple PSDGROUPs: PSDGROUP ID Must be INCRCONS, RETRBASE or NONRBASE Use of "*" for repeated values not meaningful for Source defined as both particulate and gaseous Number of Events Specified Exceeds Max: NEVE = Filename specified is too long. Maximum length = Format specified is too long. Maximum length = User-specified met data format not used; use FREE PERIOD and ANNUAL averages are both selected for Invalid Averaging Period Specified for SCREEN Mode Averaging Period .NE. 1-Hr for TOXXFILE Option Aver. Period must be .LE. 24 for EVENT Processing Error Allocating Storage for Setup Arrays! Error Allocating Storage for Result Arrays! Specified SRCID Has Not Been Defined Yet: KEYWORD= Urban Area ID Has Not Been Defined. URBID = Following SRCID Included in Multiple Urban Areas: Urban ID has already been defined. URBID = Attempt to Define Duplicate LOCATION Card for SRC: Attempt to Define Duplicate EVENTPER card for Attempt to Define Duplicate SRCPARAM Card for SRC: Specified SRCID not included in any PSD/SRCGROUP: No Sources Defined for Urban Area. URBID = No Sources Included in Specified Source Group: Input Parameter May Be Out-of-Range for Parameter Release Height Exceeds Effective Depth for OPENPIT Release Height Exceeds 3000 Meters for SRCID: Negative Exit Velocity (Set=1.0E-5) for SRCID: Mass Fraction Parameters Do Not Sum to 1 . for Src Mass Fraction Parameter Out-of-Range for Source Particle Density Out-of-Range for Source Particle Diameter Out-of-Range for Source Hint OLM or PVMRM NO2 SRCID Keyword SRCID ILEN FLD ILEN FLD Format Keyword Ihr Only EVNAME Error Code Error Code Keyword Urban ID SRCID Urban ID SRCID EVNAME SRCID SRCID Urban ID GRPID Parameter SRCID SRCID SRCID SRCID SRCID SRCID SRCID C-5 ------- Message Number 336 338 340 342 344 345 346 350 352 353 354 360 363 365 370 380 381 382 384 386 387 391 392 395 397 398 399 400 405 406 410 413 420 430 432 Error/Warning Message NO2RATIO Invalid or Not Specified for PVMRM Source Neg Emis Rate Cannot be Used with OLM/PVMRM. Src: Possible Error in PROFBASE Input: Value is < 0 Src ID Mismatch in Hourly Emissions File for ID = Missing HOUREMIS fields; EmisRate set = 0. KURDAT= Problem processing the HOUREMIS file. KURD AT = Too many fields for HOUREMIS file . KURDAT = Julian Day Out Of Range at Missing Field on MULTYEAR Keyword for Urban Roughness Length (m) May Be Out-of-Range: High-8th-High Only Required for PM-2.5 24-hr Ave 2-Digit Year Specified: Valid for Range 1950-2049 24HR and ANNUAL Averages Only for PM-2.5 NAAQS Year Input is Greater Than 2147 Invalid Date: 2/29 In a Non-leap Year. This Input Variable is Out-of-Range: Latitude in Surface File Is Not Valid: Error Decoding Latitude: Not enough fields specified for HOUREMIS; KURDAT = PARTDIAM and METHOD_2 specified for same SRCID: METHOD_2 option already specified for this SRCID: Aspect ratio (LAV) of area source greater than 100 Aspect ratio (LAV) of open pit is greater than 10 Met. Data Error; Incompatible Version of AERMET: SCREEN option used without use of SCREEN Met Data SCREEN met used without specifying SCREEN option EXP format specified with no applicable file types Output values exceed format limit; use OU FILEFORM Value of PHEE Exceeds 1.0 on KURDAT = Numer of Vertices Exceeds Maximum (NVMAX) for Src: Wind Direction Out-of-Range. KURDAT = Number of Threshold Events > 9999 for Ave Period Wind Speed Out-of-Range. KURDAT = Ambient Temperature Data Out-of-Range. KURDAT = Friction Velocity Out-of-Range. KURDAT = Hint SRCID SRCID Keyword SRCID YYMMDDHH YYMMDDHH YYMMDDHH JDAY H6H URBZO RECTABLE Keyword AVERTIME Keyword or YR = YR = Lat Lat YYMMDDHH SRCID SRCID SRCID SRCID Version Date Version Date FILEFORM = EXP YYMMDDHH SRCID YYMMDDHH YYMMDDHH YYMMDDHH YYMMDDHH C-6 ------- Message Number 435 438 439 440 441 442 450 455 456 457 458 459 460 465 470 474 475 480 481 483 484 485 486 488 489 490 495 496 499 500 510 520 Error/Warning Message Surface Roughness Length Out-of-Range. KURD AT = Convective Velocity Data Out-of-Range. KURD AT = Monin-Obukhov Length Out-of-Range. KURD AT = Calm Hour Identified in Meteorology Data File at Vert Pot Temp Grad abv ZI set to min .005, KURDAT= Vert Pot Temp Grad abv ZI exceeds 0.1 K/m, KURDAT= Record Out of Sequence in Meteorological File at: Date/time Mismatch: Hourly Emission File. KURD AT = Date/time Mismatch on Surface & Profile. KURD AT = Date/time Mismatch on OZONEFIL Data. KURD AT = Substitution made for missing ozone data. KURD AT = Missing ozone data; Full conversion used. KURD AT = Missing Hour Identified in Meteor. Data File at Number of Profile Levels Exceeds Max: MXPLVL = Mixing Height Value is < or = 0.0. KURD AT = WS RefHt invalid (<0.001); Not msg or elm: KURDAT= WS reference height is higher than 100m. KURD AT = Less Than 1 Year Found for ANNUAL Averages Data Remaining After End of Year. Number of Hours= User Start Date is Earlier Than Start of Met File Restart Date < STARTEND date or start of met file MULTYR DataGap; Restart Date < STARTEND or MetFile MULTYR Date Overlap; STARTEND Date < Restart Date First met HR.ne. 1; ST results may not be valid for First met HR.ne. 1; EV results may not be valid for Problem reading SURFFILE date for EVENTS; MNDYHR = Surface met file does not include enough variables Total precipitation in SURFFILE is zero (0.0) with PRIME plume rise error; check stack parameters for Fatal Error Occurs Opening the Data File of Fatal Error Occurs During Reading of the File of Fatal Error Occurs During Writing to the File of Hint YYMMDDHH YYMMDDHH YYMMDDHH YYMMDDHH YYMMDDHH YYMMDDHH YYMMDDHH YYMMDDHH YYMMDDHH YYMMDDHH YYMMDDHH YYMMDDHH YYMMDDHH YYMMDDHH YYMMDDHH YYMMDDHH YYMMDDHH STARTEND or YYMMDDHH STARTEND or YYMMDDHH YYMMDDHH First Day MnDyHr MnDyHr with-DEP or Non-DEP WetDepos SRCID FileName or FileType FileType FileType and Unit C-7 ------- Message Number 530 540 550 555 560 565 570 580 585 590 Error/Warning Message CAUTION! Met Station ID Mismatch with SURFFILE for No RECTABLE/MAXTABLE/DAYTABLE for Average Period File Unit/Name Conflict for the Output Option: File Unit/Name conflict across options: GRP# AVE User Specified File Unit .LT. 30 for OU Keyword: Possible Conflict With Dynamically Allocated FUNIT Problem Reading Temporary Event File for Event: End of File Reached Trying to Read the File of Output data file for INITFILE option was not found The INITFILE filename matches a SAVEFILE filename Hint Keyword nnn-HR Keyword GrpNum and AvePer Keyword Keyword EVNAME Keyword File Type ------- APPENDIX D EPA Model Clearinghouse Memorandum Dated July 9,1993 D-l ------- July 9, 1993 MEMORANDUM SUBJECT: Proposal for Calculating Plume Rise for Stacks with Horizontal Releases or Rain Caps for Cookson Pigment, Newark, New Jersey FROM: Joseph A. Tikvart, Chief Source Receptor Analysis Branch, TSD (MD-14) TO: Ken Eng, Chief Air Compliance Branch, Region II In response to your request, the Model Clearinghouse has reviewed your proposal for treating horizontal and capped stacks at Cookson Pigment so that the model (SCREEN or ISC2) will properly treat plume rise from the Cookson Pigment stacks. We concur in principle with the approach, with some relatively minor changes. First, the analysis provided by New Jersey Department of Environmental Protection is technically correct. We suggest, however, that the exit velocity for horizontal and capped stacks be set to a lower figure than 0.1 m/s. A 0.1 m/s exit velocity may still result in significant momentum plume rise being calculated, even though these kinds of sources should have zero momentum rise. We therefore suggest setting the stack exit velocity to a lower value, such as .001. For horizontal stacks that are not capped, we suggest turning stack tip downwash off, whether there are buildings or not. Stack tip downwash calculations are inappropriate for horizontal stacks. For vertical stacks that are capped, turn stack tip downwash off and reduce the stack height by three times the actual stack diameter. The cap will probably force stack tip downwash most of the time. The maximum amount of the stack tip downwash (as calculated in ISC2) is three times the stack diameter. Reducing the stack height by this amount, while turning off the stack tip downwash option, causes the maximum stack tip downwash effect. The resulting concentrations may err slightly on the high side. For stacks with small diameters, such as those at Cookson Pigment, the error should be quite small. Note, however, that this approach may not be valid for large diameter stacks (say, several meters). cc: A. Colecchia D. Wilson D-2 ------- |