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,