wolfhece 2.2.38__py3-none-any.whl → 2.2.39__py3-none-any.whl

wolfhece/hydrology/SubBasin.py

@@ -28,6 +28,7 @@ from dbfread import DBF
 import ctypes as ct
 import pandas as pd
 from pathlib import Path
+import ctypes as ct
 
 from ..PyTranslate import _
 from . import plot_hydrology as ph
@@ -80,7 +81,8 @@ class SubBasin:
     #                myHydro[i][1] : interflow
     #                myHydro[i][2] : baseflow
     #  @unit $[\si{m}^3/\si{s}}]$
-    myHydro:np.ndarray              # [m^3/s] Hydro of the subbasin only -> UPDATE: [mm/h]!!!!!
+    myHydro:np.ndarray              # [mm/h] Fortran-like array that is pointed to by ptr_q_all with Fortran DLL (old version in [m^3/s])
+    iv_saved:np.ndarray             # [mm/h] or [mm] or [-] depending on the variable
 
     intletsObj:dict
     inlets:np.ndarray
@@ -123,7 +125,8 @@ class SubBasin:
     peakTime:datetime.datetime      # datetime of the peak for total outflow -> time delay is already applied
 
     # Hello!
-    ptr_q_all:np.ndarray            # Pointer to Fortran for the result matrix
+    ptr_q_all:ct._Pointer|None      # Pointer to Fortran for the result matrix
+    ptr_iv_saved:ct._Pointer|None   # Pointer to Fortran for the saved internal variables
 
     transferParam:Wolf_Param        # Parameter file with the type of transfer applied -> e.g. timeDelay
 
@@ -171,7 +174,8 @@ class SubBasin:
         #                myHydro[i][1] : interflow
         #                myHydro[i][2] : baseflow
         #  @unit $[\si{m}^3/\si{s}}]$
-        self.myHydro = []                           # [m^3/s] Hydro of the subbasin only -> UPDATE: [mm/h]!!!!!
+        self.myHydro = np.empty((0,0),dtype=ct.c_double, order='F')     # [mm/h] hydrograph of the subbasin (old version in [m^3/s]) (potentially with Fortran pointer to ptr_q_all)
+        self.saved_iv = np.empty((0,0),dtype=ct.c_double, order='F')    # [mm/h] or [mm] or [-] depending on the variable (potentially with Fortran pointer to ptr_iv_saved)
 
         self.intletsObj = {}
         self.inlets = None
@@ -216,6 +220,7 @@ class SubBasin:
 
         # Hello!
         self.ptr_q_all = None
+        self.ptr_iv_saved = None
 
         # version of the Python WOLFHydro
         self._version = version
@@ -225,10 +230,12 @@ class SubBasin:
         # Load all the hydrographs of the sub-basins
         if(self.model==cst.measures or self.model==cst.compare_opti):
             readHydro=False
+
         # Get the main characteristics of the subbasin. If the hydro can be read, so be the main characteristics
         if(readHydro):
             self.read_myMainCharacteristics(_workingDir)
             self.init_timeDelay()
+
         if(readHydro):
             timeTest, self.myHydro = self.get_hydro(self.iDSorted, _workingDir, tzDelta=datetime.timedelta(hours=self.tz))
             if(self.time is None):
@@ -1012,6 +1019,10 @@ class SubBasin:
                     print("ERROR: Expecting only 1 file name for UH model!")
                     sys.exit()
 
+            if not Path(fileName).exists():
+                logging.error("The file " + fileName + " does not exist!")
+                return None, None
+
             with open(fileName, newline = '') as fileID:
                 data_reader = csv.reader(fileID, delimiter='\t')
                 list_data = []
@@ -1244,7 +1255,7 @@ class SubBasin:
             # Real hydrograph
             self._outFlow[nameOutFlow]["Net"]  = self.inlets + tmpHydro
         elif(self.model==cst.tom_VHM or self.model==cst.tom_2layers_linIF or self.model==cst.tom_2layers_UH or
-             self.model==cst.tom_HBV, self.model==cst.tom_SAC_SMA, self.model==cst.tom_NAM):
+             self.model==cst.tom_HBV or self.model==cst.tom_SAC_SMA or self.model==cst.tom_NAM or self.model==cst.tom_SAC_SMA_LROF):
             tmpOutFlow = np.sum(self.myHydro,1)*self.surfaceDrained/3.6
 
             # Raw hydrograph
@@ -1260,7 +1271,6 @@ class SubBasin:
             #     self.outFlow[i] = self.inlets[i] + np.sum(self.myHydro[i])
 
 
-
     def sum_inlets(self):
         """ Sum all the inlet hydrographs of a subbasin. Return an array of zeros otherwise.
 
@@ -1281,10 +1291,11 @@ class SubBasin:
                 self.inlets += curObj.get_outFlow(typeOutFlow="Net", whichOutFlow=nameInlet, lag=deltaTr)
                 self.inletsRaw += curObj.get_outFlow(typeOutFlow="Raw", whichOutFlow=nameInlet, lag=deltaTr)
         else:
-            self.inlets = np.zeros(len(self.myHydro), dtype=ct.c_double, order='F')
-            self.inletsRaw = np.zeros(len(self.myHydro), dtype=ct.c_double, order='F')
-
-
+            if self.myHydro is not None:
+                self.inlets = np.zeros(len(self.myHydro), dtype=ct.c_double, order='F')
+                self.inletsRaw = np.zeros(len(self.myHydro), dtype=ct.c_double, order='F')
+            else:
+                logging.error("ERROR: No hydrograph computed for this subbasin!")
 
 
     def add_rain(self, workingDir, fileName=None, tzDelta=datetime.timedelta(hours=0)):
@@ -1298,6 +1309,10 @@ class SubBasin:
         if(fileName is None):
             fileName = 'Subbasin_'+str(self.iDSorted)+'/simul_lumped_rain.txt'
 
+        if not Path(os.path.join(workingDir,fileName)).exists():
+            logging.error("The file " + os.path.join(workingDir,fileName) + " does not exist!")
+            return None
+
         with open(os.path.join(workingDir,fileName), newline = '') as fileID2:
             data_reader = csv.reader(fileID2, delimiter='\t')
             list_data = []
@@ -1552,7 +1567,7 @@ class SubBasin:
         """
 
         if(self.model==cst.tom_UH or self.model==cst.tom_2layers_linIF or self.model==cst.tom_2layers_UH or self.model==cst.tom_GR4 or
-           self.model==cst.tom_HBV, self.model==cst.tom_SAC_SMA, self.model==cst.tom_NAM):
+           self.model==cst.tom_HBV or self.model==cst.tom_SAC_SMA or self.model==cst.tom_NAM or self.model==cst.tom_SAC_SMA_LROF):
             # x = self.time/3600.0
             if(axis=="Hours"):
                 x = (self.time[:]-self.time[0])/3600.0
@@ -1771,7 +1786,8 @@ class SubBasin:
                 y[:,i] = self.myHydro[:,i]*self.surfaceDrained/3.6
                 tmpSum += self.myHydro[:,i]*self.surfaceDrained/3.6
             y[:,-1] = tmpSum
-        elif(self.model==cst.tom_GR4 or self.model==cst.tom_HBV, self.model==cst.tom_SAC_SMA, self.model==cst.tom_NAM):
+        elif(self.model==cst.tom_GR4 or self.model==cst.tom_HBV or
+             self.model==cst.tom_SAC_SMA or self.model==cst.tom_NAM or self.model==cst.tom_SAC_SMA_LROF):
             y[:,0] = self.myHydro[:,0]*self.surfaceDrained/3.6
         else:
             print("ERROR: this model was not implemented yet!")
@@ -1984,7 +2000,7 @@ class SubBasin:
         # y = np.zeros((len(self.outFlow)-lastElement,nbElements))
         y = []
         if(self.model==cst.tom_UH or self.model==cst.tom_2layers_linIF or self.model==cst.tom_2layers_UH or
-           cst.tom_GR4 or cst.tom_HBV or cst.tom_SAC_SMA or cst.tom_NAM):
+           cst.tom_GR4 or cst.tom_HBV or cst.tom_SAC_SMA or cst.tom_NAM or self.model==cst.tom_SAC_SMA_LROF):
             if(withDelay):
                 # y[:,0] = self.outFlow[:-1]
                 y.append(myOutFlow[:])
@@ -2053,7 +2069,7 @@ class SubBasin:
             # y_titles.append("Total")
             y_titles.append(_('Débits simulés'))
         elif(self.model==cst.tom_UH or self.model==cst.tom_2layers_linIF or self.model==cst.tom_2layers_UH or
-             self.model==cst.tom_HBV or self.model==cst.tom_SAC_SMA or self.model==cst.tom_NAM):
+             self.model==cst.tom_HBV or self.model==cst.tom_SAC_SMA or self.model==cst.tom_NAM, self.model==cst.tom_SAC_SMA_LROF):
             if(ylabel==[]):
                 y_titles.append(_('Débits simulés'))
                 # y_titles.append('Avec reconstruction Qout B. Vesdre')
@@ -2200,127 +2216,138 @@ class SubBasin:
             fileName = fileNameList[0]
             fileNameHydro = fileNameList[1]
 
-        with open(os.path.join(workingDir,fileName), newline = '') as fileID2:
-            data_reader = csv.reader(fileID2, delimiter='\t')
-            list_data = []
-            i=0
-            for raw in data_reader:
-                if i>0:
-                    list_data.append(raw)
-                i += 1
+        filename = Path(workingDir) / fileName
+        if filename.exists():
 
-        tmp = ''
-        for i in range(len(list_data)):
-            if(list_data[i][0][:4]=="Area"):
-                tmp = list_data[i][0].split()
-                self.surfaceDrained = float(tmp[1].split("[")[0])
-                self.mainCharactDict["Area"] = {}
-                self.mainCharactDict["Area"]["value"] = self.surfaceDrained
-                self.mainCharactDict["Area"]["unit"] = "["+tmp[1].split("[")[1]
-            elif(list_data[i][0][:9]=="Perimeter"):
-                tmp = list_data[i][0].split()
-                self.mainCharactDict["Perimeter"] = {}
-                self.mainCharactDict["Perimeter"]["value"] = float(tmp[1].split("[")[0])
-                self.mainCharactDict["Perimeter"]["unit"] = "["+tmp[1].split("[")[1]
-            elif(list_data[i][0][:13]=="Average slope"):
-                tmp = list_data[i][0].split()
-                self.mainCharactDict["Average slope"] = {}
-                self.mainCharactDict["Average slope"]["value"] = float(tmp[2].split("[")[0])
-                self.mainCharactDict["Average slope"]["unit"] = "["+tmp[2].split("[")[1]
-            elif(list_data[i][0][:35]=="Compactness coefficient (Gravelius)"):
-                tmp = list_data[i][0].split()
-                self.mainCharactDict["Compactness coefficient (Gravelius)"] = {}
-                self.mainCharactDict["Compactness coefficient (Gravelius)"]["value"] = float(tmp[3].split("[")[0])
-                self.mainCharactDict["Compactness coefficient (Gravelius)"]["unit"] = "[-]"
-            elif(list_data[i][0][:12]=="Max lag time"):
-                tmp = list_data[i][0].split()
-                self.mainCharactDict["Max lag time"] = {}
-                self.mainCharactDict["Max lag time"]["value"] = float(tmp[3].split("[")[0])
-                self.mainCharactDict["Max lag time"]["unit"] = "["+tmp[3].split("[")[1]
-            elif(list_data[i][0][:12] == "Min lag time"):
-                tmp = list_data[i][0].split()
-                self.timeDelay = float(tmp[3].split("[")[0])
-                self.mainCharactDict["Min lag time"] = {}
-                self.mainCharactDict["Min lag time"]["value"] = float(tmp[3].split("[")[0])
-                self.mainCharactDict["Min lag time"]["unit"] = "["+tmp[3].split("[")[1]
-            elif(list_data[i][0][:12]=="Max altitude"):
-                tmp = list_data[i][0].split()
-                self.mainCharactDict["Max altitude"] = {}
-                self.mainCharactDict["Max altitude"]["value"] = float(tmp[2].split("[")[0])
-                self.mainCharactDict["Max altitude"]["unit"] = "["+tmp[2].split("[")[1]
-            elif(list_data[i][0][:12]=="Min altitude"):
-                tmp = list_data[i][0].split()
-                self.mainCharactDict["Min altitude"] = {}
-                self.mainCharactDict["Min altitude"]["value"] = float(tmp[2].split("[")[0])
-                self.mainCharactDict["Min altitude"]["unit"] = "["+tmp[2].split("[")[1]
-            elif(list_data[i][0][:21]=="Fraction of landuse n"):
-                tmp = list_data[i][0].split()
-                self.mainCharactDict["Fraction of landuse n "+tmp[4]] = {}
-                self.mainCharactDict["Fraction of landuse n "+tmp[4]]["value"] = float(tmp[5].split("[")[0])
-                self.mainCharactDict["Fraction of landuse n "+tmp[4]]["unit"] = "["+tmp[5].split("[")[1]
+            with open(filename, newline = '') as fileID2:
+                data_reader = csv.reader(fileID2, delimiter='\t')
+                list_data = []
+                i=0
+                for raw in data_reader:
+                    if i>0:
+                        list_data.append(raw)
+                    i += 1
+
+            tmp = ''
+            for i in range(len(list_data)):
+                if(list_data[i][0][:4]=="Area"):
+                    tmp = list_data[i][0].split()
+                    self.surfaceDrained = float(tmp[1].split("[")[0])
+                    self.mainCharactDict["Area"] = {}
+                    self.mainCharactDict["Area"]["value"] = self.surfaceDrained
+                    self.mainCharactDict["Area"]["unit"] = "["+tmp[1].split("[")[1]
+                elif(list_data[i][0][:9]=="Perimeter"):
+                    tmp = list_data[i][0].split()
+                    self.mainCharactDict["Perimeter"] = {}
+                    self.mainCharactDict["Perimeter"]["value"] = float(tmp[1].split("[")[0])
+                    self.mainCharactDict["Perimeter"]["unit"] = "["+tmp[1].split("[")[1]
+                elif(list_data[i][0][:13]=="Average slope"):
+                    tmp = list_data[i][0].split()
+                    self.mainCharactDict["Average slope"] = {}
+                    self.mainCharactDict["Average slope"]["value"] = float(tmp[2].split("[")[0])
+                    self.mainCharactDict["Average slope"]["unit"] = "["+tmp[2].split("[")[1]
+                elif(list_data[i][0][:35]=="Compactness coefficient (Gravelius)"):
+                    tmp = list_data[i][0].split()
+                    self.mainCharactDict["Compactness coefficient (Gravelius)"] = {}
+                    self.mainCharactDict["Compactness coefficient (Gravelius)"]["value"] = float(tmp[3].split("[")[0])
+                    self.mainCharactDict["Compactness coefficient (Gravelius)"]["unit"] = "[-]"
+                elif(list_data[i][0][:12]=="Max lag time"):
+                    tmp = list_data[i][0].split()
+                    self.mainCharactDict["Max lag time"] = {}
+                    self.mainCharactDict["Max lag time"]["value"] = float(tmp[3].split("[")[0])
+                    self.mainCharactDict["Max lag time"]["unit"] = "["+tmp[3].split("[")[1]
+                elif(list_data[i][0][:12] == "Min lag time"):
+                    tmp = list_data[i][0].split()
+                    self.timeDelay = float(tmp[3].split("[")[0])
+                    self.mainCharactDict["Min lag time"] = {}
+                    self.mainCharactDict["Min lag time"]["value"] = float(tmp[3].split("[")[0])
+                    self.mainCharactDict["Min lag time"]["unit"] = "["+tmp[3].split("[")[1]
+                elif(list_data[i][0][:12]=="Max altitude"):
+                    tmp = list_data[i][0].split()
+                    self.mainCharactDict["Max altitude"] = {}
+                    self.mainCharactDict["Max altitude"]["value"] = float(tmp[2].split("[")[0])
+                    self.mainCharactDict["Max altitude"]["unit"] = "["+tmp[2].split("[")[1]
+                elif(list_data[i][0][:12]=="Min altitude"):
+                    tmp = list_data[i][0].split()
+                    self.mainCharactDict["Min altitude"] = {}
+                    self.mainCharactDict["Min altitude"]["value"] = float(tmp[2].split("[")[0])
+                    self.mainCharactDict["Min altitude"]["unit"] = "["+tmp[2].split("[")[1]
+                elif(list_data[i][0][:21]=="Fraction of landuse n"):
+                    tmp = list_data[i][0].split()
+                    self.mainCharactDict["Fraction of landuse n "+tmp[4]] = {}
+                    self.mainCharactDict["Fraction of landuse n "+tmp[4]]["value"] = float(tmp[5].split("[")[0])
+                    self.mainCharactDict["Fraction of landuse n "+tmp[4]]["unit"] = "["+tmp[5].split("[")[1]
+        else:
+            logging.error(f"File {filename} does not exist. Cannot read averaged characteristics.")
 
 
         data_reader = None
         list_data = []
         # fileName  = "/Subbasin_" + str(self.iDSorted) + "/" + "simul_subbasin.avrg_caractWholeHydroSubBasin"
-        with open(os.path.join(workingDir,fileNameHydro), newline = '') as fileID2:
-            data_reader = csv.reader(fileID2, delimiter='\t')
-            list_data = []
-            i=0
-            for raw in data_reader:
-                if i>0:
-                    list_data.append(raw)
-                i += 1
 
-        tmp = ''
-        for i in range(len(list_data)):
-            if(list_data[i][0][:4]=="Area"):
-                tmp = list_data[i][0].split()
-                self.surfaceDrainedHydro = float(tmp[1].split("[")[0])
-                self.mainCharactDictWholeHydro["Area"] = {}
-                self.mainCharactDictWholeHydro["Area"]["value"] = self.surfaceDrainedHydro
-                self.mainCharactDictWholeHydro["Area"]["unit"] = "["+tmp[1].split("[")[1]
-            elif(list_data[i][0][:9]=="Perimeter"):
-                tmp = list_data[i][0].split()
-                self.mainCharactDictWholeHydro["Perimeter"] = {}
-                self.mainCharactDictWholeHydro["Perimeter"]["value"] = float(tmp[1].split("[")[0])
-                self.mainCharactDictWholeHydro["Perimeter"]["unit"] = "["+tmp[1].split("[")[1]
-            elif(list_data[i][0][:13]=="Average slope"):
-                tmp = list_data[i][0].split()
-                self.mainCharactDictWholeHydro["Average slope"] = {}
-                self.mainCharactDictWholeHydro["Average slope"]["value"] = float(tmp[2].split("[")[0])
-                self.mainCharactDictWholeHydro["Average slope"]["unit"] = "["+tmp[2].split("[")[1]
-            elif(list_data[i][0][:35]=="Compactness coefficient (Gravelius)"):
-                tmp = list_data[i][0].split()
-                self.mainCharactDictWholeHydro["Compactness coefficient (Gravelius)"] = {}
-                self.mainCharactDictWholeHydro["Compactness coefficient (Gravelius)"]["value"] = float(tmp[3].split("[")[0])
-                self.mainCharactDictWholeHydro["Compactness coefficient (Gravelius)"]["unit"] = "[-]"
-            elif(list_data[i][0][:12]=="Max lag time"):
-                tmp = list_data[i][0].split()
-                self.mainCharactDictWholeHydro["Max lag time"] = {}
-                self.mainCharactDictWholeHydro["Max lag time"]["value"] = float(tmp[3].split("[")[0])
-                self.mainCharactDictWholeHydro["Max lag time"]["unit"] = "["+tmp[3].split("[")[1]
-            elif(list_data[i][0][:12] == "Min lag time"):
-                tmp = list_data[i][0].split()
-                self.timeDelay = float(tmp[3].split("[")[0])
-                self.mainCharactDictWholeHydro["Min lag time"] = {}
-                self.mainCharactDictWholeHydro["Min lag time"]["value"] = float(tmp[3].split("[")[0])
-                self.mainCharactDictWholeHydro["Min lag time"]["unit"] = "["+tmp[3].split("[")[1]
-            elif(list_data[i][0][:12]=="Max altitude"):
-                tmp = list_data[i][0].split()
-                self.mainCharactDictWholeHydro["Max altitude"] = {}
-                self.mainCharactDictWholeHydro["Max altitude"]["value"] = float(tmp[2].split("[")[0])
-                self.mainCharactDictWholeHydro["Max altitude"]["unit"] = "["+tmp[2].split("[")[1]
-            elif(list_data[i][0][:12]=="Min altitude"):
-                tmp = list_data[i][0].split()
-                self.mainCharactDictWholeHydro["Min altitude"] = {}
-                self.mainCharactDictWholeHydro["Min altitude"]["value"] = float(tmp[2].split("[")[0])
-                self.mainCharactDictWholeHydro["Min altitude"]["unit"] = "["+tmp[2].split("[")[1]
-            elif(list_data[i][0][:21]=="Fraction of landuse n"):
-                tmp = list_data[i][0].split()
-                self.mainCharactDictWholeHydro["Fraction of landuse n "+tmp[4]] = {}
-                self.mainCharactDictWholeHydro["Fraction of landuse n "+tmp[4]]["value"] = float(tmp[5].split("[")[0])
-                self.mainCharactDictWholeHydro["Fraction of landuse n "+tmp[4]]["unit"] = "["+tmp[5].split("[")[1]
+        fileNameHydro = Path(workingDir) / fileNameHydro
+        if fileNameHydro.exists():
+            with open(fileNameHydro, newline = '') as fileID2:
+                data_reader = csv.reader(fileID2, delimiter='\t')
+                list_data = []
+                i=0
+                for raw in data_reader:
+                    if i>0:
+                        list_data.append(raw)
+                    i += 1
+
+            tmp = ''
+            for i in range(len(list_data)):
+                if(list_data[i][0][:4]=="Area"):
+                    tmp = list_data[i][0].split()
+                    self.surfaceDrainedHydro = float(tmp[1].split("[")[0])
+                    self.mainCharactDictWholeHydro["Area"] = {}
+                    self.mainCharactDictWholeHydro["Area"]["value"] = self.surfaceDrainedHydro
+                    self.mainCharactDictWholeHydro["Area"]["unit"] = "["+tmp[1].split("[")[1]
+                elif(list_data[i][0][:9]=="Perimeter"):
+                    tmp = list_data[i][0].split()
+                    self.mainCharactDictWholeHydro["Perimeter"] = {}
+                    self.mainCharactDictWholeHydro["Perimeter"]["value"] = float(tmp[1].split("[")[0])
+                    self.mainCharactDictWholeHydro["Perimeter"]["unit"] = "["+tmp[1].split("[")[1]
+                elif(list_data[i][0][:13]=="Average slope"):
+                    tmp = list_data[i][0].split()
+                    self.mainCharactDictWholeHydro["Average slope"] = {}
+                    self.mainCharactDictWholeHydro["Average slope"]["value"] = float(tmp[2].split("[")[0])
+                    self.mainCharactDictWholeHydro["Average slope"]["unit"] = "["+tmp[2].split("[")[1]
+                elif(list_data[i][0][:35]=="Compactness coefficient (Gravelius)"):
+                    tmp = list_data[i][0].split()
+                    self.mainCharactDictWholeHydro["Compactness coefficient (Gravelius)"] = {}
+                    self.mainCharactDictWholeHydro["Compactness coefficient (Gravelius)"]["value"] = float(tmp[3].split("[")[0])
+                    self.mainCharactDictWholeHydro["Compactness coefficient (Gravelius)"]["unit"] = "[-]"
+                elif(list_data[i][0][:12]=="Max lag time"):
+                    tmp = list_data[i][0].split()
+                    self.mainCharactDictWholeHydro["Max lag time"] = {}
+                    self.mainCharactDictWholeHydro["Max lag time"]["value"] = float(tmp[3].split("[")[0])
+                    self.mainCharactDictWholeHydro["Max lag time"]["unit"] = "["+tmp[3].split("[")[1]
+                elif(list_data[i][0][:12] == "Min lag time"):
+                    tmp = list_data[i][0].split()
+                    self.timeDelay = float(tmp[3].split("[")[0])
+                    self.mainCharactDictWholeHydro["Min lag time"] = {}
+                    self.mainCharactDictWholeHydro["Min lag time"]["value"] = float(tmp[3].split("[")[0])
+                    self.mainCharactDictWholeHydro["Min lag time"]["unit"] = "["+tmp[3].split("[")[1]
+                elif(list_data[i][0][:12]=="Max altitude"):
+                    tmp = list_data[i][0].split()
+                    self.mainCharactDictWholeHydro["Max altitude"] = {}
+                    self.mainCharactDictWholeHydro["Max altitude"]["value"] = float(tmp[2].split("[")[0])
+                    self.mainCharactDictWholeHydro["Max altitude"]["unit"] = "["+tmp[2].split("[")[1]
+                elif(list_data[i][0][:12]=="Min altitude"):
+                    tmp = list_data[i][0].split()
+                    self.mainCharactDictWholeHydro["Min altitude"] = {}
+                    self.mainCharactDictWholeHydro["Min altitude"]["value"] = float(tmp[2].split("[")[0])
+                    self.mainCharactDictWholeHydro["Min altitude"]["unit"] = "["+tmp[2].split("[")[1]
+                elif(list_data[i][0][:21]=="Fraction of landuse n"):
+                    tmp = list_data[i][0].split()
+                    self.mainCharactDictWholeHydro["Fraction of landuse n "+tmp[4]] = {}
+                    self.mainCharactDictWholeHydro["Fraction of landuse n "+tmp[4]]["value"] = float(tmp[5].split("[")[0])
+                    self.mainCharactDictWholeHydro["Fraction of landuse n "+tmp[4]]["unit"] = "["+tmp[5].split("[")[1]
+
+        else:
+            logging.error(f"File {fileNameHydro} does not exist. Cannot read averaged characteristics of the whole hydro subbasin.")
 
 
     def get_flood(self, path='', check_coherence=False):
@@ -2921,13 +2948,15 @@ class SubBasin:
     # It will then fill the given or not dict of landuse with a name property
     def get_landuse_index_transform_default(self, landuseDict:dict={}) -> dict:
 
-        if landuseDict==None:
+        if landuseDict is None:
             landuseDict = {}
 
         for key in cst.DEFAULT_LANDUSE:
             if not key in landuseDict:
                 landuseDict[key] = {}
             landuseDict[key]["name"] = cst.DEFAULT_LANDUSE[key]
+            landuseDict[key]["surface"] = 0.0
+            landuseDict[key]["unit"] = "[-]"
 
         return landuseDict
 
@@ -2966,6 +2995,8 @@ class SubBasin:
                 if not key in landuseDict:
                     landuseDict[key] = {}
                 landuseDict[key]["name"] = line[1].strip()
+                landuseDict[key]["surface"] = 0.0
+                landuseDict[key]["unit"] = "[-]"
 
         return landuseDict
 
@@ -2986,13 +3017,19 @@ class SubBasin:
         # Get the names of all landuses
         landuseDict = {}
         landuseDict = self.get_landuse_index_transform(landuse_index_transform, landuseDict=landuseDict)
-        if landuseDict == None:
+        if landuseDict is None:
             landuseDict = self.get_landuse_index_transform_default(landuseDict=landuseDict)
 
         # Get the surfaces and units of all landuses
         for key in landuseDict:
-            landuseDict[key]["surface"] = avrgDict["Fraction of landuse n "+str(key)]["value"]
-            landuseDict[key]["unit"] = avrgDict["Fraction of landuse n "+str(key)]["unit"]
+            associated_key = "Fraction of landuse n "+str(key)
+            if associated_key in avrgDict:
+                # If the lanuse type is not present, it will be set to 0.0
+                landuseDict[key]["surface"] = avrgDict[associated_key]["value"]
+                landuseDict[key]["unit"] = avrgDict[associated_key]["unit"]
+            else:
+                landuseDict[key]["surface"] = 0.0
+                landuseDict[key]["unit"] = "[-]"
 
         if toSave:
             if onlySub:
@@ -3024,9 +3061,9 @@ class SubBasin:
         # Creation of data and legend
         data = []
         names = []
-        for key in landuseDict:
-            data.append(landuseDict[key]["surface"])
-            names.append(landuseDict[key]["name"])
+        for landuse in landuseDict.values():
+            data.append(landuse["surface"])
+            names.append(landuse["name"])
 
         # Creation of colors
         colorDict = {}
@@ -3069,7 +3106,7 @@ class SubBasin:
 
     def change_version(self, newVersion=None):
 
-        if newVersion == None:
+        if newVersion is None:
             self._version = float(cst.VERSION_WOLFHYDRO)
         elif type(newVersion) == str:
             self._version = float(newVersion)
@@ -3300,7 +3337,7 @@ class SubBasin:
         # elif self.model == cst.tom_2layers_linIF:
         #     self.activate_all_iv_2layers()
         cur_dir = os.path.join(self.fileNameRead, "Subbasin_"+str(self.iDSorted))
-
+        mc.MODELS_VAR[self.model].deactivate_all(directory=cur_dir, prefix_file='simul')
         mc.MODELS_VAR[self.model].activate_all(directory=cur_dir, prefix_file='simul', type_of_var=iv.FINAL_OUT_VAR)
         mc.MODELS_VAR[self.model].activate_all(directory=cur_dir, prefix_file='simul', type_of_var=iv.IV_VAR)
 
@@ -3716,7 +3753,7 @@ class SubBasin:
         Returns:
         - dict: A dictionary containing the summary of the fractions of the current module.
         """
-
+        print(f"Computing flow fractions for station :{self.name}")
         if all_f == {}:
             cur_dir = os.path.join(self.fileNameRead, "Subbasin_"+str(self.iDSorted))
             all_qin = mc.MODELS_VAR[self.model].get_all_iv_timeseries(directory=cur_dir,
@@ -3725,8 +3762,16 @@ class SubBasin:
                                             prefix_file='simul', type_of_var=iv.DEFAULT_VAR)
             all_f.update(all_qin)
 
-        q_simul = self.get_outFlow(unit='mm/h')
-        all_r = {"%"+key: val/q_simul * 100.0 for key, val in all_f.items()}
+        q_simul = self.get_myHydro(unit='mm/h')
+        for key, val in all_f.items():
+            if len(val) != len(q_simul):
+                logging.warning(f"Length of {key} does not match with the length of the simulated flow!")
+                continue
+        # all_r = {"%"+key: val/q_simul * 100.0 for key, val in all_f.items()}
+        all_r = {
+            f"%{key}": np.where(np.isfinite(val / q_simul), (val / q_simul) * 100.0, np.nan)
+            for key, val in all_f.items()
+        }
 
         return self._operation_on_ts(all_r, summary=summary, interval=interval)
 
@@ -3749,7 +3794,7 @@ class SubBasin:
             all_iv = mc.MODELS_VAR[self.model].get_all_iv_timeseries(directory=cur_dir,
                                             prefix_file='simul', type_of_var=iv.IV_VAR)
         if max_params != {}:
-            out_dict = {key: all_iv[key]/cur_max for key, cur_max in max_params.items()}
+            out_dict = {key: all_iv[key]/cur_max*100 for key, cur_max in max_params.items()}
         else:
             out_dict = all_iv