Update storage_all.jl

Update storage_all.jl with "VirtualChargeDischargeCost setting parameter to account for cost of virtual charging and discharging

src/model/resources/storage/storage_all.jl

@@ -58,6 +58,20 @@ function storage_all!(EP::Model, inputs::Dict, setup::Dict)
 	@expression(EP, eTotalCVarIn, sum(eTotalCVarInT[t] for t in 1:T))
 	add_to_expression!(EP[:eObj], eTotalCVarIn)
 
+	if CapacityReserveMargin > 0
+		#Variable costs of "virtual charging" for technologies "y" during hour "t" in zone "z"
+		@expression(EP, eCVar_in_virtual[y in STOR_ALL,t=1:T], inputs["omega"][t]*setup["VirtualChargeDischargeCost"]*vCAPRES_charge[y,t])
+		@expression(EP, eTotalCVarInT_virtual[t=1:T], sum(eCVar_in_virtual[y,t] for y in STOR_ALL))
+		@expression(EP, eTotalCVarIn_virtual, sum(eTotalCVarInT_virtual[t] for t in 1:T))
+		EP[:eObj] += eTotalCVarIn_virtual
+
+		#Variable costs of "virtual discharging" for technologies "y" during hour "t" in zone "z"
+		@expression(EP, eCVar_out_virtual[y in STOR_ALL,t=1:T], inputs["omega"][t]*setup["VirtualChargeDischargeCost"]*vCAPRES_discharge[y,t])
+		@expression(EP, eTotalCVarOutT_virtual[t=1:T], sum(eCVar_out_virtual[y,t] for y in STOR_ALL))
+		@expression(EP, eTotalCVarOut_virtual, sum(eTotalCVarOutT_virtual[t] for t in 1:T))
+		EP[:eObj] += eTotalCVarOut_virtual
+	end
+
 	## Power Balance Expressions ##
 
 	# Term to represent net dispatch from storage in any period