feat(options): ✨ porosity/permeability thresholds

Optionally provide arbitrary minima
to cut off smallest porosities/permeabilities
when upscaling.

Author: djmaxus

src/Options.m

@@ -4,11 +4,21 @@
         sat_tol (1,1) double {mustBeNonnegative} = 1e-4
         hydrostatic_correction (1,1) logical = false
         m_save_mip_step (1,1) {mustBeOfClass(m_save_mip_step,'logical')} = false;
+        m_perm_threshold_mD (1,3) double {mustBeNonnegative} = [0,0,0];
+        m_poro_threshold    (1,1) double {mustBeNonnegative} = 0;
     end
 
     methods % builders
         function self = save_mip_step(self,flag)
             self.m_save_mip_step = flag;
         end
+
+        function self = perm_threshold_mD(self,value)
+            self.m_perm_threshold_mD = value;
+        end
+
+        function self = poro_threshold(self,value)
+            self.m_poro_threshold = value;
+        end
     end
 end

src/upscale.m

@@ -7,6 +7,14 @@
 
 mip(1:length(saturations)) = struct('sw',nan,'sub_sw',[]);
 
+% apply thresholds before proceeding
+small_poro = porosities < options.m_poro_threshold;
+porosities = porosities .* (~small_poro) + options.m_poro_threshold .* small_poro;
+
+perm_threshold = reshape(options.m_perm_threshold_mD,[1,1,1,3])*milli*darcy;
+small_Kabs = permeabilities < perm_threshold;
+permeabilities = permeabilities .* (~small_Kabs) + perm_threshold .* small_Kabs;
+
 if isscalar(porosities)
     poro_upscaled = porosities(1);
     perm_upscaled = reshape(permeabilities(:,:,1,:),1,3);