def export_zmap_ascii(self, mfile):
"""Export to ZMAP ascii format (non-rotated)."""
# zmap can only deal with non-rotated formats; hence make a copy
# of the instance and derotate that prior to export, so that the
# original instance is unchanged
fout = xtgeo._XTGeoCFile(mfile, mode="wb")
scopy = self.copy()
if abs(scopy.rotation) > 1.0e-20:
scopy.unrotate()
zmin = scopy.values.min()
zmax = scopy.values.max()
yinc = scopy._yinc * scopy._yflip
vals = scopy.get_values1d(order="F",
asmasked=False,
fill_value=xtgeo.UNDEF)
ier = _cxtgeo.surf_export_zmap_ascii(
fout.fhandle,
scopy._ncol,
scopy._nrow,
scopy._xori,
scopy._yori,
scopy._xinc,
yinc,
vals,
zmin,
zmax,
0,
)
if ier != 0:
raise RuntimeError("Export to ZMAP Ascii went wrong, "
"code is {}".format(ier))
del scopy
fout.close()
def _export_zmap_ascii(self, mfile):
"""Export to ZMAP ascii format (non-rotated)."""
scopy = self.copy()
if abs(scopy.rotation) > 1.0e-20:
scopy.unrotate()
zmin = scopy.values.min()
zmax = scopy.values.max()
yinc = scopy._yinc * scopy._yflip
vals = scopy.get_values1d(order="C",
asmasked=False,
fill_value=xtgeo.UNDEF)
ier = _cxtgeo.surf_export_zmap_ascii(
mfile.get_cfhandle(),
scopy._ncol,
scopy._nrow,
scopy._xori,
scopy._yori,
scopy._xinc,
yinc,
vals,
zmin,
zmax,
0,
)
if ier != 0:
raise RuntimeError("Export to ZMAP Ascii went wrong, "
"code is {}".format(ier))
del scopy
mfile.cfclose()