Added a new result vector for "r_qResMissing" to keep track of any instances of "vq_resMissing".

defOutput/resultSymbols.inc

@@ -150,6 +150,7 @@ r_qGen
 r_gnTotalqGen
 r_gTotalqGen
 r_qResDemand
+r_qResMissing
 r_nTotalqResDemand
 r_solveStatus
 

inc/1d_results.gms

@@ -151,7 +151,8 @@ Parameters
     r_qGen(inc_dec, grid, node, f, t) "Dummy energy generation (increase) or consumption (generation decrease) to ensure equation feasibility (MW)"
     r_gnTotalqGen(inc_dec, grid, node) "Total dummy energy generation/consumption in gn over the simulation (MWh)."
     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_qResDemand(restype, up_down, node, f, t) "Dummy to decrease demand for a reserve (MW) before reserve commitment"
+    r_qResMissing(restype, up_down, node, f, t) "Dummy to decrease demand for a reserve (MW) after reserve commitment"
     r_nTotalqResDemand(restype, up_down, node) "Total dummy reserve provisions in n over the simulation"
     r_solveStatus(t, solve_info) "Information about the solve"
 

inc/4a_outputVariant.gms

@@ -129,13 +129,22 @@ r_investTransfer(grid, node, node_, t_invest(t))${ p_gnn(grid, node, node_, 'tra
 * --- Feasibility results -----------------------------------------------------
 
 // Dummy generation & consumption
-r_qGen(inc_dec, gn(grid, node), ft_realized(f, t))$[ord(t) > mSettings(mSolve, 't_start') + mSettings(mSolve, 't_initializationPeriod')]
+r_qGen(inc_dec, gn(grid, node), ft_realized(f, t))
+    ${  ord(t) > mSettings(mSolve, 't_start') + mSettings(mSolve, 't_initializationPeriod')
+        }
     = vq_gen.l(inc_dec, grid, node, f, t)
 ;
 // Dummy reserve demand changes
-r_qResDemand(restypeDirectionNode(restype, up_down, node), ft_realized(f, t))$[ord(t) > mSettings(mSolve, 't_start') + mSettings(mSolve, 't_initializationPeriod')]
+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)
 ;
+r_qResMissing(restypeDirectionNode(restype, up_down, node), ft_realized(f, t))
+    ${  ord(t) > mSettings(mSolve, 't_start') + mSettings(mSolve, 't_initializationPeriod')
+        }
+    = vq_resMissing.l(restype, up_down, node, f, t)
+;
 
 * --- Diagnostics Results -----------------------------------------------------