Fixed a severe bug with probability calculation.

Parameter `p_msft_probability` was not updated after scenario reduction.

CHANGELOG.md

@@ -2,6 +2,9 @@
 All notable changes to this project will be documented in this file.
 
 ## [Unreleased]
+### Fixed
+- Probabilities were not updated after using scenario reduction
+
 ### Added
 - Enable long-term samples that extend several years by using planning horizon 
   which is longer than one scenario (e.g. 3 years). Note: Cannot use all data for 

inc/3b_periodicLoop.gms

@@ -551,15 +551,6 @@ loop(uft_online(unit, f, t)${ p_u_shutdownTimeIntervals(unit) },
         = yes;
 ); // END loop(uf_online)
 
-* -----------------------------------------------------------------------------
-* --- Probabilities -----------------------------------------------------------
-* -----------------------------------------------------------------------------
-
-// Update probabilities
-Option clear = p_msft_probability;
-p_msft_probability(msft(mSolve, s, f, t))
-    = p_mfProbability(mSolve, f) / sum(f_${ft(f_, t)}, p_mfProbability(mSolve, f_)) * p_msProbability(mSolve, s);
-
 * -----------------------------------------------------------------------------
 * --- Displacements for start-up and shutdown decisions -----------------------
 * -----------------------------------------------------------------------------

inc/3c_inputsLoop.gms

@@ -431,6 +431,14 @@ if(active(mSolve, 'scenred'),
     $$include 'inc/scenred.gms'
 );
 
+* --- Update probabilities ----------------------------------------------------
+Option clear = p_msft_probability;
+p_msft_probability(msft(mSolve, s, f, t))
+    = p_mfProbability(mSolve, f)
+        / sum(f_${ft(f_, t)},
+              p_mfProbability(mSolve, f_)) * p_msProbability(mSolve, s);
+
+
 * --- Smooting of stochastic scenarios ----------------------------------------
 $ontext
 First calculate standard deviation for over all samples, then smoothen the scenarios