bt-star_main_03.21.11_log.txt

Author: Abt Associates Inc. Created: 2/5/2010 Updated: 3/22/2011 Purpose: This is a log file from the SAS program which calculates non-use benefits using a benefits transfer approach based on the Johnston et al. (2009)

1 The SAS System 12:32 Tuesday, March 22, 2011

NOTE: Copyright (c) 2002-2008 by SAS Institute Inc., Cary, NC, USA. NOTE: SAS (r) Proprietary Software 9.2 (TS2M0) Licensed to ABT ASSOCIATES INC, Site 70068845. NOTE: This session is executing on the XP_PRO platform.

NOTE: SAS initialization used: real time 0.43 seconds cpu time 0.23 seconds

1 *Author: Ryan Stapler, Abt Associates Inc. 2 *Created: 2/5/2010; 3 *Updated: 3/22/2011; 4 *Purpose: This is a log file from the SAS program which calculates non-use benefits using a 4 ! benefits transfer approach based on the Johnston et al. (2009)
5 6 7 *As of 2/09/10, seperate compliance schedules are implementd for facilities assigned IM and 7 ! CT based on "ie_year" 8 *As of 2/09/10, it includes a household weighting structure using proportion of total 8 ! habitat acres assigned to a state. 9 *Benefits are reported in $1,000. 10 11 *The location of the input and output files is here; 12 %let d=h; *(C or H drive); 13 %let mainloc=&d.:\ERD\TYCHE Contract\Assignments\WA 2-15 (316b Benefits & SP Study)\Task 5 13 ! - Benefits Analytic Support\IE_032111\BT-STAR\; 14 %let inloc=&mainloc.Inputs\; 15 %let outloc=&mainloc.Outputs\; 16 17 libname outloc "&outloc"; NOTE: Libref OUTLOC was successfully assigned as follows: Engine: V9 Physical Name: h:\ERD\TYCHE Contract\Assignments\WA 2-15 (316b Benefits & SP Study)\Task 5 - Benefits Analytic Support\IE_032111\BT-STAR\Outputs 18 19 *input assumptions; 20 %let t0=2012; *(define time zero); 21 %let yrsc=50; *(define the annualization period); 22 %let yrba=2012; *(define the base year for discounting); 23 %let ie_year=2017; *(define the first year of I&E for towers compliance, all year prior 23 ! just impingement only); 24 %let af_3=(0.03*(1.03**50))/((1.03**50)-1); *(define the 3% Page 1bt-star_main_03.21.11_log.txt

annualization factor by 24 ! changing the number compliance years); 25 %let af_7=(0.07*(1.07**50))/((1.07**50)-1); *(define the 7% annualization factor by 25 ! changing the number compliance years); 26 27 *define species that is being analyzed in the north atlantic and mid-atlantic regions; 28 %let species_1="winter.flounder"; 29 *species_#a is the name each species as it will be presented in final tables; 30 %let species_1a="Winter Flounder"; 31 32 *The name of the input file is; 2 The SAS System 12:32 Tuesday, March 22, 2011

33 %let input=BT_STAR_Inputs_03.21.11.xls; 34 %let input2=IE_outputs_22MAR11.xls; 35 *The name of the output files are; 36 %let output=BT_BenOutputs_&sysdate..xls; 37 %let tables=BT_BenTables_&sysdate..xls; 38 %let SumTables=BT_BenTablesSum_&sysdate..xls; 39 %let printout=&outloc.BT Ben Tables &sysdate..rtf; 40 %let printout3=&outloc.BT Ben Tables &sysdate..csv; 41 42 *This is a macro that can be used to import data from excel; 43 %macro import(infile,insheet,out); 44 PROC IMPORT OUT= WORK.&out 45 DATAFILE= "&infile" 46 DBMS=EXCEL REPLACE; 47 SHEET="&insheet.$"; 48 GETNAMES=YES; 49 MIXED=NO; 50 SCANTEXT=YES; 51 USEDATE=YES; 52 SCANTIME=YES; 53 RUN; 54 %mend; 55 *This is a macro that can be used to export data to excel; 56 %macro export(datafile,outfile,sheet); 57 PROC EXPORT DATA= &datafile 58 OUTFILE= &outfile 59 DBMS=EXCEL replace; 60 SHEET=&sheet; 61 RUN; 62 %mend; 63 64 *import data from excel; 65 66 *import facility data; 67 %import(&inloc.&input,facility_data,facility_data);

NOTE: WORK.FACILITY_DATA data set was successfully created. NOTE: PROCEDURE IMPORT used (Total process time): real time 0.46 seconds cpu time 0.14 seconds

68 *import household data; 69 %import(&inloc.&input,HH_data_reg,HH_1);

NOTE: WORK.HH_1 data set was successfully created. Page 2bt-star_main_03.21.11_log.txt

NOTE: PROCEDURE IMPORT used (Total process time): real time 0.48 seconds cpu time 0.06 seconds

70 *import population change data; 71 %import(&inloc.&input,abun_change,abun_change);

3 The SAS System 12:32 Tuesday, March 22, 2011

NOTE: WORK.ABUN_CHANGE data set was successfully created. NOTE: PROCEDURE IMPORT used (Total process time): real time 1.46 seconds cpu time 0.11 seconds

72 *import I&E data; 73 %import(&inloc.&input2,BNB_IE,ie);

NOTE: WORK.IE data set was successfully created. NOTE: PROCEDURE IMPORT used (Total process time): real time 5.75 seconds cpu time 1.29 seconds

74 *import impingement data; 75 %import(&inloc.&input2,BNB_IMP,i_only);

NOTE: WORK.I_ONLY data set was successfully created. NOTE: PROCEDURE IMPORT used (Total process time): real time 3.28 seconds cpu time 1.34 seconds

76 *import entrainment data; 77 %import(&inloc.&input2,BNB_ENT,ent_only);

NOTE: WORK.ENT_ONLY data set was successfully created. NOTE: PROCEDURE IMPORT used (Total process time): real time 3.25 seconds cpu time 1.12 seconds

78 *import the worksheet with the ccr.fraction data; 79 %import(&inloc.&input2,REDUCTIONS_DETAIL,reductions_detail);

NOTE: WORK.REDUCTIONS_DETAIL data set was successfully created. NOTE: PROCEDURE IMPORT used (Total process time): real time 2.32 seconds cpu time 0.29 seconds

80

************************************************************************** 81 *this section reorganizes the data in order to calculate "impingment fraction" I&E; 82 DATA reductions_detail; 83 SET reductions_detail; 84 if region="" then delete; 85 DROP regional_reduction_imp regional_reduction_ent ccr_fraction_ent 85 ! benefit_value_regional_imp; 86 RUN; Page 3bt-star_main_03.21.11_log.txt

WARNING: The variable regional_reduction_imp in the DROP, KEEP, or RENAME list has never been

referenced. WARNING: The variable regional_reduction_ent in the DROP, KEEP, or RENAME list has never been 4 The SAS System 12:32 Tuesday, March 22, 2011

referenced. WARNING: The variable ccr_fraction_ent in the DROP, KEEP, or RENAME list has never been referenced. WARNING: The variable benefit_value_regional_imp in the DROP, KEEP, or RENAME list has never been

referenced. NOTE: There were 28 observations read from the data set WORK.REDUCTIONS_DETAIL. NOTE: The data set WORK.REDUCTIONS_DETAIL has 28 observations and 3 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

87 PROC SORT data=reductions_detail; 88 By option_name region; 89 RUN;

NOTE: There were 28 observations read from the data set WORK.REDUCTIONS_DETAIL. NOTE: The data set WORK.REDUCTIONS_DETAIL has 28 observations and 3 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

90 91 Data i_only_1; 92 set i_only; 93 DO imp_benefit_value_regional=benefit_value_regional; 94 END; 95 IF species_econ~="zz.all.species" then delete; 96 IF metric~="ae1" then delete; 97 if comrec_type~="all.types" then delete; 98 IF region="zz.allregions" then delete; 99 DROP lossmode baseline_value baseline_value_regional benefit_value benefit_value_regional 99 ! regional_reduction ccr_fraction; 100 RUN;

WARNING: The variable baseline_value in the DROP, KEEP, or RENAME list has never been referenced. WARNING: The variable benefit_value in the DROP, KEEP, or RENAME list has never been referenced. WARNING: The variable regional_reduction in the DROP, KEEP, or RENAME list has never been referenced. WARNING: The variable ccr_fraction in the DROP, KEEP, or RENAME list has never been referenced. NOTE: There were 23712 observations read from the data set WORK.I_ONLY. NOTE: The data set WORK.I_ONLY_1 has 28 observations and 7 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

101 PROC SORT data=i_only_1; Page 4bt-star_main_03.21.11_log.txt

102 By option_name region value_basis species_econ metric comrec_type; 103 RUN;

NOTE: There were 28 observations read from the data set WORK.I_ONLY_1. NOTE: The data set WORK.I_ONLY_1 has 28 observations and 7 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds 5 The SAS System 12:32 Tuesday, March 22, 2011

cpu time 0.00 seconds

104 Data ent_only_1; 105 set ent_only; 106 DO ent_benefit_value_regional=benefit_value_regional; 107 END; 108 IF species_econ~="zz.all.species" then delete; 109 IF metric~="ae1" then delete; 110 if comrec_type~="all.types" then delete; 111 IF region="zz.allregions" then delete; 112 DROP lossmode baseline_value baseline_value_regional benefit_value benefit_value_regional 112 ! regional_reduction ccr_fraction; 113 RUN;

WARNING: The variable baseline_value in the DROP, KEEP, or RENAME list has never been referenced. WARNING: The variable benefit_value in the DROP, KEEP, or RENAME list has never been referenced. WARNING: The variable regional_reduction in the DROP, KEEP, or RENAME list has never been referenced. WARNING: The variable ccr_fraction in the DROP, KEEP, or RENAME list has never been referenced. NOTE: There were 19224 observations read from the data set WORK.ENT_ONLY. NOTE: The data set WORK.ENT_ONLY_1 has 28 observations and 7 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

114 PROC SORT data=ent_only_1; 115 by option_name region value_basis species_econ metric comrec_type; 116 RUN;

NOTE: There were 28 observations read from the data set WORK.ENT_ONLY_1. NOTE: The data set WORK.ENT_ONLY_1 has 28 observations and 7 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

117 Data ie_rd1; 118 set ie; 119 DO ie_benefit_value_regional=benefit_value_regional; 120 END; 121 IF species_econ~="zz.all.species" then delete; 122 IF metric~="ae1" then delete; 123 if comrec_type~="all.types" then delete; 124 IF region="zz.allregions" then delete; 125 DROP regional_reduction ccr_fraction; 126 if region="zz.allregions" then delete; 127 RUN; Page 5bt-star_main_03.21.11_log.txt

WARNING: The variable regional_reduction in the DROP, KEEP, or RENAME list has never been referenced. WARNING: The variable ccr_fraction in the DROP, KEEP, or RENAME list has never been referenced. NOTE: There were 24640 observations read from the data set WORK.IE. NOTE: The data set WORK.IE_RD1 has 28 observations and 10 variables. NOTE: DATA statement used (Total process time): 6 The SAS System 12:32 Tuesday, March 22, 2011

real time 0.01 seconds cpu time 0.00 seconds

128 PROC SORT data=ie_rd1; 129 by option_name region value_basis species_econ metric comrec_type; 130 RUN;

NOTE: There were 28 observations read from the data set WORK.IE_RD1. NOTE: The data set WORK.IE_RD1 has 28 observations and 10 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

131 DATA ie_rd2; 132 merge ie_rd1 i_only_1 ent_only_1; 133 by option_name region value_basis species_econ metric comrec_type; 134 RUN;

NOTE: There were 28 observations read from the data set WORK.IE_RD1. NOTE: There were 28 observations read from the data set WORK.I_ONLY_1. NOTE: There were 28 observations read from the data set WORK.ENT_ONLY_1. NOTE: The data set WORK.IE_RD2 has 28 observations and 12 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

135 Data ie_rd3; 136 *length impingement_fraction $50.0; 137 merge ie_rd2 reductions_detail; 138 by option_name region; 139 IF ent_benefit_value_regional="." then do ent_benefit_value_regional=0; 140 END; 141 IF imp_benefit_value_regional="." then do imp_benefit_value_regional=0; 142 END; 143 imp_frac_imp=(1-ccr_fraction_imp); 144 IF ie_benefit_value_regional~=0 then do 144 ! impingement_fraction=((imp_frac_imp*imp_benefit_value_regional)/ie_benefit_value_reg ional); 145 END; 146 IF ie_benefit_value_regional=0 then do impingement_fraction=""; 147 END; 148 DROP ie_benefit_value_regional imp_benefit_value_regional ent_benefit_value_regional 148 ! ccr_fraction_imp imp_frac_imp 149 value_basis species_econ lossmode metric comrec_type baseline_value 149 ! baseline_value_regional benefit_value benefit_value_regional; Page 6bt-star_main_03.21.11_log.txt

150 RUN;

NOTE: Character values have been converted to numeric values at the places given by: (Line):(Column). 139:32 141:32 146:62 WARNING: The variable baseline_value in the DROP, KEEP, or RENAME list has never been referenced. WARNING: The variable benefit_value in the DROP, KEEP, or RENAME list has never been referenced. NOTE: There were 28 observations read from the data set WORK.IE_RD2. 7 The SAS System 12:32 Tuesday, March 22, 2011

NOTE: There were 28 observations read from the data set WORK.REDUCTIONS_DETAIL. NOTE: The data set WORK.IE_RD3 has 28 observations and 3 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

151 Proc Sort data=ie out=ie; 152 by option_name region; 153 RUN;

NOTE: There were 24640 observations read from the data set WORK.IE. NOTE: The data set WORK.IE has 24640 observations and 9 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.04 seconds cpu time 0.03 seconds

154 PROC Sort data=ie_rd3 out=ie_rd3; 155 by option_name region; 156 RUN;

NOTE: There were 28 observations read from the data set WORK.IE_RD3. NOTE: The data set WORK.IE_RD3 has 28 observations and 3 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

157 Data ie; 158 merge ie ie_rd3; 159 by option_name region; 160 RUN;

NOTE: There were 24640 observations read from the data set WORK.IE. NOTE: There were 28 observations read from the data set WORK.IE_RD3. NOTE: The data set WORK.IE has 24640 observations and 10 variables. NOTE: DATA statement used (Total process time): real time 0.03 seconds cpu time 0.01 seconds

161 ***************************************************************** 162 *edit region names; 163 DATA facility_data; 164 Set facility_data; 165 IF Region="CN" or Region ="CS" or Region ="HI" then reg="CAL"; 166 IF Region="INL" then reg="INL"; 167 IF Region="GL" then reg="GL"; 168 IF Region="MA" then reg="MA"; Page 7bt-star_main_03.21.11_log.txt

169 IF Region="NA" then reg="NA"; 170 IF Region="SA" then reg="SA"; 171 IF Region="Mex" then reg="Mex"; 172 IF Region="VI" then reg="VI"; 8 The SAS System 12:32 Tuesday, March 22, 2011

173 IF Region="PNW" then reg="PNW"; 174 RUN;

NOTE: There were 892 observations read from the data set WORK.FACILITY_DATA. NOTE: The data set WORK.FACILITY_DATA has 892 observations and 6 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

175 176 *Under the current model framework, all species have populations which cover the North 176 ! Atlantic and Mid-Atlantic regions; 177 178 PROC Sort data=abun_change out=abun_change; 179 by reg; 180 RUN;

NOTE: There were 3 observations read from the data set WORK.ABUN_CHANGE. NOTE: The data set WORK.ABUN_CHANGE has 3 observations and 13 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

181 182

************************************************************************************

* 183

************************************************************************************

* 184 generate regional annual WTP values for mid-atlantic and north atlantic regions based on 184 ! estimated 185 increases in abundance; 186 *reports total for the NORTHEAST (North Atlantic + Mid-Atlantic); 187 Data reg_compile; 188 Merge HH_1 abun_change; 189 by reg; 190 array tmn(5) NE_baseline NE_tmn1 NE_tmn2 NE_tmn3 NE_tmn4; 191 array WTPhh(5) WTPhh_baseline WTPhh_tmn1 WTPhh_tmn2 WTPhh_tmn3 WTPhh_tmn4; 192 array WTPreg(5) WTPreg_baseline WTPreg_tmn1 WTPreg_tmn2 WTPreg_tmn3 WTPreg_tmn4; 193 do i=1 to 5; 194 *WTPhh(i) is the household WTP per year for the baseline percentage changes and the change 194 ! under each option; 195 WTPhh(i)=tmn(i)*WTPper; 196 *WTPreg(i) is the product of household WTP per and the the total number of households in 196 ! each region; 197 *Reginal WTP is expressed in 1,000s of dollars; 198 WTPreg(i)=WTPhh(i)*numHH/1000; Page 8bt-star_main_03.21.11_log.txt

199 END; 200 DROP baseline tmn1 tmn2 tmn3 tmn4 i; 201 RUN;

NOTE: There were 3 observations read from the data set WORK.HH_1. NOTE: There were 3 observations read from the data set WORK.ABUN_CHANGE. NOTE: The data set WORK.REG_COMPILE has 3 observations and 19 variables. NOTE: DATA statement used (Total process time): 9 The SAS System 12:32 Tuesday, March 22, 2011

real time 0.09 seconds cpu time 0.01 seconds

202

************************************************************************************

** 203

************************************************************************************

** 204 ALL REMAINING STEPS ARE RELATED TO DISCOUNTING/COMPLIANCE SCHEDULE!; 205 *This step trims the I&E input data set and limits it to age 1 equivalents for the species 205 ! of interest "species_1"; 206 *Table ie_2 sums up the total baselines I&E and reductions for each species over the 206 ! entirety of the Northeast (North Atlantic plus Mid-Atlantic); 207 *The purpose of the following two data steps is to manipulate the data so the compliance 207 ! schedule can be properly applied, 208 especially considering the impingement fraction; 209 data ie; 210 set ie; 211 if species_econ~="winter.flounder" then delete; 212 if metric~="ae1" then delete; 213 if comrec_type~="all.types" then delete; 214 if region="zz.allregions" then delete; 215 RUN;

NOTE: There were 24640 observations read from the data set WORK.IE. NOTE: The data set WORK.IE has 8 observations and 10 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

216 217 218 Proc sql; 219 create table ie_2 as 220 select option_name, 221 region, 222 species_econ, 223 lossmode, 224 metric, 225 comrec_type, 226 baseline_value_regional, 227 benefit_value_regional, 228 impingement_fraction, 229 sum(baseline_value_regional) as baseline_sum, 230 sum(benefit_value_regional) as benefit_sum 231 From ie group by option_name,species_econ; Page 9bt-star_main_03.21.11_log.txt

NOTE: The query requires remerging summary statistics back with the original data. NOTE: Table WORK.IE_2 created, with 8 rows and 11 columns.

232 Quit; NOTE: PROCEDURE SQL used (Total process time): real time 0.04 seconds cpu time 0.01 seconds

10 The SAS System 12:32 Tuesday, March 22, 2011

233 Data ie_3; 234 Set ie_2; 235 IF species_econ~="winter.flounder" then delete; 236 baseline_frac=baseline_value_regional/baseline_sum; 237 benefit_frac=benefit_value_regional/benefit_sum; 238 IF option_name="option.tmn1" then option="tmn1"; 239 IF option_name="option.tmn2" then option="tmn2"; 240 IF option_name="option.tmn3" then option="tmn3"; 241 IF option_name="option.tmn4" then option="tmn4"; 242 DROP option_name; 243 reg=region; 244 DROP region; 245 RUN;

NOTE: There were 8 observations read from the data set WORK.IE_2. NOTE: The data set WORK.IE_3 has 8 observations and 13 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

246 *need to input the compliance schedule; 247 248 *Import facility data with both facility id's; 249 %import(&inloc.&input,IDMAP,f1);

NOTE: WORK.F1 data set was successfully created. NOTE: PROCEDURE IMPORT used (Total process time): real time 0.62 seconds cpu time 0.17 seconds

250 *Import data with weights; 251 %import(&inloc.&input,ID_AND_WEIGHT,ID_AND_WEIGHT);

NOTE: WORK.ID_AND_WEIGHT data set was successfully created. NOTE: PROCEDURE IMPORT used (Total process time): real time 0.45 seconds cpu time 0.10 seconds

252 *Import data with region and flow; 253 %import(&inloc.&input,flowdata,f3);

NOTE: WORK.F3 data set was successfully created. NOTE: PROCEDURE IMPORT used (Total process time): real time 0.54 seconds cpu time 0.18 seconds

254 *Import data with compliance dates; Page 10bt-star_main_03.21.11_log.txt

255 %import(&inloc.&input,complianceyear,cydata);

NOTE: WORK.CYDATA data set was successfully created. 11 The SAS System 12:32 Tuesday, March 22, 2011

NOTE: PROCEDURE IMPORT used (Total process time): real time 0.40 seconds cpu time 0.12 seconds

256 257 %macro cyr; 258 cydata.tmn1 as cyr_base, 259 cydata.tmn1 as cyr_tmn1, 260 cydata.tmn2 as cyr_tmn2, 261 cydata.tmn3 as cyr_tmn3, 262 cydata.tmn4 as cyr_tmn4 263 %mend; 264 265 *Get ID for joining datasets; 266 proc sql; 267 create table f4 as select 268 id.plntcode, 269 w.id as facility, 270 w.ps_wgt as Weighting_Factor 271 from f1 as id right join ID_AND_WEIGHT as w on id.surveyid=w.id; NOTE: Table WORK.F4 created, with 892 rows and 3 columns.

272 quit; NOTE: PROCEDURE SQL used (Total process time): real time 0.03 seconds cpu time 0.01 seconds

273 274 *Drop double counted plants; 275 proc sort data=f4 out=f5; by plntcode; run;

NOTE: There were 892 observations read from the data set WORK.F4. NOTE: The data set WORK.F5 has 892 observations and 3 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

276 data f6; 277 set f5; 278 by plntcode; 279 if (not first.plntcode | not last.plntcode) and Weighting_Factor=0 then delete; 280 run;

NOTE: There were 892 observations read from the data set WORK.F5. NOTE: The data set WORK.F6 has 892 observations and 3 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

281 12 The SAS System 12:32 Tuesday, March 22, 2011 Page 11bt-star_main_03.21.11_log.txt

282 *join flow and weight data; 283 proc sql; 284 create table f7 as select 285 f6.*, 286 f3.Benefits_Region,f3.AverageFlow 287 from f6 full join f3 on f6.plntcode=f3.EIA_Code where f3.EIA_Code~=.; NOTE: Table WORK.F7 created, with 892 rows and 5 columns.

288 quit; NOTE: PROCEDURE SQL used (Total process time): real time 0.03 seconds cpu time 0.01 seconds

289 290 *recode region variable; 291 *delete facilities for other regions; 292 data f8 (keep=region AverageFlow Weighting_Factor PLNTCODE facility); 293 set f7; 294 format region $14.; 295 if Benefits_Region="CN" | Benefits_Region="CS" | Benefits_Region="CAL" | 295 ! Benefits_Region="HI" then delete; 296 else if Benefits_Region="GL" then delete; 297 else if Benefits_Region="INL" then delete; 298 else if Benefits_Region="MA" then region="Mid-Atlantic";

299 else if Benefits_Region="Mex" then delete; 300 else if Benefits_Region="NA" then region="North Atlantic"; 301 else if Benefits_Region="SA" then delete; 302 else if Benefits_Region="VI" then delete; 303 else if Benefits_Region="PNW" then delete; 304 run;

NOTE: There were 892 observations read from the data set WORK.F7. NOTE: The data set WORK.F8 has 80 observations and 5 variables. NOTE: DATA statement used (Total process time): real time 0.03 seconds cpu time 0.00 seconds

305 306 *combine and format data with average flow, region, wgt, and compliance

year; 307 proc sql; 308 create table s2a as select 309 f8.plntcode,f8.AverageFlow as avgflow,f8.Region,f8.Weighting_Factor as newweight, 310 %cyr 311 from f8 full join cydata on cydata.id=f8.facility; NOTE: Table WORK.S2A created, with 542 rows and 9 columns.

Page 12bt-star_main_03.21.11_log.txt

312 quit; NOTE: PROCEDURE SQL used (Total process time): real time 0.03 seconds cpu time 0.00 seconds

13 The SAS System 12:32 Tuesday, March 22, 2011

313 RUN; 314 DATA s2a; 315 SET s2a; 316 if region ="Inland" | region="California" | region="South Atlantic" | region="South 316 ! Atlantic" | region="Great Lakes" | region="Gulf of Mexico" then delete; 317 if cyr_base="." then delete; 318 if PLNTCODE="." then delete; 319 RUN;

NOTE: Character values have been converted to numeric values at the places given by: (Line):(Column). 317:14 318:14 NOTE: There were 542 observations read from the data set WORK.S2A. NOTE: The data set WORK.S2A has 36 observations and 9 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

320 321 322 *5j This macro calculates benefits schedule for impingement only facilities; 323 *assumes no biological lag; 324 *this step is implemented seperately for facilities getting towers; 325 %macro step5schedule(opt,cyr); 326 * 5j1 generate benefits schedule; 327 data s3; *s3 contains plntcodes weights and compliance years along with the flow arrays 327 ! (flow arrays are equal to AIF) 328 and includes the biological fraction as well. It multiplies these to get the biological 328 ! schedule; 329 set s2a; 330 where region~="Virgin Islands"; 331 array flowy_(60) flowy_1-flowy_60; 332 array bio(60) bio_1-bio_60; 333 array bschd(60) bschd_1-bschd_60; 334 do i=1 to 60; 335 336 if "&opt"="baseline" then do; 337 flowy_(i)=avgflow; 338 if region="Inland" | region="Great Lakes" | region="Gulf of Mexico" then do; 338 ! *benefits begin in first year for these regions; 339 *this assumes no biological lag; 340 if i=1 then bio(i)=1; 341 if i=2 then bio(i)=1; 342 if i=3 then bio(i)=1; Page 13bt-star_main_03.21.11_log.txt

343 if i=4 then bio(i)=1; 344 if i=5 then bio(i)=1; 345 if i=6 then bio(i)=1; 346 if i>6 and i<=45 then bio(i)=1; 347 if i=46 then bio(i)=.83; 348 if i=47 then bio(i)=.67; 349 if i=48 then bio(i)=.5; 350 if i=49 then bio(i)=.33; 351 if i=50 then bio(i)=.17; 352 if i>50 then bio(i)=0; 14 The SAS System 12:32 Tuesday, March 22, 2011

353 end; 354 else do; *benefits also begin in the first year for these regions; 355 if i=1 then bio(i)=1; 356 if i=2 then bio(i)=1; 357 if i=3 then bio(i)=1; 358 if i=4 then bio(i)=1; 359 if i=5 then bio(i)=1; 360 if i=6 then bio(i)=1; 361 if i=7 then bio(i)=1; 362 if i>7 and i<=45 then bio(i)=1; 363 if i=46 then bio(i)=.83; 364 if i=47 then bio(i)=.67; 365 if i=48 then bio(i)=.5; 366 if i=49 then bio(i)=.33; 367 if i=50 then bio(i)=.17; 368 if i>50 then bio(i)=0; 369 end; 370 end; 371 372 else do; 373 if &cyr.>=&ie_year. then delete; 374 if region="Inland" | region="Great Lakes" | region="Gulf of Mexico" then do; 375 *benefits Page 14bt-star_main_03.21.11_log.txt

begin in first year for these regions; 376 x=i-&cyr.+&t0; *e.g., if i=1 and cyr = 2011 and first year of rule is 2011 then x=1 376 ! in 2011; 377 if x<=0 then flowy_(i)=0; *x=1 in first year of compliance; 378 else flowy_(i)=avgflow; 379 if _n_=1 then put "Not BASELINE"; 380 381 if x<=0 then bio(i)=0; 382 else if x=1 then bio(i)=.2; 383 else if x=2 then bio(i)=.4; 384 else if x=3 then bio(i)=.6; 385 else if x=4 then bio(i)=.8; 386 else if x=5 then bio(i)=1; 387 else if x=6 then bio(i)=1; 388 else if x=7 then bio(i)=1; 389 else if x=8 then bio(i)=1; 390 else if x=9 then bio(i)=1; 391 else if x>=9

and i<=45 then bio(i)=1; 392 else if i=46 then bio(i)=.83; 393 else if i=47 then bio(i)=.67; 394 else if i=48 then bio(i)=.5; 395 else if i=49 then bio(i)=.33; 396 else if i=50 then bio(i)=0.17; 397 else if i>50 then bio(i)=0; 398 end; 399 else do; 400 *benefits also begin in first year for these regions; 401 x=i-&cyr.+&t0; *e.g., if i=1 and cyr = 2011 and first year of rule is 2011 then x=1 401 ! in 2012; 402 if x<0 then flowy_(i)=0;*x=1 in second year of compliance; 403 else flowy_(i)=avgflow; 404 if _n_=1 then put "Not BASELINE"; 15 The SAS System 12:32 Tuesday, March 22, 2011

405 406 if Page 15bt-star_main_03.21.11_log.txt

x<=0 then bio(i)=0; 407 else if x=1 then bio(i)=.2; 408 else if x=2 then bio(i)=.4; 409 else if x=3 then bio(i)=.6; 410 else if x=4 then bio(i)=.8; 411 else if x=5 then bio(i)=1; 412 else if x=6 then bio(i)=1; 413 else if x=7 then bio(i)=1; 414 else if x=8 then bio(i)=1; 415 else if x=9 then bio(i)=1; 416 else if x=10

then bio(i)=1; 417 else if x>=10 and i<=45 then bio(i)=1; 418 else if i=46 then bio(i)=.83; 419 else if i=47 then bio(i)=.67; 420 else if i=48 then bio(i)=.5; 421 else if i=49 then bio(i)=.33; 422 else if i=50 then bio(i)=0.17; 423 else if i>50 then bio(i)=0; 424 end; 425 end; 426 *the fractions here only reflect the fact that the intial compliance year for each facility 426 ! reflects a five year average, 427 in the inputs dataset, these values are converted to the first year by subtracting two, 428 the steps of 0.2 reflect staggered implementation over this five year period, 429 there is not biological lag incorporated in the BT model, this is strictly related to 429 ! compliance; 430

431 bschd(i)=(flowy_(i)*bio(i)); 432 end; 433 format reg $30.; 434 if region="California"

then reg="calif"; 435 else if region="Great Lakes" then reg="greatlakes"; 436 else if region="Gulf of Mexico" then reg="gulfmexico"; 437 else if region="Inland"

then reg="inland"; 438 else if region="Mid-Atlantic" then reg="midatlantic"; 439 else if region="North Atlantic" Page 16bt-star_main_03.21.11_log.txt

then reg="northatlantic"; 440 else if region="South Atlantic" then reg="southatlantic"; 441 442 run; 443 444 * 5j2 calculate total average flow by region; 445 proc sql; 446 create table s5 as select 447 s3.*, 448 sum(avgflow*newweight) as totalavgflow /*group flag*/ 449 from s3 group by region; 450 quit; 451 452 * 5j3 generate flow and benefits schedule variables (for each facility) as percentages to 452 ! total regional average flow; 453 data s6 (drop=bio_1-bio_60 flowy_1-flowy_60); 454 set s5; 455 array flowy_(60) flowy_1-flowy_60; 16 The SAS System 12:32 Tuesday, March 22, 2011

456 array bschd(60) bschd_1-bschd_60; 457 do i=1 to 60; 458 IF totalavgflow~=0 then flowy_(i)=newweight*flowy_(i)/totalavgflow; 459 IF totalavgflow=0 then flowy_(i)=0; 460 IF totalavgflow~=0 then bschd(i)=newweight*bschd(i)/totalavgflow; 461 IF totalavgflow=0 then bschd(i)=0; 462 end; 463 run; 464 465 * 5j4 sort data by region; 466 proc sort data=s6 out=s7; by region; run; 467 * 5j6 generate benefits schedule variable for region, e.g., 1 = full benefits, .5 = 50% 467 ! of full benefits; 468 data s8_&opt (keep=option region reg rbschd_1-rbschd_60); 469 retain option region reg rbschd_1-rbschd_60; 470 set ; 471 by region; 472 option="&opt"; 473 array bschd(60) bschd_1-bschd_60; 474 array rbschd(60) rbschd_1-rbschd_60; 475 do i=1 to 60; 476 if first.region then rbschd(i)=bschd(i); 477 else rbschd(i)+bschd(i); 478 end; 479 if last.region; 480 run; 481 %mend; 482 *This step generates the benefits schedule for facilities assigned
towers; 483 *assumes no biological lag; 484 %macro step6schedule(opt,cyr); 485 data s4; *s3 contains plntcodes weights and compliance years along with the flow arrays 485 ! (flow arrays are equal to AIF) 486 and includes the biological fraction as well. It multiplies these to get the biological Page 17bt-star_main_03.21.11_log.txt

486 ! schedule; 487 set s2a; 488 where region~="Virgin Islands"; 489 array flowy_(60) flowy_1-flowy_60; 490 array bio(60) bio_1-bio_60; 491 array bschd(60) bschd_1-bschd_60; 492 do i=1 to 60; 493 494 if "&opt"="baseline" then do; 495 flowy_(i)=avgflow; 496 if region="Inland" | region="Great Lakes" | region="Gulf of Mexico" then do; 497 if i=1 then bio(i)=1; 498 if i=2 then bio(i)=1; 499 if i=3 then bio(i)=1; 500 if i=4 then bio(i)=1; 501 if i=5 then bio(i)=1; 502 if i=6 then bio(i)=1; 503 if i>6 and i<=45 then bio(i)=1; 504 if i=46 then bio(i)=.83; 505 if i=47 then bio(i)=.67; 506 if i=48 then bio(i)=.5; 17 The SAS System 12:32 Tuesday, March 22, 2011

507 if i=49 then bio(i)=.33; 508 if i=50 then bio(i)=.17; 509 if i>50 then bio(i)=0;*benefits begin in first year for these regions; 510 end; 511 else do; 512 if i=1 then bio(i)=1; 513 if i=2 then bio(i)=1; 514 if i=3 then bio(i)=1; 515 if i=4 then bio(i)=1; 516 if i=5 then bio(i)=1; 517 if i=6 then bio(i)=1; 518 if i=7 then bio(i)=1; 519 if i>7 and i<=45 then bio(i)=1; 520 if i=46 then bio(i)=.83; 521 if i=47 then bio(i)=.67; Page 18bt-star_main_03.21.11_log.txt

522 if i=48 then bio(i)=.5; 523 if i=49 then bio(i)=.33; 524 if i=50 then bio(i)=.17; 525 if i>50 then bio(i)=0;*benefits also begin in the first for these regions; 526 end; 527 end; 528 529 else do; 530 if &cyr.<&ie_year. then delete; 531 if region="Inland" | region="Great Lakes" | region="Gulf of Mexico" then do; 532 *benefits begin in first year for these regions; 533 x=i-&cyr.+&t0; *e.g., if i=1 and cyr = 2011 and first year of rule is 2011 then x=1 533 ! in 2011; 534 if x<=0 then flowy_(i)=0; *x=1 in first year of compliance; 535 else flowy_(i)=avgflow; 536 if _n_=1 then put "Not BASELINE"; 537 538 if x<=0 then bio(i)=0; 539 else if x=1 then bio(i)=.2; 540 else if x=2 then bio(i)=.4; 541 else if x=3 then bio(i)=.6; 542 else if x=4 then bio(i)=.8; 543 else if x=5 then bio(i)=1; 544 else if x=6 then bio(i)=1; 545 else if x=7 then bio(i)=1; 546 else if x=8 then bio(i)=1; 547 else if x=9 then bio(i)=1; 548 else if x>=9

and i<=45 then bio(i)=1; 549 else if i=46 then bio(i)=.83; 550 else if i=47 then bio(i)=.67; 551 else if i=48 then bio(i)=.5; 552 else if i=49 then bio(i)=.33; 553 else if i=50 then bio(i)=0.17; 554 else if i>50 then bio(i)=0; 555 end; Page 19bt-star_main_03.21.11_log.txt

556 else do; 557 *benefits also begin in first year for these regions; 558 x=i-&cyr.+&t0; *e.g., if i=1 and cyr = 2011 and first year of rule is 2011 then x=1 558 ! in 2012; 18 The SAS System 12:32 Tuesday, March 22, 2011

559 if x<0 then flowy_(i)=0;*x=1 in second year of compliance; 560 else flowy_(i)=avgflow; 561 if _n_=1 then put "Not BASELINE"; 562 563 if x<=0 then bio(i)=0; 564 else if x=1 then bio(i)=.2; 565 else if x=2 then bio(i)=.4; 566 else if x=3 then bio(i)=.6; 567 else if x=4 then bio(i)=.8; 568 else if x=5 then bio(i)=1; 569 else if x=6 then bio(i)=1; 570 else if x=7 then bio(i)=1; 571 else if x=8 then bio(i)=1; 572 else if x=9 then bio(i)=1; 573 else if x=10

then bio(i)=1; 574 else if x>=10 and i<=45 then bio(i)=1; 575 else if i=46 then bio(i)=.83; 576 else if i=47 then bio(i)=.67; 577 else if i=48 then bio(i)=.5; 578 else if i=49 then bio(i)=.33; 579 else if i=50 then bio(i)=0.17; 580 else if i>50 then bio(i)=0; 581 end; 582 end; 583 *the fractions here only reflect the fact that the intial compliance year for each facility 583 ! reflects a five year average, 584 in the inputs dataset, these values are converted to the first year by subtracting two, 585 the steps of 0.2 reflect staggered implementation over this five year period, 586 there is not biological lag incorporated in the BT model, this is strictly related to Page 20bt-star_main_03.21.11_log.txt

586 ! compliance; 587 588 bschd(i)=(flowy_(i)*bio(i)); 589 end; 590 format reg $30.; 591 if region="California"

then reg="calif"; 592 else if region="Great Lakes" then reg="greatlakes"; 593 else if region="Gulf of Mexico" then reg="gulfmexico"; 594 else if region="Inland"

then reg="inland"; 595 else if region="Mid-Atlantic" then reg="midatlantic"; 596 else if region="North Atlantic" then reg="northatlantic"; 597 else if region="South Atlantic" then reg="southatlantic"; 598 599 run; 600 601 *5j2 calculate total average flow by region; 602 proc sql; 603 create table s5 as select 604 s4.*, 605 sum(avgflow*newweight) as totalavgflow /*group flag*/ 606 from s4 group by region; 607 quit; 608 609 * 5j3 generate flow and benefits schedule variables (for each facility) as percentages to 609 ! total regional average flow; 19 The SAS System 12:32 Tuesday, March 22, 2011

610 data s6 (drop=bio_1-bio_60 flowy_1-flowy_60); 611 set s5; 612 array flowy_(60) flowy_1-flowy_60; 613 array bschd(60) bschd_1-bschd_60; 614 do i=1 to 60; 615 IF totalavgflow~=0 then flowy_(i)=newweight*flowy_(i)/totalavgflow; 616 IF totalavgflow=0 then flowy_(i)=0; 617 IF totalavgflow~=0 then bschd(i)=newweight*bschd(i)/totalavgflow; 618 IF totalavgflow=0 then bschd(i)=0; 619 end; 620 run; 621 622 * 5j4 sort data by region; 623 proc sort data=s6 out=s25; by region; run; 624 * 5j6 generate benefits schedule variable for region, e.g., 1 = full benefits, .5 = 50% 624 ! of full benefits; 625 data s9_&opt (keep=option region reg rbschd_1-rbschd_60); 626 retain option region reg rbschd_1-rbschd_60; 627 set s25; 628 by region; 629 option="&opt"; 630 array bschd(60) bschd_1-bschd_60; 631 array rbschd(60) rbschd_1-rbschd_60; 632 do i=1 to 60; Page 21bt-star_main_03.21.11_log.txt

633 if first.region then rbschd(i)=bschd(i); 634 else rbschd(i)+bschd(i); 635 end; 636 if last.region; 637 run; 638 %mend; 639 * 5j This macro calculates benefits schedule, incorporating the biological lag and the lag 639 ! in compliance; 640 * 5k Generate dataset of Benefits Schedule By Option; 641 %macro optionmac00; 642 %step5schedule(baseline,cyr_base); 643 %step5schedule(tmn1,cyr_tmn1); 644 %step5schedule(tmn2,cyr_tmn2); 645 %step5schedule(tmn3,cyr_tmn3); 646 %step5schedule(tmn4,cyr_tmn4); 647 %mend; 648 649 %macro optionmac11; 650 %step6schedule(baseline,cyr_base); 651 %step6schedule(tmn1,cyr_tmn1); 652 %step6schedule(tmn2,cyr_tmn2); 653 %step6schedule(tmn3,cyr_tmn3); 654 %step6schedule(tmn4,cyr_tmn4); 655 %mend; 656 %optionmac00;

NOTE: There were 36 observations read from the data set WORK.S2A. WHERE region not = 'Virgin Islands'; NOTE: The data set WORK.S3 has 36 observations and 192 variables. NOTE: DATA statement used (Total process time): 20 The SAS System 12:32 Tuesday, March 22, 2011

real time 0.03 seconds cpu time 0.01 seconds

NOTE: The query requires remerging summary statistics back with the original data. NOTE: Table WORK.S5 created, with 36 rows and 193 columns.

NOTE: PROCEDURE SQL used (Total process time): real time 0.59 seconds cpu time 0.01 seconds

NOTE: There were 36 observations read from the data set WORK.S5. NOTE: The data set WORK.S6 has 36 observations and 73 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: There were 36 observations read from the data set WORK.S6. NOTE: The data set WORK.S7 has 36 observations and 73 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

Page 22bt-star_main_03.21.11_log.txt

NOTE: There were 36 observations read from the data set WORK.S7. NOTE: The data set WORK.S8_BASELINE has 2 observations and 63 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE 21 The SAS System 12:32 Tuesday, March 22, 2011

Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE Page 23bt-star_main_03.21.11_log.txt

Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE NOTE: There were 36 observations read from the data set WORK.S2A. WHERE region not = 'Virgin Islands'; NOTE: The data set WORK.S3 has 36 observations and 192 variables. NOTE: DATA statement used (Total process time): real time 0.03 seconds cpu time 0.01 seconds

NOTE: The query requires remerging summary statistics back with the original data. NOTE: Table WORK.S5 created, with 36 rows and 193 columns.

22 The SAS System 12:32 Tuesday, March 22, 2011

NOTE: PROCEDURE SQL used (Total process time): real time 0.03 seconds cpu time 0.01 seconds

NOTE: There were 36 observations read from the data set WORK.S5. NOTE: The data set WORK.S6 has 36 observations and 73 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

NOTE: There were 36 observations read from the data set WORK.S6. NOTE: The data set WORK.S7 has 36 observations and 73 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.03 seconds cpu time 0.01 seconds

NOTE: There were 36 observations read from the data set WORK.S7. NOTE: The data set WORK.S8_TMN1 has 2 observations and 63 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE Page 24bt-star_main_03.21.11_log.txt

Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE 23 The SAS System 12:32 Tuesday, March 22, 2011

Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE NOTE: There were 36 observations read from the data set WORK.S2A. WHERE region not = 'Virgin Islands'; NOTE: The data set WORK.S3 has 16 observations and 192 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: The query requires remerging summary statistics back with the original data. NOTE: Table WORK.S5 created, with 16 rows and 193 columns. Page 25bt-star_main_03.21.11_log.txt

NOTE: PROCEDURE SQL used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: There were 16 observations read from the data set WORK.S5. 24 The SAS System 12:32 Tuesday, March 22, 2011

NOTE: The data set WORK.S6 has 16 observations and 73 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: There were 16 observations read from the data set WORK.S6. NOTE: The data set WORK.S7 has 16 observations and 73 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

NOTE: There were 16 observations read from the data set WORK.S7. NOTE: The data set WORK.S8_TMN2 has 2 observations and 63 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

NOTE: There were 36 observations read from the data set WORK.S2A. WHERE region not = 'Virgin Islands'; NOTE: The data set WORK.S3 has 1 observations and 192 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

NOTE: The query requires remerging summary statistics back with the original data. NOTE: Table WORK.S5 created, with 1 rows and 193 columns.

NOTE: PROCEDURE SQL used (Total process time): real time 0.03 seconds cpu time 0.04 seconds

NOTE: There were 1 observations read from the data set WORK.S5. NOTE: The data set WORK.S6 has 1 observations and 73 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

NOTE: There were 1 observations read from the data set WORK.S6. NOTE: The data set WORK.S7 has 1 observations and 73 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.00 seconds Page 26bt-star_main_03.21.11_log.txt

cpu time 0.00 seconds

25 The SAS System 12:32 Tuesday, March 22, 2011

NOTE: There were 1 observations read from the data set WORK.S7. NOTE: The data set WORK.S8_TMN3 has 1 observations and 63 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE 26 The SAS System 12:32 Tuesday, March 22, 2011

Not BASELINE Page 27bt-star_main_03.21.11_log.txt

Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE NOTE: There were 36 observations read from the data set WORK.S2A. WHERE region not = 'Virgin Islands'; NOTE: The data set WORK.S3 has 36 observations and 192 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: The query requires remerging summary statistics back with the original data. NOTE: Table WORK.S5 created, with 36 rows and 193 columns.

NOTE: PROCEDURE SQL used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: There were 36 observations read from the data set WORK.S5. NOTE: The data set WORK.S6 has 36 observations and 73 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

NOTE: There were 36 observations read from the data set WORK.S6. NOTE: The data set WORK.S7 has 36 observations and 73 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: There were 36 observations read from the data set WORK.S7. NOTE: The data set WORK.S8_TMN4 has 2 observations and 63 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds 27 The SAS System 12:32 Tuesday, March 22, 2011

657 %optionmac11;

NOTE: There were 36 observations read from the data set WORK.S2A. WHERE region not = 'Virgin Islands'; NOTE: The data set WORK.S4 has 36 observations and 192 variables. Page 28bt-star_main_03.21.11_log.txt

NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: The query requires remerging summary statistics back with the original data. NOTE: Table WORK.S5 created, with 36 rows and 193 columns.

NOTE: PROCEDURE SQL used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: There were 36 observations read from the data set WORK.S5. NOTE: The data set WORK.S6 has 36 observations and 73 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: There were 36 observations read from the data set WORK.S6. NOTE: The data set WORK.S25 has 36 observations and 73 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

NOTE: There were 36 observations read from the data set WORK.S25. NOTE: The data set WORK.S9_BASELINE has 2 observations and 63 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

NOTE: There were 36 observations read from the data set WORK.S2A. WHERE region not = 'Virgin Islands'; NOTE: The data set WORK.S4 has 0 observations and 192 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: The query requires remerging summary statistics back with the original data. 28 The SAS System 12:32 Tuesday, March 22, 2011

NOTE: Table WORK.S5 created, with 0 rows and 193 columns.

NOTE: PROCEDURE SQL used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: There were 0 observations read from the data set WORK.S5. NOTE: The data set WORK.S6 has 0 observations and 73 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 secondsPage 29bt-star_main_03.21.11_log.txt

NOTE: Input data set is empty. NOTE: The data set WORK.S25 has 0 observations and 73 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

NOTE: There were 0 observations read from the data set WORK.S25. NOTE: The data set WORK.S9_TMN1 has 0 observations and 63 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

NOTE: There were 36 observations read from the data set WORK.S2A. WHERE region not = 'Virgin Islands'; NOTE: The data set WORK.S4 has 20 observations and 192 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: The query requires remerging summary statistics back with the original data. NOTE: Table WORK.S5 created, with 20 rows and 193 columns.

NOTE: PROCEDURE SQL used (Total process time): real time 0.03 seconds cpu time 0.03 seconds

NOTE: There were 20 observations read from the data set WORK.S5. NOTE: The data set WORK.S6 has 20 observations and 73 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds 29 The SAS System 12:32 Tuesday, March 22, 2011

NOTE: There were 20 observations read from the data set WORK.S6. NOTE: The data set WORK.S25 has 20 observations and 73 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.03 seconds cpu time 0.00 seconds

NOTE: There were 20 observations read from the data set WORK.S25. NOTE: The data set WORK.S9_TMN2 has 2 observations and 63 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

Page 30bt-star_main_03.21.11_log.txt

Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE 30 The SAS System 12:32 Tuesday, March 22, 2011

Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE
Not BASELINE Page 31bt-star_main_03.21.11_log.txt

NOTE: There were 36 observations read from the data set WORK.S2A. WHERE region not = 'Virgin Islands'; NOTE: The data set WORK.S4 has 35 observations and 192 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

NOTE: The query requires remerging summary statistics back with the original data. NOTE: Table WORK.S5 created, with 35 rows and 193 columns.

NOTE: PROCEDURE SQL used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

NOTE: There were 35 observations read from the data set WORK.S5. NOTE: The data set WORK.S6 has 35 observations and 73 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

NOTE: There were 35 observations read from the data set WORK.S6. NOTE: The data set WORK.S25 has 35 observations and 73 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.03 seconds 31 The SAS System 12:32 Tuesday, March 22, 2011

cpu time 0.00 seconds

NOTE: There were 35 observations read from the data set WORK.S25. NOTE: The data set WORK.S9_TMN3 has 2 observations and 63 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

NOTE: There were 36 observations read from the data set WORK.S2A. WHERE region not = 'Virgin Islands'; NOTE: The data set WORK.S4 has 0 observations and 192 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: The query requires remerging summary statistics back with the original data. NOTE: Table WORK.S5 created, with 0 rows and 193 columns.

NOTE: PROCEDURE SQL used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: There were 0 observations read from the data set WORK.S5. NOTE: The data set WORK.S6 has 0 observations and 73 variables. Page 32bt-star_main_03.21.11_log.txt

NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

NOTE: Input data set is empty. NOTE: The data set WORK.S25 has 0 observations and 73 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

NOTE: There were 0 observations read from the data set WORK.S25. NOTE: The data set WORK.S9_TMN4 has 0 observations and 63 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

658 *step5data contains the benefits schedules for all options for impingement only; 659 data Step5data; 32 The SAS System 12:32 Tuesday, March 22, 2011

660 length option $50.0; 661 set s8_baseline s8_tmn1 s8_tmn2 s8_tmn3 s8_tmn4; 662 IF region="" THEN DELETE; 663 DROP region; 664 array rbschd_(60) rbschd_1-rbschd_60; 665 do i=1 to 60; 666 *know that tmn2 does not have any impingement only, towers everywhere; 667 IF option="" then rbschd_(i)=0; 668 END; 669 run;

NOTE: There were 2 observations read from the data set WORK.S8_BASELINE. NOTE: There were 2 observations read from the data set WORK.S8_TMN1. NOTE: There were 2 observations read from the data set WORK.S8_TMN2. NOTE: There were 1 observations read from the data set WORK.S8_TMN3. NOTE: There were 2 observations read from the data set WORK.S8_TMN4. NOTE: The data set WORK.STEP5DATA has 9 observations and 63 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

670 proc sort data=step5data out=step5data_1; 671 by reg; 672 RUN;

NOTE: There were 9 observations read from the data set WORK.STEP5DATA. NOTE: The data set WORK.STEP5DATA_1 has 9 observations and 63 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

673 *step20data contains the benefits schedules for all options for towers; 674 data Step20data; Page 33bt-star_main_03.21.11_log.txt

675 length option $50.0; 676 set s9_baseline s9_tmn1 s9_tmn2 s9_tmn3 s9_tmn4; 677 IF region="" THEN DELETE; 678 DROP region; 679 array rbschd_(60) rbschd_1-rbschd_60; 680 array sbschd_(60) sbschd_1-sbschd_60; 681 Do i=1 to 60; 682 sbschd_(i)=rbschd_(i); 683 DROP rbschd_1-rbschd_60; 684 END; 685 run;

NOTE: There were 2 observations read from the data set WORK.S9_BASELINE. NOTE: There were 0 observations read from the data set WORK.S9_TMN1. NOTE: There were 2 observations read from the data set WORK.S9_TMN2. NOTE: There were 2 observations read from the data set WORK.S9_TMN3. NOTE: There were 0 observations read from the data set WORK.S9_TMN4. NOTE: The data set WORK.STEP20DATA has 6 observations and 63 variables. NOTE: DATA statement used (Total process time): 33 The SAS System 12:32 Tuesday, March 22, 2011

real time 0.00 seconds cpu time 0.00 seconds

686 proc sort data=step20data out=step20data_1; 687 by reg; 688 RUN;

NOTE: There were 6 observations read from the data set WORK.STEP20DATA. NOTE: The data set WORK.STEP20DATA_1 has 6 observations and 63 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

689 *comp_calcs contains the seperate benefit schedules for each for impingement only 689 ! (step5data) and towers (step20data); 690 data comp_calcs; 691 merge step5data_1 step20data_1; 692 by reg option; 693 run;

NOTE: There were 9 observations read from the data set WORK.STEP5DATA_1. NOTE: There were 6 observations read from the data set WORK.STEP20DATA_1. NOTE: The data set WORK.COMP_CALCS has 10 observations and 123 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

694 Proc sort data=comp_calcs out=comp_calcs; 695 by reg option; 696 Run;

NOTE: There were 10 observations read from the data set WORK.COMP_CALCS. NOTE: The data set WORK.COMP_CALCS has 10 observations and 123 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

Page 34bt-star_main_03.21.11_log.txt

697 *ie_imp brings in the impingement fraction which needs to be applied to the comp_calcs 697 ! sheet; 698 Data ie_imp; 699 LENGTH option $20; 700 SET ie (KEEP = option_name region metric comrec_type species_econ impingement_fraction); 701 If option_name="option.tmn1" then option="tmn1"; 702 else if option_name="option.tmn2" then option="tmn2"; 703 else if option_name="option.tmn3" then option="tmn3"; 704 else if option_name="option.tmn4" then option="tmn4"; 705 if species_econ~="winter.flounder" then delete; 706 if region="zz.allregions" then delete; 707 if metric~="ae1" then delete; 708 if comrec_type~="all.types" then delete; 34 The SAS System 12:32 Tuesday, March 22, 2011

709 drop option_name; 710 reg=region; 711 drop region; 712 drop species_econ; 713 drop metric; 714 drop comrec_type; 715 RUN;

NOTE: There were 8 observations read from the data set WORK.IE. NOTE: The data set WORK.IE_IMP has 8 observations and 3 variables. NOTE: DATA statement used (Total process time): real time 0.03 seconds cpu time 0.00 seconds

716 717 Proc sort data=ie_imp out=ie_imp; 718 by reg option; 719 Run;

NOTE: There were 8 observations read from the data set WORK.IE_IMP. NOTE: The data set WORK.IE_IMP has 8 observations and 3 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.03 seconds cpu time 0.00 seconds

720 *comp_imp merges the impingement only and towers benefits schedules with the schedules the 720 ! impingement fraction; 721 data comp_imp; 722 merge ie_imp comp_calcs; 723 by reg option; 724 if option="baseline" then impingement_fraction=0; 725 run;

WARNING: Multiple lengths were specified for the BY variable reg by input data sets. This may cause

unexpected results. WARNING: Multiple lengths were specified for the BY variable option by input data sets. This may

cause unexpected results. NOTE: There were 8 observations read from the data set WORK.IE_IMP. NOTE: There were 10 observations read from the data set WORK.COMP_CALCS. Page 35bt-star_main_03.21.11_log.txt

NOTE: The data set WORK.COMP_IMP has 10 observations and 124 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

726 *comp_calcs_1 calculates the final benefit schedule which incorporates the benefit fraction 726 ! ; 727 data comp_calcs_1; 728 set comp_imp; 729 array rbschd_(60) rbschd_1-rbschd_60; 730 array sbschd_(60) sbschd_1-sbschd_60; 731 array fschd_(60) fschd_1-fschd_60; 35 The SAS System 12:32 Tuesday, March 22, 2011

732 Do i=1 to 60; 733 IF rbschd_(i)="." then rbschd_(i)=0; 734 IF sbschd_(i)="." then sbschd_(i)=0; 735 fschd_(i)=impingement_fraction*rbschd_(i)+(1-impingement_fraction)*sbschd_(i); 736 max_fschd=round(max(of fschd_1-fschd_60),.000000000000001);*included to account for 736 ! rounding error. 737 *max(of fschd_1-fschd_60) ~= 1 then put _n_ "error flag 1 - fschd does not sum to 1"; 738 End; 739 Run;

NOTE: Character values have been converted to numeric values at the places given by: (Line):(Column). 733:16 734:16 NOTE: There were 10 observations read from the data set WORK.COMP_IMP. NOTE: The data set WORK.COMP_CALCS_1 has 10 observations and 185 variables. NOTE: DATA statement used (Total process time): real time 0.07 seconds cpu time 0.00 seconds

740 *data checks; 741 data comp_calcs_1; 742 set comp_calcs_1; 743 if max_fschd >1 then put _n_ "error flag 1 - fschd is greater than 1 " max_fschd; 744 if max_fschd <1 then put _n_ "error flag 2 - fschd is less than 1 " max_fschd; 745 RUN;

NOTE: There were 10 observations read from the data set WORK.COMP_CALCS_1. NOTE: The data set WORK.COMP_CALCS_1 has 10 observations and 185 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

746

****************************************************************************** 747 comp_calcs_2 contains the final compliance schedule for each option 748 the final biological schedules is listed vertically; 749 DATA comp_calcs_2; 750 SET comp_calcs_1; 751 Drop rbschd_1-rbschd_60; Page 36bt-star_main_03.21.11_log.txt

752 Drop sbschd_1-sbschd_60; 753 Drop impingement_fraction; 754 Drop i max_fschd; 755 RUN;

NOTE: There were 10 observations read from the data set WORK.COMP_CALCS_1. NOTE: The data set WORK.COMP_CALCS_2 has 10 observations and 62 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

756 757

****************************************************************************** 758 Start of repowering calculations; 36 The SAS System 12:32 Tuesday, March 22, 2011

759 760 %import(&inloc.&input2,rep_adj,rep_adj);

NOTE: WORK.REP_ADJ data set was successfully created. NOTE: PROCEDURE IMPORT used (Total process time): real time 5.34 seconds cpu time 0.26 seconds

761 %import(&inloc.&input,rpschd,rpschd);

NOTE: WORK.RPSCHD data set was successfully created. NOTE: PROCEDURE IMPORT used (Total process time): real time 0.43 seconds cpu time 0.09 seconds

762 763 Proc sort data=rep_adj out=rep_adj; 764 by option; 765 RUN;

NOTE: There were 240 observations read from the data set WORK.REP_ADJ. NOTE: The data set WORK.REP_ADJ has 240 observations and 4 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

766 Proc sort data=rpschd out=rpschd; 767 by option; 768 RUN;

NOTE: There were 10 observations read from the data set WORK.RPSCHD. NOTE: The data set WORK.RPSCHD has 10 observations and 61 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

769 770 Data repow; 771 merge rep_adj rpschd; 772 by option; Page 37bt-star_main_03.21.11_log.txt

773 adj=1/rp_adj; 774 drop rp_adj; 775 RUN;

NOTE: There were 240 observations read from the data set WORK.REP_ADJ. NOTE: There were 10 observations read from the data set WORK.RPSCHD. NOTE: The data set WORK.REPOW has 240 observations and 64 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds 37 The SAS System 12:32 Tuesday, March 22, 2011

776 777 Data repow_1; 778 set repow; 779 array rpschd_(60) rpschd_1-rpschd_60; 780 array rpschd_adj_(60) rpschd_adj_1-rpschd_adj_60; 781 do i=1 to 60; 782 rpschd_adj_(i)=rpschd_(i)*(1-adj); 783 if rpschd_(i)=0 then rpschd_adj_(i)=0; 784 end; 785 drop rp_adj rpschd_1-rpschd_60 i; 786 if lossmode="imp" then delete; 787 if lossmode="ent" then delete; 788 if region="zz.allregions" | region="calif" | region="southatlantic" | region="greatlakes" 788 ! | region="gulfmexico" | region="inland" then delete; 789 if option= "baseline_0" | option="tmn1_0" | option="tmn2_0" | option="tmn3_0" | 789 ! option="tmn4_0" then delete; 790 reg=region; 791 drop region lossmode; *all values only refer to IE combined; 792 RUN;

WARNING: The variable rp_adj in the DROP, KEEP, or RENAME list has never been referenced. NOTE: There were 240 observations read from the data set WORK.REPOW. NOTE: The data set WORK.REPOW_1 has 10 observations and 63 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

793 794 Proc Sort data=repow_1 out=repow_1; 795 by option reg; 796 RUN;

NOTE: There were 10 observations read from the data set WORK.REPOW_1. NOTE: The data set WORK.REPOW_1 has 10 observations and 63 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

797 798 Proc Sort data=comp_calcs_2 out=comp_calcs_2; 799 by option reg; 800 RUN;

Page 38bt-star_main_03.21.11_log.txt

NOTE: There were 10 observations read from the data set WORK.COMP_CALCS_2. NOTE: The data set WORK.COMP_CALCS_2 has 10 observations and 62 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

38 The SAS System 12:32 Tuesday, March 22, 2011

801 802 Data comp_calcs_2; 803 merge repow_1 comp_calcs_2; 804 by option reg; 805 RUN;

WARNING: Multiple lengths were specified for the BY variable option by input data sets. This may

cause unexpected results. NOTE: There were 10 observations read from the data set WORK.REPOW_1. NOTE: There were 10 observations read from the data set WORK.COMP_CALCS_2. NOTE: The data set WORK.COMP_CALCS_2 has 10 observations and 123 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

806 807 Data comp_calcs_2; 808 set comp_calcs_2; 809 array fschd_(60) fschd_1-fschd_60; 810 array rpschd_adj_(60) rpschd_adj_1-rpschd_adj_60; 811 array fschd_adj_(60) fschd_adj_1-fschd_adj_60; 812 do i = 1 to 60; 813 fschd_adj_(i)=fschd_(i)*adj + rpschd_adj_(i); 814 fschd_(i) = fschd_adj_(i); 815 END; 816 drop i rpschd_adj_1-rpschd_adj_60 fschd_adj_1-fschd_adj_60 adj; 817 RUN;

NOTE: There were 10 observations read from the data set WORK.COMP_CALCS_2. NOTE: The data set WORK.COMP_CALCS_2 has 10 observations and 62 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

818 819 *end of repowering calculations 820 *************************************************************************; 821 822 PROC SORT data=ie_3 out=ie_3; 823 by reg; 824 RUN;

NOTE: There were 8 observations read from the data set WORK.IE_3. NOTE: The data set WORK.IE_3 has 8 observations and 13 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.03 seconds cpu time 0.00 seconds

Page 39bt-star_main_03.21.11_log.txt

825 *make changes to ie_3 to delete extraneous variables; 826 Data ie_3; 827 LENGTH option $20; 39 The SAS System 12:32 Tuesday, March 22, 2011

828 Set ie_3; 829 DROP lossmode metric comrec_type baseline_value_regional benefit_value_regional 829 ! impingement_fraction baseline_sum benefit_sum; 830 DROP lossmode metric comrec_type baseline_value_regional benefit_value_regional 830 ! baseline_sum benefit_sum; 831 RUN;

NOTE: There were 8 observations read from the data set WORK.IE_3. NOTE: The data set WORK.IE_3 has 8 observations and 5 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

832 *this data step generates baseline option for inclusion in the ie data
set; 833 DATA ie_baseline; 834 Set ie_3; 835 IF option~="tmn1" then delete; 836 benefit_frac=baseline_frac; 837 option="baseline"; 838 RUN;

NOTE: There were 8 observations read from the data set WORK.IE_3. NOTE: The data set WORK.IE_BASELINE has 2 observations and 5 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

839 *this data step merges the baseline option back with the rest of the ie data set; 840 *baseline_frac can be dropped because benefit_frac and baseline_frac have been equalized 840 ! for the baseline option; 841 PROC SORT data=ie_baseline out=ie_baseline; 842 by option species_econ reg; 843 RUN;

NOTE: There were 2 observations read from the data set WORK.IE_BASELINE. NOTE: The data set WORK.IE_BASELINE has 2 observations and 5 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.03 seconds cpu time 0.01 seconds

844 PROC SORT data=ie_3 out=ie_3; 845 by option species_econ reg; 846 RUN;

NOTE: There were 8 observations read from the data set WORK.IE_3. NOTE: The data set WORK.IE_3 has 8 observations and 5 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.03 seconds cpu time 0.00 secondsPage 40bt-star_main_03.21.11_log.txt

40 The SAS System 12:32 Tuesday, March 22, 2011

847 DATA ie_all; 848 MERGE ie_baseline ie_3; 849 by option species_econ reg; 850 DROP baseline_frac; 851 RUN;

NOTE: There were 2 observations read from the data set WORK.IE_BASELINE. NOTE: There were 8 observations read from the data set WORK.IE_3. NOTE: The data set WORK.IE_ALL has 10 observations and 4 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

852 *comp_calcs_2 contains the final benefits shedule for the mid-atlantic and north atlantic 852 ! regions 853 these need to be combined in order to treat the northeast (north atlantic & mid-atlantic) 853 ! as a single population; 854 *the final schedules presented in comp_calcs_2 are generalized for all facilities and do 854 ! not yet reflect species-specific adjustments 855 related to the location of population losses (north atlantic vs mid-atlantic); 856 857 *sort comp_calcs_2 and ie_all before merging with ie_data; 858 PROC SORT data=comp_calcs_2 out=comp_calcs_2; 859 by option reg; 860 RUN;

NOTE: There were 10 observations read from the data set WORK.COMP_CALCS_2. NOTE: The data set WORK.COMP_CALCS_2 has 10 observations and 62 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.03 seconds cpu time 0.00 seconds

861 862 PROC SORT data=ie_all out=ie_all; 863 by option reg species_econ; 864 RUN;

NOTE: There were 10 observations read from the data set WORK.IE_ALL. NOTE: The data set WORK.IE_ALL has 10 observations and 4 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.03 seconds cpu time 0.00 seconds

865 866 *comp_calcs_4 contains each option/species combination along the final benefits schedule 866 ! (impingement fraction adjusted) 867 for each; 868 DATA comp_calcs_4; 869 LENGTH option $20; 870 Merge ie_all comp_calcs_2; Page 41bt-star_main_03.21.11_log.txt

871 by option reg; 41 The SAS System 12:32 Tuesday, March 22, 2011

872 DROP baseline_frac; 873 RUN;

WARNING: The variable baseline_frac in the DROP, KEEP, or RENAME list has never been referenced. NOTE: There were 10 observations read from the data set WORK.IE_ALL. NOTE: There were 10 observations read from the data set WORK.COMP_CALCS_2. NOTE: The data set WORK.COMP_CALCS_4 has 10 observations and 64 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

874 *the next few stesps breaks comp_calcs_4 into two separate data set for the north_atlantic 874 ! and mid_atlantic which will be merged back together; 875 *comp_calcs_5 is restricted to the mid-atlantic schedule; 876 DATA comp_calcs_5; 877 Set comp_calcs_4; 878 IF reg="northatlantic" then delete; 879 array ma_schd_(60) ma_schd_1-ma_schd_60; 880 array fschd_(60) fschd_1-fschd_60; 881 Do i=1 to 60; 882 ma_schd_(i)=fschd_(i)*benefit_frac; 883 END; 884 DROP fschd_1-fschd_60 benefit_frac reg i; 885 RUN;

NOTE: There were 10 observations read from the data set WORK.COMP_CALCS_4. NOTE: The data set WORK.COMP_CALCS_5 has 5 observations and 62 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

886 *comp_calcs_6 is restricted to the north atlantic schedule; 887 DATA comp_calcs_6; 888 Set comp_calcs_4; 889 IF reg="midatlantic" then delete; 890 array na_schd_(60) na_schd_1-na_schd_60; 891 array fschd_(60) fschd_1-fschd_60; 892 Do i=1 to 60; 893 na_schd_(i)=fschd_(i)*benefit_frac; 894 END; 895 DROP fschd_1-fschd_60 benefit_frac reg i; 896 RUN;

NOTE: There were 10 observations read from the data set WORK.COMP_CALCS_4. NOTE: The data set WORK.COMP_CALCS_6 has 5 observations and 62 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

897 PROC SORT data=comp_calcs_5 out=comp_calcs_5; 898 by option species_econ; 42 The SAS System 12:32 Tuesday, March 22, 2011

Page 42bt-star_main_03.21.11_log.txt

899 RUN;

NOTE: There were 5 observations read from the data set WORK.COMP_CALCS_5. NOTE: The data set WORK.COMP_CALCS_5 has 5 observations and 62 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

900 PROC SORT data=comp_calcs_6 out=comp_calcs_6; 901 by option species_econ; 902 RUN;

NOTE: There were 5 observations read from the data set WORK.COMP_CALCS_6. NOTE: The data set WORK.COMP_CALCS_6 has 5 observations and 62 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.03 seconds cpu time 0.00 seconds

903 904 *comp_calcs_7 merges the adjusted schedules for the mid-atlantic and north atlantic regions 904 ! and calculates 905 a weighted northeast schedule (ne_schd); 906 Data comp_calcs_7; 907 Merge comp_calcs_5 comp_calcs_6; 908 by option species_econ; 909 array ma_schd_(60) ma_schd_1-ma_schd_60; 910 array na_schd_(60) na_schd_1-na_schd_60; 911 array ne_schd_(60) ne_schd_1-ne_schd_60; 912 Do i=1 to 60; 913 IF na_schd_(i)="." then na_schd_(i)=0; 914 ne_schd_(i)=ma_schd_(i) + na_schd_(i); 915 END; 916 RUN;

NOTE: Character values have been converted to numeric values at the places given by: (Line):(Column). 913:17 NOTE: There were 5 observations read from the data set WORK.COMP_CALCS_5. NOTE: There were 5 observations read from the data set WORK.COMP_CALCS_6. NOTE: The data set WORK.COMP_CALCS_7 has 5 observations and 183 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

917 *This step modifies the Reg_compile dataset so it can be merged together with the final 917 ! compliance schedule; 918 DATA Reg_compile_1; 919 Set reg_compile; 920 DROP tmn1 tmn2 tmn3 tmn4 baseline; 921 DROP WTPper; 922 DROP WTPhh_tmn1 WTPhh_tmn2 WTPhh_tmn3 WTPhh_tmn4 WTPhh_baseline; 923 RUN; 43 The SAS System 12:32 Tuesday, March 22, 2011

WARNING: The variable tmn1 in the DROP, KEEP, or RENAME list has never been referenced. Page 43bt-star_main_03.21.11_log.txt

WARNING: The variable tmn2 in the DROP, KEEP, or RENAME list has never been referenced. WARNING: The variable tmn3 in the DROP, KEEP, or RENAME list has never been referenced. WARNING: The variable tmn4 in the DROP, KEEP, or RENAME list has never been referenced. WARNING: The variable baseline in the DROP, KEEP, or RENAME list has never been referenced. NOTE: There were 3 observations read from the data set WORK.REG_COMPILE. NOTE: The data set WORK.REG_COMPILE_1 has 3 observations and 13 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

924 *reg_compile_2 lists each species under evaulation in the north atlantic and mid-atlantic 924 ! regions along with 925 the reg benefits under each option for each species; 926 DATA reg_compile_2; 927 Set reg_compile_1; 928 array WTPreg_(5) WTPreg_tmn1 WTPreg_tmn2 WTPreg_tmn3 WTPreg_tmn4 WTPreg_baseline; 929 array tmn(5) tmn1 tmn2 tmn3 tmn4 baseline; 930 do i=1 to 5; 931 tmn(i)=WTPreg_(i); 932 END; 933 DROP WTPreg_tmn1 WTPreg_tmn2 WTPreg_tmn3 WTPreg_tmn4 WTPreg_baseline NE_baseline NE_tmn1 933 ! NE_tmn2 NE_tmn3 NE_tmn4; 934 DROP regHH; 935 DROP i; 936 RUN;

WARNING: The variable regHH in the DROP, KEEP, or RENAME list has never been referenced. NOTE: There were 3 observations read from the data set WORK.REG_COMPILE_1. NOTE: The data set WORK.REG_COMPILE_2 has 3 observations and 8 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

937 *The following three data/proc steps reorganize the data so that the WTP results can be 937 ! merged with the 938 weighted northeast compliance schedule; 939 PROC sort data=reg_compile_2 out=reg_compile_2; 940 by reg species_econ; 941 RUN;

NOTE: There were 3 observations read from the data set WORK.REG_COMPILE_2. NOTE: The data set WORK.REG_COMPILE_2 has 3 observations and 8 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.03 seconds cpu time 0.00 seconds

942 *reg_compile_2 should be exported and saved as it contains the annual regional benefits for 942 ! each species/option 943 for both the north atlantic and mid-atlantic regions; 44 The SAS System 12:32 Page 44bt-star_main_03.21.11_log.txt

Tuesday, March 22, 2011

944 PROC Transpose data=reg_compile_2 out=reg_compile_3; 945 By reg species_econ; 946 RUN;

NOTE: There were 3 observations read from the data set WORK.REG_COMPILE_2. NOTE: The data set WORK.REG_COMPILE_3 has 18 observations and 5 variables. NOTE: PROCEDURE TRANSPOSE used (Total process time): real time 0.09 seconds cpu time 0.00 seconds

947 DATA reg_compile_3; 948 SET reg_compile_3 (rename=(_NAME_=option COl1=WTPreg)); 949 RUN;

NOTE: There were 18 observations read from the data set WORK.REG_COMPILE_3. NOTE: The data set WORK.REG_COMPILE_3 has 18 observations and 5 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

950 PROC SORT data=reg_compile_3 out=reg_compile_3; 951 by option species_econ reg; 952 RUN;

NOTE: There were 18 observations read from the data set WORK.REG_COMPILE_3. NOTE: The data set WORK.REG_COMPILE_3 has 18 observations and 5 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

953

************************************************************************************

**** 954 *comp_calcs_8 strips out the seperate schedules for the mid-atlantic and north atlantic 955 to estimate the net compliance shedule for northeast waters; 956 Data comp_calcs_8; 957 SET comp_calcs_7; 958 DROP ma_schd_1-ma_schd_60 na_schd_1-na_schd_60 i; 959 RUN;

NOTE: There were 5 observations read from the data set WORK.COMP_CALCS_7. NOTE: The data set WORK.COMP_CALCS_8 has 5 observations and 62 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.01 seconds

960 961 PROC SORT data=comp_calcs_8 out=comp_calcs_8; 962 by option species_econ; 963 RUN;

NOTE: There were 5 observations read from the data set WORK.COMP_CALCS_8. 45 The SAS System 12:32 Tuesday, March 22, 2011

NOTE: The data set WORK.COMP_CALCS_8 has 5 observations and 62 variables. Page 45bt-star_main_03.21.11_log.txt

NOTE: PROCEDURE SORT used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

964 *merge together the annual regional WTP values and the weighted northeast compliance 964 ! schedule for each species/option; 965 DATA reg_compile_4; 966 merge reg_compile_3 comp_calcs_8; 967 by option species_econ; 968 RUN;

WARNING: Multiple lengths were specified for the BY variable option by input data sets. This may

cause unexpected results. WARNING: Multiple lengths were specified for the BY variable species_econ by input data sets. This

may cause unexpected results. NOTE: There were 18 observations read from the data set WORK.REG_COMPILE_3. NOTE: There were 5 observations read from the data set WORK.COMP_CALCS_8. NOTE: The data set WORK.REG_COMPILE_4 has 18 observations and 65 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

969

************************************************************************************

******* 969 ! ****************** 970 *incorporate year into the data set for discounting and implement all discounting and 970 ! annualization calcs; 971 *reg_compile_5 includes raw results prior to sorting; 972 DATA reg_compile_5; 973 Set reg_compile_4; 974 array yr_(60) yr_1-yr_60; 975 array ne_schd_(60) ne_schd_1-ne_schd_60; 976 array df3_(60) df3_1-df3_60; 977 array df7_(60) df7_1-df7_60; 978 array WTP3_(60) WTP3_1-WTP3_60; 979 array WTP7_(60) WTP7_1-WTP7_60; 980 do i=1 to 60; 981 yr_(i)=i; 982 df3_(i)=1/(1.03**yr_(i)); 983 df7_(i)=1/(1.07**yr_(i)); 984 WTP3_(i)=WTPreg*ne_schd_(i)*df3_(i); 985 WTP7_(i)=WTPreg*ne_schd_(i)*df7_(i); 986 END; 987 *calculate NPV for 3% discount rate and 7% discount rate; 988 WTP_NPV_3=sum(of WTP3_1-WTP3_60); 989 WTP_NPV_7=sum(of WTP7_1-WTP7_60); 990 Ann_WTP_3=WTP_NPV_3*&af_3.; 991 Ann_WTP_7=WTP_NPV_7*&af_7.; 992 RUN;

NOTE: Missing values were generated as a result of performing an operation on missing values. Each place is given by: (Number of times) at (Line):(Column). 46 The SAS System 12:32 Tuesday, March 22, 2011

Page 46bt-star_main_03.21.11_log.txt

180 at 984:17 180 at 984:29 180 at 985:17 180 at 985:29 3 at 988:12 3 at 989:12 3 at 990:1 3 at 991:1 NOTE: There were 18 observations read from the data set WORK.REG_COMPILE_4. NOTE: The data set WORK.REG_COMPILE_5 has 18 observations and 370 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

993 Data discounting_out; 994 set reg_compile_5; 995 if reg~="northeast" then delete; 996 if _LABEL_="numHH" then delete; 997 drop _LABEL_ species_econ ne_schd_1-ne_schd_60 yr_1-yr_60 df3_1-df3_60 df7_1-df7_60 i 997 ! WTP_NPV_3 WTP_NPV_7 Ann_WTP_3 Ann_WTP_7; 998 RUN;

NOTE: There were 18 observations read from the data set WORK.REG_COMPILE_5. NOTE: The data set WORK.DISCOUNTING_OUT has 5 observations and 123 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.00 seconds

999 1000 *reg_compile_6 contains organized results with extra information excluded; 1001 Data reg_compile_6; 1002 Set reg_compile_5 (Keep = reg species_econ option WTPreg Ann_WTP_3 Ann_WTP_7); 1003 if option ="numHH" then delete; 1004 if reg ="northeast" then delete; 1005 RUN;

NOTE: There were 18 observations read from the data set WORK.REG_COMPILE_5. NOTE: The data set WORK.REG_COMPILE_6 has 10 observations and 6 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.00 seconds

1006 1007 %export(reg_compile_6,"&outloc.&output",Results_summary);

NOTE: File "h:\ERD\TYCHE Contract\Assignments\WA 2-15 (316b Benefits & SP Study)\Task 5 - Benefits Analytic Support\IE_032111\BT-STAR\Outputs\BT_BenOutputs_22MAR11.xls" will be created if the export process succeeds. NOTE: "RESULTS_SUMMARY" range/sheet was successfully created. NOTE: PROCEDURE EXPORT used (Total process time): real time 0.54 seconds cpu time 0.07 seconds

1008 %export(discounting_out,"&outloc.&output", Discounting_out);

NOTE: "DISCOUNTING_OUT" range/sheet was successfully created. 47 The SAS System 12:32 Tuesday, March 22, 2011

Page 47bt-star_main_03.21.11_log.txt

NOTE: PROCEDURE EXPORT used (Total process time): real time 0.29 seconds cpu time 0.09 seconds

1009

NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414 NOTE: The SAS System used: real time 29.00 seconds cpu time 6.74 seconds

Page 48