def test_slice(self):
reg = ecl.EclRegion(self.grid, False)
reg.select_islice(0, 5)
OK = True
for gi in reg.global_list:
(i, j, k) = self.grid.get_ijk(global_index=gi)
if i > 5:
OK = False
self.assertTrue(OK)
self.assertTrue(self.grid.ny * self.grid.nz *
6 == len(reg.global_list))
reg.select_jslice(7, 8, intersect=True)
OK = True
for gi in reg.global_list:
(i, j, k) = self.grid.get_ijk(global_index=gi)
if i > 5:
OK = False
if j < 7 or j > 8:
OK = False
self.assertTrue(OK)
self.assertTrue(2 * self.grid.nz * 6 == len(reg.global_list))
reg2 = ecl.EclRegion(self.grid, False)
reg2.select_kslice(3, 5)
reg &= reg2
OK = True
for gi in reg.global_list:
(i, j, k) = self.grid.get_ijk(global_index=gi)
if i > 5:
OK = False
if j < 7 or j > 8:
OK = False
if k < 3 or k > 5:
OK = False
self.assertTrue(OK)
self.assertTrue(2 * 3 * 6 == len(reg.global_list))
def test_kw_isub(self):
P = self.rst_file["PRESSURE"][5]
fipnum = self.init_file["FIPNUM"][0]
fipnum_copy = fipnum.deep_copy()
reg = ecl.EclRegion(self.grid, False)
reg.select_more(P, 260)
fipnum.sub(1, mask=reg)
self.assertFalse(fipnum.equal(fipnum_copy))
fipnum.add(1, mask=reg)
self.assertTrue(fipnum.equal(fipnum_copy))
#!/usr/bin/env python
import ert.ecl.ecl as ecl
g = ecl.EclGrid("data/eclipse/case/ECLIPSE")
r = ecl.EclRegion(g, False)
ll = g.get_xyz(ijk=(0, 0, 0))
ur = g.get_xyz(ijk=(39, 63, 0))
dx = 0.5 * (ur[0] - ll[0])
dy = 0.5 * (ur[1] - ll[1])
x0 = ll[0] + dx * 0.5
y0 = ll[1] + dy * 0.5
x1 = x0 + dx
y1 = y0 + dy
points = [(x0, y0), (x1, y0), (x1, y1), (x0, y1)]
print " X ----------------------------------- (%10.3f , %10.3f)" % (
x1 - x0, y1 - y0)
print " | |"
print " | |"
print " | |"
print " | |"
print " | |"
print " | |"
print " | |"
print " | |"
print " | |"
print "(%10.3f , %10.3f) ---------------------------------- X\n\n" % (x0 - x0,