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 Corporation 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 Source Potential-to-Emit Calculations, Maximum 24-

hr Actual Emissions Rates, and 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 Areas within 300km of Chemical Products Corporation 

 
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 Page 2 

DRAFT

1. Introduction 

Chemical Products Corporation (“CPC”) operates an inorganic chemical manufacturing plant 

in Cartersville, Georgia and primarily processes barite ore (barium sulfate) and celestite ore 

(strontium sulfate) into finished products such 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 reduction kilns and various types of 

dryers. CPC also owns a furnace that has been used 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 source with respect to the Part 70 (Title V) 

operating permit program and operates under Georgia 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 address 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 one 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-eligible sources” if is it determined that the 

source emits any air pollutant which may reasonably 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 to determine 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 conditions 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 processing facility, were in existence on August 7, 

1977, and began operation after August 7, 1962 and have the potential-to emit greater than 

250 tpy of at least one visibility-impairing pollutant. 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 contribute to visibility impairment in a Class I Area. 

Periodic monitoring, recordkeeping and reporting requirements required to demonstrate 
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 Page 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 modeling 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 finer 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 initial 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 the 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 contributes 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.doc Page 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 cause 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 tiny pollution particles in the air. Some 

light is absorbed by particles and other light 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, particularly during humid conditions. The 

predicted change in haze is measured by calculating natural light extinction (bext) 

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

is calculated using the following equation, 

, 10ln10 ⎟

⎟

⎠

⎞

⎜

⎜

⎝

⎛ = bext 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

Organic PM

Ammonium Nitrates Condesable

34

Ammonium Sulfates 

ext 4 2 4 OC4NONHRH3fSONHRH3fb bEC10PMC0.6Soil

2.5 2.510

= + + + + ++

−

14 244 4 344 14 244 4 344 321 142 32143

The concentrations, in square brackets, are in micrograms per cubic meter 

(ug/m3

) and bext is in units of inverse megameters (Mm-1). The Rayleigh 

scattering term (bray) has a default value of 10 Mm-1, as recommended in EPA 

guidance for tracking reasonable progress. 

2. Chemical Products Corporation BART-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 units that comprise the BART-eligible source at 
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 Page 5 

DRAFT

Chemical Products Corporation. Potential emissions for NOx, SO2 and PM10 are 

determined 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 parameters 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 the 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.424 
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 Page 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 modeling domains are strategically designed to 

cover all potential BART eligible sources within VISTAS states and all Class I Areas 

within 300 km of those sources. The extents of the 4-km sub-regional domains are 

shown in Figure 3-1. 

Figure 3-1 VISTAS 4-km Sub-Regional Domains
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 Page 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) files were used by VISTAS to generate the 

land use data at 4-km resolution for input to the 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 chemical 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-urban monitors within and just outside the 

computational grid shall be used as input 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 output per 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-1 for 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 10eb bSoil Mm 10Soil

2.5

− = +=+=

[ ] [ ] ( ) 10 1 HI

3 back 10eSoils Mm 01

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.doc Page 8 

DRAFT

Table 3-1 Natural Background Conditions for Class I Areas 

Class I Area bext

(Mm-1) 

Annual 

Average 

Haze Index 

(dv) 

20% Best 

Days Haze 

Index 

(dv) 

20% Worst 

Days Haze 

Index 

(dv) 

Soil bback

Annual Avg. 

(μg/m3

) 

Soil bback

20% Best 

Days Avg. 

(μg/m3

) 

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 amplify the visibility impairment of hygroscopic 

sulfates and nitrates were characterized using “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 CALPUFF as “the best modeling application 

available for predicting a single source’s contribution to visibility impairment.” As a 
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 Page 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 recirculation 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 Summary Report and Recommendations for 

Modeling Long Range Transport Impacts (EPA-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 selection 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 shall 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 procedures 
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 Page 10 

DRAFT

proposed to be used determining the change 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 be used 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 background 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 default 
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 Page 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 consultant, Dr. Ivar Tombach. The new IMPROVE 

algorithm takes into account site-specific Rayleigh 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 TBD 
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

DRAFT

Appendix A 

Chemical Products Corporation BART-Eligible Source Potential-to-Emit 

Calculations, Maximum 24-hr Actual Emissions Rates, and Modeled 

Emissions Inventory 
Chemical Products Corporation - Best Available Retrofit Technology Eligible Emissions Units (BART-Eligible Source)

Table A-1. Chemical Products Corporation BART-Eligible Source Potential-to-Emit Coke

Stack ID 

No. a 

Emissions 

Unit ID No. a Emission Unit Description a Construction Date a

Maximum 

Rated Heat 

Input, 

(mmBtu/hr) a

Fuel Oil 

Consumption, 

(gph) a

Coke 

Consumption, 

(tph) a

Design 

Capacity, (tph) a

LCC

Easting, 

(km) b

LCC

Northing, 

(km) b

Stack

Height, 

(ft) c

Exhaust

Temp., t 

(

oF) c

Exit

Velocity, (fps) c

Stack

Diameter, 

(ft) c

Georgia Rule (d)

PM10 Emissions

Permit Limit, (lb/hr)

Georgia Rule (e)

PM10 Emissions

Permit Limit, (lb/hr)

Georgia Rule (g)

Fuel Sulfur

Permit Limit, (%)

Coke Sulfur 

Content, (%) d

SO2 Emissions

Factor, Mass 

Balance, 

(lb/ton) e

NOx Emissions

Factor, AP-42 1.4-1/-2, 

(lb/mmBtu)

SO2 Emiss

Factor, AP-42 

(lb/mmB

Chemical 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.0005BDZ1 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.0005BDZ3 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.0005BDZ4 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.0005BKZ1 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.0005GBZ1 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.0005

Notes: 

a. Stack and emissions unit identification, construction date, maximum rated heat input, fuel oil consumption rates, petroleum coke consumption and design capacity determined from Title V Permit No. 2819-015-0008-V-02-0/-1 and Title V Application No. 16069. 

b. 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 performed by VISTAS. d. Highest petroleum coke sulfur content based on historical data for calendar years 2001, 2002, and 2003 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 natural gas, propane, or fuel oil. SO2 emissions rate for the East Rotary Kiln includes estimated SO2 from processing coke. 

Table A-2. Chemical Products Corporation BART-Eligible Source Maximum 24-hr Actual Emissions Rates (Fuel Oil Combustion, 0.5% Sulfur Content, AP-42) Coke

Stack ID 

No. a 

Emissions 

Unit ID No. a Emission Unit Description a Construction Date a

Maximum 

Rated Heat 

Input, 

(mmBtu/hr) a

Fuel Oil 

Consumption, 

(gph) a

Coke 

Consumption, 

(tph) a

Design 

Capacity, (tph) a

LCC

Easting, 

(km) b

LCC

Northing, 

(km) b

Stack

Height, 

(ft) c

Exhaust

Temp., t 

(

oF) c

Exit

Velocity, (fps) c

Stack

Diameter, 

(ft) c

Georgia Rule (d)

PM10 Emissions

Permit Limit, (lb/hr)

Georgia Rule (e)

PM10 Emissions

Permit Limit, (lb/hr)

Maximum Actual

Fuel Sulfur

Permit Limit, (%)

Maximum 

Actual Coke 

Sulfur Content, 

(%) d

SO2 Emissions

Factor, Mass 

Balance, 

(lb/ton) e

NOx Emissions

Factor, AP-42 1.4-1/-2, 

(lb/mmBtu)

SO2 Emiss

Factor, AP-42 

(lb/mmB

Chemical 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 

BDZ1 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 

BDZ3 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 

BDZ4 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 

BKZ1 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.1500 0.000 0.000 

GBZ1 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 

Notes: 

a. Stack and emissions unit identification, construction date, maximum rated heat input, fuel oil consumption rates, petroleum coke consumption and design capacity determined from Title V Permit No. 2819-015-0008-V-02-0/-1 and Title V Application No. 16069. 

b. 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 performed by 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 Specication

Chemical Products Corporation, Title V Permit No. 2819-015-0008-V-02-0/-1 

BCZ1 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.000 0 100 BDZ1 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 BDZ3 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 BDZ4 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 BKZ1 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 GBZ1 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 

Stack ID 

No.

Emissions 

Unit ID No. Emission Unit Description

Base 

Elevation, (ft) H2SO4, (lb/hr) NOx, (lb/hr) PM10, (lb/hr)

PM2.5, (lb/hr) NH3, (lb/hr)

PM10 Particle Speciation (Unknow

Filterable, (%) Condesable

Exhaust

Temp., t 

(

oF)

Exit

Velocity, 

(fps)

Stack

Diameter, (ft) SO2, (lb/hr)

Location

Latitude Longitude Zone

Location LCC

Easting, (km) Northing, (km)

Stack

Height, (ft)

Natural Gas Natural Gas 
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

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 1 
Appendix~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 4 
INPUT 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 a1 lowed? (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 5 
Aqueous 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 

CGRUP name wi 7 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 12 
INPUT GROUP: 5 -- Output Options -------------- 4 

FILE DEFAULT VALUE ---- ------------- 

Concentrations (ICON) 

Dry Fluxes (IDRY) 

wet luxes (IWET) 

20 Temperature (ITZD) 

2D Density (IRHO) 

Re1 ative ~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 14 
Met. 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 17 
NOX = 

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 22 
AppendixBcohutta2001~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 23 
INPUT 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 = 1 b/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 = 1 b/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 1di ng 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 

INP 
Data 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 36 
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

DRAFT

Appendix C 

Example of Proposed POSTUTIL Input 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~. 

wh: it is the inout source of Ammonia? 

(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 8 
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

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 on1 y 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~~~HRCDAT 

Page 1 
visibility plot DAILY-VISIB-VUNAM.DAT " VUNAM = % 

(~ai 1 yPeak 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 2 
Species 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 informati on .................... 

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 ! 

(-1 OR 1 <= IEGRID <= NX) 

Page 3 
Y 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 will be processed if it 7. 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 4 
Modeled 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/~m per 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 igno red) 

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 7 
Appendix~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 (%) 

nddi tional Debug Output ....................... 

Page 9 
Output 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 10 
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

DRAFT

Appendix E 

Example of CALPOST List File Output using the Existing and New 

IMPROVE Algorithm for Calculating Visibility Impairment at Class I 

Areas within 300km of Chemical Products Corporation
CALPOST-IMPROVE Processor for the East 

0 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 0

6. Enter desirratio (default 

New 

eigh & Sea Salt" 

g column 

CALPOST 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 calcbelow agains"Ranked Dail

YEAR DAY HR RECEPTOR COORDINATES (km) TYPE BEXT(Model) BEXT(BKG BEXT(Total)%CHANGE F(RH) bxSO4 bxNO3 bxOC bxEC bxPMC bxPMF Rank

dv(total) d

#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 "Raylworksheet

4. (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) ---------------------------

YEAR DAY HR RECEPTOR COORDINATES (km) TYPE BEXT(Source BEXT(BKG BEXT(Total)%CHANGE RH(%) bsSO4 bsNO3 bsOC bsEC bsPMC bsPMF baNO2Rank dv(total) d

0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A 0 0 0 0 0 #N/A #N/A #N/A #N/A #N/A #N/A #N/A 0 0 0 0 0 1 #N/A