Itemize costs incurred by the objective function. Use real dollars instead of NPV when reporting future energy costs per MWh. Update the output files for the two examples that include all output files.

Author: josiahjohnston

examples/copperplate0/outputs/costs_itemized.csv

@@ -0,0 +1,4 @@
+PERIOD,Component,AnnualCost_NPV,AnnualCost_Real,Component_type
+2020,TotalGenFixedCosts,11213737.1287,1601083.86234,annual
+2020,GenVariableOMCostsInTP,1608195.12658,229616.160528,timepoint
+2020,FuelCostsPerTP,5682691.43038,811367.828532,timepoint

examples/copperplate0/outputs/electricity_cost.csv

@@ -1,2 +1,2 @@
-PERIOD,EnergyCost_per_MWh,SystemCostPerPeriod
-2020,49.6695083562,18504623.6857
+PERIOD,EnergyCostReal_per_MWh,SystemCostPerPeriod_NPV,SystemCostPerPeriod_Real,SystemDemand_MWh
+2020,7.09175249668,18504623.6857,2642067.8514,372555.0

examples/copperplate0/outputs/results.pickle

@@ -0,0 +1,5 @@
+PERIOD,Component,AnnualCost_NPV,AnnualCost_Real,Component_type
+2020,TotalGenFixedCosts,16582355.87,2367608.77112,annual
+2020,StorageEnergyInstallCosts,35397.2613454,5053.97828224,annual
+2020,GenVariableOMCostsInTP,7867155.78173,1123263.01393,timepoint
+2020,FuelCostsPerTP,0.0,0.0,timepoint

examples/storage/outputs/costs_itemized.csv

@@ -0,0 +1,2 @@
+PERIOD,EnergyCostReal_per_MWh,SystemCostPerPeriod_NPV,SystemCostPerPeriod_Real,SystemDemand_MWh
+2020,9.3836501009,24484908.913,3495925.76334,372555.0

examples/storage/outputs/electricity_cost.csv

@@ -308,17 +308,53 @@ def load_inputs(mod, switch_data, inputs_dir):
     switch_data.load(filename=os.path.join(inputs_dir, 'financials.dat'))
 
 def post_solve(instance, outdir):
-    mod = instance
+    m = instance
+    # Overall electricity costs
     normalized_dat = [
         {
         	"PERIOD": p,
-        	"SystemCostPerPeriod": value(mod.SystemCostPerPeriod[p]),
-        	"EnergyCost_per_MWh": value(
-        	    mod.SystemCostPerPeriod[p] / 
-        	    sum(mod.zone_total_demand_in_period_mwh[z,p] for z in mod.LOAD_ZONES)
-        	)
-        } for p in mod.PERIODS
+        	"SystemCostPerPeriod_NPV": value(m.SystemCostPerPeriod[p]),
+        	"SystemCostPerPeriod_Real": value(
+        	    m.SystemCostPerPeriod[p] / m.bring_annual_costs_to_base_year[p]
+        	),
+        	"EnergyCostReal_per_MWh": value(
+        	    m.SystemCostPerPeriod[p] / m.bring_annual_costs_to_base_year[p] /
+        	    sum(m.zone_total_demand_in_period_mwh[z,p] for z in m.LOAD_ZONES)
+        	),
+        	"SystemDemand_MWh": value(sum(
+        	    m.zone_total_demand_in_period_mwh[z,p] for z in m.LOAD_ZONES
+        	))
+        } for p in m.PERIODS
     ]
     df = pd.DataFrame(normalized_dat)
     df.set_index(["PERIOD"], inplace=True)
     df.to_csv(os.path.join(outdir, "electricity_cost.csv"))
+    # Itemized annual costs
+    annualized_costs = [
+        {
+        	"PERIOD": p,
+        	"Component": annual_cost,
+        	"Component_type": "annual",
+        	"AnnualCost_NPV": value(
+        	    getattr(m, annual_cost)[p] * m.bring_annual_costs_to_base_year[p]
+        	),
+        	"AnnualCost_Real": value(getattr(m, annual_cost)[p])
+        } for p in m.PERIODS for annual_cost in m.Cost_Components_Per_Period
+    ] + [
+        {
+        	"PERIOD": p,
+        	"Component": tp_cost,
+        	"Component_type": "timepoint",
+        	"AnnualCost_NPV": value(sum(
+        	    getattr(m, tp_cost)[t] * m.tp_weight_in_year[t]
+        	    for t in m.TPS_IN_PERIOD[p]
+        	) * m.bring_annual_costs_to_base_year[p]),
+        	"AnnualCost_Real": value(sum(
+        	    getattr(m, tp_cost)[t] * m.tp_weight_in_year[t]
+        	    for t in m.TPS_IN_PERIOD[p]
+        	))
+        } for p in m.PERIODS for tp_cost in m.Cost_Components_Per_TP
+    ]
+    df = pd.DataFrame(annualized_costs)
+    df.set_index(["PERIOD", "Component"], inplace=True)
+    df.to_csv(os.path.join(outdir, "costs_itemized.csv"))