[Bug] Climate Adjustment Factors to Hydropower not Reflected in PRAS data

[jvcarag]

Summary

The hydropower capacity factor (CF) data ingested by PRAS does not include climate-based adjustments applied to hydropower when GSw_ClimateHydro == 1. This is because PRAS uses the monthly-indexed hydcf.csv input file for hydropower CFs, which is a precursor to the rep/stress-period-indexed cf_hyd.csv file that the adjustment factors are applied to:

PRAS ingesting hydcf.csv:

ReEDS/reeds/resource_adequacy/reeds2pras/src/utils/reeds_input_parsing.jl

Lines 340 to 342 in 62f6381

	function get_hydro_data(data::ReEDSdatapaths)
	filepath_cf = joinpath(data.ReEDSfilepath, "inputs_case", "hydcf.csv")
	hydcf = DataFrames.DataFrame(CSV.File(filepath_cf))

Place where climate adjustments to hydropower CF are applied in ReEDS:

ReEDS/reeds/core/solve/2_temporal_params.gms

Lines 462 to 496 in 62f6381

	* Written by reeds/input_processing/hourly_writetimeseries.py
	$onempty
	parameter cf_hyd(i,allszn,r,allt) "--fraction-- hydro capacity factors by season and year"
	/
	$offlisting
	$ondelim
	$include inputs_case%ds%%temporal_inputs%%ds%cf_hyd.csv
	$include inputs_case%ds%stress%stress_year%%ds%cf_hyd.csv
	$offdelim
	$onlisting
	/ ;
	$offempty
	$ifthen.climatehydro %GSw_ClimateHydro% == 1
	* Written by climateprep.py
	table climate_hydro_seasonal(r,allszn,allt) "annual/seasonal nondispatchable hydropower availability"
	$offlisting
	$ondelim
	$include inputs_case%ds%%temporal_inputs%%ds%climate_hydadjsea.csv
	$offdelim
	$onlisting
	;
	* adjust cf_hyd based on annual/seasonal climate multipliers
	* non-dispatchable hydro gets new seasonal profiles as well as annually-varying CF
	* dispatchable hydro keeps the original seasonal profiles; only annual CF changes. Reflects the assumption
	* that reservoirs will be utilized in the same seasonal pattern even if seasonal inflows change.
	cf_hyd(i,szn,r,t)$[hydro_nd(i)$(yeart(t)>=Sw_ClimateStartYear)] =
	sum{allt$att(allt,t), cf_hyd(i,szn,r,t) * climate_hydro_seasonal(r,szn,allt) } ;
	cf_hyd(i,szn,r,t)$[hydro_d(i)$(yeart(t)>=Sw_ClimateStartYear)] =
	sum{allt$att(allt,t), cf_hyd(i,szn,r,t) * climate_hydro_annual(r,allt) } ;
	$endif.climatehydro

Suggested Fix

Since climate adjustment factors are first read into the model at monthly resolution, they can be applied to hydcf.csv in the hydcf.py input processing script. PR #129 will address this Issue.