Reid Gardner Generating Station Modeling files for EPA SIP action

Naming convention for input and output files

The name "rgset2" for "RGGS set #2" is associated with the

main CALPUFF simulations. Sensitivity simulations did not follow

this convention; see below.

--------------------

CALPUFF and POSTUTIL

•CC_cNN_YYYY.EXT

•CC = cf for CALPUFF, and pu for POSTUTIL

•cNN = control scenario, per table below

•YYYY = meteorology year 2001, 2002, or 2003

•.EXT = .inp for control file input,

••.lst for list file output,

••.con for concentration file output

Example: pu_c15_2001.inp

Since there are 8 controls scenarios (omitting c14) and 3 years,

there are 8*4 = 24 CALPUFF and POSTUTIL runs and corresponding

.inp, .lst, and .con files.

--------------------

CALPOST

•ct_cNN_vmNX_AREA_YYYY.EXT

•ct for "CALPOST"

•vmNX = visibility method (MVISBK) and background, see table

•AREA = Class I area abbreviation, see table

•.EXT = .inp or .lst only

Example: ct_c15_vm8b_grca_2001.inp

Since there are 4 visibility methods and 5 areas, there are

4*5 = 20 CALPOST runs for each of the 24 scenario/year

combinations, or a total of 480 CALPOST .inp and .lst files.

--------------------

Extracted visibility impacts

Visibility impacts were extracted from the CASPOST.LST files,

and combined into comma-separated value files. The "p98"

and "table" files were imported into Excel to create some

of the impact tables in the TSD and the accompanying

spreadsheet "RGGS_TSD_CALPUFF_tables.xls"

•METRIC_SUBSET_rgset2_ct6.221.csv

•METRIC = bext for extinction, or ddv for delta deciviews

•subset = empty for all daily values,

••p98 for 98th percentile for all years together,

••p98yrs for 98th percentile for individual years,

••table for 98th percentile and/or exceedances

•••of 5% and 10% bext or 0.5 and 1.0 dv thresholds

•rgset2 = RGGS set #2, the full set of CALPUFF runs

•ct6.221 = CALPOST version 6.221

•.csv = extension indicating comma-separated values

File names:

•bext_rgset2_ct6.221.csv

•bext_p98_rgset2_ct6.221.csv

•bext_p98yrs_rgset2_ct6.221.csv

•bext_table_rgset2_ct6.221.csv

•ddv_rgset2_ct6.221.csv
•ddv_p98_rgset2_ct6.221.csv

•ddv_p98yrs_rgset2_ct6.221.csv

•ddv_table_rgset2_ct6.221.csv

(.CSV files without the "_ct6.221" used a different version

of CALPOST. CALPOST 5.6394 gave identical results to

CALPOST 6.221 for visibility methods vm6a and vm6b.)

-------------------------------------------------------

Control scenarios

cNN•scenario name•Corresponding NDEP NOx case

c00•WRAP •WRAP baseline

c02•WRAP NOx •NA

c04•Baseline NOx•"Baseline"

c11•Enh. LNB+OFA•1 or "LNB+OFA"

c12•SNCR •2

c13•ROFA+Rotamix•3

c14•SCR+ROFA •4

c15•SCR •5

c16•SCR 0.06 lb/MMBtu•NA

All scenarios except "WRAP baseline" use NDEP BART determination emission rate for SO2 

and PM emissions.

SCR scenario c14 was not modeled since it has the same emissions as c15.

-------------------------------------------------------

Visibility Methods Used

vmNX•MVISBK•natural background concentrations

vm6a•6•annual average

vm6b•6•best 20% of days

vm8a•8•annual average

vm8b•8•best 20% of days

-------------------------------------------------------

Class I Areas within 300 km of Reid Gardner Generating Station

AREA•Name •dist•nr •ST•Agency

brca•Bryce Canyon • 226•213•UT•NPS

grca•Grand Canyon • 85•791•AZ•NPS

jotr•Joshua Tree • 292•815•CA•NPS

syca•Sycamore Canyon• 288•316•AZ•USFS

zion•Zion • 148•224•UT•NPS

=======================================================

Sensitivity simulations

There were three sets of sensitivity simulations: Multiple units,

WRAP comparison, and Gridded WRAP comparison

--------------------

Multiple units

These simulations were to assess the affect of modeling the

three RGGS stacks separately as opposed to together.

The naming convention was the same as for the main simulations

above, except that "_u1", "_u2", or "_u3", was inserted after

the control scenario abbreviation. Only control scenario c00,

visibility method 6a and year 2001, was run, so there are

1*1*3 = 3 CALPUFF and POSUTIL runs for the three units, and

for each of these 1*5 CALPOST runs for a total of 45.

Examples: cf_c00_u2_2001.inp, ct_c00_u2_vm6a_jotr_2001.inp

Only a single extraction file was created: ddv_p98_c00_units.csv
The results from this are in a table in the TSD.

--------------------

WRAP comparison

These simulations are a spectrum between the rgset2 main

simulations above, and the WRAP's, a sequence to check the

effect of gradual changes to the inputs. Only control scenario

c00, visibility method 6a, and year 2001 was run. Only the

98th percentile delta deciviews was used in comparing results.

The naming convention is similar to the main one, but the

initial "CC" indicator was usually changed to indicate the

version of the executable, and subdirectory names were used

to remove ambiguity and ensure that various outputs remained

separate (e.g. CALPOST TUNAM outputs, which have limited

scope for variation, like rank(8)_no3_24hr_conc_XXXXXXXX.grd).

"EPA input" means the file was edited from the rgset2 main

simulation input files above.

"WRAP input" means the file was edited from the WRAP's own

input files, downloaded the WRAP Regional Modeling Center web

site at http://pah.cert.ucr.edu/aqm/308/bart.shtml

File names (only .inp and only Grand Canyon shown)

EPA input, but with WRAP emissions and PM speciation

•rgset2_wrapem\cf_cwe_2001.inp

•rgset2_wrapem\ct_cwe_brca_2001.inp

•rgset2_wrapem\pu_cwe_2001.inp

•rgset2_wrapem\ ddv_p98_wrap0_cwe.csv

WRAP input, but CALPUFF 5.8 and CALPOST 5.6394

•WRAP0_cf5.8_ct5.6\cf58_WRAP0_2001.inp

•WRAP0_cf5.8_ct5.6\pu15_cf58_WRAP00_2001.inp

•WRAP0_cf5.8_ct5.6\ct56_cf58_wrap000_grca_2001.inp

•WRAP0_cf5.8_ct5.6\ddv_p98_wrap000_cf58_ct56.csv

WRAP input, but CALPUFF 5.8 and CALPOST 6.131

•WRAP0_cf5.8_ct6.1\cf58_WRAP0_2001.inp

•WRAP0_cf5.8_ct6.1\pu15_cf58_WRAP00_2001.inp

•WRAP0_cf5.8_ct6.1\ct61_cf58_WRAP000_grca_2001.inp

•WRAP0_cf5.8_ct6.1\ddv_p98_wrap000_cf58_ct61.csv

WRAP input, but CALPUFF 6.112 and CALPOST 6.131

•WRAP0_cf6.1_ct6.1\cf61_WRAP0_2001.inp

•WRAP0_cf6.1_ct6.1\pu15_WRAP00_2001.inp

•WRAP0_cf6.1_ct6.1\ct61_WRAP000_grca_2001.inp

•WRAP0_cf6.1_ct6.1\ddv_p98_wrap000_cf61_cf61.csv

(The "WRAP0", "WRAP00" and "WRAP000" names reflect the original

idea, which was to integers incremented for each change

in CALPUFF, POSTUTIL, and CALPOST, respectively. This did

not work in practice because of the initial "cf", "pu", or "ct"

assumed in certain scripts used to call the model executables.)

--------------------

Gridded WRAP comparison

These simulations were the same as for the WRAP comparison,

except they considered gridded receptors over the whole

computational domain instead of just discrete receptors at

Class I areas. The naming convention is also similar,

except that "g" is inserted after the first zero, or at the

end, e.g. "wrap0g00" or "cweg", and in place of visibility
method and area, gridXXX was used, where XXX = "NO3", "SO4",

or "visib: for visibility (extinction).

Example for EPA input, but with WRAP emissions and PM speciation

(CAPUFF, POSUTIL, .lst and .con not shown)

•rgset2_wrapem\ct_cweg_gridNO3_2001.inp

•rgset2_wrapem\ct_cweg_gridSO4_2001.inp

•rgset2_wrapem\ct_cweg_gridvis_2001.inp

Examples for corresponding TUNAM and VUNAM ranked output

(These were used in creating isopleth plots)

•rank(0)_no3_8760hr_conc_w000gd01.grd

•rank(0)_so4_8760hr_conc_w000gd01.grd

•rank(0)_visib_8760hr_conc_w000gd01.grd