xtgeo 4.6.0__cp39-cp39-win_amd64.whl → 4.7.0__cp39-cp39-win_amd64.whl

cxtgeo.py

@@ -538,9 +538,6 @@ def grd3d_get_all_corners(nx, ny, nz, swig_np_dbl_in_v1, swig_np_dbl_in_v2, swig
 def grd3d_well_ijk(nx, ny, nz, swig_np_dbl_in_v1, swig_np_dbl_in_v2, swig_np_int_in_v1, swig_np_dbl_in_v3, swig_np_int_in_v2, nval, p_utme_v, p_utmn_v, p_tvds_v, ivector, jvector, kvector, iflag):
     return _cxtgeo.grd3d_well_ijk(nx, ny, nz, swig_np_dbl_in_v1, swig_np_dbl_in_v2, swig_np_int_in_v1, swig_np_dbl_in_v3, swig_np_int_in_v2, nval, p_utme_v, p_utmn_v, p_tvds_v, ivector, jvector, kvector, iflag)
 
-def grdcp3d_from_cube(ncol, nrow, nlay, swig_np_dbl_inplaceflat_v1, swig_np_flt_inplaceflat_v1, swig_np_int_inplaceflat_v1, xori, yori, zori, xinc, yinc, zinc, rotation, yflip, option):
-    return _cxtgeo.grdcp3d_from_cube(ncol, nrow, nlay, swig_np_dbl_inplaceflat_v1, swig_np_flt_inplaceflat_v1, swig_np_int_inplaceflat_v1, xori, yori, zori, xinc, yinc, zinc, rotation, yflip, option)
-
 def grd3cp3d_xtgformat1to2_geom(ncol, nrow, nlay, swig_np_dbl_inplaceflat_v1, swig_np_dbl_inplaceflat_v2, swig_np_dbl_inplaceflat_v3, swig_np_flt_inplaceflat_v1, swig_np_int_inplaceflat_v1, swig_np_int_inplaceflat_v2):
     return _cxtgeo.grd3cp3d_xtgformat1to2_geom(ncol, nrow, nlay, swig_np_dbl_inplaceflat_v1, swig_np_dbl_inplaceflat_v2, swig_np_dbl_inplaceflat_v3, swig_np_flt_inplaceflat_v1, swig_np_int_inplaceflat_v1, swig_np_int_inplaceflat_v2)
 
@@ -559,9 +556,6 @@ def grdcp3d_get_vtk_esg_geometry_data(ncol, nrow, nlay, swig_np_dbl_inplaceflat_
 def grdcp3d_get_vtk_grid_arrays(ncol, nrow, nlay, swig_np_dbl_inplaceflat_v1, swig_np_flt_inplaceflat_v1, swig_np_dbl_aout_v1, swig_np_dbl_aout_v2, swig_np_dbl_aout_v3):
     return _cxtgeo.grdcp3d_get_vtk_grid_arrays(ncol, nrow, nlay, swig_np_dbl_inplaceflat_v1, swig_np_flt_inplaceflat_v1, swig_np_dbl_aout_v1, swig_np_dbl_aout_v2, swig_np_dbl_aout_v3)
 
-def grdcp3d_conv_grid_roxapi(ncol, nrow, nlay, swig_np_dbl_inplaceflat_v1, swig_np_flt_inplaceflat_v1, swig_np_dbl_aout_v1, swig_np_dbl_aout_v2, swig_np_dbl_aout_v3):
-    return _cxtgeo.grdcp3d_conv_grid_roxapi(ncol, nrow, nlay, swig_np_dbl_inplaceflat_v1, swig_np_flt_inplaceflat_v1, swig_np_dbl_aout_v1, swig_np_dbl_aout_v2, swig_np_dbl_aout_v3)
-
 def grdcp3d_quality_indicators(ncol, nrow, nlay, swig_np_dbl_inplaceflat_v1, swig_np_flt_inplaceflat_v1, swig_np_int_inplaceflat_v1, swig_np_flt_inplaceflat_v2):
     return _cxtgeo.grdcp3d_quality_indicators(ncol, nrow, nlay, swig_np_dbl_inplaceflat_v1, swig_np_flt_inplaceflat_v1, swig_np_int_inplaceflat_v1, swig_np_flt_inplaceflat_v2)
 

cxtgeoPYTHON_wrap.c

@@ -18687,160 +18687,6 @@ fail:
 }
 
 
-SWIGINTERN PyObject *_wrap_grdcp3d_from_cube(PyObject *self, PyObject *args) {
-  PyObject *resultobj = 0;
-  int arg1 ;
-  int arg2 ;
-  int arg3 ;
-  double *arg4 = (double *) 0 ;
-  long arg5 ;
-  float *arg6 = (float *) 0 ;
-  long arg7 ;
-  int *arg8 = (int *) 0 ;
-  long arg9 ;
-  double arg10 ;
-  double arg11 ;
-  double arg12 ;
-  double arg13 ;
-  double arg14 ;
-  double arg15 ;
-  double arg16 ;
-  int arg17 ;
-  int arg18 ;
-  int val1 ;
-  int ecode1 = 0 ;
-  int val2 ;
-  int ecode2 = 0 ;
-  int val3 ;
-  int ecode3 = 0 ;
-  PyArrayObject *array4 = NULL ;
-  int i4 = 1 ;
-  PyArrayObject *array6 = NULL ;
-  int i6 = 1 ;
-  PyArrayObject *array8 = NULL ;
-  int i8 = 1 ;
-  double val10 ;
-  int ecode10 = 0 ;
-  double val11 ;
-  int ecode11 = 0 ;
-  double val12 ;
-  int ecode12 = 0 ;
-  double val13 ;
-  int ecode13 = 0 ;
-  double val14 ;
-  int ecode14 = 0 ;
-  double val15 ;
-  int ecode15 = 0 ;
-  double val16 ;
-  int ecode16 = 0 ;
-  int val17 ;
-  int ecode17 = 0 ;
-  int val18 ;
-  int ecode18 = 0 ;
-  PyObject *swig_obj[15] ;
-  
-  (void)self;
-  if (!SWIG_Python_UnpackTuple(args, "grdcp3d_from_cube", 15, 15, swig_obj)) SWIG_fail;
-  ecode1 = SWIG_AsVal_int(swig_obj[0], &val1);
-  if (!SWIG_IsOK(ecode1)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode1), "in method '" "grdcp3d_from_cube" "', argument " "1"" of type '" "int""'");
-  } 
-  arg1 = (int)(val1);
-  ecode2 = SWIG_AsVal_int(swig_obj[1], &val2);
-  if (!SWIG_IsOK(ecode2)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "grdcp3d_from_cube" "', argument " "2"" of type '" "int""'");
-  } 
-  arg2 = (int)(val2);
-  ecode3 = SWIG_AsVal_int(swig_obj[2], &val3);
-  if (!SWIG_IsOK(ecode3)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode3), "in method '" "grdcp3d_from_cube" "', argument " "3"" of type '" "int""'");
-  } 
-  arg3 = (int)(val3);
-  {
-    array4 = obj_to_array_no_conversion(swig_obj[3], NPY_DOUBLE);
-    if (!array4 || !require_c_or_f_contiguous(array4)
-      || !require_native(array4)) SWIG_fail;
-    arg4 = (double*) array_data(array4);
-    arg5 = 1;
-    for (i4=0; i4 < array_numdims(array4); ++i4) arg5 *= array_size(array4,i4);
-  }
-  {
-    array6 = obj_to_array_no_conversion(swig_obj[4], NPY_FLOAT);
-    if (!array6 || !require_c_or_f_contiguous(array6)
-      || !require_native(array6)) SWIG_fail;
-    arg6 = (float*) array_data(array6);
-    arg7 = 1;
-    for (i6=0; i6 < array_numdims(array6); ++i6) arg7 *= array_size(array6,i6);
-  }
-  {
-    array8 = obj_to_array_no_conversion(swig_obj[5], NPY_INT);
-    if (!array8 || !require_c_or_f_contiguous(array8)
-      || !require_native(array8)) SWIG_fail;
-    arg8 = (int*) array_data(array8);
-    arg9 = 1;
-    for (i8=0; i8 < array_numdims(array8); ++i8) arg9 *= array_size(array8,i8);
-  }
-  ecode10 = SWIG_AsVal_double(swig_obj[6], &val10);
-  if (!SWIG_IsOK(ecode10)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode10), "in method '" "grdcp3d_from_cube" "', argument " "10"" of type '" "double""'");
-  } 
-  arg10 = (double)(val10);
-  ecode11 = SWIG_AsVal_double(swig_obj[7], &val11);
-  if (!SWIG_IsOK(ecode11)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode11), "in method '" "grdcp3d_from_cube" "', argument " "11"" of type '" "double""'");
-  } 
-  arg11 = (double)(val11);
-  ecode12 = SWIG_AsVal_double(swig_obj[8], &val12);
-  if (!SWIG_IsOK(ecode12)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode12), "in method '" "grdcp3d_from_cube" "', argument " "12"" of type '" "double""'");
-  } 
-  arg12 = (double)(val12);
-  ecode13 = SWIG_AsVal_double(swig_obj[9], &val13);
-  if (!SWIG_IsOK(ecode13)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode13), "in method '" "grdcp3d_from_cube" "', argument " "13"" of type '" "double""'");
-  } 
-  arg13 = (double)(val13);
-  ecode14 = SWIG_AsVal_double(swig_obj[10], &val14);
-  if (!SWIG_IsOK(ecode14)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode14), "in method '" "grdcp3d_from_cube" "', argument " "14"" of type '" "double""'");
-  } 
-  arg14 = (double)(val14);
-  ecode15 = SWIG_AsVal_double(swig_obj[11], &val15);
-  if (!SWIG_IsOK(ecode15)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode15), "in method '" "grdcp3d_from_cube" "', argument " "15"" of type '" "double""'");
-  } 
-  arg15 = (double)(val15);
-  ecode16 = SWIG_AsVal_double(swig_obj[12], &val16);
-  if (!SWIG_IsOK(ecode16)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode16), "in method '" "grdcp3d_from_cube" "', argument " "16"" of type '" "double""'");
-  } 
-  arg16 = (double)(val16);
-  ecode17 = SWIG_AsVal_int(swig_obj[13], &val17);
-  if (!SWIG_IsOK(ecode17)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode17), "in method '" "grdcp3d_from_cube" "', argument " "17"" of type '" "int""'");
-  } 
-  arg17 = (int)(val17);
-  ecode18 = SWIG_AsVal_int(swig_obj[14], &val18);
-  if (!SWIG_IsOK(ecode18)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode18), "in method '" "grdcp3d_from_cube" "', argument " "18"" of type '" "int""'");
-  } 
-  arg18 = (int)(val18);
-  {
-    char *err;
-    clear_exception();
-    grdcp3d_from_cube(arg1,arg2,arg3,arg4,arg5,arg6,arg7,arg8,arg9,arg10,arg11,arg12,arg13,arg14,arg15,arg16,arg17,arg18);
-    if ((err = check_exception())) {
-      PyErr_SetString(PY_XTGeoCLibError, err);
-      return NULL;
-    }
-  }
-  resultobj = SWIG_Py_Void();
-  return resultobj;
-fail:
-  return NULL;
-}
-
-
 SWIGINTERN PyObject *_wrap_grd3cp3d_xtgformat1to2_geom(PyObject *self, PyObject *args) {
   PyObject *resultobj = 0;
   long arg1 ;
@@ -19490,143 +19336,6 @@ fail:
 }
 
 
-SWIGINTERN PyObject *_wrap_grdcp3d_conv_grid_roxapi(PyObject *self, PyObject *args) {
-  PyObject *resultobj = 0;
-  long arg1 ;
-  long arg2 ;
-  long arg3 ;
-  double *arg4 = (double *) 0 ;
-  long arg5 ;
-  float *arg6 = (float *) 0 ;
-  long arg7 ;
-  double *arg8 = (double *) 0 ;
-  long arg9 ;
-  double *arg10 = (double *) 0 ;
-  long arg11 ;
-  double *arg12 = (double *) 0 ;
-  long arg13 ;
-  long val1 ;
-  int ecode1 = 0 ;
-  long val2 ;
-  int ecode2 = 0 ;
-  long val3 ;
-  int ecode3 = 0 ;
-  PyArrayObject *array4 = NULL ;
-  int i4 = 1 ;
-  PyArrayObject *array6 = NULL ;
-  int i6 = 1 ;
-  PyObject *array8 = NULL ;
-  PyObject *array10 = NULL ;
-  PyObject *array12 = NULL ;
-  PyObject *swig_obj[8] ;
-  int result;
-  
-  (void)self;
-  if (!SWIG_Python_UnpackTuple(args, "grdcp3d_conv_grid_roxapi", 8, 8, swig_obj)) SWIG_fail;
-  ecode1 = SWIG_AsVal_long(swig_obj[0], &val1);
-  if (!SWIG_IsOK(ecode1)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode1), "in method '" "grdcp3d_conv_grid_roxapi" "', argument " "1"" of type '" "long""'");
-  } 
-  arg1 = (long)(val1);
-  ecode2 = SWIG_AsVal_long(swig_obj[1], &val2);
-  if (!SWIG_IsOK(ecode2)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "grdcp3d_conv_grid_roxapi" "', argument " "2"" of type '" "long""'");
-  } 
-  arg2 = (long)(val2);
-  ecode3 = SWIG_AsVal_long(swig_obj[2], &val3);
-  if (!SWIG_IsOK(ecode3)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode3), "in method '" "grdcp3d_conv_grid_roxapi" "', argument " "3"" of type '" "long""'");
-  } 
-  arg3 = (long)(val3);
-  {
-    array4 = obj_to_array_no_conversion(swig_obj[3], NPY_DOUBLE);
-    if (!array4 || !require_c_or_f_contiguous(array4)
-      || !require_native(array4)) SWIG_fail;
-    arg4 = (double*) array_data(array4);
-    arg5 = 1;
-    for (i4=0; i4 < array_numdims(array4); ++i4) arg5 *= array_size(array4,i4);
-  }
-  {
-    array6 = obj_to_array_no_conversion(swig_obj[4], NPY_FLOAT);
-    if (!array6 || !require_c_or_f_contiguous(array6)
-      || !require_native(array6)) SWIG_fail;
-    arg6 = (float*) array_data(array6);
-    arg7 = 1;
-    for (i6=0; i6 < array_numdims(array6); ++i6) arg7 *= array_size(array6,i6);
-  }
-  {
-    npy_intp dims[1];
-    if (!PyInt_Check(swig_obj[5]))
-    {
-      const char* typestring = pytype_string(swig_obj[5]);
-      PyErr_Format(PyExc_TypeError,
-        "Int dimension expected.  '%s' given.",
-        typestring);
-      SWIG_fail;
-    }
-    arg9 = (long) PyInt_AsLong(swig_obj[5]);
-    dims[0] = (npy_intp) arg9;
-    array8 = PyArray_SimpleNew(1, dims, NPY_DOUBLE);
-    if (!array8) SWIG_fail;
-    arg8 = (double*) array_data(array8);
-  }
-  {
-    npy_intp dims[1];
-    if (!PyInt_Check(swig_obj[6]))
-    {
-      const char* typestring = pytype_string(swig_obj[6]);
-      PyErr_Format(PyExc_TypeError,
-        "Int dimension expected.  '%s' given.",
-        typestring);
-      SWIG_fail;
-    }
-    arg11 = (long) PyInt_AsLong(swig_obj[6]);
-    dims[0] = (npy_intp) arg11;
-    array10 = PyArray_SimpleNew(1, dims, NPY_DOUBLE);
-    if (!array10) SWIG_fail;
-    arg10 = (double*) array_data(array10);
-  }
-  {
-    npy_intp dims[1];
-    if (!PyInt_Check(swig_obj[7]))
-    {
-      const char* typestring = pytype_string(swig_obj[7]);
-      PyErr_Format(PyExc_TypeError,
-        "Int dimension expected.  '%s' given.",
-        typestring);
-      SWIG_fail;
-    }
-    arg13 = (long) PyInt_AsLong(swig_obj[7]);
-    dims[0] = (npy_intp) arg13;
-    array12 = PyArray_SimpleNew(1, dims, NPY_DOUBLE);
-    if (!array12) SWIG_fail;
-    arg12 = (double*) array_data(array12);
-  }
-  {
-    char *err;
-    clear_exception();
-    result = (int)grdcp3d_conv_grid_roxapi(arg1,arg2,arg3,arg4,arg5,arg6,arg7,arg8,arg9,arg10,arg11,arg12,arg13);
-    if ((err = check_exception())) {
-      PyErr_SetString(PY_XTGeoCLibError, err);
-      return NULL;
-    }
-  }
-  resultobj = SWIG_From_int((int)(result));
-  {
-    resultobj = SWIG_Python_AppendOutput(resultobj,(PyObject*)array8);
-  }
-  {
-    resultobj = SWIG_Python_AppendOutput(resultobj,(PyObject*)array10);
-  }
-  {
-    resultobj = SWIG_Python_AppendOutput(resultobj,(PyObject*)array12);
-  }
-  return resultobj;
-fail:
-  return NULL;
-}
-
-
 SWIGINTERN PyObject *_wrap_grdcp3d_quality_indicators(PyObject *self, PyObject *args) {
   PyObject *resultobj = 0;
   long arg1 ;
@@ -20627,14 +20336,12 @@ static PyMethodDef SwigMethods[] = {
 	 { "grd3d_zminmax", _wrap_grd3d_zminmax, METH_VARARGS, NULL},
 	 { "grd3d_get_all_corners", _wrap_grd3d_get_all_corners, METH_VARARGS, NULL},
 	 { "grd3d_well_ijk", _wrap_grd3d_well_ijk, METH_VARARGS, NULL},
-	 { "grdcp3d_from_cube", _wrap_grdcp3d_from_cube, METH_VARARGS, NULL},
 	 { "grd3cp3d_xtgformat1to2_geom", _wrap_grd3cp3d_xtgformat1to2_geom, METH_VARARGS, NULL},
 	 { "grd3cp3d_xtgformat2to1_geom", _wrap_grd3cp3d_xtgformat2to1_geom, METH_VARARGS, NULL},
 	 { "grdcp3d_process_edges", _wrap_grdcp3d_process_edges, METH_VARARGS, NULL},
 	 { "grdcp3d_corners", _wrap_grdcp3d_corners, METH_VARARGS, NULL},
 	 { "grdcp3d_get_vtk_esg_geometry_data", _wrap_grdcp3d_get_vtk_esg_geometry_data, METH_VARARGS, NULL},
 	 { "grdcp3d_get_vtk_grid_arrays", _wrap_grdcp3d_get_vtk_grid_arrays, METH_VARARGS, NULL},
-	 { "grdcp3d_conv_grid_roxapi", _wrap_grdcp3d_conv_grid_roxapi, METH_VARARGS, NULL},
 	 { "grdcp3d_quality_indicators", _wrap_grdcp3d_quality_indicators, METH_VARARGS, NULL},
 	 { "well_geometrics", _wrap_well_geometrics, METH_VARARGS, NULL},
 	 { "well_trunc_parallel", _wrap_well_trunc_parallel, METH_VARARGS, NULL},

xtgeo/__init__.py

@@ -142,6 +142,7 @@ from xtgeo.grid3d.grid import (
     grid_from_cube,
     grid_from_file,
     grid_from_roxar,
+    grid_from_surfaces,
 )
 from xtgeo.xyz.polygons import (
     polygons_from_file,
@@ -210,6 +211,7 @@ __all__ = [
     "grid_from_cube",
     "grid_from_file",
     "grid_from_roxar",
+    "grid_from_surfaces",
     "grid_properties",
     "grid_properties",
     "grid_property",

xtgeo/common/_xyz_enum.py

@@ -14,6 +14,7 @@ class _AttrName(ExtendedEnum):
     XNAME = "X_UTME"
     YNAME = "Y_UTMN"
     ZNAME = "Z_TVDSS"
+    PNAME = "POLY_ID"
     M_MD_NAME = "M_MDEPTH"
     Q_MD_NAME = "Q_MDEPTH"
     M_AZI_NAME = "M_AZI"
@@ -24,6 +25,12 @@ class _AttrName(ExtendedEnum):
     J_INDEX = "J_INDEX"
     K_INDEX = "K_INDEX"
     R_HLEN_NAME = "R_HLEN"
+    HNAME = "H_CUMLEN"
+    DHNAME = "H_DELTALEN"
+    TNAME = "T_CUMLEN"
+    DTNAME = "T_DELTALEN"
+    WELLNAME = "WELLNAME"
+    TRAJECTORY = "TRAJECTORY"
 
 
 @unique
@@ -38,6 +45,6 @@ class _AttrType(ExtendedEnum):
 class _XYZType(ExtendedEnum):
     """Enumerate type of context"""
 
-    POINTS = 1
-    POLYGONS = 2  # ie. same here as PolyLines
-    WELL = 3
+    POINTS = "POINTS"
+    POLYGONS = "POLYGONS"
+    WELL = "WELL"

xtgeo/common/version.py

@@ -17,5 +17,5 @@ __version__: str
 __version_tuple__: VERSION_TUPLE
 version_tuple: VERSION_TUPLE
 
-__version__ = version = '4.6.0'
-__version_tuple__ = version_tuple = (4, 6, 0)
+__version__ = version = '4.7.0'
+__version_tuple__ = version_tuple = (4, 7, 0)

xtgeo/grid3d/_ecl_inte_head.py

@@ -144,5 +144,5 @@ class InteHead:
         return (
             None
             if value is None
-            else TypeOfGrid.alternate_code(cast(Literal[0, 1, 2, 3], value))
+            else TypeOfGrid.alternate_code(cast("Literal[0, 1, 2, 3]", value))
         )

xtgeo/grid3d/_find_gridprop_in_eclrun.py

@@ -90,7 +90,7 @@ def remainder_saturations(saturations):
     rest = sum(saturations.values())
     if len(saturations) == 2 or np.allclose(rest, 1.0):
         missing = set(sat_keys).difference(set(saturations.keys()))
-        return {m: 1.0 - rest for m in missing}
+        return dict.fromkeys(missing, 1.0 - rest)
     return {}
 
 

xtgeo/grid3d/_grid_etc1.py

@@ -16,6 +16,7 @@ from xtgeo import _cxtgeo
 from xtgeo.common.calc import find_flip
 from xtgeo.common.constants import UNDEF_INT, UNDEF_LIMIT
 from xtgeo.common.log import null_logger
+from xtgeo.common.types import Dimensions
 from xtgeo.grid3d.grid_properties import GridProperties
 from xtgeo.xyz.polygons import Polygons
 
@@ -24,10 +25,10 @@ from ._grid3d_fence import _update_tmpvars
 from .grid_property import GridProperty
 
 if TYPE_CHECKING:
-    from xtgeo.common.types import Dimensions
     from xtgeo.grid3d import Grid
     from xtgeo.grid3d.types import METRIC
     from xtgeo.surface.regular_surface import RegularSurface
+    from xtgeo.surface.surfaces import Surfaces
     from xtgeo.xyz.points import Points
 
 logger = null_logger(__name__)
@@ -42,6 +43,9 @@ def create_box(
     flip: Literal[1, -1],
 ) -> dict[str, np.ndarray]:
     """Create a shoebox grid from cubi'sh spec, xtgformat=2."""
+
+    from xtgeo.cube.cube1 import Cube
+
     ncol, nrow, nlay = dimension
     nncol = ncol + 1
     nnrow = nrow + 1
@@ -51,33 +55,100 @@ def create_box(
     zcornsv = np.zeros((nncol, nnrow, nnlay, 4), dtype=np.float32)
     actnumsv = np.zeros((ncol, nrow, nlay), dtype=np.int32)
 
-    option = 0 if not oricenter else 1
-
-    _cxtgeo.grdcp3d_from_cube(
-        ncol,
-        nrow,
-        nlay,
-        coordsv,
-        zcornsv,
-        actnumsv,
-        origin[0],
-        origin[1],
-        origin[2],
-        increment[0],
-        increment[1],
-        increment[2],
-        rotation,
-        flip,
-        option,
+    cube = Cube(
+        ncol=ncol,
+        nrow=nrow,
+        nlay=nlay,
+        xinc=increment[0],
+        yinc=increment[1],
+        zinc=increment[2],
+        xori=origin[0],
+        yori=origin[1],
+        zori=origin[2],
+        rotation=rotation,
     )
 
+    cubecpp = _internal.cube.Cube(cube)
+    coordsv, zcornsv, actnumsv = _internal.grid3d.create_grid_from_cube(
+        cubecpp, oricenter, flip
+    )
     return {
         "coordsv": coordsv,
         "zcornsv": zcornsv,
-        "actnumsv": actnumsv,
+        "actnumsv": actnumsv.astype(np.int32),
     }
 
 
+def create_grid_from_surfaces(
+    srfs: Surfaces,
+    ij_dimension: tuple[int, int] | None = None,
+    ij_origin: tuple[float, float] | None = None,
+    ij_increment: tuple[float, float] | None = None,
+    rotation: float | None = None,
+    tolerance: float = _internal.numerics.TOLERANCE,
+) -> Grid:
+    """Use a stack of surfaces to create a nonfaulted grid.
+
+    Technically, a shoebox grid is made first, then the layers are adjusted to follow
+    surfaces.
+    """
+    from xtgeo.grid3d.grid import create_box_grid
+
+    n_surfaces = len(srfs.surfaces)
+
+    # ensure that surfaces are consistent
+    if not srfs.is_depth_consistent():
+        raise ValueError(
+            "Surfaces are not depth consistent, they must not cross is depth"
+        )
+    top = srfs.surfaces[0]
+    base = srfs.surfaces[-1]
+
+    zinc = (base.values.mean() - top.values.mean()) / (n_surfaces - 1)
+    kdim: int = n_surfaces - 1
+    zori = top.values.mean()
+    ncol: int = top.ncol - 1  # since surface are nodes while grid is cell centered
+    nrow: int = top.nrow - 1
+
+    if ij_dimension:  # mypy needs this:
+        dimension = Dimensions(int(ij_dimension[0]), int(ij_dimension[1]), kdim)
+    else:
+        dimension = Dimensions(ncol, nrow, kdim)
+
+    increment = (*ij_increment, zinc) if ij_increment else (top.xinc, top.yinc, zinc)
+    origin = (*ij_origin, zori) if ij_origin else (top.xori, top.yori, zori)
+    rotation = rotation if rotation is not None else top.rotation
+
+    grd = create_box_grid(
+        dimension=dimension,
+        origin=origin,
+        increment=increment,
+        rotation=rotation,
+        oricenter=False,
+        flip=1,
+    )
+
+    # now adjust the grid to surfaces
+    surf_list = []
+    for surf in srfs.surfaces:
+        cpp_surf = _internal.regsurf.RegularSurface(surf)
+        surf_list.append(cpp_surf)
+
+    grd_cpp = _internal.grid3d.Grid(grd)
+    new_zcorns, new_actnum = grd_cpp.adjust_boxgrid_layers_from_regsurfs(
+        surf_list, tolerance
+    )
+
+    grd._zcornsv = new_zcorns
+    grd._actnumsv = new_actnum
+
+    # set the subgrid index (zones)
+    subgrids = {f"zone_{i + 1}": 1 for i in range(n_surfaces - 1)}
+    grd.set_subgrids(subgrids)
+
+    return grd
+
+
 method_factory = {
     "euclid": _cxtgeo.euclid_length,
     "horizontal": _cxtgeo.horizontal_length,
@@ -1185,7 +1256,7 @@ def reduce_to_one_layer(self: Grid) -> None:
 
 def translate_coordinates(
     self: Grid,
-    translate: tuple[int, int, int] = (0, 0, 0),
+    translate: tuple[float, float, float] = (0.0, 0.0, 0.0),
     flip: tuple[int, int, int] = (1, 1, 1),
 ) -> None:
     """Translate grid coordinates."""