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