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