def get_dz(self, name="dZ", flip=True, asmasked=True):
"""Get dZ as property"""
self._xtgformat1()
ntot = (self._ncol, self._nrow, self._nlay)
dzv = GridProperty(
ncol=self._ncol,
nrow=self._nrow,
nlay=self._nlay,
values=np.zeros(ntot, dtype=np.float64),
name=name,
discrete=False,
)
dz = np.zeros(self.ntotal, dtype=np.float64)
nflip = 1
if not flip:
nflip = -1
option = 0
if asmasked:
option = 1
_cxtgeo.grd3d_calc_dz(
self._ncol,
self._nrow,
self._nlay,
self._zcornsv,
self._actnumsv,
dz,
nflip,
option,
)
dzv.values = np.ma.masked_greater(dz, xtgeo.UNDEF_LIMIT)
# return the property object
logger.info("DZ mean value: %s", dzv.values.mean())
return dzv
def get_dz(self, name="dZ", flip=True, asmasked=True):
"""Get dZ as property"""
ntot = (self._ncol, self._nrow, self._nlay)
dzv = GridProperty(
ncol=self._ncol,
nrow=self._nrow,
nlay=self._nlay,
values=np.zeros(ntot, dtype=np.float64),
name=name,
discrete=False,
)
ptr_dz_v = _cxtgeo.new_doublearray(self.ntotal)
nflip = 1
if not flip:
nflip = -1
option = 0
if asmasked:
option = 1
_cxtgeo.grd3d_calc_dz(
self._ncol,
self._nrow,
self._nlay,
self._p_zcorn_v,
self._p_actnum_v,
ptr_dz_v,
nflip,
option,
XTGDEBUG,
)
_gridprop_lowlevel.update_values_from_carray(dzv, ptr_dz_v, np.float64, delete=True)
# return the property object
logger.info("DZ mean value: %s", dzv.values.mean())
return dzv