Modify write outputs scripts for CCS

src/write_outputs/energy_share_requirement/write_esr_revenue.jl

@@ -61,5 +61,4 @@ function write_esr_revenue(path::AbstractString, inputs::Dict, setup::Dict, dfPo
 	dfESRRev.Total = sum(eachcol(dfESRRev[:, 6:nESR + 5]))
 	CSV.write(joinpath(path, "ESR_Revenue.csv"), dfESRRev)
 	return dfESRRev
-end
-
+end

src/write_outputs/write_co2.jl

@@ -37,14 +37,15 @@ end
 function write_co2_capture_plant(path::AbstractString, inputs::Dict, setup::Dict, EP::Model)
     dfGen = inputs["dfGen"]
     G = inputs["G"]     # Number of resources (generators, storage, DR, and DERs)
+    CCS = inputs["CCS"]
     T = inputs["T"]     # Number of time steps (hours)
     Z = inputs["Z"]     # Number of zones
 
-    dfCapturedEmissions_plant = DataFrame(Resource=inputs["RESOURCES"], Zone=dfGen[!, :Zone], AnnualSum=zeros(G))
-    if any(dfGen.CO2_Capture_Fraction .!= 0)
+    dfCapturedEmissions_plant = DataFrame(Resource=inputs["RESOURCES"][CCS], Zone=dfGen[CCS, :Zone], AnnualSum=zeros(length(CCS)))
+    if !isempty(CCS)
         # Captured CO2 emissions by plant
-        emissions_captured_plant = zeros(G, T)
-        emissions_captured_plant = (value.(EP[:eEmissionsCaptureByPlant]))
+        emissions_captured_plant = zeros(length(CCS), T)
+        emissions_captured_plant = (value.(EP[:eEmissionsCaptureByPlant]).data)
         if setup["ParameterScale"] == 1
             emissions_captured_plant *= ModelScalingFactor
         end

src/write_outputs/write_costs.jl

@@ -108,6 +108,7 @@ function write_costs(path::AbstractString, inputs::Dict, setup::Dict, EP::Model)
 		FLEX_ZONE = intersect(inputs["FLEX"], Y_ZONE)
 		COMMIT_ZONE = intersect(inputs["COMMIT"], Y_ZONE)
 		ELECTROLYZERS_ZONE = intersect(inputs["ELECTROLYZER"], Y_ZONE)
+		CCS_ZONE = intersect(inputs["CCS"], Y_ZONE)
 
 		eCFix = sum(value.(EP[:eCFix][Y_ZONE]))
 		tempCFix += eCFix
@@ -218,8 +219,9 @@ function write_costs(path::AbstractString, inputs::Dict, setup::Dict, EP::Model)
 		tempCNSE = sum(value.(EP[:eCNSE][:,:,z]))
 		tempCTotal += tempCNSE
 
-		if any(dfGen.CO2_Capture_Fraction .!=0)
-			tempCCO2 = sum(value.(EP[:ePlantCCO2Sequestration][Y_ZONE,:]))
+		# if any(dfGen.CO2_Capture_Fraction .!=0)
+		if !isempty(CCS_ZONE) 
+			tempCCO2 = sum(value.(EP[:ePlantCCO2Sequestration][CCS_ZONE]))
 			tempCTotal += tempCCO2		
 		end
 

src/write_outputs/write_net_revenue.jl

@@ -12,6 +12,7 @@ function write_net_revenue(path::AbstractString, inputs::Dict, setup::Dict, EP::
 	STOR_ALL = inputs["STOR_ALL"]
 	VRE_STOR = inputs["VRE_STOR"]
 	dfVRE_STOR = inputs["dfVRE_STOR"]
+	CCS = inputs["CCS"]
 	if !isempty(VRE_STOR)
 		VRE_STOR_LENGTH = size(inputs["VRE_STOR"])[1]
 		SOLAR = inputs["VS_SOLAR"]
@@ -124,7 +125,8 @@ function write_net_revenue(path::AbstractString, inputs::Dict, setup::Dict, EP::
 	# Add CO2 releated sequestration cost or credit (e.g. 45 Q) to the dataframe
 	dfNetRevenue.CO2SequestrationCost = zeros(nrow(dfNetRevenue))
 	if any(dfGen.CO2_Capture_Fraction .!= 0)
-		dfNetRevenue.CO2SequestrationCost .= value.(EP[:ePlantCCO2Sequestration])
+		dfNetRevenue.CO2SequestrationCost = zeros(G)
+		dfNetRevenue[CCS, :CO2SequestrationCost] = value.(EP[:ePlantCCO2Sequestration]).data
  	end
 	if setup["ParameterScale"] == 1
 		dfNetRevenue.CO2SequestrationCost *= ModelScalingFactor^2 # converting Million US$ to US$
@@ -196,7 +198,7 @@ function write_net_revenue(path::AbstractString, inputs::Dict, setup::Dict, EP::
 					  .+ dfNetRevenue.Charge_cost 
 					  .+ dfNetRevenue.EmissionsCost 
 					  .+ dfNetRevenue.StartCost 
-					  .+dfNetRevenue.CO2SequestrationCost)
+					  .+ dfNetRevenue.CO2SequestrationCost)
 	dfNetRevenue.Profit = dfNetRevenue.Revenue .- dfNetRevenue.Cost
 
 	CSV.write(joinpath(path, "NetRevenue.csv"), dfNetRevenue)