Activating two equations in the schedule model and adding the dummy capacity...

Activating two equations in the schedule model and adding the dummy capacity to ensure capacity margin equation feasibility to the results file.

defModels/schedule.gms

@@ -61,8 +61,8 @@ Model schedule /
     q_boundCyclic
 
     // Policy
-*    q_inertiaMin
-*    q_instantaneousShareMax
+    q_inertiaMin
+    q_instantaneousShareMax
     q_constrainedOnlineMultiUnit
 *    q_capacityMargin
 *    q_constrainedCapMultiUnit

defOutput/resultSymbols.inc

@@ -151,6 +151,7 @@ r_gnTotalqGen
 r_gTotalqGen
 r_qResDemand
 r_nTotalqResDemand
+r_qCapacity
 r_solveStatus
 
 // Diagnostics

inc/1d_results.gms

@@ -153,6 +153,7 @@ Parameters
     r_gTotalqGen(inc_dec, grid) "Total dummy energy generation/consumption in g over the simulation (MWh)."
     r_qResDemand(restype, up_down, node, f, t) "Dummy to decrease demand for a reserve (MW)"
     r_nTotalqResDemand(restype, up_down, node) "Total dummy reserve provisions in n over the simulation"
+    r_qCapacity(grid, node, f, t) "Dummy capacity to ensure capacity margin equation feasibility (MW)"
     r_solveStatus(t, solve_info) "Information about the solve"
 
 ; // END PARAMETER DECLARATION

inc/4a_outputVariant.gms

@@ -157,6 +157,10 @@ r_qGen(inc_dec, gn(grid, node), ft_realized(f, t))$[ord(t) > mSettings(mSolve, '
 r_qResDemand(restypeDirectionNode(restype, up_down, node), ft_realized(f, t))$[ord(t) > mSettings(mSolve, 't_start') + mSettings(mSolve, 't_initializationPeriod')]
     = vq_resDemand.l(restype, up_down, node, f, t)
 ;
+// Dummy capacity
+r_qCapacity(gn(grid, node), ft_realized(f, t))$[ord(t) > mSettings(mSolve, 't_start') + mSettings(mSolve, 't_initializationPeriod')]
+    = vq_capacity.l(grid, node, f, t)
+;
 
 * --- Diagnostics Results -----------------------------------------------------