diff --git a/run_model_on_snellius/exe/STEMMUS_SCOPE b/run_model_on_snellius/exe/STEMMUS_SCOPE
index e85f8560..5bf5133e 100755
Binary files a/run_model_on_snellius/exe/STEMMUS_SCOPE and b/run_model_on_snellius/exe/STEMMUS_SCOPE differ
diff --git a/src/+groundwater/calculateGroundwaterRecharge.m b/src/+groundwater/calculateGroundwaterRecharge.m
index 2435756e..dd12c596 100644
--- a/src/+groundwater/calculateGroundwaterRecharge.m
+++ b/src/+groundwater/calculateGroundwaterRecharge.m
@@ -103,13 +103,14 @@ soilThick               cumulative soil layers thickness (from top to bottom)
         warning('Index of the bottom layer boundary that includes groundwater head ~= index of the recharge layer + 1!');
     end
 
-    % Outputs to be exported in csv
+    % Outputs to be exported in csv or exposed to BMI
     gwfluxes.QLh = QLh_flip(indxRchrg);
     gwfluxes.QLT = QLT_flip(indxRchrg);
     gwfluxes.QLa = QLa_flip(indxRchrg);
     gwfluxes.QVH = QVH_flip(indxRchrg);
     gwfluxes.QVT = QVT_flip(indxRchrg);
     gwfluxes.QVa = QVa_flip(indxRchrg);
+    gwfluxes.indxRchrg = indxRchrg;
     % gwfluxes.recharge_init = recharge_init;
     % gwfluxes.S = S;
     % gwfluxes.sy = sy;
diff --git a/src/STEMMUS_SCOPE.m b/src/STEMMUS_SCOPE.m
index 351cbd24..8c6e155a 100644
--- a/src/STEMMUS_SCOPE.m
+++ b/src/STEMMUS_SCOPE.m
@@ -783,6 +783,7 @@
             end
         else
             gwfluxes.recharge = 0;
+            gwfluxes.indxRchrg = NaN;
         end
 
         % set SoilVariables for the rest of the loop