wolfhece 2.2.40__py3-none-any.whl → 2.2.42__py3-none-any.whl

wolfhece/wolf_array.py

@@ -116,6 +116,30 @@ try:
 except ImportError as e:
     raise Exception(_('Error importing pydownloader module'))
 
+
+import math
+from fractions import Fraction
+
+def pgcd_decimal(a, b):
+    # Convertir les décimaux en fractions exactes
+    fa = Fraction(a).limit_denominator()
+    fb = Fraction(b).limit_denominator()
+
+    # Extraire les numérateurs et dénominateurs
+    num_a, den_a = fa.numerator, fa.denominator
+    num_b, den_b = fb.numerator, fb.denominator
+
+    # Mettre au même dénominateur
+    lcm_den = math.lcm(den_a, den_b)
+    a_int = num_a * (lcm_den // den_a)
+    b_int = num_b * (lcm_den // den_b)
+
+    # PGCD des entiers
+    pgcd_int = math.gcd(a_int, b_int)
+
+    # Revenir à l’échelle décimale
+    return pgcd_int / lcm_den
+
 WOLF_ARRAY_HILLSHAPE = -1
 WOLF_ARRAY_FULL_SINGLE = 1
 WOLF_ARRAY_FULL_DOUBLE = 2
@@ -7288,6 +7312,8 @@ class WolfArray(Element_To_Draw, header_wolf):
         else:
             return
 
+        need_wx = False
+
         raster:gdal.Dataset
         if crop is not None :
             if not os.path.exists(fn):
@@ -7297,44 +7323,47 @@ class WolfArray(Element_To_Draw, header_wolf):
             tmpdx = self.dx
 
             fn_crop = fn + '_crop.tif'
-            if type(crop) is np.ndarray:
-                pass
-            elif type(crop) is list:
-                pass
-            else:
-                if not self.wx_exists:
-                    logging.error(_('Crop must be a list or a numpy array with 4 values - xmin, xmax, ymin, ymax'))
-                    return
 
-                raster_in:gdal.Dataset
-                raster_in = gdal.Open(fn)
+            raster_in:gdal.Dataset
+            raster_in = gdal.Open(fn)
 
-                proj = raster_in.GetProjection()
-                # check if proj is epsg:31370
-                if 'Lambert 72' not in proj and "31370" not in proj:
-                    logging.error(_('{} projection is not Lambert 72 - EPSG:31370').format(fn))
+            proj = raster_in.GetProjection()
+            # check if proj is epsg:31370
+            if 'Lambert 72' not in proj and "31370" not in proj:
+                logging.error(_('{} projection is not Lambert 72 - EPSG:31370').format(fn))
 
-                    raster_in2 = gdal.Dataset()
-                    gdal.Warp(raster_in2, raster_in, dstSRS='EPSG:31370')
+                raster_in2 = gdal.Dataset()
+                gdal.Warp(raster_in2, raster_in, dstSRS='EPSG:31370')
 
-                    # Close the raster_in
-                    raster_in.FlushCache()
-                    raster_in = None
+                # Close the raster_in
+                raster_in.FlushCache()
+                raster_in = None
 
-                    raster_in = raster_in2
+                raster_in = raster_in2
 
-                    if raster_in is None:
-                        logging.error(_('Could not reproject the file {fn}'))
-                        return
+                if raster_in is None:
+                    logging.error(_('Could not reproject the file {fn}'))
+                    return
+
+            geotr = raster_in.GetGeoTransform()
+            self.dx = geotr[1]
+            self.dy = abs(geotr[5])
+
+            ulx = geotr[0]
+            uly = geotr[3]
+            llx = ulx + geotr[1] * raster_in.RasterXSize
+            lly = uly + geotr[5] * raster_in.RasterYSize
 
-                geotr = raster_in.GetGeoTransform()
-                self.dx = geotr[1]
-                self.dy = abs(geotr[5])
+            if type(crop) is np.ndarray:
+                pass
+            elif type(crop) in [list, tuple]:
+                pass
+            else:
+                if not self.wx_exists:
+                    logging.error(_('Crop must be a list or a numpy array with 4 values - xmin, xmax, ymin, ymax'))
+                    return
 
-                ulx = geotr[0]
-                uly = geotr[3]
-                llx = ulx + geotr[1] * raster_in.RasterXSize
-                lly = uly + geotr[5] * raster_in.RasterYSize
+                need_wx = True
 
                 newcrop = CropDialog(None, self.get_mapviewer())
 
@@ -7375,7 +7404,7 @@ class WolfArray(Element_To_Draw, header_wolf):
 
             xmin, xmax, ymin, ymax = crop
 
-            if self.wx_exists:
+            if need_wx:
 
                 with wx.FileDialog(None, _('Save the cropped file for later'), wildcard="Tiff files (*.tif)|*.tif",
                                     style=wx.FD_SAVE | wx.FD_OVERWRITE_PROMPT) as fileDialog:
@@ -7393,15 +7422,13 @@ class WolfArray(Element_To_Draw, header_wolf):
 
                 fn  = fn_crop
             else:
-                from tempfile import NamedTemporaryFile
-                tmpfile = NamedTemporaryFile(suffix='.tif')
+                fn_tmp = '__tmp_wolf.tif'
 
                 if uly > lly:
-                    raster = gdal.Translate(str(tmpfile.name), fn, projWin=[xmin, ymax, xmax, ymin])
+                    raster = gdal.Translate(fn_tmp, fn, projWin=[xmin, ymax, xmax, ymin])
                 else:
-                    raster = gdal.Translate(str(tmpfile.name), fn, projWin=[xmin, ymin, xmax, ymax])
-
-                fn = str(tmpfile.name)
+                    raster = gdal.Translate(fn_tmp, fn, projWin=[xmin, ymin, xmax, ymax])
+                fn = fn_tmp
         else:
             raster = gdal.Open(fn)
 
@@ -7501,6 +7528,12 @@ class WolfArray(Element_To_Draw, header_wolf):
         raster.FlushCache()
         raster = None
 
+        try:
+            if fn_tmp == '__tmp_wolf.tif':
+                os.remove(fn_tmp)
+        except:
+            pass
+
     def add_ops_sel(self):
         """
         Adding selection manager and operations array
@@ -9235,7 +9268,13 @@ class WolfArray(Element_To_Draw, header_wolf):
             newArray.origz = self.origz
             newArray.translz = self.translz
 
-        if type(other) == float:
+        if type(other) in [float, int]:
+
+            if type(other) == int and self.wolftype in [WOLF_ARRAY_FULL_DOUBLE, WOLF_ARRAY_FULL_SINGLE]:
+                other = float(other)
+            elif type(other) == float and self.wolftype in [WOLF_ARRAY_FULL_INTEGER, WOLF_ARRAY_FULL_INTEGER16, WOLF_ARRAY_FULL_INTEGER16_2, WOLF_ARRAY_FULL_INTEGER8, WOLF_ARRAY_FULL_UINTEGER8]:
+                other = int(other)
+
             if other != 0.:
                 newArray.array = np.ma.masked_array(self.array + other, self.array.mask, dtype=self.array.dtype)
         else:
@@ -9264,7 +9303,13 @@ class WolfArray(Element_To_Draw, header_wolf):
             newArray.origz = self.origz
             newArray.translz = self.translz
 
-        if type(other) == float:
+        if type(other) in [float, int]:
+
+            if type(other) == int and self.wolftype in [WOLF_ARRAY_FULL_DOUBLE, WOLF_ARRAY_FULL_SINGLE]:
+                other = float(other)
+            elif type(other) == float and self.wolftype in [WOLF_ARRAY_FULL_INTEGER, WOLF_ARRAY_FULL_INTEGER16, WOLF_ARRAY_FULL_INTEGER16_2, WOLF_ARRAY_FULL_INTEGER8, WOLF_ARRAY_FULL_UINTEGER8]:
+                other = int(other)
+
             if other != 0.:
                 newArray.array = np.ma.masked_array(self.array * other, self.array.mask, dtype=self.array.dtype)
             else:
@@ -9296,7 +9341,13 @@ class WolfArray(Element_To_Draw, header_wolf):
             newArray.origz = self.origz
             newArray.translz = self.translz
 
-        if type(other) == float:
+        if type(other) in [float, int]:
+
+            if type(other) == int and self.wolftype in [WOLF_ARRAY_FULL_DOUBLE, WOLF_ARRAY_FULL_SINGLE]:
+                other = float(other)
+            elif type(other) == float and self.wolftype in [WOLF_ARRAY_FULL_INTEGER, WOLF_ARRAY_FULL_INTEGER16, WOLF_ARRAY_FULL_INTEGER16_2, WOLF_ARRAY_FULL_INTEGER8, WOLF_ARRAY_FULL_UINTEGER8]:
+                other = int(other)
+
             if other != 0.:
                 newArray.array = np.ma.masked_array(self.array - other, self.array.mask, dtype=self.array.dtype)
         else:
@@ -9325,7 +9376,18 @@ class WolfArray(Element_To_Draw, header_wolf):
             newArray.origz = self.origz
             newArray.translz = self.translz
 
-        newArray.array = np.ma.masked_array(self.array ** other, self.array.mask, dtype=self.array.dtype)
+        if type(other) in [float, int]:
+
+            if type(other) == int and self.wolftype in [WOLF_ARRAY_FULL_DOUBLE, WOLF_ARRAY_FULL_SINGLE]:
+                other = float(other)
+            elif type(other) == float and self.wolftype in [WOLF_ARRAY_FULL_INTEGER, WOLF_ARRAY_FULL_INTEGER16, WOLF_ARRAY_FULL_INTEGER16_2, WOLF_ARRAY_FULL_INTEGER8, WOLF_ARRAY_FULL_UINTEGER8]:
+                other = int(other)
+
+            newArray.array = np.ma.masked_array(self.array ** other, self.array.mask, dtype=self.array.dtype)
+        else:
+            logging.error(_('Power operator only implemented for scalars'))
+            newArray.array[:,:] = 0
+
         newArray.count()
 
         assert newArray.array.dtype == self.array.dtype, _('Bad dtype')
@@ -9350,9 +9412,17 @@ class WolfArray(Element_To_Draw, header_wolf):
             newArray.origz = self.origz
             newArray.translz = self.translz
 
-        if type(other) == float:
+        if type(other) in [float, int]:
+
+            if type(other) == int and self.wolftype in [WOLF_ARRAY_FULL_DOUBLE, WOLF_ARRAY_FULL_SINGLE]:
+                other = float(other)
+            elif type(other) == float and self.wolftype in [WOLF_ARRAY_FULL_INTEGER, WOLF_ARRAY_FULL_INTEGER16, WOLF_ARRAY_FULL_INTEGER16_2, WOLF_ARRAY_FULL_INTEGER8, WOLF_ARRAY_FULL_UINTEGER8]:
+                other = int(other)
+
             if other != 0.:
                 newArray.array = np.ma.masked_array(self.array / other, self.array.mask, dtype=self.array.dtype)
+            else:
+                logging.error(_('Division by 0 -- Nothing to do !'))
         else:
             newArray.array = np.ma.masked_array(np.where(other == 0., 0., self.array / other.array), self.array.mask, dtype=self.array.dtype)
         newArray.count()
@@ -10217,7 +10287,7 @@ class WolfArray(Element_To_Draw, header_wolf):
             self.write_txt_header()
             self.write_array()
 
-    def get_rebin_shape_size(self, factor:float) -> tuple[tuple[int, int], tuple[float, float]]:
+    def get_rebin_shape_size(self, factor:float, reshape_array_if_necessary:bool = False) -> tuple[tuple[int, int], tuple[float, float]]:
         """
         Return the new shape after rebinning.
 
@@ -10237,13 +10307,44 @@ class WolfArray(Element_To_Draw, header_wolf):
 
         newnbx = self.nbx
         newnby = self.nby
-        if np.mod(self.nbx,factor) != 0 or np.mod(self.nby,factor) != 0 :
-            newnbx = self.nbx
-            newnby = self.nby
-            if np.mod(self.nbx,factor) !=0:
-                newnbx = self.nbx + factor - np.mod(self.nbx,factor)
-            if np.mod(self.nby,factor) !=0:
-                newnby = self.nby + factor - np.mod(self.nby,factor)
+
+        if factor >=1.0:
+            # Decrease resolution
+            to_reshape = False
+            if np.mod(self.nbx,factor) != 0 or np.mod(self.nby,factor) != 0 :
+                newnbx = self.nbx
+                newnby = self.nby
+                if np.mod(self.nbx,factor) !=0:
+                    newnbx = int(self.nbx + factor - np.mod(self.nbx,factor))
+                    to_reshape = True
+                if np.mod(self.nby,factor) !=0:
+                    newnby = int(self.nby + factor - np.mod(self.nby,factor))
+                    to_reshape = True
+
+                if to_reshape and reshape_array_if_necessary:
+                    logging.warning(f"The shape ({self.nbx}, {self.nby}) is not a multiple of the factor {factor}")
+                    new_array = ma.masked_array(np.ones((newnbx, newnby), dtype= self.dtype) * self.nullvalue, dtype= self.dtype)
+                    new_array[:self.nbx, :self.nby] = self.array
+                    self.array = new_array
+        else:
+            # Increase resolution
+            if abs(int(1/factor) - 1/factor) > 1e-6:
+                # find the greatest common divisor of between self.dx and old_dx
+                gcd = pgcd_decimal(self.dx, self.dx * factor)
+                intermediate_factor = gcd / self.dx
+
+                [newnbx,newnby], [dx,dy] = self.get_rebin_shape_size(intermediate_factor)
+
+                factor = factor / intermediate_factor
+
+                if np.mod(newnbx,factor) != 0 or np.mod(newnby,factor) != 0 :
+                    if np.mod(newnbx,factor) !=0:
+                        newnbx = int(newnbx + factor - np.mod(newnbx,factor))
+                    if np.mod(newnby,factor) !=0:
+                        newnby = int(newnby + factor - np.mod(newnby,factor))
+
+                newdx = dx * factor
+                newdy = dy * factor
 
         newnbx = int(newnbx / factor)
         newnby = int(newnby / factor)
@@ -10290,6 +10391,9 @@ class WolfArray(Element_To_Draw, header_wolf):
 
         """
 
+        op_str = operation
+        (testnbx, testnby), __ = self.get_rebin_shape_size(factor)
+
         if operation_matrix is not None:
             tmp_header = self.get_rebin_header(factor)
             if not operation_matrix.is_like(tmp_header):
@@ -10311,36 +10415,11 @@ class WolfArray(Element_To_Draw, header_wolf):
                 logging.error(_("Operator not supported -- Must be a string in ['sum', 'mean', 'min', 'max', 'median'] or a Rebin_Ops Enum"))
 
 
-        if np.mod(self.nbx,factor) != 0 or np.mod(self.nby,factor) != 0 :
-            # The shape is adjusted to be a multiple of the factor.
-            # Fill the array with nullvalue
-            newnbx = self.nbx
-            newnby = self.nby
-            if np.mod(self.nbx,factor) !=0:
-                newnbx = int(self.nbx + factor - np.mod(self.nbx,factor))
-            if np.mod(self.nby,factor) !=0:
-                newnby = int(self.nby + factor - np.mod(self.nby,factor))
-
-            # We ensure that the null_value is converted to the type of the array.
-            # Not doing so leads to a change in the array type because
-            # of numpy2 type promotion.
-            newarray = np.ma.masked_array( np.full((newnbx,newnby), self.nullvalue, dtype = self.dtype) )
-            #newarray = np.ma.ones((newnbx,newnby), dtype = self.dtype) * self.nullvalue
-            newarray[:self.nbx,:self.nby] = self.array
-            newarray.mask[:self.nbx,:self.nby] = self.array.mask
-            self.array = newarray
-
-            self.nbx = newnbx
-            self.nby = newnby
-
-        self.nbx = int(self.nbx / factor)
-        self.nby = int(self.nby / factor)
-
-        self.dx = self.dx * float(factor)
-        self.dy = self.dy * float(factor)
-        new_shape = (self.nbx, self.nby)
-
-        if factor>1.:
+        new_shape, (dx, dy) = self.get_rebin_shape_size(factor, reshape_array_if_necessary= True)
+
+        if factor >=1.0:
+            # Decrease resolution
+
             if operation_matrix is not None:
                 # Reshape the input array to split it into blocks of size f x f
                 reshaped_a = self.array.reshape(new_shape[0], int(factor), new_shape[1], int(factor))
@@ -10370,15 +10449,43 @@ class WolfArray(Element_To_Draw, header_wolf):
             self.set_nullvalue_in_mask()
         else:
 
-            if (int(1/factor) - 1/factor) > 1e-6:
+            # Increase resolution
+
+            if abs(int(1/factor) - 1/factor) > 1e-6:
                 logging.warning(f"The factor {factor} doesn't lead to an integer dimension for the Kronecker product")
+                logging.warning("The array will be rebinned firstly using the most common factor and then using the operation")
+
+                # find the greatest common divisor of between self.dx and old_dx
+                gcd  = pgcd_decimal(self.dx, self.dx * factor)
+                intermediate_factor = gcd / self.dx
+
+                logging.warning(f"The greatest common divisor is {gcd}")
 
-            ones = np.ones( (int(1/factor), int(1/factor)), dtype=int)
-            self.array = np.kron(self.array, ones).astype(self.array.dtype)
+                tmp = WolfArray(mold=self, whichtype=self.wolftype)
+                tmp.rebin(intermediate_factor)
 
+                if op_str is None:
+                    op_str = 'min'
+
+                tmp.rebin(factor / intermediate_factor, operation= op_str)
+                self.array = tmp.array
+                tmp.array = None
+                del tmp
+            else:
+
+                ones = np.ones( (int(1/factor), int(1/factor)), dtype=int)
+                self.array = np.kron(self.array, ones).astype(self.array.dtype)
 
             self.mask_data(self.nullvalue)
 
+        self.nbx = new_shape[0]
+        self.nby = new_shape[1]
+        self.dx = dx
+        self.dy = dy
+
+        assert isinstance(self.array, ma.MaskedArray), _('The array must be a masked array')
+        assert self.array.shape == (testnbx, testnby), _(f'Bad shape: {self.array.shape} != {(testnbx, testnby)}')
+
         self.count()
 
         # rebin must not change the type of the array

wolfhece/wolfresults_2D.py

@@ -5431,7 +5431,7 @@ class Wolfresults_2D(Element_To_Draw):
                     head[ij]= z[ij] + v[ij]**2./2/9.81
 
                     # Fill the time of arrival
-                    danger_map_matrix_toa[curblock].array[(danger_map_matrix_toa.array == DEFAULT_TOA) & (wd.array > hmin)] = cur_time
+                    danger_map_matrix_toa[curblock].array[(danger_map_matrix_toa.array == DEFAULT_TOA) & (wd.array > hmin) & (~wd.array.mask)] = cur_time
 
                     # Fill the time of maximum
                     #  Searching where wd > h_max
@@ -5466,7 +5466,7 @@ class Wolfresults_2D(Element_To_Draw):
                 head[ij]= z[ij] + v[ij]**2./2/9.81
 
                 # Fill the time of arrival
-                danger_map_matrix_toa.array[(danger_map_matrix_toa.array == DEFAULT_TOA) & (wd.array > hmin)] = cur_time
+                danger_map_matrix_toa.array[(danger_map_matrix_toa.array == DEFAULT_TOA) & (wd.array > hmin)  & (~wd.array.mask)] = cur_time
 
                 # Fill the time of maximum
                 #  Searching where wd > h_max
@@ -5616,7 +5616,7 @@ class Wolfresults_2D(Element_To_Draw):
                     danger_toa = danger_map_matrix_toa[i].array
                     res_toa = result[5][i].array
 
-                    ij = np.where((danger_toa.data == DEFAULT_TOA) & (res_toa.data > 0.))
+                    ij = np.where(((danger_toa.data == DEFAULT_TOA) & (res_toa.data > 0.) & (~res_toa.mask)) | (danger_toa.data > res_toa.data))
                     danger_toa.data[ij] = res_toa.data[ij]
                     danger_toa.mask[ij] = False
 
@@ -5657,7 +5657,7 @@ class Wolfresults_2D(Element_To_Draw):
                 danger_toa = danger_map_matrix_toa.array
                 res_toa = result[5].array
 
-                ij = np.where((danger_toa.data == DEFAULT_TOA) & (res_toa.data > 0.))
+                ij = np.where(((danger_toa.data == DEFAULT_TOA) & (res_toa.data > DEFAULT_TOA) & (~res_toa.mask)) | (danger_toa.data > res_toa.data))
                 danger_toa.data[ij] = res_toa.data[ij]
                 danger_toa.mask[ij] = False
 

{wolfhece-2.2.40.dist-info → wolfhece-2.2.42.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: wolfhece
-Version: 2.2.40
+Version: 2.2.42
 Author-email: Pierre Archambeau <pierre.archambeau@uliege.be>
 Project-URL: Homepage, https://uee.uliege.be/hece
 Project-URL: Issues, https://uee.uliege.be/hece

{wolfhece-2.2.40.dist-info → wolfhece-2.2.42.dist-info}/RECORD

@@ -9,7 +9,7 @@ wolfhece/MulticriteriAnalysis.py,sha256=vGmkzYagZohNe0XjwGJ6VUXcDPjOt80lNFthXpzx
 wolfhece/PandasGrid.py,sha256=etfVhIHzja4Z1EUY6BcDOKX-w7V-Xou1yaf0NMqmclo,4599
 wolfhece/PyConfig.py,sha256=13DDWjJdohYHwn1uRVHB0s8Jcwq_b9pwcwbAr8NlZyc,19667
 wolfhece/PyCrosssections.py,sha256=6K9xhupnIVXzLdsLTPN9e1v0xck9BdWZd0-EgpWg9us,133277
-wolfhece/PyDraw.py,sha256=XjIi5yDcl56C5fFHQUFdJvyU7r8jJ9z7aRVPOZijQcI,727297
+wolfhece/PyDraw.py,sha256=jCXBonWGmqTQD9AXsii-yalFpUFLKTrfPP2wzh-57NE,741859
 wolfhece/PyGui.py,sha256=GpVRxNpR8WNDFyHnDvhtHFFsq_cZZlyVgSkFiS-ARYI,185342
 wolfhece/PyGuiHydrology.py,sha256=dmBlRO8AljsvCPH6eVt0l9ZLx7g5j7Ubl9Srk7ECwyA,34693
 wolfhece/PyHydrographs.py,sha256=1P5XAURNqCvtSsMQXhOn1ihjTpr725sRsZdlCEhhk6M,3730
@@ -56,24 +56,24 @@ wolfhece/pyGui1D.py,sha256=9g7OS3YiKsqy--6y0cBD7x2gaqTTYFXWkxImpgnTA20,121937
 wolfhece/pyLandUseFlanders.py,sha256=4R4J4BgP1-W0sebnb-TBreStZx3DKOh1zYj_dg23Ac4,8657
 wolfhece/pybridges.py,sha256=bFAqjL4ColeJtwvyCPGQ8VllWoq1RbVWXxFrdfrvqm8,65954
 wolfhece/pydike.py,sha256=dRb6qGkqoTXjf107KcajcIk1F_FuMPaOZLSwixT3dgA,11196
-wolfhece/pydownloader.py,sha256=EpHYZwEBRt3cyxRd_a7E7lx1VlwlYeQYNOAKVacJO4c,11518
+wolfhece/pydownloader.py,sha256=jdyOpIZPGgHlo6MakklLHIqsTmDMrbDVeKchSpij_K4,16245
 wolfhece/pylogging.py,sha256=4TI8hgBB65z-zpvU5Rfa2jkPXPhJaqXjHVPwbcdzTNc,4528
 wolfhece/pypolygons_scen.py,sha256=NWaNeK0RSUeOkgukeogK9FLmQiDjGZ9yhqs9208fojM,46237
 wolfhece/pyshields.py,sha256=KMtUO5kD0lisKnJD1NsDz-qaY5DpFcmS4O3WkXtUSmo,27898
 wolfhece/pyviews.py,sha256=zuZjWUptRDm1MTE1PN4Xj_qSITnojgDMG0LlFIBH3SE,13739
-wolfhece/pywalous.py,sha256=6FTDh0dzKg7H6F7_YXEKT2gmE0OM_a8hvOLhEu4kRpY,26729
+wolfhece/pywalous.py,sha256=Vl0WH0OZLOguV4JWhfJNwo0CgGneyD3DznFbyfAME-w,32221
 wolfhece/rain_SPWMI.py,sha256=qCfcmF7LajloOaCwnTrrSMzyME03YyilmRUOqrPrv3U,13846
 wolfhece/textpillow.py,sha256=7hgfsLYAaE_rNKD-g8xsON8sdWvoV8vbqnGGxIayShE,14137
 wolfhece/tools2d_dll.py,sha256=TfvvmyZUqEZIH0uHwUCJf0bdmCks_AiidDt23Unsp5w,13550
 wolfhece/tools_mpl.py,sha256=gQ3Jg1iuZiecmMqa5Eli2ZLSkttu68VXL8YmMDBaEYU,564
 wolfhece/toolshydrology_dll.py,sha256=cIGyhxV8H5f7GXhDqAamM7uC0W0hQTou3eTkqZdnqBE,5656
-wolfhece/wolf_array.py,sha256=vhsxF79jcdK1tY1GFjCWx8YupC8D9r-C7nZ9Is2pMO4,574421
+wolfhece/wolf_array.py,sha256=39O7kYjksEbCC2yn0mu6AfabmVrBWrmoybbfjw4gfIQ,578890
 wolfhece/wolf_hist.py,sha256=fTEb60Q4TEwobdZsRU4CFXAId1eOKdWAqF8lnF1xEWc,3590
 wolfhece/wolf_texture.py,sha256=Pt1j_lX74p70Fj3y3qYxYMuN8gghVd8_ih1vFhTIdkA,23884
 wolfhece/wolf_tiles.py,sha256=v-HohqaWuMYdn75XLnA22dlloAG90iwnIqrgnB0ASQ4,10488
 wolfhece/wolf_vrt.py,sha256=wbxXVN7TL9zgdyF79S-4e3pje6wJEAgBEfF_Y8kkzxs,14271
 wolfhece/wolf_zi_db.py,sha256=kz1TofsPfuBQqwYEM32wZl4b-jqgAygorZSm2JnGJ_k,42450
-wolfhece/wolfresults_2D.py,sha256=0c3akPDFETE0Uq0JiZ0FZEx6Rxiy9JB5s9P6sDei0E0,245735
+wolfhece/wolfresults_2D.py,sha256=BTewbH0zX55iLLbY3KuHVMAOCsvif2FMvF1SQMJFgf0,245893
 wolfhece/xyz_file.py,sha256=1pzLFmmdHca4yBVR9Jitic6N82rY28mRytGC1zMbY28,6615
 wolfhece/acceptability/Parallels.py,sha256=2wVkfJYor4yl7VYiAZiGGTFwtAab2z66ZfRtBliVweE,4088
 wolfhece/acceptability/__init__.py,sha256=hfgoPKLDpX7drN1Vpvux-_5Lfyc_7feT2C2zQr5v-Os,258
@@ -95,7 +95,7 @@ wolfhece/apps/curvedigitizer.py,sha256=lEJJwgAfulrrWQc-U6ij6sj59hWN3SZl4Yu1kQxVz
 wolfhece/apps/hydrometry.py,sha256=lhhJsFeb4zGL4bNQTs0co85OQ_6ssL1Oy0OUJCzhfYE,656
 wolfhece/apps/isocurrent.py,sha256=dagmGR8ja9QQ1gwz_8fU-N052hIw-W0mWGVkzLu6C7I,4247
 wolfhece/apps/splashscreen.py,sha256=EdGDN9NhudIiP7c3gVqj7dp4MWFB8ySizM_tpMnsgpE,3091
-wolfhece/apps/version.py,sha256=1RyOO8Q7UPJwuggHsbQr86F755yJaRu6lbSflfgskcg,388
+wolfhece/apps/version.py,sha256=aOHEtZqRpvSjcC4VtpzMUp8extpQ1AUSw4MRuICsvz0,388
 wolfhece/apps/wolf.py,sha256=mRnjYsUu4KIsRuamdQWAINFMuwN4eJgMo9erG-hkZ70,729
 wolfhece/apps/wolf2D.py,sha256=4z_OPQ3IgaLtjexjMKX9ppvqEYyjFLt1hcfFABy3-jU,703
 wolfhece/apps/wolf_logo.bmp,sha256=ruJ4MA51CpGO_AYUp_dB4SWKHelvhOvd7Q8NrVOjDJk,3126
@@ -131,7 +131,7 @@ wolfhece/eva/bootstrap.py,sha256=Ys4xTDIvG_QtxCKWLYzb3_XAZU441jGX7fHIbd9Mvr0,840
 wolfhece/eva/hydrogramme_mono.py,sha256=uZFIgJJ-JogMFzt7D7OnyVaHvgxCQJPZz9W9FgnuthA,8138
 wolfhece/eva/joint_models.py,sha256=RqPLb0g775zJIe3X3nwaf1g1dqnEIUbJjPMae_J_nmM,5863
 wolfhece/eva/mixture_models.py,sha256=Eq9fU7GjbSLHqCE5mdkIC5oAJhv6aNxY9kwnqlcOmc8,16221
-wolfhece/eva/pyseries.py,sha256=Oz2QCPGrbebAwQL1Vl-kc4t1OkzS-pCuj5rZwDGEJ-I,118536
+wolfhece/eva/pyseries.py,sha256=iUTM31qZefT-iKRvMEmVfWXEZMr3_h0ypPqwFQ1Zm9E,119590
 wolfhece/fonts/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 wolfhece/fonts/arial.ttf,sha256=lXZrWPfYabD6LPbm_rJsGyHN8mMfHFhj_JvSBtXG6O4,915212
 wolfhece/fonts/helvetica.ttf,sha256=X4Zd3zdUmuRGMLE6UB-BMIbirpdK3Ia5czfNnuSx5P8,317968
@@ -161,8 +161,8 @@ wolfhece/hydrology/read.py,sha256=H-sK4hU23snrZC3BDLUHJqUTUBwvQ48ZFgxE-6PxDYg,11
 wolfhece/hydrology/slope_manager.py,sha256=vlek0z8qcqB61eleiksyOe3QR1vpbtwfeowy6ms7_Fg,5580
 wolfhece/hydrology/wolfMap_treatment.py,sha256=eAxr24zJGwmDof1aZpcxewVvv_bWDvoO8t9Wwf99Mlo,10606
 wolfhece/hydrometry/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-wolfhece/hydrometry/kiwis.py,sha256=rhfvRKlKEKIDiwnGfHnT5oDa-I-RHU2ikPqD25urdHw,65463
-wolfhece/hydrometry/kiwis_gui.py,sha256=s4xCsqkKvI1YcIZ4X2A7hV11zyJFI8Q_4i82jOO6_ys,25060
+wolfhece/hydrometry/kiwis.py,sha256=aXf9QpbDld2K3jWOYLVcGj5Y_GBs2wFCUv7HGVHTK3w,65557
+wolfhece/hydrometry/kiwis_gui.py,sha256=bXUfIRlP0Cmcoe-sw-0KaFLoAUfQ8yCVPBX-xLSfuBs,25113
 wolfhece/hydrometry/kiwis_wolfgui.py,sha256=GPa5YNp2V2ZlxYQjPiCXERgPuWB_zij4ec2yHpRvoFA,4027
 wolfhece/icons/folder_open.png,sha256=ykBlU2FtGcpjAu1YO4_eGlDg2svgs_IAs-4d5O2e3AM,2329
 wolfhece/icons/folder_plus.png,sha256=Z1Jcs87Z6pNcudxGDTNhC_hgolSSaXe6_fmQE_n-ero,1963
@@ -235,7 +235,7 @@ wolfhece/mar/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 wolfhece/mar/commontools.py,sha256=SiSxpv5BFYEBCEydsE4ZmBuot3KTy0UYMx2aa-4fbuQ,52549
 wolfhece/mar/interface_MAR_WOLF.py,sha256=MWeXaHLDT4Eo9jZOAvz013lmpgGYT1v9VUYGAgBgSRU,21454
 wolfhece/math_parser/__init__.py,sha256=Mi7YTrlJtcflyrRdZHHgE-uNPUFfOWmsf8FsOwKBRPI,27961
-wolfhece/math_parser/calculator.py,sha256=Avyx-ikIwsficqd-Q4erX6U0zGQ3cAWfADyq6_uHpB0,7493
+wolfhece/math_parser/calculator.py,sha256=NYPcNWMBhbnYTK3tqEnbn1vyAcyw_ogNDki1SN7Do_4,7497
 wolfhece/mesh2d/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 wolfhece/mesh2d/bc_manager.py,sha256=BmNT-M0uWe5--AhxiS4EWziHad4RL-D7wUySqXTICsY,54237
 wolfhece/mesh2d/cell_tracker.py,sha256=mPmnD5lEf3gLPuLqtAIo-Gp-ipAwQdPxzjWOGt0b7jM,8958
@@ -243,7 +243,7 @@ wolfhece/mesh2d/config_manager.py,sha256=DcdxCIIs_dyC6ayJOBULeY364LONogL9PBaqBtC
 wolfhece/mesh2d/cst_2D_boundary_conditions.py,sha256=r43pHHdCtmNp5R2zh1Ckb7EzwQDf6C4YMLwRFTl4KMc,5006
 wolfhece/mesh2d/gpu_2d.py,sha256=xyZst3ZSYmRp9G0kxlUvTahAiMC9sNEH0MRUlEjqFZI,25999
 wolfhece/mesh2d/simple_2d.py,sha256=wqENJwpUPxKQcpGIcymQXUj2KgkGWCVH6cs4Os9h9Gs,112581
-wolfhece/mesh2d/wolf2dprev.py,sha256=oK9r94CYJaeOBSv4s4mkuhLA8rCsnrkjG2sGW_zR1x0,493127
+wolfhece/mesh2d/wolf2dprev.py,sha256=823RckzGGCQ00h0uti-nL7JoVsAZ_Hv9C02lnvizDIg,495518
 wolfhece/models/5_coul.pal,sha256=OI1UqcNIDBpJn2k_VDel__r-hKjjvdob0eqinGCI3QY,160
 wolfhece/models/6_coul.pal,sha256=z7NK2dg0tAQBUweRQV54dIwJbPM1U5y1AR2LLw19Idw,148
 wolfhece/models/7_coul.pal,sha256=XTnnUyCE8ONokScB2YzYDnSTft7E6sppmr7P-XwMsCE,205
@@ -265,6 +265,7 @@ wolfhece/models/shields_cst.pal,sha256=zUGFI6HiL0bsHeOzcWNih3F9cxXKXLLZYA5rtqRbz
 wolfhece/models/vulnerability.pal,sha256=Fevrc_9owywLhbPMBunXDcrGXPJhARo9iSV1eOq3roA,106
 wolfhece/models/walous_niv1.pal,sha256=mHMjCB-ja47mV3ZsvDOhS2CEK8YN6ewOkf1W7l7JQ8k,138
 wolfhece/models/walous_niv2.pal,sha256=B5wt5-O88dpaiA6yR4uriOy0tfDVU_cms0Xr6Dw9ZGg,377
+wolfhece/models/walous_ocs.pal,sha256=DlicMO2qlH3_WPm_bwVBK6lTuUn2MXRcfIpkjl8D2oo,227
 wolfhece/models/waterdepths.pal,sha256=8rcQfuZOeLKzYv5sARPkhpvZYc1OToj3ZukcbuRUgIY,136
 wolfhece/models/white_black.pal,sha256=BKjVN1kuk1OtJB9S3uM5fxteWbggI3Pb8p9WdDptT4c,51
 wolfhece/opengl/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
@@ -290,7 +291,7 @@ wolfhece/report/tools.py,sha256=ZA4PfHQXgvDtTnrSZh-oSwhQT_cQtb49LuZB8VI7nyI,9685
 wolfhece/report/wolf_report.png,sha256=NoSV58LSwb-oxCcZScRiJno-kxDwRdm_bK-fiMsKJdA,592485
 wolfhece/scenario/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 wolfhece/scenario/check_scenario.py,sha256=d-LWa_FxmPxTSc_H1lDHwqLB6TCqj1IUrRJhatfPMMA,5623
-wolfhece/scenario/config_manager.py,sha256=qosuoNOVbBuDuBUjupX5MFygYG2YJyr2AFgxxFqqdLo,126593
+wolfhece/scenario/config_manager.py,sha256=K4KLL2ZFhlZl88uxFL4A4t_W6qan2zeBehlnjvi709w,126969
 wolfhece/scenario/imposebc_void.py,sha256=PqA_99hKcaqK5zsK6IRIc5Exgg3WVpgWU8xpwNL49zQ,5571
 wolfhece/scenario/update_void.py,sha256=Yb7TMIUx9Gzm9_6qRMJnF39Uqi17dIkMmscSXo2WaTs,10033
 wolfhece/shaders/fragment_shader_texture.glsl,sha256=w6h8d5mJqFaGbao0LGmjRcFFdcEQ3ICIl9JpuT71K5k,177
@@ -319,8 +320,8 @@ wolfhece/ui/wolf_multiselection_collapsiblepane.py,sha256=u4C7CXe_bUyGKx7c_Bi0x9
 wolfhece/ui/wolf_times_selection_comparison_models.py,sha256=ORy7fz4dcp691qKzaOZHrRLZ0uXNhL-LIHxmpDGL6BI,5007
 wolfhece/wintab/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 wolfhece/wintab/wintab.py,sha256=8A-JNONV6ujgsgG3lM5Uw-pVgglPATwKs86oBzzljoc,7179
-wolfhece-2.2.40.dist-info/METADATA,sha256=XQvIHP3e-urKdNHAi5G8nvpKNuFSyBxxYuvePwSh59w,2785
-wolfhece-2.2.40.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-wolfhece-2.2.40.dist-info/entry_points.txt,sha256=Jr187pyvA3EeJiQLjZK9yo6mJX7IAn6ygZU9T8qF_gQ,658
-wolfhece-2.2.40.dist-info/top_level.txt,sha256=EfqZXMVCn7eILUzx9xsEu2oBbSo9liWPFWjIHik0iCI,9
-wolfhece-2.2.40.dist-info/RECORD,,
+wolfhece-2.2.42.dist-info/METADATA,sha256=qjWwy3fkLflM8KkQ6xXqvlK-GY598fnWupRva1TyNPM,2785
+wolfhece-2.2.42.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+wolfhece-2.2.42.dist-info/entry_points.txt,sha256=Jr187pyvA3EeJiQLjZK9yo6mJX7IAn6ygZU9T8qF_gQ,658
+wolfhece-2.2.42.dist-info/top_level.txt,sha256=EfqZXMVCn7eILUzx9xsEu2oBbSo9liWPFWjIHik0iCI,9
+wolfhece-2.2.42.dist-info/RECORD,,