BUG,ENH: add algorithm 3 for cube window

jcrivenaes · Sep 14, 2023 · e33e29b · e33e29b
1 parent 3e10946
commit e33e29b
Show file tree

Hide file tree

Showing 4 changed files with 592 additions and 9 deletions.
diff --git a/src/clib/xtg/surf_cube_attr_intv.c b/src/clib/xtg/surf_cube_attr_intv.c
@@ -145,7 +145,7 @@ pillar_attrs(double *pillar, int np, double *pattr)
 static void
 pillar_attrs_sums(double *pillar, int np, double *pattr)
 {
-    /* SUM attributes are only menaingful on a discrete level*/
+    /* SUM attributes are only meaningful on a discrete level*/
 
     int n;
     for (n = 11; n <= 13; n++)

diff --git a/src/xtgeo/surface/_regsurf_cube_window.py b/src/xtgeo/surface/_regsurf_cube_window.py
@@ -4,9 +4,11 @@
 
 import numpy as np
 import numpy.ma as ma
-from xtgeo.common import XTGeoDialog
-from xtgeo.common import XTGShowProgress
-from . import _regsurf_cube_window_v2 as cwv
+
+from xtgeo.common import XTGeoDialog, XTGShowProgress
+
+from . import _regsurf_cube_window_v2 as cwv2
+from . import _regsurf_cube_window_v3 as cwv3
 
 xtg = XTGeoDialog()
 
@@ -69,8 +71,25 @@ def slice_cube_window(
             deadtraces=deadtraces,
         )
 
+    elif algorithm == 2:
+        attrs = cwv2.slice_cube_window(
+            self,
+            cube,
+            zsurf=zsurf,
+            other=other,
+            other_position=other_position,
+            sampling=sampling,
+            mask=mask,
+            zrange=zrange,
+            ndiv=ndiv,
+            attribute=attribute,
+            maskthreshold=maskthreshold,
+            snapxy=snapxy,
+            showprogress=showprogress,
+            deadtraces=deadtraces,
+        )
     else:
-        attrs = cwv.slice_cube_window(
+        attrs = cwv3.slice_cube_window(
             self,
             cube,
             zsurf=zsurf,