def test_report_list_file(self, fname, rlist0):
rlist = ecl.EclFile.file_report_list(fname)
self.assertTrue(approx_equalv(rlist, rlist0))
f = ecl.EclFile(fname)
rlist = f.report_list
self.assertTrue(approx_equalv(rlist, rlist0))
def testRect(self):
a1 = 1.0
a2 = 2.0
a3 = 3.0
grid = ecl.EclGrid.create_rectangular((9,9,9) , (a1,a2,a3))
grid.save_EGRID( "/tmp/rect.EGRID" )
grid2 = ecl.EclGrid( "/tmp/rect.EGRID")
self.addFile( "/tmp/rect.EGRID" )
self.assertTrue( grid )
self.assertTrue( grid2 )
(x,y,z) = grid.get_xyz( ijk=(4,4,4) )
self.assertTrue( approx_equalv( [x,y,z],[4.5 * a1,4.5*a2,4.5*a3] ))
v = grid.cell_volume( ijk=(4,4,4) )
self.assertTrue( approx_equal( v , a1*a2*a3 ))
z = grid.depth( ijk = (4,4,4 ))
self.assertTrue( approx_equal( z , 4.5*a3 ))
g1 = grid.global_index( ijk = (2,2,2) )
g2 = grid.global_index( ijk = (4,4,4) )
(dx,dy,dz) = grid.distance( g2 , g1 )
self.assertTrue( approx_equalv([dx,dy,dz],[2*a1,2*a2,2*a3] ) )
self.assertTrue( grid.cell_contains(2.5*a1 , 2.5*a2, 2.5*a3 , ijk=(2,2,2)))
ijk = grid.find_cell( 1.5*a1 , 2.5*a2, 3.5*a3 )
self.assertTrue( approx_equalv( ijk , [1 , 2 , 3]))
def testRect(self):
a1 = 1.0
a2 = 2.0
a3 = 3.0
grid = ecl.EclGrid.create_rectangular((9, 9, 9), (a1, a2, a3))
grid.save_EGRID("/tmp/rect.EGRID")
grid2 = ecl.EclGrid("/tmp/rect.EGRID")
self.addFile("/tmp/rect.EGRID")
self.assertTrue(grid)
self.assertTrue(grid2)
(x, y, z) = grid.get_xyz(ijk=(4, 4, 4))
self.assertTrue(
approx_equalv([x, y, z], [4.5 * a1, 4.5 * a2, 4.5 * a3]))
v = grid.cell_volume(ijk=(4, 4, 4))
self.assertTrue(approx_equal(v, a1 * a2 * a3))
z = grid.depth(ijk=(4, 4, 4))
self.assertTrue(approx_equal(z, 4.5 * a3))
g1 = grid.global_index(ijk=(2, 2, 2))
g2 = grid.global_index(ijk=(4, 4, 4))
(dx, dy, dz) = grid.distance(g2, g1)
self.assertTrue(approx_equalv([dx, dy, dz], [2 * a1, 2 * a2, 2 * a3]))
self.assertTrue(
grid.cell_contains(2.5 * a1, 2.5 * a2, 2.5 * a3, ijk=(2, 2, 2)))
ijk = grid.find_cell(1.5 * a1, 2.5 * a2, 3.5 * a3)
self.assertTrue(approx_equalv(ijk, [1, 2, 3]))
def test_report_list_file(self , fname , rlist0):
rlist = ecl.EclFile.file_report_list( fname )
self.assertTrue( approx_equalv( rlist , rlist0 ))
f = ecl.EclFile( fname )
rlist = f.report_list
self.assertTrue( approx_equalv( rlist , rlist0 ))
def test_interp(self):
sum = self.sum
self.assertTrue( approx_equal( sum.get_interp( "WWCT:OP_3" , days = 750 ) , 0.11719122))
self.assertTrue( approx_equal( sum.get_interp( "WWCT:OP_3" , date = datetime.date( 2004,1,1)) , 0.603358387947 ))
v = sum.get_interp_vector( "WOPT:OP_1" , days_list = [100 , 200 , 400 , 800] )
self.assertTrue( approx_equalv( [805817.11875 , 1614955.34677419 , 3289267.67857143 , 6493021.6218035 ] , v))
v = sum.get_interp_vector( "WGPT:OP_2" , date_list = [datetime.date(2002,1,1) , datetime.date(2003,1,1) , datetime.date(2004 , 1 , 1)])
self.assertTrue( approx_equalv( v , [ 8.20773632e+08 , 9.68444032e+08 , 1.02515213e+09]) )
def test_wells(self):
sum = self.sum
wells = sum.wells()
wells.sort()
self.assertTrue( approx_equalv( wells , ["OP_1" , "OP_2" , "OP_3" , "OP_4" , "OP_5" , "WI_1" , "WI_2" , "WI_3"]))
wells = sum.wells( pattern = "*_3")
wells.sort()
self.assertTrue( approx_equalv( wells , ["OP_3" , "WI_3"]))
groups = sum.groups()
groups.sort()
self.assertTrue( approx_equalv( groups , ['GMWIN','OP','WI']))
def test_wells(self):
sum = self.sum
wells = sum.wells()
wells.sort()
self.assertTrue(
approx_equalv(wells, [
"OP_1", "OP_2", "OP_3", "OP_4", "OP_5", "WI_1", "WI_2", "WI_3"
]))
wells = sum.wells(pattern="*_3")
wells.sort()
self.assertTrue(approx_equalv(wells, ["OP_3", "WI_3"]))
groups = sum.groups()
groups.sort()
self.assertTrue(approx_equalv(groups, ['GMWIN', 'OP', 'WI']))
def test_stringlist_reference(self):
sum = ecl.EclSum(case)
wells = sum.wells()
self.assertTrue(
approx_equalv(wells, [
'OP_1', 'OP_2', 'OP_3', 'OP_4', 'OP_5', 'WI_1', 'WI_2', 'WI_3'
]))
self.assertTrue(isinstance(wells, StringList))
def test_interp(self):
sum = self.sum
self.assertTrue(
approx_equal(sum.get_interp("WWCT:OP_3", days=750), 0.11719122))
self.assertTrue(
approx_equal(
sum.get_interp("WWCT:OP_3", date=datetime.date(2004, 1, 1)),
0.603358387947))
v = sum.get_interp_vector("WOPT:OP_1", days_list=[100, 200, 400, 800])
self.assertTrue(
approx_equalv([
805817.11875, 1614955.34677419, 3289267.67857143,
6493021.6218035
], v))
v = sum.get_interp_vector("WGPT:OP_2",
date_list=[
datetime.date(2002, 1, 1),
datetime.date(2003, 1, 1),
datetime.date(2004, 1, 1)
])
self.assertTrue(
approx_equalv(v, [8.20773632e+08, 9.68444032e+08, 1.02515213e+09]))
def test_stringlist_reference(self):
sum = ecl.EclSum( case )
wells = sum.wells()
self.assertTrue( approx_equalv( wells , ['OP_1','OP_2','OP_3','OP_4','OP_5','WI_1','WI_2','WI_3']))
self.assertTrue( isinstance( wells , StringList ))
