## Atlanta, Georgia Prepared by 

# July 2008 

# Chemical Products Corporation BART-Eligible Source

DRAFT

Cartersville, Georgia (Bartow County) 

# CALPUFF BART-Exclusion Modeling Protocol (4-km Grid)Chemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol

P:\Chemical Products\BART\2008 BART\VISTAS\Chemical Products VISTAS BART Protocol.DRAFT.GAEPD.Rev2.doc 

# DRAFT

Chemical Products Corporation VISTAS BART CALPUFF (4-km) Modeling Protocol

Table of Contents Page

1. Introduction 2 1.1 Objectives 3 1.2 Location of Source vs. Relevant Class I Areas 3 1.3 Source Impact Evaluation Criteria 4


### 2. Chemical Products Corporatio n BART-Eligible Source Description 4 2.1 CPC BART-Eligible Source Emissions Units 4 2.2 CPC BART-Eligible Source Model Emissions Inventory 5


### 3. Geophysical and Meteorological Data 6 3.1 Modeling Domain, Terrain, and Meteorological Dataset 6 3.2 Air Quality Database 7 3.3 Natural Conditions at Class I Areas 7


### 4. Air Quality Modeling Methodology 8 4.1 Plume Model Selection 8 4.2 CALPUFF Computational Domain and Receptors 9 4.3 CALPUFF Model Options Selection 9 4.4 POSTUTIL and CALPOST Processing for Light Extinction and Haze Impact Calculations 9


### 5. BART Exemption Modeling Results Review 11 (To be completed after approval of modeling protocol)


### Appendix A Chemical Products Corporation BART-Eligible So urce Potential-to-Emit Calculations, Maximum 24- hr Actual Emissions Rates, an d Modeled Emissions Inventory


### Appendix B Example of Proposed CALPUFF Input File


### Appendix C Example of Proposed POSTUTIL Input File



### Appendix D Example of Proposed CALPOST Input File



### Appendix E Example of CALPOST List File Output using the Existing and New IMPROVE Algorithm for Calculating Visibility Impairment at Class I Area s within 300km of Chemical Products CorporationChemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol

P:\Chemical Products\BART\2008 BART\VISTAS\Chemical Products VISTAS BART Protocol.DRAFT.GAEPD.Rev2.docPage 2 

# DRAFT

1. Introduction

Chemical Products Corporation (“CPC”) operates an inorganic chemical manufacturing plant in Cartersville, Georgia and prim arily processes barite ore (barium sulfate) and celestite ore (strontium sulfate) into finished products su ch as barium carbonate, barium chloride, strontium carbonate and strontium nitrate. These chemicals are produced using a number of different types of heated process equipment such as reduct ion kilns and various types of dryers. CPC also owns a furnace that has beenused in the production of sodium silicate glass. Along with the primary products, CPC produces several secondary products such as barium metaborate, elemental sulfur, sodium hydrosulfide, sodium sulfide, ammonium sulfide and dry calcium carbonates. CPC is a major so urce with respect to the Part 70 (Title V) operating permit program and operates underGeorgia Title V Permit No. 2819-015-0008-V02-0, as amended May 1, 2007 (V-02-1).

In 1999, the US EPA published a final rule to a ddress a type of visibility impairment known as “regional haze” (64 FR 35714) in 156 Federally-protected national parks and wilderness areas (“mandatory Class I Federal Areas”, or just “Class I Areas”). The Regional Haze Rule requires the evaluation of the best available retrofit technology (BART) for certain major stationary sources defined as “BART-eligible” which are those sources which have the potential-to-emit 250 tons per year (tpy) of on e or more visibility-impairing air pollutants, were “in existence” before August 7, 1977 but began operation after August 7, 1962, and whose operations fall within one or more of 26 specific source categories; chemical process plants are one of the listed source categories, and include process units and support facilities that are considered to be part of major Standard Industrial Classification 28 – Chemicals and Allied Products. BART is required for “BART-elig ible sources” if is it determined that the source emits any air pollutant which may reasonab ly be anticipated to cause or contribute to any impairment of visibility in a Class I Area.

According to US EPA guidelines for BART determinations under the Regional Haze Rule (70 FR 39104), dispersion modeling can be used todetermine if an individual source can reasonably be anticipated to cause or contribute to visibility impairment in a Class I area and is thus not subject to BART. A BART-eligible source that is responsible for a 0.5 deciview (dv) change or more from natural background conditions is considered to “contribute” to visibility impairment in a Class I area while a BART-eligible source that is responsible for a 1.0 dv change or more natural background co nditions is considered to “cause” a visibility impairment in a Class I area. The Visibility Improvement State and Tribal Association of the Southeast (VISTAS) requires the application of the CALPUFF modeling system for carrying out air quality modeling in support of BART determinations using the maximum actual 24-hr emissions rates of NOx, SO2, and PM10. The VISTAS application of CALPUFF for BART determinations uses a two-tier approach: (1 ) initial exemption modeling using 12-km grid CALMET data and (2) finer resolution exclusion modeling using 4-km grid or smaller CALMET data.

Chemical Products Corporation has six emission units that comprise the BART-eligible source because the units operate at a chemical processi ng facility, were in existence on August 7, 1977, and began operation after August 7, 1962 an d have the potential-to emit greater than 250 tpy of at least one visibility-impairing polluta nt. To avoid the requirements of BART, CPC has a federally enforceable permit limitation on SO2 emissions of less than 250 tpy during any 12 consecutive month period from the six BART-eligible emissions units unless exemption modeling conducted in accordance with 40 CFR 51 .303(a)(2) demonstrates that it cannot be reasonably be anticipated to cause or contributeto visibility impairment in a Class I Area. Periodic monitoring, recordkeeping and report ing requirements required to demonstrateChemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol

P:\Chemical Products\BART\2008 BART\VISTAS\Chemical Products VISTAS BART Protocol.DRAFT.GAEPD.Rev2.docPage 3 

# DRAFT

compliance with the 250 tpy SO2 BART avoidance limitation are allowed to be voided if Chemical Products Corporation can demonstrate that it is exempt from the requirements of BART and is thus the purpose of this protocol.

Initial 12-km grid BART exemption modelin g was conducted on behalf of CPC through funding provided by the Visibility Improvement State and Tribal Association of the Southeast (VISTAS). Now, CPC is electing to conduct fin er 4-km grid BART exclusion modeling on its own, in conjunction with Smith Aldridge, Inc. The source configuration, pollutant speciation, Class I Area receptors, ozone data sets and CALPUFF model options for finer grid exclusion modeling are essentially the same as the initia l exemption modeling. In addition to the better grid resolution and the introduction of MM5 observational data, several other modeling refinements enhance the accuracy of the finer 4-km grid modeling, including the use of the higher resolution terrain data in defining th e terrain fields and application of the Ammonia Limiting Method in the POSTUTIL post-processor.

1.1. Objectives . The purpose of this document is to(a) establish a modeling protocol that describes the procedure that will be used for the application of the CALPUFF modeling system for finer grid (4-km) BART exclusion modeling for the Chemical Products Corporation BART-eligible source and (b) to determine if the Chemical Products Corporation BART-eligible source emits any air pollutant(s) in an amount which may reasonably be anticipated to contribute to a visibility impairment in a Class I Area. The following protocol is based on VISTAS BART common modeling protocol, revised August 31, 2006.

1.2. Location of Source vs. Relevant Class I Areas . The VISTAS BART Modeling Protocol specifies that all Class I Areas within 300 km of a BART-eligible source must be initially evaluated to determine whether the source co ntributes to visibility impairment. Figure 1-1 and Table 1-1 summarize the distances separating the CPC BART-eligible source from mandatory Class I Federal Areas within 300 km.

300 km

Joyce-Kilmer Wilderness

Cohutta Wilderness Shining Rock Wilderness

Great Smokey Mountains National Park

Sipsey Wilderness

Chemical Products Corporation Bartow County, Georgia

Figure 1-1 Class I Areas within 300 km of Chemical Products Corporation (Latitude 34 °09’ 03” N, Longitude 84 °46’ 09” W)Chemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol

P:\Chemical Products\BART\2008 BART\VISTAS\Chemical Products VISTAS BART Protocol.DRAFT.GAEPD.Rev2.docPage 4 

# DRAFT

Table 1-1 Distances Separating CPC and Class I Areas within 300 km

Cohutta Wilderness Area Great Smokey Mountains National Park Joyce-Kilmer Wilderness Area Sipsey Wilderness Area Shining Rock Wilderness Area

78.9 km 167.0 km 153.4 km 237.3 km 217.4 km

Distances determined from Chemical Products Corporation (Latitude 34 °09’ 03” N, Longitude 84 °46’ 09” W) to closest receptor in each Class I Area.

1.3. Source Impact Evaluation Criteria

.

1.3.1. BART Exclusion Modeling (4-km Grid) . For 4-km grid BART exclusion modeling, a BART-eligible source is considered to ca use or contribute to visibility impairment in a Class I Area if the 98th percentile 24-hr average predicted change in the haze index over a three-year period exceeds 0.5 dv from natural background conditions. The 98th percentile 24-hr average predicted impact is the greatest of the three annual highest 8th-high or highest 22nd-high predicted change over the entire three year period, whichever is greater, per EPA BART guidelines and the VISTAS BART common modeling protocol.

Haze is caused when sunlight encounters ti ny pollution particles in the air. Some light is absorbed by particles and other lig ht is scattered away before it reaches an observer. More pollutants mean more absorption and scattering of light, which reduce the clarity and color of what we see. Some types of particles such as sulfates, scatter more light, partic ularly during humid conditions. The predicted change in haze is measured bycalculating natural light extinction ( b ext

) which is used to represent the haze index ( HI ), expressed in dv. The haze index is calculated using the following equation,

, 10 ln10 ⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛

= ext b HI



## ( )( ) The change in light extinction is affected by various chemical species and the Rayleigh scattering phenomenon and can be calculated using the IMPROVE methodology as shown in the following equation,



## [ ] ( )

## [ ] [ ] [ ] [ ] [ ]

{ { Scattering Rayleigh ray

Carbon Elemental Filterable

PM PMCoarse Filterable

PM PMFine Filterable PM Organic Condesable Nitrates Ammonium 34

Sulfates Ammonium 4 2 4 ext bEC10PMC0.6Soil OC4NONHRH3fSONHRH3fb

2.510 2.5 ++ + + + + =

− 32143 42 1 321

44 4 3

44 4 2 1

44 4 3

44 4 2 1

The concentrations, in square brackets, are in micrograms per cubic meter (ug/m 3 ) and b ext

is in units of inverse megameters (Mm -1 ). The Rayleigh scattering term ( b ray

) has a default value of 10 Mm -1 , as recommended in EPA guidance for tracking reasonable progress.



### 2. Chemical Products Corporation BA RT-Eligible Source Description

2.1. CPC BART-Eligible Source Emissions Units . Table 2-1 provided a summary of the applicability criteria for the six emissions unitsthat comprise the BART-eligible source atChemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol

P:\Chemical Products\BART\2008 BART\VISTAS\Chemical Products VISTAS BART Protocol.DRAFT.GAEPD.Rev2.docPage 5 

# DRAFT

Chemical Products Corporation. Potential emissions for NOx, SO2 and PM10 aredetermined in data provided in Table A-1 of Appendix A to this document, using worstcase fuel for each unit, Georgia Rule (g) for fuel oil sulfur content, and Georgia Rules (d) and (e) for particulate matter.

Table 2-1 CPC BART-Eligible Source Emissions Units BART Applicability Summary

Emissions Unit Source Code Construction Date Latitude Longitude Potential SO2, (tpy) Potential NOx, (tpy) Potential PM10, (tpy)
South Rotary Calciner BC01 1968 34 °9’ 3.40”N 84 °47’ 10.29 W 278.33 16.08 56.27

Rotary Dryer BD01 1964 34 °9’ 3.37”N 84 °47’ 9.52 W 102.62 5.96 45.46

South Spray Dryer BD03 1968 34 °9’ 2.87”N 84 °47’ 10.25 W 68.42 3.92 41.57

Barium Chloride Rotary Dryer BD04 1976 34 °9’ 3.10”N 84 °47’ 10.61 W 52.87 3.02 37.49

East Rotary Kiln BK02 1966 34 °9’ 2.32”N 84 °47’ 7.33 W 1444.66 18.31 52.79

B&W Process Steam Boiler GBX1 1973 34 °9’ 1.78”N 84 °47’ 6.62 W 1391.95 90.95 69.25

2.2. CPC BART-Eligible Source Model Emissions Inventory . An emissions inventory of modeled source paramete rs is provided in Tables A-2 and A-3 of Appendix A to this document. All modeled source parameters are consistent with those previously modeled by VISTAS, excluding the maximum 24-hr emissions rates modeled. The maximum actual 24-hr emissions used for the 4-km finer grid exclusion modeling for NOx and SO2 were determined based on the combustion of fuel oil with 0.5% sulfur by weight and use of emissions factors published by US EPA in AP-42 ( http://www.epa.gov/ttn/chief/ap42/ch01/final/c01s03.pdf ). The maximum actual 24-hr emission rate of SO2 from the East Rotary Kiln includes an emissions factor derived for the processing of petroleum coke at one ton per hour with a maximum 7.50% sulfur content by weight assuming a conversion factor of sulfur to SO2 of two during processing. Particulate matter emissions rates were determined from Georgia Rules (d) and (e), as applicable, and used as a surrogate for PM10. For PM10 speciation, all PM10 was assumed to be condensable organic particulate matter in the submicron particle size categories (PM081 and PM056) consistent with the guidance provided by VISTAS for BART modeling. Table 2-2 summarizes th e maximum 24-hr actual emissions rates proposed to be modeled.

Table 2-2 CPC BART-Eligible Source Maximum Actual 24-hr Emissions Rates

Emissions Unit Source Code SO2 (lb/hr) NOx (lb/hr) PM10 (lb/hr) PM081 (lb/hr) PM056 (lb/hr)

South Rotary Calciner BC01 12.709 3.580 12.848 6.424 6.424Chemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol

P:\Chemical Products\BART\2008 BART\VISTAS\Chemical Products VISTAS BART Protocol.DRAFT.GAEPD.Rev2.docPage 6 

# DRAFT

Table 2-2 CPC BART-Eligible Source Maximum Actual 24-hr Emissions Rates

Emissions Unit Source Code SO2 (lb/hr) NOx (lb/hr) PM10 (lb/hr) PM081 (lb/hr) PM056 (lb/hr)

Rotary Dryer BD01 4. 686 1.320 10.379 5.190 5.190

South Spray Dryer BD03 3.124 0.880 9.491 4.745 4.745

Barium Chloride Rotary Dryer BD04 2.414 0.680 8.560 4.280 4.280

East Rotary Kiln BK02 306.97 4.180 12.053 6.026 6.026

B&W Process Steam Boiler GBX1 52.966 14.92 15.811 7.906 7.906



### 3. Geophysical and Meteorological Data

3.1. Modeling Domain, Terrain, and Meteorological Data Sets . VISTAS has developed 12-km coarse domain and five sub-regional 4-km CALMET meteorological databases for three years (2001-2003). The sub-regional modelin g domains are strategically designed to cover all potential BART eligible sources wi thin VISTAS states and all Class I Areas within 300 km of those sources. The exte nts of the 4-km sub-regional domains are shown in Figure 3-1.

Figure 3-1 VISTAS 4-km Sub-Regional DomainsChemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol

P:\Chemical Products\BART\2008 BART\VISTAS\Chemical Products VISTAS BART Protocol.DRAFT.GAEPD.Rev2.docPage 7 

# DRAFT

USGS 90-meter Digital Elevation Model (DEM) files were used by VISTAS to generate the terrain data at 4-km resolution for input to the 4-km sub-regional CALMET run. Likewise, USGS 90-meter Composite Theme Grid (CTG) file s were used by VISTAS to generate the land use data at 4-km resolution for input tothe 4-km sub-regional CALMET run. Three years of MM5 data (2001-2003) were used by VISTAS to generate the 4-km sub-regional meteorological datasets.

3.2. Air Quality Database . The CALPUFF model is capable of simulating linear chemical transformation effects by using pseudo-first-order chemical reaction mechanisms for the conversions of SO2 to SO4, and NOx, which consists of nitrogen oxide (NO) and nitrogen dioxide (NO2), to nitrate (NO3) and nitric acid (HNO3). In this application, chemical transformations involving five species (SO2, SO4, NOx, HNO3, and NO3) will be modeled using the MESOPUFF II chem ical transformation scheme. Ambient concentrations of ozone and ammonia concentrations as represented in the model affect the MESOPUFF II chemical transformation simulation.

3.2.1. Ozone Concentrations . Ambient ozone concentrations are used in the CALPUFF “MESOPUFF II” chemical transformation scheme as a surrogate for OH radicals throughout the daylights hours during SO2 and NOx oxidation. Hourly measurements of ozone from all non-urba n monitors within and just outside the computational grid shall be used as inpu t to CALPUFF. The model default of 80 ppb was used for the background ozone concentration in the instance when hourly data was missing for each station. The standard ozone data files for 2001-2003 will be provided by Peter Courtney of the GA EPD Data and Modeling Unit.

3.2.2. Ammonia Concentrations . For the background ammonia value, a constant background value of 0.5 ppb shall be used rather than using ammonia data derived from CMAQ model outputper recommendation by VISTAS.

3.3. Natural Conditions at Class I Areas . For the Class I Areas within 300 km of the CPC BART-eligible source, Table 3- 1 summarizes the default natural background conditions as tabulated in Appendix B of US EPA’s Guidance for Estimating Natural Visibility Conditions under the Regional Haze Rule. To represent natural background conditions in the absence of anthropogenic sources of sulfates and nitrates, the VISTAS BART Modeling Protocol expresses monthly background extinction coefficients in terms of Rayleigh scattering (default value of 10 Mm -1for Class I Areas) and scattering due to soils (i.e., fine particles) only and is calculated from the tabulated haze index (HI) values from the following equation,



## ( ) 

## [ ]

{ { 

## [ ] 1

Scattering Rayleigh ray

PM PMFine Fiterable 10 HI
back Mm 10Soil bSoil 10eb

2.5 − +=+= =



## [ ] [ ] ( ) 1 10 HI

3 back Mm 01 10eSoils
m μ g b

−

− == ⎟ ⎠ ⎞ ⎜ ⎝ ⎛Chemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol

P:\Chemical Products\BART\2008 BART\VISTAS\Chemical Products VISTAS BART Protocol.DRAFT.GAEPD.Rev2.docPage 8 

# DRAFT

Table 3-1 Natural Background Conditions for Class I Areas

Class I Area b ext (Mm -1 ) Annual Average Haze Index (dv) 20% Best Days Haze

Index (dv) 20% Worst Days Haze

Index (dv) Soil b back Annual Avg.

( μ g/m 3 ) Soil b back 20% Best Days Avg. ( μ g/m 3 )

Cohutta Wilderness Area 21.39 7.60 3.76 11.44 11.38 4.56

Great Smokey Mountains National Park 21.39 7.60 3.76 11.44 11.38 4.56

Joyce-Kilmer Wilderness Area 21.40 7.61 3.77 11.45 11.40 4.58

Sipsey Wilderness Area 21.28 7.55 3.71 11.39 11.28 4.49

Shining Rock Wilderness Area 21.40 7.61 3.77 11.45 11.40 4.58

The effects of relative humidity to amp lify the visibility impairment of hygroscopic sulfates and nitrates were characterized usin g “Method 6,” which computes the change in light extinction using a monthly average relative humidity adjustment particular to each Class I Area applied to background and modeled sulfate and nitrate. Table 3-2 summarizes the monthly average humidity values that were applied for the Class I Areas considered in this analysis, as tabulated in Table A-3 of US EPA’s Guidance for Estimating Natural Visibility Conditions under the Regional Haze Rule.

Table 3-2 Monthly Site-Specific Average 

### f (RH) for Class I Areas

Class I Area Jan Feb Mar April May June July Aug Sept Oct Nov Dec

Cohutta Wilderness Area 3.3 3.1 3.0 2.8 3.4 3.8 4.0 4.2 4.2 3.8 3.4 3.5

Great Smokey Mountains National Park 3.3 3.0 2.9 2.7 3.2 3.9 3.8 4.0 4.2 3.8 3.3 3.4

Joyce-Kilmer Wilderness Area 3.3 3.1 2.9 2.7 3.3 3.8 4.0 4.2 4.2 3.8 3.3 3.5

Sipsey Wilderness Area 3.4 3.1 2.9 2.8 3.3 3.7 3.9 3.9 3.9 3.6 3.3 3.4

Shining Rock Wilderness Area 3.3 3.0 2.9 2.7 3.4 3.9 4.1 4.5 4.4 3.8 3.3 3.4



### 4. Air Quality Modeling Methodology

4.1. Plume Model Selection . CALPUFF has been adopted by the U.S. EPA as a “Guideline” model for source-receptor distances greater than 50 km (i.e. long range transport), and for use on a case-by-case basis in complex flow situations for shorter distances. CALPUFF is also recommended for Class I impact assessments by FLAG and IWAQM. The final BART guidance recommends CALP UFF as “the best modeling application available for predicting a single source’s co ntribution to visibility impairment.” As aChemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol

P:\Chemical Products\BART\2008 BART\VISTAS\Chemical Products VISTAS BART Protocol.DRAFT.GAEPD.Rev2.docPage 9 

# DRAFT

result of these recommendations, the VISTAS BART Modeling Protocol is based on the use of CALPUFF modeling system (Version 5.754) for its BART determinations.

4.2. CALPUFF Computational Domain and Receptors . The CPC BART-eligible source modeling shall use the meteorological dataset representing VISTAS sub-regional domain 4 and shall be provided by GA EPD. The computational domain shall be a subset of the VISTAS regional domain with the size of the domain selected to encompass the CPC BART-eligible source and each Class I Area, plus a minimum 50 km buffer in all directions to allow for the puff recirculat ion beyond CPC and the Class I Areas being evaluated. The receptor grid for each Class I Area shall be those receptor locations and elevations provided by the National Park Service (NPS) Air Resources Division (ARD) available at http://www.nature.nps.gov/air/maps/receptors/index.cfm . Table 4-1 summarizes the extents of the computational domains that shall be used for each Class I Area.

Table 4-1 Computational Domains used for CALPUFF for each Class I Area
CALPUFF Computation Grid Variables Class I Area IBCOMP JBCOMP IECOMP JECOMP Approximate Model Buffer

Cohutta Wilderness Area 76 136 127 208 + 100 km

Great Smokey Mountains National Park 76 149 136 230 + 100 km

Joyce-Kilmer Wilderness Area 76 140 136 224 + 100 km

Sipsey Wilderness Area 15 125 133 186 + 100 km

Shining Rock Wilderness Area 76 165 136 227 + 100 km

4.3. CALPUFF Modeling Options Selection . The CALPUFF modeling shall be conducting in accordance with the IWAQM Phase 2 Su mmary Report and Recommendations for Modeling Long Range Transport Impacts (E PA-454/R-98-019, December 1998), except as modified by the VISTAS BART Modeling Protocol (Section 4.3.3, Revision 3.2, August 2006). An example of the proposed CALPUFF input file is provided in Appendix B describing the CALPUFF model options se lection for Cohutta Wilderness Area for meteorological year 2001.

4.3.1. Building Downwash. Building downwash shall be ignored since the CPC BARTeligible source is located at a distance greater than 50 km from the nearest Class I Area.

4.3.2. Ozone Dataset . Observed hourly ozone data from CASTNet and AIRS stations shall be used. A default of 80 ppb sha ll be used for the background ozone concentration in the instance when hourly data was missing for each station.

4.3.3. Background Ammonia Concentration . A constant value of 0.5 ppb shall be used for background ammonia concentration.

4.4. POSTUTIL and CALPOST Processing for Light Extinction and Haze Impact Calculations

. Examples of the proposed POSTUTIL and CALPOST input files are provided in Appendices C and D, respectively, which describe the post-processing proceduresChemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol

P:\Chemical Products\BART\2008 BART\VISTAS\Chemical Products VISTAS BART Protocol.DRAFT.GAEPD.Rev2.docPage 10 

# DRAFT

proposed to be used determining the chan ge in light extinction as measured at receptors in the Cohutta Wilderness Area for meteorological year 2001.

4.4.1. Species Considered . Species considered shall include sulfur dioxide, oxides of nitrogen and particulate matter with an aerodynamic diameter less than 10 microns (SO2, NOx, and PM10). Speciation of PM10 shall be performed in accordance with “Procedure for Speciation of Emissions for VISTAS BART Modeling”, July 18, 2006. In the absence of continuous emissions monitoring data and specific particulate matter AP-42 factors relevant to the source classification code for each emissions units, Georgia Rules (d) and (e) particulate matter permitted emissions limits shall beused as a surrogate for PM10 which are in turn speciated evenly into condensable organic (OC) particulate matter in the submicron size categories PM081 (PM0.625-1) and PM056 (PM0.5-0.625). PM081 and PM056 are modeled explicitly in CALPUFF and then summed in POSTUTIL to represent secondary organic aerosols using the following input file modification;

****** INPUT GROUP: 1 – General run control parameters

“Number of species to compute from those modeled”,

! NSPECCMP=1 !,

SUBGROUP (2b)

“The following NSPECOUT species will be written”,

! ASPECO = PM081 ! ! ASPECO = PM056 ! ! ASPECO = SOA !

SUBGROUP (2c)

“The following NSPECCMP species will be computed by scaling and summing on or more of the processed input species. Identify the name(s) of the computed species and provide the scaling factors for each of the NSPECINP input species.”

! CSPECCMP = SOA ! ! PM081 = 1.0 ! ! PM056 = 1.0 ! ! END ! ******

4.4.2. Nitrate Partitioning . The Ammonia Limiting Method may be used in POSTUTIL to repartition nitrate using the default ba ckground ammonia concentration of 0.5 ppb.

4.4.3. Background Light Extinction Method . A monthly average natural background shall be calculated using Visibility Method 6, based on the annual average defaultChemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol

P:\Chemical Products\BART\2008 BART\VISTAS\Chemical Products VISTAS BART Protocol.DRAFT.GAEPD.Rev2.docPage 11 

# DRAFT

natural conditions specified in Table 3-1 and the monthly average 

### f (RH) values specified in Table 3-2 for the centroid of each Class I Area.

4.4.4. Light Extinction Efficiencies. The US EPA (2003a) default IMROVE algorithm light extinction efficiencies shall be used.



### 5. BART Exemption Modeling Results Review

To be determined after GA EPD acceptance of modeling protocol.

Results shall be provided in Table 5-1 using the existing IMPROVE algorithm. Results shall be provided in Appendix E using both the existing and new IMPROVE algorithm. The evaluation of impacts using the new IMPROVE formula shall be performed using the guidance and calculation spreadsheet provided by VISTAS technical cons ultant, Dr. Ivar Tombach. The new IMPROVE algorithm takes into account site-specific Rayle igh scattering and a Class I-specific sea salt concentration. Where site-specific data were not available, the default EPA Rayleigh scattering value was used with a sea salt concentration of zero.

Table 5-1 Summary of CALPOST IMPROVE Change in Light Extinction Results Number of Days > 0.5 dv in Class I Area Class I Area Distance (km) 2001 2002 2003 Highest 8th-High

(Max) Highest 22nd-High

3 Years

Cohutta Wilderness Area 78.9 TBD TBD TBD TBD TBD

Great Smokey Mountains National Park 167.0 TBD TBD TBD TBD TBD

Joyce-Kilmer Wilderness Area 153.4 TBD TBD TBD TBD TBD

Sipsey Wilderness Area 237.3 TBD TBD TBD TBD TBD

Shining Rock Wilderness Area 217.4 TBD TBD TBD TBD TBDChemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol 

# DRAFT

Appendix A
Chemical Products Corporation BART-Eligible Source Potential-to-Emit Calculations, Maximum 24-hr ActualEmissions Rates, and Modeled Emissions InventoryChemical Products Corporation - Best Available Retrofit Technology Eligible Emissions Units (BART-Eligible Source)Table A-1. Chemical Products Corporation BART-Eligible Source Potential-to-EmitCokeStack IDNo.aEmissionsUnit ID No.aEmission Unit DescriptionaConstructionDateaMaximumRated HeatInput,(mmBtu/hr)aFuel OilConsumption,(gph)aCokeConsumption,(tph)aDesignCapacity,(tph)aLCCEasting,(km)bLCCNorthing,(km)bStackHeight,(ft)cExhaustTemp., t(oF)cExitVelocity, (fps)cStackDiameter,(ft)cGeorgia Rule (d)PM10 EmissionsPermit Limit, (lb/hr)Georgia Rule (e)PM10 EmissionsPermit Limit, (lb/hr)Georgia Rule (g)Fuel SulfurPermit Limit, (%)Coke SulfurContent, (%)dSO2 EmissionsFactor, MassBalance,(lb/ton)eNOx EmissionsFactor, AP-42 1.4-1/-2,(lb/mmBtu)SO2 EmissionsFactor, AP-42 1.4-1/-2,(lb/mmBtu)NOx EmissionsFactor, AP-42 1.5-1,(lb/mmBtu)SO2 EmissionsFactor, AP-42 1.5-1,(lb/mmBtu)Fuel OilTypeNOx EmissionsFactor, AP-42 1.3-1,(lb/103gal)SO2 EmissionsFactor, AP-42 1.3-1,(lb/103gal)PM10, PotentialEmissions Rate,(lb/hr)fNOx, PotentialEmissions Rate,(lb/hr)gSO2, PotentialEmissions Rate,(lb/hr)gChemical Products Corporation, Title V Permit No. 2819-015-0008-V-02-0/-1BCZ1 BC01 Carbonate Product South Rotary Calciner 1969 24.00 179 -- 5.50 1118.276 -572.125 45 400 50.0 2.26 -- 12.85 2.5 -- -- 0.098 0.00059 0.153 0.000218 No. 2 20.00 355.00 12.848 3.672 63.545BDZ1 BD01 Carbonate Product Rotary Dryer 1964 8.89 66 -- 4.00 1118.296 -572.124 48 300 13.0 2.00 -- 10.38 2.5 -- -- 0.098 0.00059 0.153 0.000218 No. 220.00 355.00 10.379 1.360 23.430BDZ3 BD03 Carbonate Product South Spray Dryer 1968 5.85 44 -- 3.50 1118.279 -572.141 48 300 35.0 1.35 -- 9.49 2.5 -- -- 0.098 0.00059 0.153 0.000218 No. 2 20.00 355.00 9.491 0.895 15.620BDZ4 BD04 Barium Chloride Rotary Dryer 1976 4.50 34 -- 3.00 1118.269 -572.136 30 250 12.0 1.67 -- 8.56 2.5 -- -- 0.098 0.00059 0.153 0.000218 No. 2 20.00 355.00 8.560 0.689 12.070BKZ1 BK02 Ore Reduction East Rotary Kiln 1966 10.20 76 1.00 5.00 1118.356 -572.148 195 400 54.0 4.00 -- 12.05 2.5 7.50 0.1500 0.098 0.00059 0.153 0.000218 No. 6 55.00 392.50 12.053 4.180 329.830GBZ1 GBX1 B&W Process Steam Boiler 1973 100.00 746 -- -- 1118.376 -572.162 55 500 38.0 2.99 15.81 -- 3.0 -- -- 0.098 0.00059 0.208 0.000262 No. 2 20.00426.00 15.811 20.765 317.796Notes:PTE (tpy) PM NOx SO2a. Stack and emissions unit identification, construction date, maximum rated heat input, fuel oil consumption rates, petroleumcoke consumption and design capacity determined from Title V Permit No. 2819-015-0008-V-02-0/-1 and Title V Application No. 16069.South Rotary Calciner 56.27 16.08 278.33b. Lambert Conformal Conic (LCC) coordinates determined from conversion from latitude and longitude coordinates determined from Google Earth and stack locations provided by Chemical Products Corporation.Rotary Dryer 45.46 5.96 102.62c. Stack characteristics determined from CALPUFF.INP files provided by GA EPD from 2006 12km BART exemption modeling performedby VISTAS.South Spray Dryer 41.57 3.92 68.42d. Highest petroleum coke sulfur content based on historical data for calendar years 2001, 2002, and 2003 provided by ChemicalProducts Corporation.Barium Chloride Rotary Dryer 37.49 3.02 52.87e. SO2 emissions factor for processing of petroleum coke determined from twice the highest sulfur content in coke processed based on historical data and 100% conversion to SO2; (2)*(2000 lb/hr coke processed)*(6.48 lb sulfur per 100 lb coke processed)/(2000 lb/ton).East Rotary Kiln 52.79 18.31 1444.66f. PM10 potential-to-emit determined from Georgia Rules (d) and (e).B&W Boiler 69.25 90.95 1391.95g. NOx and SO2 potential-to-emit determined from maximum emissions rates resulting from combustion of natural gas, propane, orfuel oil. SO2 emissions rate for the East Rotary Kiln includes estimated SO2 from processing coke.BART-Eligible Source Total302.84 138.24 3338.83Table A-2. Chemical Products Corporation BART-Eligible Source Maximum 24-hr Actual Emissions Rates (Fuel Oil Combustion, 0.5%Sulfur Content, AP-42)CokeStack IDNo.aEmissionsUnit ID No.aEmission Unit DescriptionaConstructionDateaMaximumRated HeatInput,(mmBtu/hr)aFuel OilConsumption,(gph)aCokeConsumption,(tph)aDesignCapacity,(tph)aLCCEasting,(km)bLCCNorthing,(km)bStackHeight,(ft)cExhaustTemp., t(oF)cExitVelocity, (fps)cStackDiameter,(ft)cGeorgia Rule (d)PM10 EmissionsPermit Limit, (lb/hr)Georgia Rule (e)PM10 EmissionsPermit Limit, (lb/hr)Maximum ActualFuel SulfurPermit Limit, (%)MaximumActual CokeSulfur Content,(%)dSO2 EmissionsFactor, MassBalance,(lb/ton)eNOx EmissionsFactor, AP-42 1.4-1/-2,(lb/mmBtu)SO2 EmissionsFactor, AP-42 1.4-1/-2,(lb/mmBtu)NOx EmissionsFactor, AP-42 1.5-1,(lb/mmBtu)SO2 EmissionsFactor, AP-42 1.5-1,(lb/mmBtu)Fuel OilTypeNOx EmissionsFactor, AP-42 1.3-1,(lb/103gal)SO2 EmissionsFactor, AP-42 1.3-1,(lb/103gal)PM10, PotentialEmissions Rate,(lb/hr)fNOx, PotentialEmissions Rate,(lb/hr)gSO2, PotentialEmissions Rate,(lb/hr)gChemical Products Corporation, Title V Permit No. 2819-015-0008-V-02-0/-1BCZ1 BC01 Carbonate Product South Rotary Calciner 1969 24.00 179 -- 5.50 1118.276 -572.125 45 400 50.0 2.26 -- 12.85 0.50 -- --0.000 0.000 0.000 0.000No. 2 20.00 71.00 12.848 3.580 12.709BDZ1 BD01 Carbonate Product Rotary Dryer 1964 8.89 66 -- 4.00 1118.296 -572.124 48 300 13.0 2.00 -- 10.38 0.50 -- --0.000 0.000 0.000 0.000No. 2 20.00 71.00 10.379 1.320 4.686BDZ3 BD03 Carbonate Product South Spray Dryer 1968 5.85 44 -- 3.50 1118.279 -572.141 48 300 35.0 1.35 -- 9.49 0.50 -- --0.000 0.000 0.000 0.000No. 2 20.00 71.00 9.491 0.880 3.124BDZ4 BD04 Barium Chloride Rotary Dryer 1976 4.50 34 -- 3.00 1118.269 -572.136 30 250 12.0 1.67 -- 8.56 0.50 -- --0.000 0.000 0.000 0.000No. 2 20.00 71.00 8.560 0.680 2.414BKZ1 BK02 Ore Reduction East Rotary Kiln 1966 10.20 76 1.00 5.00 1118.356 -572.148 195 400 54.0 4.00 -- 12.05 0.50 7.50 0.15000.000 0.000 0.000 0.000No. 6 55.00 78.50 12.053 4.180 305.966GBZ1 GBX1 B&W Process Steam Boiler 1973 100.00 746 -- -- 1118.376 -572.162 55 500 38.0 2.99 15.81 -- 0.50 -- --0.000 0.000 0.000 0.000No. 2 20.00 71.00 15.811 14.920 52.966Notes:PTE (tpy) PM NOx SO2a. Stack and emissions unit identification, construction date, maximum rated heat input, fuel oil consumption rates, petroleumcoke consumption and design capacity determined from Title V Permit No. 2819-015-0008-V-02-0/-1 and Title V Application No. 16069.BART-Eligible Source Total302.84 111.95 1672.57b. Lambert Conformal Conic (LCC) coordinates determined from conversion from latitude and longitude coordinates determined from Google Earth and stack locations provided by Chemical Products Corporation.c. Stack characteristics determined from CALPUFF.INP files provided by GA EPD from 2006 12km BART exemption modeling performedby VISTAS.d. Highest expected petroleum coke sulfur content based on historical data provided by Chemical Products Corporation.e. SO2 emissions factor for processing of petroleum coke determined from twice the highest sulfur content in coke processed based on historical data and 100% conversion to SO2; (2)*(2000 lb/hr coke processed)*(6.48 lb sulfur per 100 lb coke processed)/(2000 lb/ton).f. PM10 potential-to-emit determined from Georgia Rules (d) and (e).g. NOx and SO2 potential-to-emit determined from maximum emissions rates resulting from combustion of fuel oil with 0.50% sulfur content. SO2 emissions rate for the East Rotary Kiln includes estimated SO2 from processing coke.Table A-3. Chemical Products Corporation BART-Eligible Source CALPUFF Modeling Input Parameters and PM10 SpecicationChemical Products Corporation, Title V Permit No. 2819-015-0008-V-02-0/-1PM081 PM056BCZ1 BC01 Carbonate Product South Rotary Calciner 34.1509 84.7862 16 1118.276 -572.125 45 750 400 50.0 2.26 12.709 0.000 3.580 12.848 0.000 0.0000 100 50 50 6.424 6.424BDZ1 BD01 Carbonate Product Rotary Dryer 34.1509 84.7860 16 1118.296 -572.124 48 753 300 13.0 2.00 4.686 0.000 1.320 10.379 0.000 0.000 0 100 50 50 5.190 5.190BDZ3 BD03 Carbonate Product South Spray Dryer 34.1508 84.7862 16 1118.279 -572.141 48 753 300 35.0 1.35 3.124 0.000 0.880 9.491 0.000 0.000 0 100 5050 4.745 4.745BDZ4 BD04 Barium Chloride Rotary Dryer 34.1509 84.7863 16 1118.269 -572.136 30 735 250 12.0 1.67 2.414 0.000 0.680 8.560 0.000 0.000 0 100 50 50 4.2804.280BKZ1 BK02 Ore Reduction East Rotary Kiln 34.1506 84.7854 16 1118.356 -572.148 195 900 400 54.0 4.00 305.966 0.000 4.180 12.053 0.000 0.000 0 100 50 50 6.026 6.026GBZ1 GBX1 B&W Process Steam Boiler 34.1505 84.7852 16 1118.376 -572.162 55 760 500 38.0 2.99 52.966 0.000 14.920 15.811 0.000 0.000 0 100 50 50 7.906 7.906Stack IDNo.EmissionsUnit ID No.Emission Unit DescriptionOrganic Condensable (OC)Assumed Fraction of PM100.625 - 1.0μm, (%)0.5 - 0.625 mm, (%)BaseElevation, (ft)H2SO4, (lb/hr) NOx, (lb/hr) PM10, (lb/hr)Organic Condensable (OC)0.625 - 1.0μm, (lb/hr)0.5 - 0.625 mm, (lb/hr)PM2.5, (lb/hr) NH3, (lb/hr)PM10 Particle Speciation (Unknown)Filterable, (%) Condesable, (%)ExhaustTemp., t(oF)ExitVelocity,(fps)StackDiameter, (ft)SO2, (lb/hr)Fuel OilFuel OilLocationLatitude Longitude ZoneLocation LCCEasting, (km) Northing, (km)StackHeight, (ft)Natural Gas PropaneNatural Gas PropaneP:\Chemical Products\BART\2008 BART\VISTAS\Chemical Products VISTAS BART.AppA.Rev1.xlsAppendix A ChemProdBARTSource Page 1 of 1Chemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol 

# DRAFT

Appendix B Example of Proposed CALPUFF Input File for Cohutta Wilderness Area (2001)AppendixBCohutta2001~~LPuFF.1~~ VISTAS BART - chem Products Corp. 

## - 4km Domain NO. 

### 4 

## - cohutta 

## - 2001 Computational Grid X (Easterly): From 76 To 127 Computational Grid Y (~ortherly): From 

### 136 To 208


### ---------------- Run 

### title 

### (3 lines)



### ..........................................

INPUT GROUP: 0 

## -- Input 

### and output File Names



## -------------- oefault Name Type File Name



### ------------ 

## ---- - - - - - - - - - CALMET.DAT input 

# " METOAT = ..\CALMET.DAT 

## "

or 1SCMET.DAT input 

# " ISCDAT = 

## *

or PLMMET.DAT 

### input 

## * PLMDAT = 

### +
or PROFILE.DAT input 

# " PRFDAT = SURFACE.DAT 

### inuut 

## * SFCOAT = RESTART6



## . OAT input 

## * RSTARTB= 

### + 

### -------------------------------________________________________________-------________________________________________---------------------------------------- CALPUFF.LST Output 

### ! PUFLST =CALPUFF.LST ! CONC



# .DAT Output ! CONDAT =CONC.DAT DFLX. OAT Output 

### ! DFDAT =DFLX.OAT I WFLX.DAT Output ! WFDAT =WFLX.OAT 

### !

VISB



# .DAT Output 

### ! VISDAT =VISB.DAT I TKZD. DAT OUtDUt 

## " TZDDAT = 

### + RHOZD. DAT output 

# " RHODAT= RESTARTE.DAT OUtDUt 

## * RSTARTE=

Emission Fi 

### 1

es


### -------------- PTEMARB.DAT input 

# " PTDAT = 

### + VOLEMARB.DAT 

### input 

# " VOLDAT = 5; BAEMARB.DAT input 

# " ARDAT = & LNEMARB.DAT input 

## * LNDAT = 1


### other Files

OZONE



## . DAT 

### input ! OZDAT = 

### F:\VISTASmet\~ART4km~4\03data\ozone2OOldom4.dat I VD



# . DAT input "VDDAT = 2% CHEM



# .DAT 

### input 

## * CHEMOAT= f H202. DAT 

### input 

# " H202DAT= % HILL. DAT 

### input 

## * HILDAT= ? HILLRCT.DAT input "RCTDAT= 4 COASTLN.DAT input 

# " CSTDAT= 57 FLUXBDY .DAT 

### input 

## * BDYDAT= % BCON



## . DAT 

### input 

# " BCNDAT= 

### + DEBUG. DAT Output 

### ! DEBUG = DEBUG.DAT 

### ! MASSFLX.DAT Output 

# " FLXDAT= 

### + MASSBAL.DAT Output 

# " BALDAT= 

### + FOG. DAT Output 

# " FOGDAT= 

## *

~ll file names 

### will be converted to 

### lower case 

### if LCFILES = T 

### otherwise, 

### if LCFILES = F, 

### file names 

### will be converted to UPPER CASE T = lower 

### case ! LCFILES = F 

### ! F = UPPER CASE Page 

### 1Appendix~Cohutta2OOlCALPUFF.INP NOTE: 

### (1) file/path names can be up to 70 characters in length

provision for multiple input files

Number of CALMET.DAT files for run (NMETDAT) Default: 

### 1

Number of PTEMARB.DAT files for run (NPTDAT) Default: 0

Number of BAEMARB. DAT files for run (NARDAT) Default: 0

Number of VOLEMARB.DAT files for run (NVOLDAT) Default: 0

!END!



### ------------- subgroup (Oa) ! NMETDAT = 

### 36 !

! NPTDAT = 0 

### !

! NARDAT 

### = 0 

### !

! NVOLDAT = 0 

### !

The following CALMET.DAT filenames are processed in sequence 

### if NMETDAT>~

none input ! METOAT= 

### ~:\~1STASmet\B~~~4kmD4\2001-4~~\met2OOl-dom4-Ola.dat !

none input 

### ! METDAT= 

### F:\VISTASmet\BART4kmD4\2001-4KM\met2001-dom4-Olb.dat !

none input 

### ! METDAT= 

### ~:\~1STASmet\~~~~4kmD4\2001-4~~\met2OOl-dom4-Olc.dat !

none input 

### ! METDAT= 

### ~:\~1ST~Smet\~~~~4kmD4\2001-4~~\met2001-dom4-02a.dat !

none input ! METDAT= 

### F:\V1STASmet\BART4kmD4\2001-4~~\met2001-dom4-02b.dat !

none input ! METDAT= 

### F:\~1STASmet\B~~~4kmD4\2001-4~~\met2001-dom4-02c.dat !

none input 

### ! METDAT= ~:\~1STASmet\~~~~4kmD4\2001-4~~\met2001-dom4-03a.dat 

### !

none input ! METDAT= 

### F:\V1STASmet\BART4kmD4\2001-4K~\met2001-dom4-03b.dat !

none input ! METDAT= 

### F:\~1STASmet\B~~~4kmD4\2001-4~~\met2001-dom4-03c.dat !

none input 

### ! METDAT= ~:\~15TASmet\~~~~4kmD4\2001-4~~\met2001-dom4-04a.dat 

### !

none input '! METDAT~ 

### ~:\~1STASmet\B~~~4kmD4\2001-4~~\met2001-dom4-04b.dat !

none input 

### ! METDAT= 

### F:\VISTASmet\BART4kmD4\2001-4KM\met2OOl-dom4-04c.dat !

none input .! METDAT~ 

### ~:\~1STASmet\B~~~4kmD4\2001-4~~\met2OOl-dom4-05a.dat !

none input 

### ! METDAT= ~:\~1STASmet\B~~~4kmD4\2001-4~~\met2OOl-dom4-05b.dat !

none input 

### ! METDAT= 

### F:\V1STASmet\BART4kmD4\2001-4KM\met2OOl-dom4-05c.dat !

none input 

### ! METDAT= 

### ~:\~1STASmet\~~~~4kmD4\2001-4~~\met2OOl-dom4-06a.dat !

none input ! METDAT= 

### ~:\V1STASmet\~~~~4kmD4\2001-4~~\met2OOl-dom4-06b.dat !

none input 

### ! METDAT= Page 

### 2 !END!

!END!

!END!

!END!

!END!

!END!

!END 

### !

!END!

!END!

! END !

!END!

!END!

!END!

!END 

### !

!END!

!END!

!END!

# . INP

~~~ - -- - - - . -

... . . . -

. -- - 

### . none 'input 

### '! METDAT= 

### F:\VIST~Smet\~ART4kmD4\2001-4KM\met2OOl-dom4-07a.dat !

none input ! METDAT= 

### F:\VISTASmet\BART4kmD4\2001-4KM\met2001-dom4-07b.dat !

-~ ~~~-~-

~ none input 

### '! METDAT~ 

### F:\VIsTAsmet\~ART4kmD4\2001-4KM\met2OOl-dom4-07c.dat !

none input ! METDAT= 

### F:\V1S~~Smet\BART4kmD4\2001-4KM\met2OOl-dom4-08a.dat !

none input 

### ! METDAT= 

### ~:\~1~~~~met\B~~~4kmD4\2001-4~~\met2OOl-dom4-08b.dat !

none input ! METDAT= 

### F:\~1S~~Smet\BART4kmD4\2001-4K~\met2OOl-dom4-08c.dat !

none input ! METDAT= 

### ~:\~15~~~met\B~~~4kmD4\2001-4K~\met2OOl-dom4-09a.dat !

none input 

### ! METDAT= 

### ~:\~1s~~~met\~~~~4kmD4\2001-4~~\met2001-dom4-09b.dat !

none input 

### ! METDAT= 

### ~:\~1S~~Smet\BART4kmD4\2001-4K~\met2001-dom4-09c.dat !

none input ! METDAT= 

### ~:\~1~~~~met\B~~~4kmD4\2001-4~~\met2OOl-dom4-lOa.dat !

none input !METDAT= 

### ~:\~1~~~~met\~~~~4kmD4\2001-4~~\met2OOl-dom4-lOb.dat !

none input 

### ! METDAT= 

### F:\VIS~~Smet\BART4kmD4\2001-4KM\met2OOl-dom4-lOc.dat !

none input !METDAT= 

### F:\V1S~~Smet\BART4kmD4\2001-4KM\met2OOl-dom4-lla.dat !

none input ! METDAT= 

### F:\V1S~~Smet\BART4kmD4\2001-4KM\met2OOl-dom4-llb.dat !

none input ! METDAT= 

### ~:\~1S~~~met\B~~~4kmD4\2001-4~~\met2OOl-dom4-llc.dat !

none input ! METDAT= 

### ~:\~1S~~~met\B~~~4kmD4\2001-4~~\met2OOl-dom4-12a.dat !

none input ! METDAT= 

### ~:\~1S~~Smet\B~~~4kmD4\2001-4~~\met2OOl-dom4-12b.dat !

none input 

### ! METDAT= 

### ~:\~1~~~~met\~~~~4kmD4\2001-4~~\met2OOl-dom4-12c.dat ! END !

END!

END !

END !

INPUT GROUP: 

### 1 

# -- General run control parameters



### --------------

Option to run all periods found in the met. file (METRUN) oefault: 0 

### ! METRUN = 0 

### !

METRUN = 0 

## - Run period explicitly defined below METRUN = 

### 1 

## - Run all periods in met. file
Starting date: Year (IBYR) 

## -- No default ! IBYR = 2001 

### ! (used only 

### if Month (IBMD) 

### -- No default ! IBMO = 

### 1 !

METRUN = 0) Day (IBDY) 

## -- NO default 

### ! IBDY = 

### 1 ! Hour (IBHR) 

## -- No default 

### ! IBHR = 

### 1 !

Note: IBHR is the time at the END of the first hour of the simulation (IBHR=~, the first hour of a day, runs from 00:OO to 01:OO)

Base time zone (XBTZ) 

### -- NO default ! XBTZ = 

### 5 ! The zone is the number of hours that must be ADDED to the time to obtain UTC (or GMT) Examples: PST = 8., MST = 7. CST = 

### 6., EST = 

### 5. Page 3

### Appendix~Cohutta2001~~L~u~~.1~~

Length of run (hours) (IRLG) 

## -- No default 

### ! IRLG = 8760 

### !



### umber of chemical species (NSPEC) ~efault: 

### 5 ! NSPEC = 

### 11 !

Number of chemical species to be emitted (NSE) Default: 3 ! NSE = 

### 4 !

Flag to stop run after SETUP phase (ITEST) Default: 

### 2 

### ! ITEST = 

### 2 ! (used to allow checking of the model inputs, files, etc.) ITEST = 

### 1 

## - STOPS program after SETUP phase ITEST = 

### 2 

## - Continues with execution of program after SETUP

Restart configuration :

Control flag (MRESTART) ~efault: 0 

### ! MRESTART = 0 

### !

0 = Do not read or write a restart file 

### 1 = Read a restart file at the beginning of the run 

### 2 = write a restart file during run 3 = Read a restart file at beginning of run and write a restart file during run

Number of periods in Restart output cycle (NRESPD) Default: 0 

### ! NRESPD 

### = 0 

### !

0 = File written only at last period >O = File updated every NRESPD periods
Meteorological Data Format (METFM) ~efault: 

### 1 ! METFM = 

### 1 !

METFM = 

### 1 

## - CALMET binary file (CALMET.MET) METFM = 

### 2 

## - ISC ASCII file (1SCMET.MET) MLlFM = 

### 3 

## - AUSPLUME ASCII file (PLMMET.MET) METFM = 4 

## - CTDM O~US tower file fPROFILE.DAT) and surfake parameters fije (SURFACE.DAT) METFM = 5 

## - AERMET tower file (PROFILE.DAT) and surface parameters 'file (SURFA~E.DAT)

Meteorologi ca1 Profile Data Format (MPRFFM) (used only for METFM = 

### 1, 2, 

### 3) ~efault: 

### 1 ! MPRFFM 

### = 

### 1 !

MPRFFM = 

### 1 

## - CTDM plus tower file (PROFILE.DAT) MPRFFM = 

### 2 

## - AERMET tower file (PROFILE.DAT)

PG sigma-y is adjusted by the factor (AVET/PGTIME)*"0.2 ~veraging Time (minutes) (AVET) ~efault: 60.0 

### ! AVET = 60. ! PG Averaging Time (minutes) (PGTIME) Default: 60.0 ! PGTIME = 60. 

### !

Page 

### 4INPUT GROUP: 

### 2 

## -- ~echnical options



## --------------

vertical distribution used in the
near field (MGAUSS) 

### ~efault: 1 0 = 

### uniform 

### 1 = 

### Gaussian

Terrain adjustment method (MCTADJ) 

### ~efault: 3 0 = 

### no adjustment 1= 

### ISC-tvoe of terrain adjustment 2 = 

### simple', CALPUFF-type 

## OF terrain adjustment 3 = 

### partial plume path adjustment
Subgrid-scale complex terrain flag (MCTSG) 

### Default: 0 0 = 

### not modeled 

### 1 = 

### modeled

Near-field puffs modeled as elongated slugs? (MSLUG) 0 = no 

### 1 = 

### yes (slug model used) Default: 0



### Transitional plume rise modeled ? (MTRANS) 

### oefault: 1 0 

### =no (i.e., final rise only) 1 = 

### yes 

### (i .e., transitional rise computed)

stack tip downwash? (MTIP) 

### Default: 1 0 = 

### no 

### (i .e., no stack tip downwash) 1 = 

### yes (i.e., use stack tip downwash)

Method used to simulate building downwash? (MBDW) 

### Default: 1 1 = 1% 

### method 2 = PRIME 

### method

Vertical wind shear modeled above stack top? (MSHEAR) 

### Default: 0 0 = 

### no 

### (i .e. 

## , vertical wind shear not modeled) 1 = 

### yes 

### (i .e., vertical wind shear modeled)
puff splitting a1lowed? (MSPLIT) 

### Default: 0 0 = 

### no 

### (i .e., puffs not split) 1 = 

### yes 

### (i .e., puffs are split)

chemical mechanism flag (MCHEM) 

### Default: 1 0 = 

### chemical transformation not modeled 

### 1= 

### transformation rates computed internally (MESOPUFF 

### II scheme) 2 = 

### user-specified transformation rates used 3 = 

### transformation rates computed internally (RIVAD/ARM~ 

### scheme) 4 = 

### secondary organic aerosol formation computed (MESOPUFF 11 

### scheme for OH) ! MGAUSS = 

### 1

! MCTADJ = 

### 3

! 

### MCTSG = 0



### ! MSLUG = 0

! 

### MTRANS = 

### 1

! 

### MTIP = 

### 1 !

! MBDW = 

### 1 !

! MSHEAR = 0 !



### ! MSPLIT = 0 

### !

! MCHEM = 

### 1 !

Page 5Aqueous hase transformation flag (MAQCHEM) (used 0n?y 

### if MCHEM = 

### 1, or 3) Default: 0 ! MAQCHEM = 0 ! 0 = aqueous phase transformation not modeled 

### 1 = transformation rates adjusted for aqueous phase reactions

wet removal modeled ? (MWET) ~efault: 

### 1 ! MWET= 

### 1 ! 0 = no 

### 1 = yes Dry deposition modeled ? (MDRY) oefault: 

### 1 ! MORY = 

### 1 ! 0 = no 

### 1 = yes (dry deposition method specified

for each species in Input Group 3)

Gravitational settling (plume tilt) modeled ? (MTILT) Default: 0 

### ! MTILT= 0 

### ! 0 = no 

### 1 = yes (puff center falls at the gravitational

settling velocity for 

### 1 particle species)

Restrictions:

- MDRY = 

### 1



## - NSPEC = 

### 1 (must be particle species as well)

- Sg = 0 GEOMETRIC STANDARD DEVIATION in Group 8 is set to zero for a single particle diameter

Method used to compute dispersion coefficients (MDISP) Default: 3 ! MDISP = 3 !



### 1 = disoersion coefficients comouted from measured values of turbulence, sigma v, sigma w 2 = dispersion coefficients from internally calculated sigma v, sigma w using mi crometeorological variables (u*, w* 

## , L, etc.) 3 = PG dispersion coefficients for RURAL areas (computed using the ISCST multi-seqment approximation) and MP coefficients in

-

. . urban areas 4 = same as 3 except PG coefficients computed using the MESOPUFF I1 eqns. 

### 5 = CTDM sigmas used for stable and neutral conditions. For unstable conditions, sigmas are computed as in MDISP = 3, described above. MDISP = 

### 5 assumes that measured values are read
Sigma-v/sigma-theta, sigma-w measurements used? (MTURBVW) (used only 

### if MDISP = 

### 1 or 5) ~efault: 3 ! MTURBVW = 

### 3 ! 1 = use sigma-v or sigma-theta measurements from PROFILE.DAT to Compute Sigma-y (valid for METFM = 

### 1, 2, 

### 3, 4, 

### 5) 2 = use s~gma-w measurements from PR~FILE.DAT to compute sigma-z (valid for METFM = 

### 1, 2, 3, 4, 

### 5) 3 = use both sigma-(v/theta) and sigma-w from PROFILE.DAT to compute sigma-y and sigma-z (valid for METFM = 

### 1, 2, 3, 

### 4, 5) 4 = use sigma-theta measurements from PLMMET.DAT to compute sigma-y (valid only 

### if METFM = 3) Page 6

### Back-up method used to compute dispersion when measured turbulence data are missing (MDISPZ) 

### Default: 3 ! MDISP2 = 

### 3 ! 

### (used only 

### if MDISP = 

### 1 or 5) 2 = 

### dispersion coefficients from internal1 y calculated sigma v, sigma w 

### using micrometeorological variables (Uf W* 

## , L, 

### etc.) 3 = PG 

### hispersion coefficients for RURAL areas (computed using the ISCsT multi-segment approximation) and MP 

### coefficients in urban areas 4 = 

### same as 3 except PG 

### coefficients computed using the MESOPUFF 11 

### eqns. [DIAGNOSTIC FEATURE] Method used for Lagrangian timescale for Sigma-y (used only 

### if MDISP=1,2 or MDISP2=1,2) (MTAULY) 

### Default: 0 ! 

### MTAULY = 0 ! 0 = 

### Draxler default 617.284 (s) 1= 

### Computed as Lag. Length 

### / (.75 q) 

## -- after SCIPUFF 10 < Direct user input (5) 

## -- e.g., 306.9

Method used for ~dv&tive-~eca~ timescale for Turbul (used only if 

### MDISP=2 or MDISP2=2) (MTAUADV) 

### Default: 0 0 = 

### No turbulence advection 1 = 

### computed (OPTION NOT IMPLEMENTED) 10 < Direct user input (s) 

## -- e.g., 300 ence

! MTAUADV

Method used to compute turbulence sigma-v & 

### sigma-w using mi 

### crometeorologi cal variables (used only if MDISP = 

### 2 or MDISP2 = 

### 2) (MCTURB) 

### Default: 1 ! MCTURB = 

### 1 ! 1 = 

### standard CALPUFF subroutines 

### 2 = AERMOD 

### subroutines

PG 

### sigma-y,z adj. for roughness? Default: 0 ! MROUGH = 0 ! (MROUGH) 

### 0 = 

### no 1= 

### yes

Partial plume penetration of ~efault: 1 ! MPARTL = 1 

### ! elevated inversion? (MPARTL) 

### 0 = 

### no 1 = 

### yes

Wren th of temperature inversion Default: 0 ! MTINV = 0 

### ! proviied in PROFILE.DAT extended records? (MTINV) 0 = 

### no (computed from measured/default gradients) 

### 1 = 

### yes

PDF used for dispersion under convective conditions? ~efault: 0 ! MPDF = 0 

### ! (MPDF) 0 = 

### no 1 = 

### yes
Sub-Grid TIBL module used for shore line? Page 7

### AppendixBCohutta2OOlCALPUFF.1~~ Oefault: 0 ! 

### MSGTIBL = 0 !

0 = 

### no 

### 1 = 

### yes Boundary conditions (concentration) modeled?

(MBCON) 

### Default: 0 

### ! MBCON = 0 !



### 1= 

### yes, using formatted BCON.DAT 

### file 2 = 

### yes, using unformatted CONC.OAT 

### file

Note: MBCON z 0 

### reauires that the last soecies modeled be 'BCON'



# . dass is placed in spec;es BCON 

### when generatinq boundarv condition ouffs so that clean air entering the modeling domain can be simulated in the same way as polluted air. specify zero emission of species BCON 

### for all regular sources.

Individual source contributions saved? Default: 0 ! MSOURCE = 0 ! (MSOURCE) 0 = 

### no 

### 1 = 

### yes

Analyses of fogging and icing impacts due to emissions from arrays of mechanically-forced cooling towers can be performed using CALPUFF in conjunction with a cooling tower emissions processor (CTEMISS) and its associated postprocessors. Hourly emissions of water vapor and temperature from each cooling tower cell are computed for the current cell configuration and ambient conditions by CTEMISS. CALPUFF models the disoersion of these emissions and provides cloud information in a' specialized format for further analysis. Output to FOG.DAT is orovided in either 'plume mode' or 'receptor.mode' format.
Configure for FOG Model output? Default: 0 

### ! MFOG = 0 (MFOG) 0 = 

### no 1 = 

### yes - 

### report results in PLUME 

### Mode format 2 = 

### yes - 

### report results in RECEPTOR 

### Mode format

Test options specified to see 

### if they conform to regulatory values? (MREG) 

### oefault: 1 ! MREG = 0

0 = NO 

### checks are made 1 = 

### Technical options must conform to USEPA 

### Long Range Transport (LRT) guidance METFM 

### 1 or 2 AVET 60. 

### (min) PGTIME 

### 60. (min) MGAUSS 1 

### MCTADJ 

### 3 MTRANS 

### 1 MTIP 

### 1 MCHEM 

### 1 or 3 

### (if modeling Sox, NOx) MWET 

### 1 MORY 

### 1 MDISP 

### 2 or 3 

### Page 8

### AppendixBCohutta2OOl~~~~u~~.1~~ MPDF 0 

### if MDISP=3 

### 1 

### if MDISP=2 MROUGH 0 MPARTL 

### 1 SYTDEP 550. (m) MHFTSZ 0 SVMIN 0.5 (m/s)

INPUT GROUP: 3a, 3b 

### -- species 

### list



## -------------------

The following species are modeled:



### ! CSPEC = SO2 

### ! !END! ! CSPEC = SO4 ! !END! 

### ! CSPEC = NOX 

### ! !END! ! CSPEC = HN03 ! !END! 

### ! CSPEC = NO3 

### ! !END! ! CSPEC = PM800 ! !END! 

### ! CSPEC = PM425 

### ! !END! ! CSPEC = PM187 

### ! !END! 

### ! CSPEC = PM112 ! !END! ! CSPEC = PM081 ! !END! 

### ! CSPEC = PMOS6 ! !END!

D~Y OUTPUT GROUP SPECIES MODELED EMInED DEPOSITED NUMBER

NAME (O=NO, YES) (O=NO, YES) (O=NO, (O=NONE, (Limit: 12 ~=COMPUTED-GAS l=lst CGRUP 

### , characters 2=COMPUTED-PARTICLE 2=2nd CGRUP, 

### in length) 3=USER-SPECIFIED) 3= etc.)

NOX = HN03 = ND3 = PM800 = PM425 = PM187 = PMlI.2 = PM081 = PMOS6 =

Note: 

### he last species 

### in (3a) 

### must 

### be 'BCON' 

### when using 

### the boundary condition option (MBCON 

### > 0). Species BCON 

### should typically 

### be 

### modeled as inert 

### (no chem transformation or removal). Page 9

## ------------- subgroup (3b)



## ------------- The following names are used for Species-~roups in which results for certain s ecies are combined (added) prior to output. The 

# 7 CGRUP name wi 1 be used as the species name in output files.
use this feature to model specific particle-size distributions by treating each size-range as a separate species. order must be consistent with 3(a) above.

INPUT GROUP: 4 

## -- Map Projection and Grid control parameters



### --------------
Projection for all 

### (X,Y):



## ......................... Map projection (PMAP) ~efault: UTM ! PMAP = LCC !

UTM : Universal Transverse Mercator TM : Tan ential Transverse Mercator LCC : Lam

# 8 ert Conformal Conic

PS : Polar Stereoaraohic EM : Equatorial ~ercator LAZA : Lambert Azimuthal Equal Area

False Easting and Northing (km) at the projection origin (used only 

### if PMAP= TTM, LCC, or LAZA) (FEAST) ~efaul t=0.0 

### ! FEAST = 0 

### ! (FNORTH) ~efaul t=O .O ! FNORTH = 0 

### !

UTM zone (1 to 60) (used only 

### if PMAP=UTM) (IUTMZN) NO Default ! IUTMZN = 17 !
Hemisphere for UTM projection? (used only 

### if PMAP=UTM) (UTMHEM) ~efault: N ! UTMHEM = N 

### ! N : Northern hemisphere projection 5 : southern hemisphere projection

Latitude and Longitude (decimal degrees) of projection origin (used only 

### if PMAP= TTM, LCC, PS, EM, or LAZA)

(RLATO) NO Default 

### ! RLATO = 40N ! (RLONO) NO Default 

### ! RLONO = 97w 

### !

TM : RLONO identifies central (true N/S) meridian of projection RLATO selected for convenience LCC : RLONO identifies central (true N/S) meridian of projection RLATO selected for convenience PS : RLONO identifies central (grid N/S) meridian of projection RLATO selected for convenience EM : RLONO identifies central meridian of projection RLATO is REPLACED by O.ON 

### quato tor) LAZA: 

### RLONO identifies longitude of tangent-point of mapping plane 

### RLATO identifies latitude of tangent-point of mapping plane Page 10

### Matching parallel (s) of latitude (decimal degrees) for projection (used only if PMAP= LCC 

### or PS)

(XLAT~) NO 

### Default ! XLAT~ = 

### 33N ! 

### (XLAT2) NO 

### Default ! XLATZ = 

### 45N !



### LCC : 

### Projection cone slices through Earth's surface at RLAT~ 

### and RLAT~ PS : 

### Projection plane slices through Earth at RLATl (RLAT2 is not used)

Note: ~atitudes and longitudes should be positive, and include a letter N,S,E, 

### or W 

### indicating north or south latitude, and east or west longitude. For example, 35.9 N 

### Latitude = 

### 35.9N 118.7 E 

### Longitude = 

### 118.7E

Datum- regi on

The Datum-Region for the coordinates is identified by a character string. Many ma ping roducts currently available use the model of the Earth known as tRe woryd Geodetic System 1984 (WGS-84). other local models may be in use, and their selection in CALMET will 

### make its output consistent with local mapping products. The list of Datum-Regions with official transformation parameters is provided by the National Imagery and Mapping Agency (NIMA).

NIMA 

### Datum 

## - ~egions(Examp1es)

WGS-84 WGS-84 Reference ~lliwsoid and ~eoid. Global coveraae IWGS841 NAS-C 

### NORTH AMERICAN 

### 1927 c'larke 1866 spheroid, MEAN 

### FOR 

## ON^ (~~b27) NAR-C 

### NORTH AMERICAN 

### 1983 GRS 

### 80 Spheroid, MEAN FOR CONUS 

### (NAD83) Nws-84 NWS 

### 6370~~ Radius, sphere ESR-5 ESRI 

### REFERENCE 6371KM Radius, Sphere
Datum-region for output coordinates (DATUM) 

### Default: WGS-84 ! 

### DATUM = 

### NWS-84 !

METEOROLOGICAL Grid:
Rectangular grid defined for projection PMAP, 

### with x the Easting and Y 

### the Northing coordinate

NO. 

### x grid cells (NX) 

### No default ! NX = 

### 248 !
NO. Y 

### grid cells (NY) NO 

### default ! NY 

### = 257 !
NO. 

### vertical layers (NZ) NO 

### default ! NZ 

### = 10 !



### Grid spacing (DGRIDKM) NO 

### default ! DGRIDKM = 

### 4 ! units: km

cell face heights (zFAcE(~Z+~)) NO 

### defaults units: m 

### ! ZFACE = O., 

### 20, 40, 80, 160, 320, 640, 1200, 2000, 3000, 4000 !

~eference coordinates of SOUTHWEST 

### corner of grid cell(1, 

### 1):

x coordinate (XORIGKM) NO 

### default ! XORIGKM = 

### 718.005 ! Page 11

### ~ppendixBCohutta2OOlCALPUFF.INP Y 

### coordinate (YORIGKM) NO 

### default ! YORIGKM = 

### -1214.003 ! 

### units: km

COMPUTATIONAL Grid:

The computational grid is identical to or a subset of the MET. grid. The lower left (LL) corner of the com utational grid is at grid point (IBCOMP, IBCOMP) 

### of the MET. grid. T

# R e upper right CUR) 

### corner of the computational grid is at grid point (IECOMP, IECOMP) 

### of the MET. grid. The grid spacing of the computational grid is the same as the MET. grid.



### x index of LL 

### corner (IBCOMP) NO 

### default ! IBCOMP = 

### 76 ! (1 <= IBCOMP 

### <= NX)

Y 

### index of LL 

### corner (IBCOMP) NO 

### default ! IBCOMP = 

### 136 ! (1 <= IBCOMP C= NY)

X 

### index of UR 

### corner (IECOMP) 

### No default ! IECOMP = 

### 127 ! (1 <= IECOMP 

### <= NX)

Y 

### index of UR 

### corner (IECOMP) NO 

### default ! IECOMP = 

### 208 ! 

### (1 <= IECOMP 

### <= NY)



### SAMPLING Grid (GRIDDED 

### RECEPTORS) :



### The lower left (LL) corner of the Sam ling grid is at grid point (IBSAMP, IBSAMP) 

### of the MET. grid. Tle upper right CUR) 

### corner of the sampling grid is at grid point (IESAMP, IESAMP) 

### of the MET. grid. The sampling grid must be identical to or a subset of the computational grid. ~t 

### may be a nested grid inside the computational grid. The grid spacing of the sampling grid is DGRIDKM/MESHDN. 

### Logical flag indicating 

### if gridded receptors are used (LSAMP) 

### Default: T (T=yes 

## , F=no)



### x index of LL 

### corner (IBSAMP) NO 

### default

(IBCOMP 

### <= IBSAMP <= IECOMP)

Y 

### index of LL 

### corner (IBSAMP) NO 

### default

(IBCOMP <= IBSAMP 

### <= IECOMP)

X 

### index of UR 

### corner (IESAMP) 

### No default

(IBCOMP <= IESAMP 

### <= IECOMP)

Y 

### index of UR 

### corner (IESAMP) NO 

### default

(IBCOMP <= IESAMP <= IECOMP)



### Nesting factor of the sampling

grid (MESHDN) 

### CMESHDN is an integer >= 

### 1) ! LSAMP = F !

! IBSAMP = 

### 76

! IBSAMP = 

### 136

! IESAMP = 

### 127

! IESAMP = 

### 208

Default: 1 

### ! MESHDN = 1 

### !

Page 12INPUT GROUP: 

### 5 

### -- Output Options



### --------------
4 FILE DEFAULT VALUE



## ---- -------------

Concentrations (ICON) Dry Fluxes (IDRY) wet 

### luxes (IWET) 20 Temperature (ITZD) 2D Density (IRHO) Re1ative ~umi di ty 

### (IvIs) (relative humidity file

required for visibility analvsi s)

use data cbmpression option (LCOMPRS) output file? Oefault: T 4 VALUE THIS RUN



### --------------

! ICON = 

### 1 ! ! IDRY = 

### 1 ! ! IWET = 

### 1 ! ! IT~D = 0 

### ! ! IRHO = 0 ! 

### ! IVIS = 1 !

! LCOMPRS = T !

4 0 = Do not create file, 

### 1 = create file

QA PLOT FILE OUTPUT OPTION:

Create a standard series of output files (e.g. locations of sources, receptors, grids 

# ...)
suitable for plotting? (IQAPLOT) Default: 

### 1 ! IQAPLOT = 

### 1 ! 0 = no 

### 1 = yes

DIAGNOSTIC MASS FLUX OUTPUT OPTIONS:

Mass flux across specified boundaries for selected species reported? (IMFLX) Default: 0 ! IMFLX = 0 

### ! 0 = no 

### 1 = yes (FLUXBDY.DAT and MASSFLX.DAT filenames

are specified in Input Group 0)

Mass balance for each species reported? (IMBAL) Default: 0 

### ! IMBAL = 0 ! 0 = no 

### 1 = yes (MASSBAL.DAT filename is specified in Input Group 0)

LINE PRINTER OUTPUT OPTIONS:

Print concentrations (ICPRT) Print drv fluxes (IDPRT) print wet fluxes (IWPRT~ 

### (0 = Do not print, 

### 1 = Print)

Concentration wrint interval (ICFRQ) in timesteps Dry flux wrint interval (IDFRQ) in timesteps Wet flux print interval (IWFRQ) in timesteps ~efault: 0 

### ! ICPRT = 0 

### ! ~efault: 0 ! IOPRT = 0 

### ! Default: 0 ! IWPRT = 0 !

Default: 

### 1 ! ICFRQ = 

### 1 !

~efault: 

### 1 ! IDFRQ = 

### 1 !

~efault: 

### 1 ! IWFRQ= 

### 1 !

Page 13

### Appendix~Cohutta2OOl~~~~~~~.1~~ Units for Line Printer Output (IPRTU) Default: 

### 1 ! IPRTU= 

### 1 ! for for concentration Deposition 

### 1 = g/m**3 g/m**2/s 2 = mg/m**3 mg/m"*2/5 3 = ug/m""3 ug/mr*2/s 4 = ng/mr*3 ng/m**2/s 5 = Odour units
Messages tracking progress of run written to the screen ? (IMESG) Default: 

### 2 ! IMESG = 

### 2 ! 0 = no 

### 1= yes (advection step, puff ID) 2 = yes (YYWJJJHH, # old puffs, # emitzed puffs)

SPECIES (or GROUP for combined species) LIST FOR OUTPUT OPTIONS



## ---- CONCENTRATIONS 

## ---- FLUXES 

### ------ 

## -- MASS FLUX 

## -- SPECIES /GROUP PRINTED? SAVED ON DISK? SAVED ON DISK? SAVED ON DISK?

SO2 = 0, 0 ! SO4 = 0, 0 

### ! NOX = 0, 0 

### ! HND3 = 0, 0 

### ! NO3 = 0, 0 ! PM800 = 0, 0 ! PM425 = 0, 0 

### ! PM187 = 0, 0 ! PMll2 = 0, 0 ! PM081 = 0, 0 

### ! PMOS6 = 0, 0 ! ----*- DRY FLUXES 

### ------

PRINTED? SAVED ON DISK? 

## ------ WET

Note: Species BCON (for MBCON 

### > 0) does not need to be saved on disk.

OPTIONS FOR PRINTING DEBUG QUANTITIES (much output)
Logical for debug output (LDEBUG)

First puff to track (IPFDEB) ~efault: F 

### ! LDEBUG = F !

oefault: 

### 1 ! IPFDEB = 

### 1 !

Number of puffs to track (NPFDEB) ~efault: 

### 1 ! NPFDEB = 1000 

### ! Page 14Met. period to start output (NN~) ~efault: 

### 1 ! NN~ = 

### 1 !

Met. period to end output (NNZ) Default: 10 ! NN2 = 8760 

### !

INPUT GROUP: 6a, 6b, & 6c -- subgrid scale complex terrain inputs



### .........................



### --------------- Subgroup (6a)



### --------------- Number of terrain features (NHILL) Default: 0 ! NHILL = 0 

### !

Number of special complex terrain receptors (NCTREC) ~efault: 0 

### ! NCTREC 

### = 0 

### !

~errain and CTSG Receptor data for C~SG hills input in CTDM format ? (MHILL) NO ~efaul t ! MHILL = 2 

### ! 1 = Hill and Receptor data created by CTDM processors & read from HILL.DAT and HILLRCT.DAT files 2 = Hill data created by OPTHILL & input below in Subgroup (6b); Receptor data in subgroup (6c)

Factor to convert horizontal dimensions Default: 1.0 

### ! XHILLZM = 

### 1. ! to meters (MHILL=~)

Factor to convert vertical dimensions ~efault: 1.0 

### ! ZHILLZM = 

### 1. ! to meters (MHILL=l)
X-origin of CTDM system relative to NO ~efault ! XCTDMKM = 0.0~00 ! CALPUFF coordinate system, in Kilometers (MHILL=l)
Y-origin of CTDM system relative to NO Default 

### ! YCTDMKM = 0.0~00 ! CaLPUFF coordinate system, in Kilometers (MHILL=~)

! END 

### !

1 

# ""

HILL information

Page 

### 15

### COMPLEX TERRAIN RECEPTOR INFORMATION

XRCT Y 

### RCT ZRCT (km)



### ------ Ckm)



### ----- (m)



### ------ XHH



### 1 oescription of Complex ~errain Variables: xC, YC = 

### Coordinates of center of hill THETAH = 

### Orientation of major axis of hill (clockwise from North) ZGRID = 

### eight 

### of the 0 of the grid above mean sea 1

eve1 RELIEF = 

### Height of the crest of the hill above the grid elevation EXPO 

### 1 = 

### Hill-shape exponent for the major axis EXPO 

### 2 = 

### Hill-shape exponent for the major axis SCALE 1 = 

### Horizontal length scale along the major axis SCALE 2 = 

### Horizontal len th scale along the minor axis AMAX = 

### Maximum allowe8 axis length for the major axis BMAX = 

### Maximum allowed axis length for the major axis
XRCT, YRCT = 

### coordinates of the complex terrain receptors ZRCT = 

### Height of the ground (MSL) 

### at the complex terrain ReceDtor xHH = 

### ~ill'number associated with each complex terrain receptor (NOTE: MUST BE ENTERED AS A 

### REAL NUMBER) 

## * f NOTE: DATA 

### for each hill and CTSG 

### receptor are treated as a separate input subgroup and therefore must end with an input group terminator.

INPUT GROUP: 7 

## -- Chemical parameters for dry deposition of gases


## --------------

SPECIES 

### DIFFUSIVITY ALPHA STAR REACTIVITY MESOPHYLL RESISTANCE HENRY'S 

### LAW COEFFICIENT NAME (cIII**~/s) 

### (dimensionless) (s/cm)

502 = 

### 0.1509, 1000, 8, 0.04 ! I NOX = 

### 0.1656, 1, 8, 

### 3.5 ! I 

### HN03 = 

### 0.1628, 1, BE-08 ! 

### 18,

Page 16

### AppendixBCohutta2001~~L~U~~.1~~ INPUT GROUP: 

### 8 

## -- size parameters for dry deposition of particles



## --------------

FOT SINGLE SPECIES, 

### the mean and standard deviation are used to compute a deposition velocity for NINT (see group 9) size-ranges, and these are then averaged to obtain a mean deposition velocity.

For GROUPED 

### SPECIES, the size distribution should be explicitly specified (by the 'species' in the group), and the standard deviation for each should be entered as 0. The model will 

### then use the deposition velocity for the stated mean diameter.

SPECIES 

### GEOMETRIC MASS MEAN 

### GEOMETRIC STANDARD

NAME 

### DIAMETER DEVIATION (microns) (microns)



## - - - - - - - ------------------- --*--------------- ! 

### SO4 = 

### 0.48, 2 ! NO3 = 

### 0.48, 2 ! I 

### PM800 = 

### 8, 0 ! I 

### PM425 = 

### 4.25, 0 ! ! 

### PM187 = 

### 1.875, 0 ! 

### ! PM112 = 

### 1.125, 0 ! 1 

### PM081 = 

### 0.8125, 0 ! 

### ! PMOS6 = 

### 0.5625, 0 !

!END!



### INPUT GROUP: 

### 9 

### -- Miscellaneous dry deposition parameters



## --------------

Reference cuticle resistance (s/cm) (RCUTR) 

### ~efault: 30 ! RCUTR 

### = 30.0 ! Reference ground resistance (s/cm) (RGR) 

### Default: 10 ! RGR = 

### 10.0 ! 

### Reference pollutant reactivity (REACTR) 

### Default: 8 ! 

### REACTR = 

### 8.0 !

Number of particle-size intervals used to evaluate effective particle deposition velocity (NINT) 

### Default: 9 ! 

### NINT = 

### 9 !
Vegetation state in unirrigated areas (IVEG) 

### Default: 1 ! IVEG = 

### 1 ! IVEG=l for active and unstressed vegetation IVEG=2 for active and stressed vegetation IVEG=3 for inactive vegetation

INPUT GROUP: 

### 10 

## -- wet Deposition Parameters



## ---------------

Scavenging Coefficient 

### -- units: (set)**(-1)


### Page 17NOX = HN03 = NO3 = PM800 = PM425 = PM187 = PMllZ = PM081 = PM056 = A 

### pendixBCohutta2OOlC~~~u~~.1~~ Liqui

# S Precip. Frozen Preci 

### w.

INPUT GROUP: 11 

## -- Chemistry Parameters



### ---------------

Ozone data input option (MOZ) Default: 

### 1 ! MOZ = 

### 1 ! (used only 

### if MCHEM = 1, 3, or 4) 0 = use a monthly background ozone value 

### 1 = read hourly ozone concentrations from the 0ZONE.DAT data file
Monthly ozone concentrations (used only if MCHEM = 

### 1, 

### 3, or 

### 4 and

MOZ = 0 or MOZ = 

### 1 and all hourly 03 data missing) (BCK03) in ppb Default : 12*80. 

### ! BCK03 = 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80 !
Monthly ammonia concentrations (used only 

### if MCHEM = 1, or 3) (BCKNH3) in ppb Default: 12*10. 

### ! BCKNH3 = 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5 !
Nighttime so2 loss rate (RNITE~) in percent/hour Default: 0.2 

### ! RNITE~ = 0.2 

### !
Nighttime NOX loss rate (RNITEZ) in percent/hour oefault: 2.0 ! RNITEZ = 2.0 

### !
Nighttime HN03 formation rate (RNITE3) in percent/hour oefault: 2.0 ! RNITE3 = 2.0 

### !

H202 data input option (MH202) Default: 

### 1 ! MHz02 = 

### 1 ! (used only 

### if MAQCHEM = 

### 1) 0 = use a monthly background HZ02 value 

### 1 = read hourly HZ02 concentrations from the HZO2.DAT data file
Monthly HZ02 concentrations (used only 

### if MQACHEM = 

### 1 and

MHz02 = 0 or MHz02 = 

### 1 and all hourly HZ02 data missing) (BCKH202) in ppb oefault: 12*1. ! BCKHZOZ = 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00 

### !

Page 18

### AppendixBCohutta2OOl~~~~u~~.1~~



## --- Data for SECONDARY ORGANIC 

### AEROSOL (SOA) 

### Option (used only 

### if MCHEM = 

### 4)

The 

### Son module uses monthlv values of: Fine particulate concentration in ug/mA3 (BCKPMF) Organic fraction of fine particulate (OFRAC) VOC 

### / NOX 

### ratio (after reaction) (VCNX) 

### to characterize the air mass when computing the formation of 

### son from 

### voc emissions. Typical values for several distinct air mass types are:

Month 

### 1 2 3 4 

### 5 6 

### 7 8 9 10 11 12 Jan ~eb Mar Apr May iun 1 ~ug sep Oct NOV 

### oec
Clean Continental BCKPMF1.1.1.1.1.1.1.1.1.1.1.1. 

### OFRAC 

### .I5 .I5 .20 .20 .20 .20 .20 .20 .20 .20 .20 .15 VCNX 

### 50. 50. 50. 50. 50. 50. 50. 50. 50. 50. 50. 50.

Clean Marine (surface) BCKPMF 

### .5 .5 .5 .5 .5 .5 .5 .5 

### .5 .5 .5 .5 OFRAC 

### .25 .25 .30 .30 .30 .30 .30 .30 .30 .30 .30 .25 VCNX 

### 50. 50. 50. 50. 50. 50. 50. 50. 50. 50. 50. 50.

urban 

## - low biogenic (controls present) BCKPMF 

### 30. 30. 30. 30. 30. 30. 30. 30. 30. 30. 30. 30. OFRAC 

### .20 .20 .25 .25 .25 .25 .25 .25 .20 .20 .20 .20 VCNX 

### 4.4.4.4.4.4.4.4.4.4.4.4.

urban - 

### high biogenic (controls present) BCKPMF 

### 60. 60. 60. 60. 60. 60. 60. 60. 60. 60. 60. 60. OFRAC 

### .25 .25 .30 .30 .30 .55 .55 .55 .35 .35 .35 .25 VCNX 

### 15. 15. 15. 15. 15. 15. 15. 15. 15. 15. 15. 15. Regional Plume BCKPMF 

### 20. 20. 20. 20. 20. 20. 20. 20. 20. 20. 20. 20. OFRAC 

### .20 .20 .25 .35 .25 .40 .40 .40 .30 .30 .30 .20 VCNX 

### 15. 15. 15. 15. 15. 15. 15. 15. 15. 15. 15. 15.

urban 

## - no controls present BCKPMF 

### 100. 100. 100. 100. 100. 100. 100. 100. 100. 100. 100. 100. OFRAC 

### .30 .30 .35 .35 .35 

### .55 .55 .55 .35 .35 .35 .30 VCNX 

### 2.2.2.2.2.2.2.2.2.2.2.2.

Default: Clean Continental ! BCKPMF = 

### 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00, 1.00

! 

### OFRAC = 

### 0.15, 0.15, 0.20, 0.20, 0.20, 0.20, 0.20, 0.20, 0.20, 0.20, 0.20, 0.15 !

! VCNX = 

### 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0 !

INPUT GROUP: 

### 12 

## -- Misc. oispersion and Computational Parameters

Horizontal size of puff 

### (m) beyond which time-dependent dispersion equations (Heffter) Page 

### 19

### Appendix~Cohutta2001~~~~~~~.1~~
are used to determine sigma-y and sigma-z (SYTDEP) Default: 550. ! SYTDEP = 5.5~02 

### !

Switch for using Heffter equation for sigma z as above (0 = Not use Heffter; 

### 1 = use Heffter (MHFTSZ) Default: 0 

### ! MHFTSZ = 0 !
stability class used to determine plume growth rates for puffs above the boundary layer (ISUP) oefault: 

### 5 ! ISUP= 

### 5 !

vertical dis ersion constant for stable conditions 

# &I in Eqn. 2.7-3) (coNK1) Default: 0.01 

### ! CONK1 = 0.01 

### !

vertical dispersion constant for neutral/ unstable conditions (k2 in Eqn. 2.7-4) (CONK2) Default: 0.1 ! CONK2 = 0.1 !

Factor for determining Transition- oint from Schulman-Scire to Huber-Snyder Buiyding Downwash scheme (ss used for Hs 

### < Hb 

## + TBD 

## * HL) (TBD) Default: 0.5 

### ! TBD=O.S ! TBD 

### < 0 

### ==> always use Huber-Snyder TED = 

### 1.5 ==> always use Schulman-Scire TBD = 0.5 ==z ISC Transition-point
Range of land use categories for which urban dispersion is assumed (IURB1, IURB2) Default: 10 

### ! IURBl = 10 !

site characterization parameters for single-point Met data files 

### --------- (needed for METFM = 2,3,4,5)

Land use category for modeling domain (ILANDUIN) Default: 20 

### ! ILANDUIN = 20 

### !
Roughness length (m) for modeling domain (ZOIN) Default: 0.25 ! ZOIN = 0.25 !

Leaf area index for modeling domain (XLAIN) Default: 3.0 

### ! XLAIIN = 3.0 !

Elevation above sea level 

### (m) (ELEVIN) Default: 0.0 ! ELEVIN = 0.0 !

Latitude (degrees) for met location (XLATIN) Default: -999. ! XLATIN = -999. 

### !
Longitude (degrees) for met location (XLONIN) Default: -999. 

### ! XLONIN = -999. 

### !
specialized information for interpreting single-point Met data files 

### -----
Anemometer height 

### (m) (used only 

### if METFM = 2,3) (ANEMHT) Default: 10. 

### ! ANEMHT = 10. 

### !

Form of lateral turbulance data in PROFILE.DAT file (used only 

### if METFM = 4,s or MTURBVW = 

### 1 or 3) (ISIGMAV) ~efault: 

### 1 ! ISIGMAV = 

### 1 ! 0 = read sigma-theta 

### 1 = read sigma-v
Page 20

### AppendixBCohutta2OOlCALPUFF.INP choice of mixing heights (used only 

### if METFM = 4) (IMIXCTDM) Default: 0 0 = read PREDICTED mixing hei hts 1 = read OBSERVED mixing heig 

# ! ts

Maximum length of a slug (met. grid units) (XMXLEN) Default: 1.0

Maximum travel distance of a puff/slug (in grid units) during one sampling step (XSAMLEN) Default: 1.0

Maximum Number of sluas/~uffs release from 

# -

.. one source during one time step (MXNEW) Default: 99

Maximum Number of sampling steps for
one puff/slug during one time step (MXSAM) Default: 99

Number of iterations used when computing the transport wind for a sampling step that includes gradual rise (for CALMET and PROFILE winds) (NCOUNT) Default: 2

Minimum sigma y for a new puff/slug (m) (SYMIN) Default: 1.0

Minimum sigma z for a new puff/slug 

### (m) (SZMIN) Default: 1.0 ! IMIXCTDM = 0 !

! XMXLEN = 

### 1. !

! XSAMLEN = 

### 1. !

! MXNEW = 99 

### !

! MXSAM = 99 !

! NCOUNT = 2 !



### ! SYMIN = 

### 1. !

! SZMIN = 

### 1. !

Default minimum turbulence velocities sigma-v and sigma-w for each stability class over land and over water (m/s) (SVMIN(12) and SWMIN 

### (12))

Default SVMIN : .SO, .50, .50, .50, .50, .50, .37, .37, .37, .37, .37, .37 Default SWMIN : .20, .12, .08, .06, .03, .016, .20, .12, .08, .06, .03, .016



### ! SVMIN = 0.500, 0.500, 0.500, 0.500, 0.500, 0.500, 0.370, 0.370, 0.370, 0.370, 0.370, 0.370 ! ! SWMIN = 0.200, 0.120, 0.080, 0.060, 0.030, 0.016, 0.200, 0.120, 0.080, 0.060, 0.030, 0.016 

### !
Divergence criterion for dw/dz across puff used to initiate adjustment for horizontal convergence (1/s) Partial adjustment starts at 

### cDIv(l), and full adjustment is reached at CDIV(~) (cDIv(2)) Default: 0.0,o.o ! CDIV = .0, .O !

Minimum wind speed (m/s) allowed for non-calm conditions. Also used as minimum speed returned when using power-law extrapolation toward surface (WSCALM) Default: 0.5 ! WSCALM = 0.5 !

Maximum mixing height 

### (m) (XMAXZI)Minimum mixing height (m) (XMINZI) Default: 50. 

### ! XMINZI = 50.0 !

Default wind speed classes 

## -- 5 upper bounds (m/s) are entered; the 6th class has no upper 

### limit (WSCAT(~)) Default : ISC RURAL : 1.54, 3.09, 5.14, 8.23, 10.8 (10.8+)

wind Speed class : 

### 1



### --- 2



### --- 3



## --- 4



### --- 5



## --- ! WSCAT = 1.54, 3.09, 5.14, 8.23, 10.80 !

Default wind speed profile power-law exponents for stabilities 1-6 (pLXO(6)) Default : ISC RURAl va111e<

~ - - - . . -

. . . -

. 

## -

. 

## - - - 1st RURAL : .07, .07, 

## .lo, .15, .35, .55 ISC URBAN : .15, .15, .20, .25, .30, .30 Stability Class : A B 

### c D



## --- --- E



### --- F



## - - - - - - - - - ! PLXO = 0.07, 0.07, 0.10, 0.15, 0.35, 0.55 !

Default potential temperature gradient for stable classes E, F (deg~/m) (PTGO(2)) Default: 0.020, 0.035 ! PTGO = 0.020, 0.035 !

Default lume path coefficients for each stagility class (used when option for partial plume height terrain adjustment is selected -- MCTAD7=3) (PPC(6)) Stability class : A B C D E F Default PPC: .50, .SO, .SO, .SO, .35, .35



## --- --- --- --- --- 

## --- ! PPC = 0.50, 0.50, 0.50, 0.50, 0.35, 0.35 !
Slug-to-puff transition criterion factor equal to sigma-y/length of slug (SL2PF) Default: 10. ! SLZPF = 5.0 ! puff-splitting control variables



## ........................ VERTICAL SPLIT

Number of puffs that result every time a puff is split 

## - nsplit=2 means that 

### 1 puff splits
intn 7 

# - (NSPLIT) Default: 3 ! NSPLIT = 3 ! Time(s) of a day when split puffs are eligible to be split once aqain: this is tv~icallv set once per day, around-sunset before nocturnal shear develops. 24 values: 0 is midnight (00:OO) and 23 is 

### 11 PM (23:OO) O=do not re-split l=eligible for re-split IRES SPLIT(^^)) ~efault: Hour 17 = 

### 1 IRE~PL~T=0,~,~,~,0,~,~,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,

split is allowed only if last hour's mixing Page 22AppendixBcohutta2001~LPUFF.INP height 

### (m) exceeds a minimum value (ZISPLIT) Default : 100.
split is allowed only 

### if ratio of last hour's mixing ht to the maximum mixing ht experienced by the puff is less than a maximum value (this postpones a split until a nocturnal layer develops) (ROLDMAX) Default: 0.25

HORIZONTAL SPLIT



## ----------------

Number of puffs that result every time a puff is split 

## - nsplith=5 means that 

### 1 puff splits into 

### 5 (NSPLITH) Default: 

### 5

Minimum sigma- (Grid Cells Units) of puff before it may e split (SYSPLITH) 

# i:

Default: 1.0

Minimum puff elongation rate (SYSPLITH/hr) due to wind shear, before it may be split (SHSPLITH) Default: 2.

Minimum concentration (g/mA3) of each species in puff before it may be split Enter array of NSPEC values; 

### if a single value is entered, it 

### will be used for ALL species (CNSPLITH) Default: 1.0~-07 I
Integration control variables



## ........................ Fractional convergence criterion for numerical SLUG sampling integration (EPSSLUG) Default: 1.0e-04

Fractional converaence criterion for numerical AREA source integration (EPSAREA) Default: l.Oe-OE
Trajectory step-length (m) used for numerical rise integration (DSRISE) Default: 1.0 ! ZISPLIT = 100.0 !



### ! ROLDMAX = 0.25 !



### ! NSPLITH = 5 

### !

! SYSPLITH = 1.0 

### !

! SHSPLITH = 2.0 !



### ! CNSPLITH = 1.OE-07



### ! EPSSLUG = 1.OE-04 

### !

! EPSAREA 

### = 1.OE-06 

### !

! DSRISE = 1.0 !
Boundary Condition (BC) Puff control variables



## ........................
Minimum height (m) to which BC puffs are mixed as they are emitted (MBCDN=~ ONLY). 

### AC~U~J heiaht is reset to the current mixino heiaht



### ., 2 - it the release point 

### if greater than this minimum. (HTMINBC) ~efault: 500. 

### ! HTMINBC = 500.0 

### !

search radius (Itm) about a receptor for sampling nearest BC BC puffs are typically emitted with a spacing of one grid ce 

### 1 length, so the search radius should be greater than DGRIDKM. 

# Tuff.

(RSAMPBC) Default: 10. ! RSAMPBC = 10.0 

### ! Near-Surface depletion adjustment to concentration profile used when sampling BC puffs? (MDEPBC) ~efault: 1 

### ! MDEPBC = 

### 1 ! 0 = Concentration is NOT adjusted for depletion 

### 1= Adjust concentration for depletion Page 23INPUT GROUPS: 

### 13a, 13b, 13c, 13d 

## -- Point source parameters

subgroup (13a)

Number of point sources with parameters provided below (NPT1) NO 

### default ! 

### NPTl = 

### 6 !



### units used for point source emissions below (IPTU) 

### Default: 1 

### ! IPTU = 

### 3 ! 1 = 

### g/s 2 = 3 = 

### 4 = 

### tons/yr 5 = 

### odour unit 

## * m**3/s (vol. flux of odour compound) 6 = 

### odour unit 

## * m**3/min 7 = 

### metric tons/yr

Number of source-species combinations with variable emissions scaling factors provided below in (13d) (NSPT~) 

### Default: 0 ! NSPT~ = 

### 0 !

Number of point sources with variable emission parameters provided in external file (NPT2) NO 

### default ! NPT2 = 

### 0 !


### (If NPTZ > 0, these point
source emissions are read from the file: PTEMARB.DAT)

!END!

C 

### source X Y 

### stack ~ase stack Exit Exit ~misSi0n NO. 

### Coordinate Coordinate Height Elevation Diameter Vel



# . Temp. b
Bldg.

Dwash

1 ! SRCNAM = BCZl 

### ! 1 ! X = 

### 1118.2763, -572.1253, 13.72, 228.6, 0.69, 15.24, 477.59, 

# - u., 12.709, 0, 3.58, 0, 0, 0, 0, 

### 0, 0, 6.425, 6.425! 1 ! SIGYZI = 

### 0. 0 ! 

### 1 ! ZPLTFM = 

### .O ! Page 24

### AppendixBCohutta2001C~LPUFF.INP I! FMFAC= O! !END! 

### 2 ! SRCNAM = BDZl 

### ! 2 ! X = 

### 1118.2958, -572.1236, 14.63, 229.51, 0.61, 3.96, 422.04, 

## o., 4.686, 0, 1.32, 0, 0, 0, 0, 0, 0, 5.19, 5.19! 2 ! SIGYZI = 

### 0, 0 ! 

### 2 ! ZPLTFM = .O 

### ! 2 ! FMFAC = 

### 0 ! !END! 

### 3 ! SRCNAM = 

### BDZ3 ! 3 ! X = 

### 1118.2795, -572.1414, 14.63, 229.51, 0.41, 10.67, 422.04, O., 

### 3.124, 0, 0.88, 0, 0, 0, 0, 0, 0, 4.745, 4.745! 3 ! SIGYZI = 

### 0, 0 ! 

### 3 ! ZPLTFM = .O 

### ! 3! FMFAC= O! !END! 

### 4 ! SRCNAM = 

### BDZ4 ! 4 ! X = 

### 1118.2694, -572.1356, 9.14, 224.03, 0.51, 3.66, 394.26, 

# - u., 2.414, 4 ! SIGYZI = 

### 4 ! ZPLTFM = 

### 4 ! FMFAC = 

### 5 ! SRCNAM = 

### 5! X= 

### 477.59, 

## o., 305.966, 5 ! SIGYZI = 

### 5 ! ZPLTFM = 5 ! 

### FMFAC = 

### 6 ! SRCNAM = 

### 6! X= 

### 533.15, O., 

### 52.966, 6 

### ! SIGYZI = 

### 6 ! 

### ZPLTFM = 

### 6 ! FMFAC =



### a

Data for each source are treated as a separate input subgroup and therefore must end with an input group terminator.

SRCNAM 

### is a 12-character name for a source (NO 

### default) 

### x is an array holding the source data listed by the column headings (NO 

### default) SIGYZI 

### is an array holding the initial sigma-y and sigma-z 

### (m) (Default: 0.,0.) FMFAC 

### is a vertical momentum flux factor (0. or 1.0) used to represent the effect of rain-caps or other physical configurations that reduce momentum rise associated with the actual exit velocity. (Default: 1.0 

## -- full momentum used) ZPLTFM is the platform height (m) 

### for sources influenced by an isolated structure that has a significant o en area between the surface and the bulk of the srructure, suc

# r: as an offshore oil platform. The Base Elevation is that of the surface (around or ocean). and the stack Height is the release height -above the Base (not above the platform). Building heights entered in subgroup 13c must be those of the buildings on the platform, measured from the platform deck. ZPLTFM is used only with MBDW=l (ISC 

### downwash method) for sources with building downwash. (Default: 0.0) Page 25

## - 0. = 

### No buildinq downwash modeled 1. = 

### Downwash modeled for buildings resting on the surface 2. = 

### Downwash modeled for buildings raised above the surface 

# CZPLTFM z 

### 0.) NOTE: must be entered as a REAL number (i.e., with decimal point)

C

### An emission rate must be entered for every pollutant modeled.

Enter emission rate of zero for secondary pollutants that are modeled, but not emitted. units are specified by IPTU (e.g. 1 

### for 9/51.

subgroup (13~)



## ---------------

BUILDING DIMENSION 

### DATA FOR SOURCES 

### SUBJECT TO DOWNWASH



## ....................................................... source a

NO. 

### Effective building height, width, length and X/Y 

### offset (in meters) every 10 degrees. LENGTH, XBADJ, and YBAD3 

### are only needed for MBDW=2 (PRIME 

### downwash option)



### ------ -------------------------------------+------------------------------



### a

Building height, width, length, and x/Y 

### offset from the source are treated
as a separate input subgroup for each source and therefore must end with
an input group terminator. The X/Y 

### offset is the position, relative to the stack, of the center of the upwind face of the projected building, with the x-axis pointing along the flow direction.

a POINT SOURCE: 

### VARIABLE EMISSIONS 

### DATA

use this subgroup to describe temporal variations in the emission rates given in 13b. Factors entered multiply the rates in 13b. Skip sources here that have constant emissions. For more elaborate variation in source parameters, use PTEMARB.DAT 

### and NPTZ > 0.

IVARY 

### determines the type of variation, and is source-specific: (IVARY) 

### Default: 0 0 = 

### Constant 1 = 

### Diurnal cycle (24 scaling factors: hours 1-24) 2 = 

### Monthly cycle (12 scaling factors: months 1-12) 3 = 

### Hour & 

### Season (4 groups of 24 hourly scaling factors, where first group is DEC-JAN-FEE) 4 = 

### speed & 

### Stab. (6 groups of 6 scaling factors, where first group is Stability class A, 

### and the speed classes have upper bounds (m/s) defined in Group 12 5 = 

### Temperature (12 scaling factors, where temperature classes have upper bounds 

### (c) of: 0, 5, 10, 15, 

### 20, 25, 30, 35, 40, 45, 50, 50+)

Page 26

### a Data for each species are treated as a separate input subgroup and therefore must end with an input group terminator.

INPUT GROUPS: 

### 14a, 14b, 14c, 14d 

## -- Area source parameters



### ................................

subgroup (14a)



### umber of polygon area sources with parameters specified below (NAR~) NO 

### default ! NAR~ = 0 !



### units used for area source emissions below (IARU) 

### Default: 1 ! IARU = 

### 1 ! 1 = 

### g/m2*2/s 2 = 

### kg/m**2/hr 3 = 

### 1b/me*2/hr 4 = 

### tons/m**2/yr 5 = 

### Odour unit 

## * m/s (vol. flux/m**2 of odour compound) 6 = 

### odour unit 

## * m/min 7 = 

### metric tons/m**2/yr

Number of source-species combinations with variable emissions scaling factors provided below in (14d) (NSAR~) 

### oefault: 0 

### ! NSAR~ = 0 

### !

Number of buoyant polygon area sources with variable location and emission parameters (NARZ) NO 

### default ! 

### NAR2 = 0 ! 

### (If NAR2 z 0, 

### ALL parameter data for these sources are read from the file: BAEMARB.DAT)



### a AREA SOURCE: 

### CONSTANT DATA



### a

Data for each source are treated as a separate input subgroup and therefore must end with an input group terminator. h 

## -

An emission rate must be entered for every pollutant modeled. Enter emission rate of zero for secondary pollutants that are Page 27

### Appendix~Cohutta2OOlCALPUF~.1~~ modeled, but not emitted. units are specified by IARU 

### (e.g. 1 for g/m2*2/s).

Subgroup (14c)

COORDINATES (km) FOR EACH 

### VERTEX(4) OF EACH 

### POLYGON



### ........................................................

- 

## - - - - - - - a

Data for each source are treated as a separate input subgroup and therefore must end with an input group terminator.



## --------------- subgroup (14d)

a AREA SOURCE: 

### VARIABLE EMISSIONS 

### DATA



### ......................................
Use this subgroup to describe temporal variations in the emission rates given in 14b. Factors entered multiply the rates in 14b. skip sources here that have constant emissions. For more elaborate variation in source parameters, use BAEMARB.DAT 

### and NAR2 > 0.

IVARY 

### determines the type of variation, and is source-specific: (IVARY) 

### Default: 0 0 = 

### constant 1 

### = Diurnal cycle (24 scaling factors: hours 1-24) 2 = 

### Monthly cycle (12 scaling factors: months 1-12) 3 = 

### Hour & 

### Season (4 groups of 24 hourly scaling factors, where first group is DEC-JAN-FEB) 4 = 

### speed & 

### Stab. (6 groups of 6 scaling factors, where first arouw is stabilitv class A. 

### and the speed classes hive upper' bounds (m/s) defined in Group 12 5 = 

### Temperature (12 scaling factors, where temperature classes have upper bounds 

### (c) of: 0, 5, 

### 10, 15, 20, 25, 30, 35, 

### 40, 45, 50, SO+)



### a

Data for each species are treated as a separate input subgroup and therefore must end with an input group terminator.



## --------------- subgroup (15a) Page 28

### Number of buoyant line sources with variable location and em1 ssion parameters (NLN2) No default ! NLN~ = 0 

### !


### (~f NLN2 > 0, ALL 

### parameter data for these sources are read from the file: LNEMARB.DAT)



### Number of buoyant line sources (NLINES) No default ! 

### NLINES = 0 

### !

units used for line source emissions below (ILNU) 

### Default: 1 ! 

### ILNU = 1 ! 

### 1= 

### 91s 2 = 

### Whr 3 = 

### 1b/hr 4 = 

### tons/yr 5 = 

### odour unit 

## * m**3/s (vol



## . flux of odour compound) 6 = 

### odour unit 

## * m**3/min 7 = 

### metric tons/yr

Number of source-species combinations with variable emissions scaling factors provided below in (15c) (NSLN~) 

### Default: 

### 0 ! NSLN~ = 0 !



### Maximum number of segments used to model each line (MXNSEG) 

### Default: 

### 7 ! MXNSEG = 

### 7 !



### The following variables are required only 

### if NLINES z 0. 

### They are used in the buoyant line source plume rise calculations.

Number of distances at which transitional rise is computed
Average building length (XL)

Average building height (HBL)

Average bui 1ding width (WBL)

Average line source width (WML) Default: 6 

### ! NLRISE = 6 !

NO 

### default ! XL 

### = 0. ! (in meters)

No default ! 

### HBL 

### = 0. 

### ! (in meters)

NO 

### default ! 

### WBL 

### = 0. 

### ! (in meters)

NO 

### default ! 

### WML 

### = 0. 

### ! (in meters)
Average separation between buildings (DXL) 

### No default ! 

### DxL = 0. 

### ! (in meters)
Average buoyancy parameter (FPRIMEL) No default ! FPRIMEL = 0. 

### ! 

### (in m**4/s*"3)

Subgroup (15b)

Page 29

### a Data for each source are treated as a separate input subgroup and therefore must end with an input group terminator.

b

An emission rate must be entered for every pollutant modeled. Enter emission rate of zero for secondary pollutants that are modeled, but not emitted. Units are specified by ILNTU (e.g. 1 for g/s).

subgroup (15c)



### ---------------
a BUOYANT LINE SOURCE: 

### VARIABLE EMISSIONS 

### DATA

use this subgroup to describe temporal variations in the emission rates given in 15b. Factors entered multiply the rates in 15b. skip sources here that have constant emissions.

IVARY 

### determ (IVARY) 

### 0 = 

### ines the type of variation, and is source-specific: Default: 0 constant Di~rna1 cycle (24 scaling factors: hours 1-24) Monthly cycle (12 scaling factors: months 1-12) Hour & 

### Season (4 groups of 24 hourly scaling factors, where first group is DEC-JAN-FEE) 

### speed & 

### stab. (6 groups of 6 scaling factors, where first arouD is stabilitv class A. 

### and the speed classes have upper' bounds (m/s) defined in Group 12 Temperature (12 scaling factors, where temperature classes have upper bounds (C) 

### of: 0, 5, 10, 15, 

### 20, 25, 30, 35, 40, 45, 50, 50+)

a

Data for each species are treated as a separate input subgroup and therefore must end with an input group terminator.

INPUT GROUPS: 

### 16a, 16b, 16c 

## -- volume source parameters



## ...........................

Page 30

### Number of volume sources with parameters provided in 16b,c (NVL~) NO 

### default ! 

### NVLl = 0 

### !

units used for volume source emissions below in 16b (IVLU) 

### Default: 1 ! IVLU = 

### 1 ! 

### 1 = 

### 2 = 

### 915

3 = 

### kg/hr Ib/hr 4 = 

### tons/yr 5 = 

### odour unit 

## * m4*3/s (vol



# . flux of odour compound) 6 = 

### Odour unit 

## * m**3/min 

### 7 = 

### metric tons/yr

Number of source-species combinations with variable emissions scaling factors provided below in (16c) (NSVL~) 

### Default: 0 ! NSVL~ = 0 !



### Number of volume sources with variable location and emission parameters (NVLZ) No default ! NVLZ = 0 ! 

### (~f NVLZ z 0, ALL 

### parameter data for these sources are read from the VOLEMARB.DAT 

### file(s) )



### subgroup (16b)

a

Data for each source are treated as a separate input subgroup and therefore must end with an input group terminator.

b

An emission rate must be entered for every pollutant modeled. Enter emission rate of zero for secondary pollutants that are modeled, but not emitted. units are specified by IVLU 

### (e.g. 1 for g/s).

subgroup (16c)

a VOLUME SOURCE: VARIABLE EMISSIONS 

### DATA



### ........................................

use this subgroup to describe temporal variations in the emission rates given in 16b. Factors entered multiply the rates in 16b. skip sources here that have constant emissions. For more elaborate Page 31

### AppendixBCohutta2OOlCALPUFF.1~~ variation in source parameters, use VOLEMARB.DAT and NVL2 

### > 0.

IVARY determines the type of variation, and is source-specific: (IVARY) Default: 0 0 = Constant 

### 1 = Diurnal cycle (24 scaling factors: hours 1-24) 2 = Monthly cycle (12 scaling factors: months 1-12) 3 = Hour & Season (4 groups of 24 hourly scaling factors, where first group is DEC-JAN-FEB) 4 = speed & stab. (6 groups of 6 scalin factors, where first group is stabi

# ? ity class A, and the speed classes have upper bounds (m/s) defined in Group 12 

### 5 = Temperature (12 scaling factors, where temperature classes have upper bounds 

### (c) of: 0, 

### 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 50+)



## - - - - - - - -
a

Data for each species are treated as a separate input subgroup and therefore must end with an input group terminator.

INPUT GROUPS: 17a & 17b -- Non-gridded (discrete) receptor information

subgroup (17a)

Number of non-gridded receptors (NREC) No default 

### ! NREC 

### = 220 

### !

!END!

a NON-GRIDDED (DISCRETE) RECEPTOR DATA



### x Y Ground Height b Receptor coordinate coordinate Elevation nbove Ground NO. (km) (km)

0 

### ! !END! 0 

### ! !END!

0 

### ! !END! O!! O!! O!! O!!

0 

### ! 0 

### !

0 

### ! 0 ! 0 ! END 

### ! END! END ! END ! !END! !END!

!END! !END! !END! Page 32

# . INP !END! !END! !END!

!END! !END! !END! !END! 

### ! END ! !END! !END! 

### ! END 

### ! !END! !END! END! END 

### ! !END!

: END ! !END! !END! 

### ! END 

### ! ! END ! !END! !END! 

### ! END 

### !

!END! 

### ! END !END! 

### ! END ! ! END 

### ! !END! !END! !END! !END! !END! 

### ! END ! !END! !END! 

### ! END ! !END! !END! !END! !END! 

### ! END ! !END! !END! 

### ! END ! !END! !END! ! END 

### ! ! END !END! 

### ! END 

### ! !END! !END! !END! !END! !END! !END! !END! !END 

### ! !END! !END! 

### ! END !Page 

### 35 INPData for each receptor are treated as a separate input subgroup and therefore must end with an input group terminator.

Receptor height above ground is optional. 

### If no value is entered, the receptor is placed on the ground.

Page 36Chemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol 

# DRAFT

Appendix C Example of Proposed POSTUTIL Inpu t File for Cohutta Wilderness Area (2001)

### Appendixccohutta2OOlPOS~~TI~.1~~ VISTAS BART 

## - chem Products Gorp. 

## - 4km Domain No. 

### 4 

## - COHUTTA 

## - 2001 Computational Grid 

### x (Easter1 ): From 76 To 127 Computational Grid Y (Northeryy): From 136 To 208


## ---------------- ~un title (3 lines) 

## -------------------

POSTUTIL MODEL CONTROL FILE



## ...........................

Output Files



## ------------

File Default ~ile Name



### ---- 

## ----------------- List File POSTUTIL.LST 

### ! UTLLST =POSTUTIL.LST Data File MODEL. DAT ! UTLDAT =PostUti~ConcOu

Input Files


## -----------

A time-varying file of "background" concentrations can be included when the ammonia-limiting method CALM) for setting the HN03/N03 concentration partition is accomplished in 

### 1 step. This option is selected by settin MNITRATE=~in Input Group 

### 1. Species required in t

# z
e "background" concentration file are: 504, N03, HN03 and TNH3 (total NH3 = NH3gaseous 

## + NH3particulate).

File Default File Name



### ---- 

### ----------------- BCKG File BCKGALM.DAT 

## * BCKGALM =BCKGALM.DAT 5

A number of CALPUFF data files may be processed in this application. The files may represent individual CALPUFF simulations that were made for a specific set of s ecies and/or sources. Specify the total number of CALPUFF runs you wisE to combine, and provide the filename for each in subgroup Ob.

Number of CALPUFF data files (NFILES) Default: 1 

### ! NFILES = 

### 1 !

Meteorological data files are needed for the HN03/N03 partition option. Three types of meteorological data files can be used: METFM 

### = 0 

## - CALMET.DAT METFM= 

### 1 - 1-D file with RH, Temp and Rhoair timeseries METFM= 

### 2 

## - 2-0 files with either Rh, Temp or Rhoair in each (3 2-D files are needed) 

### he default is to use CALMET.DAT files.

Page 

### 1

### ~ppendixCCohutta2OOlPoSTUTIL.1NP Multiple meteorological data files may be used in se uence to span the processing period. specify the number of time- erio% files (NMET) that you need to use, and provide a filename for eacR in subgroup Ob.



## - NMET is 0 

### if no meteorological files are provided



## - NMET is 1 

### if METFM=~ (multiple file feature is not available)



## - NMET is 

### 1 or more 

### if METFM=O or 2 (multiple CALMET files or 2DMET files)

Default: 0 

### ! NMET = 0 !



### ~ll filenames will be converted to lower case 

### if LCFILES = T otherwise, 

### if LCFILES = F, filenames will be converted to UPPER CASE

Convert filenames to lower case? Default: T ! LCFILES = T 

### ! T = lower case F = UPPER CASE

NMET CALMET Data Files (METFM=O) 

### :

~nput File oefault File Name



### ---------- ----------------- 1 MET. DAT 

# " UTLMET =CALMET.OAT 

# " "END*

NMET 

### 1-D Data Files CMETFM=~) :

NMET 2-D Data Files of Each Type (METFMd):

~nput File oefault File Name



## ---------- 

## ----------------- 1 RHUMD. DAT 

## * MZDRHU = RELHUM.DAT 

# " "END* 

### 1 TEMP. DAT 

# " MZDTMP = TEMP.DAT 

# " "END" 

### 1 RHOAIR.DAT 

## * MZDRHO = RHOAIR.DAT 

# " "END"

NFILES CALPUFF Data Files:

~nput File ~efault File Name



### ---------- ----------------- 0 CALPUFF.DAT 

### ! MODDAT = ..\Cohutta\conc.dat 

### ! !END!

Note: provide NMET lines of the form 

# " UTLMET 

### = name 

# " "END"

Or 

# " METlD 

### = name 

# " "END" Page 

### 2

### or 

## * M2DRHU = 

### name 

## * *END* 

### (and) 

## * MZDTMP = 

### name 

## * "END* 

### (and) 

## * MZDRHO = 

### name 

## * *END*



### and NFILES lines of the form 

## * MODDAT = 

### name 

## * *END*



### where the 

## * should be replaced with an exclamation point, the special delimiter character.

INPUT GROUP: 

### 1 

## -- General run control Darameters

starting date:



### umber of periods

Number of species

Number of species

Number of species Year (ISYR) 

## -- No default Month (ISMO) 

## -- No default Day (ISDY) 

## -- No default Hour (ISHR) 

## -- No default

to process (NPER) 

### -- NO 

### default

to process from CALPUFF 

### runs (NSPECINP) 

## -- NO 

### default

to write to output file (NSPECOUT) 

## -- No default

to compute from those modeled (must be no greater than NSPECOUT) (NSPECCMP) 

## -- NO 

### default ! ISYR = 

### 2001 ! 

### ! ISMO = 1 

### ! ! ISDY = 

### 1 ! ! ISHR = 

### 1 !



### ! NPER = 

### 8760 !



### ! NSPECINP = 

### 11 !

! 

### NSPECOUT = 

### 6 !

! NSPECCMP = 

### 1 !

when multiple files are used, a species name may appear in more than
one file. Data for this species will 

### be summed (appropriate 

### if the CALPUFF,runs use different source groups). 

### If this summing is not appropriate, remove duplicate species from the file(s).
stop run 

### if duplicate species names are found? (MDUPLCT) 

### Default: 0 ! MDUPLCT = 

### 0 ! 

### 0 = 

### no 

### (i .e., duplicate species are summed) 1 = 

### yes (i .e., run is halted)

Data for each species in a CALPUFF data file may also be scaled as they are read. This can be done to alter the emission rate of all soui-ces that were modeled in a particular CALPUFF ap lication. The scaling factor for each species is entered in Su group (Zd), for each file for which scaling is requested. 

# 1

Number of CALPUFF data files that will be scaled (must be no greater than NFILES) (NSCALED) 

### Default: 0 ! 

### NSCALED = 

### 0 !

~mmonia-Limiting Method Option to recompute the HN03/N03 concentration p a r t i t i o n p r i o r to performing other actions is controlled by MNITRATE. This option will NOT alter any deposition fluxes contained in the CALPUFF file(s). Three partition selections are provided. 

### he first two are typically used in sequence (POSTUTIL 

### is run more than once). The first selection (MNITRATE=l) computes the partition for the TOTAL (all sources) concentration fields (504, N03, HN03; NH3), and the second (MNITRATE=2) Page 3

### AppendixCCohutta2001POsTu~1~.I~~ uses this partition (from the previous application of POSTUTIL) 

### to compute the partition for individual source groups. 

### he third selection (MNITRATE=3) can be used instead in 

### a single POSTUTIL application 

### if a f i l e o f background concentrations is provided (BCKGALM 

### in Input Group 0).
Required information for MNITRATE=l includes: species N03, HN03, and SO4 NH~ 

### concentration(s) met. data file for RH 

### and T 

### ~equi red information for MNITRATE=2 includes : 

### species NO3 and HN03 for a source group species N03ALL and HN03ALL for all source groups, properly partitioned
Required information for MNITRATE=3 includes : 

### species N03, HN03, and SO4 for a source group species N03, HN03, 504 and TNH3 from the background BCKGALM file 

### If TNH3 is not in the background BCKGALM file, monthly TNH3 concentrations are used (BCKTNH~) TNH3= total NH3 = ~~3gaseous+~~3particulate


### Recompute the HN03/N03 partition for concentrations? (MNITRATE) 

### Default: 0 ! MNITRATE = 

### 0 ! 0 = 

### no 1 = 

### yes, for all sources combined 2 = 

### yes, for a source group 3 = 

### yes, ALM 

### application in one step

SOURCE OF AMMONIA:



### Ammonia may be available as a modeled species in the CALPUFF files, and it 

### may or may not be appropriate to use it 

### for repartitioning N03/H~03 (in Option MNITRATE=l or MNITRATE=3). Its use is contolled by NH3TYP. when NH3 is listed as a processed species in Subgroup 

### (Za), as one of the NSPECINP ASPECI 

### entries, and the right option is chosen for NH3TYP, the NH3 modeled values from the CALPUFF concentration files will 

### be used in the chemical equil 

### i brium calculation.

NH3TYP also controls when monthly background ammonia values are used. ~oth gaseous (NH3) and total (TNH3=NH3gaseous+NH3particulate) 

### ammonia can be provided monthly as BCKNH~/BCKTNH~.

it 

### is the inout source of Ammonia? wh: (NH~TYP) NO 

### Default 0 = 

### No backqround will be used. ONLY NH~ 

### from the concentration files listed in subqrouw (2a) as a processed species-will be used. (Cannot be used with MNITRATE=~)

1 = 

### NH3 Monthly averaged background (BCKNH~) 

### listed below will be added to NH3 from concentration files listed in subgroup 

### (2a)

2 = 

### NH3 from background concentration file BCKGALM will 

### be added to NH3 from concentration files listed in subgroup (2a) (ONLY possi bl e for MNITRATE=3)

3 = 

### NH3 Monthly averaged background (BCKNH3) listed below will 

### be used alone. Page 4! ASPECO = 

### ! ASPECO = ! ASPECO = ! ASPECO = 

### ! ASPECO =



### ------------- Subgroup (2c) 

### AppendixCCohutta2OOlPOS~u~1~.1~~ 504 ! !END! NOX ! !END! 

### ~~03 ! !END! 

### NO3 ! !END! SOA ! !END!



### The following NSPECCMP 

### species will 

### be computed by scaling and summing
one or more of the processed input species. Identify the name(s) of the computed species and provide the scaling factors for each of the NSPECINP 

### input species (NSPECCMP 

### groups of NSPECINPC~ 

### lines each):

! CSPECCMP = SOA 

### ! ! PM081 = 

### 1.0 ! ! PMOS6 = 

### 1.0 ! ! END 

### !



### ------------- subgroup (2d)



### -------------

Each species in NSCALED CALPUFF data files may be scaled before being processed (e.g., to change the emission rate for all sources modeled in the run that produced a data file). For each file, identify the file name and then provide the name(s) of the scaled species and the corresponding scaling factors (A,B 

### where x' = 

### Ax+B).



# --. " SO4 

## - L.,, C 

### HNO? = 

### n R 

## " NO! "END"



### Page 6

### Appendix~Cohutta20OlC~L~oST.1~~ 1) Visibility: daily visibility tabulations are always reported for the selected rece tors when ASPEC = VISIB. 

### In addition, any of tRe other tabulations listed below may be chosen to characterize the light extinction coefficients. [List file or Plot/nnalysis File]

2) Top 50 table for each averaging time selected [List file only] (LT50) 

## -- Default: T 

### ! ~~50 = F 

### !
3) Top 'N' 

### table for each averaging time selected [List file or plot file] (LTOPN) 

## -- Default: F ! LTOPN = F !



### -- Number of 'Top-N' values at each receptor selected (NTOP 

### must be 

### <= 4) (NTOP) 

### -- Default: 4 ! NTOP = 1 !



### -- Suecific ranks of 'TOW-N' 

### values reuorted (NTOP 

### values must be entered) (ITOP(4) array) 

## -- Default: ! ITOP = 

### 0 !



### 4) Threshold exceedance counts for each receptor and each averaging time selected [List file or Plot file] (LExcD) 

### -- Default: F ! 

### LEXCD = F !



## -- Identify the threshold for each averaging time by assigning a non-negative value (output units).



## -- Default: -1.0 Threshold for 1-hr averages (THRESHl) ! THRESH1 = 

### -1.0 ! ~hreshold for 3-hr averages (THRESH3) ! 

### THRESH3 = 

### -1.0 ! 

### ~hreshold for 24-hr averages (THRESH24) ! 

### THRESH24 = 

### -1.0 ! 

### ~hreshold for NAVG-hr averages (THRESHN) ! 

### THRESHN = 

### -1.0 !



## -- Counts for the shortest averaging period selected can be tallied daily, and receptors that experience more than NCOUNT counts over any NDAY 

### period will 

### be reported. This type of exceedance violation output is triggered only 

### if NDAY 

### > 0.

~ccumulation period(oays) (NDAY) 

## -- Default: 0 ! NDAY = 

### 0 ! 

### umber of exceedances allowed (NCOUNT) 

## -- Default: 

### 1 ! NCOUNT = 

### 1 !



### 5) selected day table(s)

~cho option 

## -- Many records are written each averaging period selected and output is grouped by day [List file or Plot file] (LECHO) 

## -- Default: F ! LECHO = F !



### Timeseries Option 

## -- Averages at all selected receptors for each selected averaging period are written to timeseries files. Each file contains one averaging period, and all receptors are written to a single record each averaging time. 

### [TSERIESASPEC-ttHR-CONCTSUNAM.DAT files] Page 8Chemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol 

# DRAFT

Appendix D Example of Proposed CALPOST Input File for Cohutta Wilderness Area (2001)VISTAS BART - chem products Corp. 4km Domain NO. 

### 4 - COHUTA 

## - 2001 Computational Grid X (~asterly): From 76 To 127 Computational Grid Y (Northerlv): From 136 TO 7nR

> 

# -

~ --.
. 

## - --- ---------------- Run title (3 

### lines)



### ..........................................

INPUT GROUP: 0 

## -- Input and Output File Names



### --------------

Input Files



### -----------

File



## ---- Default File Name



### ----------------- Conc/~ep Flux File MODEL. DAT 

### ! MODDAT = ..\~ost~til~onc~ut.dat ! Relative ~umidity File VISE .DAT 

# " VISDAT = ..\VISB.DAT 

# " Background Data File BACK. DAT "BACKDAT = 

## *
Transmissometer or VSRN



# .DAT "VSRDAT = 

# "
Newhelometer Data File or -. DATSAV weather Data File or Prognostic weather File
Output Files



### ------------

File



### ---- List File Default File Name



### ----------------- CALPOST.LST ! PSTLST = Cal~ost .?st !

Pathname for Timeseries ~ileS (blank) 

## " TSPATH = 

## * (activate with exclamation points on1y 

### if providing NON-BLANK character string)

pathname for Plot Files (blank) 

## * PLPATH = 

## * (activate with exclamation points only if

providing NON-BLANK character string)

user Character string 

### (u) to augment default filenames (activate with exclamation points only 

### if

providing NON-BLANK character string)

Timeseries TSERIESASPEC-ttHR-CONC-TSUNAM.DAT peak value PEAKvAL~ASPEC-~~HRCONCTSUNAM.DAT

Top ~th ~ank Plot RANK(ALL)ASPEC-ttHR-CONUUNAM.DAT

or 

### RANK(?~)-ASPEC-~~HRRCONC-TUNAM.GRD

Exceedance Plot EXCEEDASPEC-ttHR-CONC2UNAM.DAT

or EXCEEDASPEC-ttHRCONC_XUNAM.GRD

~~~~~M~~~D~~~~~oo(UTC~~~~~)~LOO~ASPEC~~~HRCON~. DAT Page 

### 1

### visibility plot DAILY-VISIB-VUNAM.DAT 

# " VUNAM = % 

### (~ai 1

y Peak summary)

nuxi 1 iary output Fi 1 es



### ......................
File oefault ~ile Name



### ---- 

## ----------------- visibility Change DELVIS.OAT ! 

### DVISDAT = 

### deciview.dat !



### ...................................................................... rill file names will be converted to lower case 

### if LCFILES = T 

### otherwise, 

### if LCFILES = F, 

### file names 

### will be converted to UPPER CASE T = 

### lower case ! LCFILES = T ! F = UPPER CASE 

### NOTE: (1) 

### file/path names can be up to 132 characters in length NOTE: 

### (2) Filenames for ALL PLOT 

### and TIMESERIES FILES 

### are constructed using a template that includes a pathname, user-supplied character(s), and context-specific strings, where ASPEC = 

### species Name

CONC = CONC Or 

### WFLX Or 

### DFLX Or TFLX tt = 

### Averaging period (e.g. 03) 

### ii = 

### Rank (e.g. 02) hh = 

### ~our(ending) in LST
szzzz = 

### LST time zone shift (EST is -0500)

yyyy = 

### Year(L5~)
mm = 

### Month(LST) dd = 

### day of month (LST) 

### are determined internally based on selections made below. 

### ~f a path or user-supplied character(s1 are supplied, each must contain at least 1 non-blank character.

!END!

INPUT GROUP: 

### 1 -- 

### General run control parameters



## --------------
option to run all periods found in the met. file(s) (METRUN) 

### Default: 0 ! METRUN = 

### 1 !



### METRUN = 

### 0 

## - Run period explicitly defined below METRUN = 

### 1 

## - Run all periods in CALPUFF data file(s)

Startina date: year (ISYR) 

### -- NO 

### default ! ISYR = 

### 1990 ! 

## - Month (ISMO) 

## -- NO 

### default ! ISMO = 

### 0 ! Day (ISDY) 

## -- NO 

### default ! ISDY = 

### 0 ! 

### Starting time: Hour (ISHR) 

### -- No default ! ISHR = 

### 0 !

Number of hours to process (NHRS) 

### -- No default ! NHRS = 

### 0 !



### Process every hour of data?(N~~P) 

### -- Default: 1 ! NREP = 

### 1 ! (1 = 

### every hour processed,

2 = 

### every 2nd hour processed, 5 = 

### every 5th hour processed, etc.)
Species & 

### Concentration/Deposition Information Page 2Species to process (ASPEC) 

## -- NO default ! ASPEC 

### = VISIB ! (ASPEC = VISIB for visibility processing)
~ayer/deposition code (ILAYER) 

## -- ~efault: 

### 1 ! ILAYER = 

### 1 ! '1' for CALPUFF concentrations, 

### '-1' for dry deposition fluxes, 

### '-2' for wet deposition fluxes, '-3' for wetidry deposition fluxes.
scaling factors of the form: -- Defaults: !A= 0 

### ! x(new) = X(o1d) 

# " A 

## + B A = 0.0 !B= 0 ! (NOT applied 

### if A = B = 0.0) B = 0.0

Add Hourly ~ackground Concentrations/~luxes? (LBACK) 

## -- Default: F ! LBACK 

### = F !

Source Information

Option to process source contributions: 0 = Process only total reported contributions 1= Sum all individual source contributions and process 

### 2 = Run in TRACEBACK mode to identify source contributions at a SINGLE receptor (MSOURCE) 

## -- Default: 0 ! MSOURCE = 0 !
Receptor information



## ....................
Gridded receptors processed? (LG) 

## -- Default: F ! LG = F ! Discrete receptors processed? (LD) 

### -- Default: F 

### ! LD = T 

### ! CTsG Complex terrain receptors processed? (LCT) 

## -- Default: F 

### ! LCT 

### = F 

### !
--Report results 

### by DISCRETE receptor RING? (only used when LD = T) (LDRING) 

## -- Default: F 

### ! LDRING = F 

### !

--select range of DISCRETE receptors (only used when LD = T):

select ALL DISCRETE receptors by setting NDRECP flag to 

### -1; OR Select SPECIFIC DISCRETE receptors by entering 

### a flag (0,l) for each 0 = discrete receptor not processed 

### 1 = discrete receptor processed using repeated value notation to select blocks of receptors: 23"1, 15*0, 12"l Flag for all receptors after the last one assigned is set to 0 (NDRECP) 

## -- Default: 

### -1 ! NDRECP = 

### -1 !

--select range of GRIDDED receptors (only used when LG = T):



### x index of LL corner (IBGRID) 

## -- Default: 

### -1 ! IBGRID = 

### -1 ! (-1 OR 

### 1 <= IBGRID 

### <= NX)

Y index of LL corner (JBGRID) -- Default: 

### -1 ! JBGRID = 

### -1 ! (-1 OR 

### 1 C= IBGRID 

### <= NY)



### x index of UR corner (IEGRID) 

## -- Default: -1 

### ! IEGRID = 

### -1 ! 

### (-1OR 1 

### <= IEGRID 

### <= NX) Page 3Y 

### index of UR 

### corner (~EGRID) 

## -- oefault: -1 ! JEGRID = -1 ! 

### (-1 OR 

### 1 <= 

### 3EGRID <= NY)



### Note: Entire grid is processed 

### if 

### IBGRID=3BGRID=IEGRID=lEGRID=-I

--Specific gridded receptors can also be excluded from CALPOST processing by filling 

### a processing grid array with 0s and 1s. 

### ~f the processing flag for rece tor index 

### (i ,j) 

### is 1 (ON), 

### that receptor 

# 7. will be processed 

### if it 

### ies within the range delineated by IBGRIo, JBGRID,IEGRID, IEGRID 

### and 

### if LG=T. 

### ~f it 

### is 

### 0 (OFF), it will 

### not be processed in the run. By default, all array values are set to 1 (ON).



### Number of gridded receptor rows provided in subgroup (la) to identify specific gridded receptors to process (NGONOFF) 

## -- Default: 0 ! 

### NGONOFF = 

### 0 !



### -------------- subgroup (la) 

### -- Specific gridded receptors included/excluded



## -------------- specific gridded receptors are excluded from CALPOST processing by filling a processing grid arra with 0s and 1s. A 

### total of NGONOFF lines are read here. Eacx line corresponds to one 'row' in the sampling grid, starting with the NORTHERNMOST 

### row that contains receptors that you wish to exclude, and finishing with
row 1 to the SOUTH (no intervening rows may be skipped). within
a row, each receptor position is assigned either a 0 or 1, starting with the westernmost receptor. 0 = 

### gridded receptor not processed 1 = 

### processed
Repeated value notation may be used to select blocks of receptors: 23"1, 15"0, 12"1

Because all values are initially set to 1, any receptors north of the first row entered, or east of the last value provided in a row, remain ON.

(NGXRECP) 

## -- Default: 1

INPUT GROUP: 

### 2 

### -- Visibility Parameters (ASPEC = ~~SIB)



### --------------
Identify the Base Time Zone for the CALPUFF simulation (BTZONE) 

## -- NO 

### default ! BTZONE = 5 !



### Particle growth curve ~(RH) 

### for hygroscopic species (MFRH) 

## -- Default: 2 ! 

### MFRH = 

### 2 !

1 = IWAQM 

### (1998) ~(RH) 

### curve (originally used with MVISBK=l) 2 = 

### FLAG (2000) ~(RH) 

### tabulation 3 = EPA 

### (2003) ~(RH) 

### tabulation

Maximum relative humidity (%) 

### used in particle growth curve (RHMAX) -- 

### Default: 98 ! RHMAX = 

### 98 ! 

### Page 4Modeled s~ecies to be included in Include Include Include Include Include Include SULFATE? NITRATE? ORGANIC CARBON? COARSE PARTICLES? FINE PARTICLES? ELEMENTAL CARBON? (LVS04) (LVN03) (LVOC) (LVPMC) (LVPMF) (LVEC) computing the 1



## -- Default: T



## -- Default: T



## -- Default: T



### -- Default: T



## -- Default: T



## -- Default: T ght extinction ! LVS04 = T 

### ! ! LVN03 = T ! 

### ! LVOC = T 

### ! ! LVPMC = F ! 

### ! LVPMF = F ! ! LVEC = F 

### !
And, when ranking for TOP-N, TOP-50, and Exceedance tables,

Include BACKGROUND? (LVBK) 

## -- Default: T 

### ! LVBK = T 

### !
species name used for particulates in MODEL.DAT file COARSE (SPECPMC) 

## -- ~efault: PMC ! SPECPMC = 

### ! FINE (SPECPMF) 

## -- Default: PMF 

### ! SPECPMF = !
Extinction Efficiency 

### (l/~mper ug/m**3)



## ........................................ MODELED particulate species : PM COARSE (EEPMC) 

## -- Default: 0.6 ! EEPMC = 0.6 PM FINE (EEPMF) 

### -- Default: 1.0 

### ! EEPMF = 1.0 BACKGROUND particulate species: PM COARSE (EEPMCBK) 

## -- Default: 0.6 ! EEPMCBK = 0.6 other species: AMMONIUM SULFATE (EESD4) 

## -- Default: 3.0 ! EES04 = 3. AMMONIUM NITRATE (EEN03) 

## -- Default: 3.0 ! EEN03 = 3. ORGANIC CARBON (EEOC) 

## -- Default: 4.0 ! EEOC = 4. SOIL (EESDIL)-- Default: 1.0 ! EESOIL = 1. ELEMENTAL CARBON (EEEC) 

## -- Default: 10. ! EEEC = 10.
Background Extinction Computation



## .................................
Method used for the 24h-average of percent change of light extinction: Hourly ratio of source light extinction 

### / background light extinction is averaged? (LAVER) 

## -- Default: F 

### ! LAVER = F 

### !

Method used for background light extinction (MVISBK) 

### -- Default: 2 

### ! MVISBK = 6 !



### 1 

### = Supply single light extinction and hygroscopic fraction



## - Hourly F(RH) adjustment applied to hygroscopic background and modeled sulfate and nitrate 

### 2 = Compute extinction from speciated PM measurements (A)



## - Hourly F(RH) adjustment applied to observed and modeled sulfate and nitrate



## - F(RH) factor is capped at F(RHMAX) 3 

### = Compute extinction from speciated PM measurements (B)

- Hourly F(RH) adjustment applied to observed and modeled sulfate and nitrate



## - Receptor-hour excluded 

### if RH>RHMAX



## - Receptor-day excluded 

### if fewer than 6 valid receptor-hours 4 

### = Read hourly transmissometer background extinction measurements



## - Hourly F(RH) adjustment applied to modeled sulfate and nitrate

- Hour excluded 

### if measurement invalid (missing, interference,
or large RH)



## - Receptor-hour excluded 

### if RHzRHMAX



## - Receptor-day excluded 

### if fewer than 6 valid receptor-hours 

### 5 

### = Read hourly nephelometer background extinction measurements



## - ~ayleigh extinction value (BEXTRAY) added to measurement

- Hourly F(RH) adjustment applied to modeled sulfate and nitrate



## - Hour excluded 

### if measurement invalid (missing, interference, Page 

### 5

### Appendix~Cohutta2OOlCALPOS~.1~~ or large RH)



## - Receptor-hour excluded 

### if RH>RHMAX



## - Receptor-day excluded 

### if fewer than 6 valid receptor-hours 6 = 

### Compute extinction from speciated PM 

### measurements

- 

### FLAG monthly RH 

### adjustment factor applied to observed and modeled sulfate And nitrate 7 = 

### use observed weather or prognostic weather information for background extinction during weather events; otherwise, use Method 2



## - Hourly F(RH) 

### adjustment applied to modeled sulfate and nitrate


# - F(RH) factor is capped at F(RHMAX)



## - During observed weather events, compute BeXt from visual range 

### if using an observed weather data file, or



## - During prognostic weather events, use Bext from the prognostic weather file



## - use Method 2 for hours without a weather event
Additional inputs used for MVISBK = 

### 1:



## ......................................

Background light extinction (l/Mm) (BEXTBK) 

## -- NO 

### default ! BEXTBK = 

### ! Percentage of particles affected by relative humidity (RHFRAC) 

## -- NO 

### default ! RHFRAC = !


### Additional inputs used for MVISBK = 

### 6:



### ......................................

Extinction coefficients for hygroscopic species (modeled and background) are computed using a monthly RH 

### adjustment factor in place of an hourly RH 

### factor (VISB.DAT 

### file is NOT 

### needed). Enter the 12 monthly factors here (RHFAC). Month 1 is January.

(RHFAC) 

## -- NO 

### default ! RHFAC = 

### 3.3, 3.1, 

### 3, 2.8, 3.4, 3.8, 4, 4.2, 4.2, 3.8, 3.4, 3.5 

### !

Additional inputs used for MVISBK = 7:



### The weather data file (DATSAV abbreviated space-delimited) that is identified as VSRN.DAT may contain data for more than one station. Identify the stations that are needed in the order in which they will be used to obtain valid weather and visual range. The first station that contains valid data for an hour will be used. Enter up to MXWSTA (set in PARAMS 

### file) integer station IDS 

### of up to 6 digits each as variable IDWSTA, 

### and enter the corresponding time zone for each, as variable TZONE (= 

### UTC-LST).

A 

### prognostic weather data file with sext for weather events may be used in place of the observed weather file. Identify this as the VSRN.DAT 

### file and use a station ID 

### of IDWSTA = 999999, 

### and TZONE = 0.



### NOTE: TZONE identifies the time zone used in the dataset. The DATSAV 

### abbreviated space-delimited data usually are prepared with UTC time rather than local time, so TZONE is typically set to zero.
(IDWSTA) -- NO 

### default 

# " IDWSTA = 

### 000000 

# " (TZONE) 

### -- NO 

### default 

# " TZONE = 0. 

# "

ndditional inputs used for MVISBK = 

### 2,3,6,7:



### ........................................

Background extinction coefficients are computed from monthly CONCENTRATIONS of ammonium sulfate (BKS04), ammonium nitrate (BKND31,
coarse particulates (BKPMC) 

## , organic carbon (BKOC) 

## , soi 1 (BKSOIL) 

## , and elemental carbon (BKEC)



## . Month 1 

### is January. (ug/m*"3) Page 6(BKS04) 

## -- No default 

### ! BKS04 = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 

### ! (BKN03) 

### -- NO default ! BKN03 = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ! (BKPMC) 

## -- NO default 

### ! BKPMC = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 

### ! (BKOC) 

## -- NO default 

### ! BKOC = 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 

### ! (BKSOIL) 

## -- NO default 

### ! BUSOIL= 11.38, 11.38, 11.38, 11.38, 11.38, 11.38, 11.38, 11.38, 11.38, 11.38, 11.38, 11.38 

### ! (BKEC) 

## -- NO default !BKEC =0,0,0,0,0,0,0,0,0,0,0,0!

~dditional inputs used for MVISBK = 2,3,5,6,7:



## ..............................................

Extinction due to ~ayleigh scattering is added 

### (I/M~) (BEXTRAY) 

### -- Default: 10.0 ! BEXTRAY 

### = 10 

### !

INPUT GROUP: 3 

### -- OUtPUt options



## --------------

Documentation



## -------------

Documentation records contained in the header of the CALPUFF output file may be written to the list file. Print documentation image? (LDOC) 

### -- Default: F ! LDOC = F Output Units



### ------------ units for 

### ~ll output (IPRTU) 

## -- Default: for for Concentration Deposition 

### 1 = g/mQ"3 g/m""2/s 2 = mg/mQQ3 mg/m:22/s 3 = ug/m"*3 ug/m*_z/s 4 = ng/mQ*3 ng/m--2/s 5 = odour units
Visibility: extinction expressed in 1/Meg 

### 1 ! IPRTU = 3

let :ers (IPRTU is ignored)

Averaging time(s) reported



## ..........................

1-hr averages (L~HR) 

## -- ~efault: T

3-hr averages (L3HR) 

### -- Default: T

24-hr averages (L~~HR) 

## -- Default: T
~un-length averages (LRUNL) 

## -- Default: T ! LlHR = F 

### !

! L3HR = F !



### ! L24HR = F 

### !

! LRUNL = F !
user-specified averaging time in hours - results for
an averaging time of NAVG hours are reported for NAVG greater than 0: (NAVG) 

### -- ~efaul t : 0 

### ! NAVG = 0 

### !

Types of tabulations reported

Page 7Appendix~Cohutta2001c~~~os~.1~~ 

### 1) Visibility: daily visibility tabulations are always reported for the selected rece tors when ASPEC = VISIB. 

### In addition, any of tEe other tabulations listed below may be chosen to characterize the light extinction coefficients. [List file or Plot/Analysis File]

2) Top 50 table for each averaging time seleczed [List file only] (LT50) 

## -- Default: T 

### ! LT50 = F ! 

### 3) Top 'N' 

### table for each averaging time selected [List file or Plot file] (LTOPN) 

## -- Default: F ! LTOPN = F 

### !



## -- Number of 'Top-N' values at each receptor selected (NTOP 

### must be <= 

### 4) (NTOP) 

## -- Default: 4 ! NTOP = 1 

### !



## -- Specific ranks of 'TOP-N' 

### values reported (NTOP 

### values must be entered) (ITOP(4) array) 

## -- oefault: ! ITOP = 

### 0 !

4) Threshold exceedance counts for each receptor and each averaging time selected [List file or Plot file] (LEXCD) 

### -- Default: F ! LEXCD = F 

### !



## -- rdentify the threshold for each averaging time by assigning a non-negative value (output units)



## .

-- Default: -1.0 Threshold for 1-hr averages (THRESH1) ! 

### THRESH1 = 

### -1.0 ! Threshold for 3-hr averages THRESH^) ! 

### THRESH3 = 

### -1.0 ! 

### Threshold for 24-hr averages THRESH^^) 

### ! THRESH24 = 

### -1.0 I 

### Threshold for NAVG-hr averages (THRESHN) ! 

### THRESHN = 

### -1.0 

# i

-- counts for the shortest averaging period selected can be tallied daily, and receptors that experience more than NCOUNT counts over any NDAY 

### period 

### will be reported. This type of exceedance violation output is triggered only 

### if NDAY > 0.



### Accumulation period(~ays) (NDAY) 

### -- Default: 

### 0 ! NDAY = 

### 0 ! Number of exceedances allowed (NCOUNT) 

## -- Default: 1 

### ! NCOUNT = 

### 1 !

5) Selected day table(s)

Echo Option 

## -- Many records are written each averaging period selected and output is grouped by day [List file or Plot file] (LECHO) 

### -- Default: F ! LECHO = F 

### !

Timeseries Option 

## -- Averages at all selected receptors for each selected averaging period are written to timeseries files. Each file contains one averaging period, and all receptors are written to a single record each averaging time. [TSERIESASPEC-tt~~-CONC-TSUNAM.DAT 

### files] Page 8

### Appendix~Cohutta2OOl~~L~OS~.1~~ (LTIME) 

## -- Default: F ! 

### LTIME = F !



### Peak value Option 

## -- Averages at all selected receptors for each selected averaging period are screened and the peak value each period is written to timeseries files. Each file contains one averaging period. [PEAKVAL~ASPEC~~~HRCONC~TSUNAM.DAT 

### files] (LPEAK) 

## -- Default: F ! LPEAK = F 

### !



## -- Days selected for output (IECHO(366)) 

## -- Default: 366"O ! IECHO = 

### 366"O ! (366 values must be entered)

Plot output options

Plot files can be created for the TOP-N, 

### Exceedance, and Echo tables selected above. Two formats for these files are available, DATA 

### and GRID. 

### In the DATA 

### format, results at all receptors are listed along with the receptor location [x, y,vall,val2,.



# . .I



# . In the GRID 

### format, results at only gridded receptors are written, using a compact representation. The gridded values are written in
rows (x varies), starting with the most southern row of the grid. The GRID 

### format is given the .GRD 

### extension, and includes headers compatible with the SURFER(R) 

### plotting software.

A 

### plotting and analysis file can also be created for the daily peak visibility summary output, in DATA 

### format only.

Generate Plot file output in addition to writing tables to List file? (LPLT) 

## -- Default: F ! 

### LPLT = F !



### Use GRID 

### format rather than DATA format, when avai 1 able? (LGRD) 

### -- Default: F ! 

### LGRD = F 

### !

Auxiliary output Files (for subsequent analyses)



## ................................................ visibility

A 

### separate output file may be requested that contains the change in visibility at each selected receptor when ASPEC = VISIB. 

### This file can be processed to construct visibility measures that are not available in CALPOST.
Output file with the visibility change at each receptor? (MDVIS) 

## -- Default: 0 ! MDVIS = 

### 0 !

0 = 

### Do Not create file 

### 1 = 

### Create file of DAILY (24 hour) ~~lta-oeciview 2 = 

### Create f3le of DAILY (24 hour) Extinction change (%) 

### 3 = 

### create file of HOURLY 

### ~elta-~eciview 4 = 

### create file of HOURLY 

### Extinction Change (%)



### ndditional Debug Output



### ....................... Page 9Output selected information to List file for debugging? (LDEBuG) 

## -- 0efault: F ! LDEBUG = F !

output hourly extinction information to REPORT.HRV?

(Vi si bi 1 

### i ty ~ethod 

### 7) (LVEXTHR) 

## -- ~efault: F ! LVEXTHR = F 

### !

Page 10Chemical Products Corporation – AIRS No. 04-13-105-00008 102 Old Mill Road, Cartersville, Georgia, 30120, Bartow County
Visibility Improvement State and Tribal Association of the Southeast Best Available Retrofit Control Technology CALPUFF Visibility Modeling Protocol 

# DRAFT

Appendix E Example of CALPOST List File Ou tput using the Existing and New IMPROVE Algorithm for Calculating Visibility Impairment at Class I Areas within 300km of Chemical Products CorporationCALPOST-IMPROVE Processor for the EastREVISED DRAFT2 August 2006000000000000000000000000000000000000000000006. Enter desired NO2/Nratio (default is 0)Neweigh & Sea Salt"g columnCALPOST Recalculation with New IMPROVE Algorithm---------------------------- INPUT from CALPOST (based on old IMPROVE algorithm) -----------1. At cell A7, import "Ranked Daily Visibility Change" (bext) table, includinheadings, from CALPOST (22 days, max)5. Check calculated valubelow against CALPOST"Ranked Daily VisibilityYEARDAYHR RECEPTORCOORDINATES (TYPEkm)BEXT(Model)BEXT(BKGBEXT(Total)%CHANGEF(RH)bxSO4bxNO3 bxOC bxEC bxPMCbxPMFRankdv(total) dv(bkg)#NUM! #NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!#NUM!3. Enter value of site-specific Rayleigh scattering coefficient, from "Raylworksheet4. (Optional) Insert annual average sea salt concentration, from "Rayleigh & Sea Salt"worksheet. Leave blank if not used, i.e. default is 0.---------------------------------- OUTPUT (based on new IMPROVE algorithm) ---------------------------YEARDAYHRRECEPTORCOORDINATES (TYPEkm)BEXT(SourceBEXT(BKGBEXT(Total)%CHANGERH(%)bsSO4bsNO3bsOCbsECbsPMCbsPMFbaNO2dv(total)dv(bkg)Rank00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A00000#N/A#N/A#N/A#N/A#N/A#N/A#N/A000001#N/A#N/A