add rechargeTemp to BMI variables

src/+groundwater/calculateGroundwaterRecharge.m

@@ -22,6 +22,7 @@ gwfluxes                structure that includes the outputs (groundwater recharg
         indxRchrg               index of the soil layer where the recharge is calculated
         recharge                groundwater recharge, which is the upper boundary flux of the top layer of the phreatic aquifer (after the correction of the specific yield)
         recharge_init           upper boundary flux into the moving balancing domain (before the correction of the specific yield)
+        rechargeTemp            temperature of the groundwater recharge
         SS                      specific storage of MODFLOW aquifers
         SY                      large-scale specific yield of the phreatic aquifer
         sy                      small-scale specific yield (dynamically changing water yield) caused by fluctuation of the water table
@@ -97,6 +98,10 @@ soilThick               cumulative soil layers thickness (from top to bottom)
         gwfluxes.recharge = 0;
     end
 
+    % Get rechrage temperature
+    TT_flip(1:1:ModelSettings.NN) = SoilVariables.TT(ModelSettings.NN:-1:1);
+    gwfluxes.rechargeTemp = TT_flip(indxRchrg);
+
     % check 1
     diff = abs(GroundwaterSettings.indxBotmLayer_R - indxRchrg);
     if diff > 1

src/STEMMUS_SCOPE.m

@@ -369,6 +369,9 @@
         if isnan(GroundwaterSettings.tempBotm)
             GroundwaterSettings.tempBotm = SoilVariables.TT(GroundwaterSettings.indxBotmLayer);
         end
+
+        % Initialize groundwater recharge temperature
+        gwfluxes.rechargeTemp = 0;
     end
 
     % Will do one timestep in "update mode", and run until the end if in "full run" mode.