mirror of
https://github.com/OPM/ResInsight.git
synced 2026-10-14 17:32:13 -05:00
Jacob Støren
f5d1d2d1f9
To include phase split well connection rates and some summary intersect support. Supposedly Final update before release.
58 lines
1.9 KiB
Python
58 lines
1.9 KiB
Python
#!/usr/bin/env python
|
|
# Copyright (C) 2015 Statoil ASA, Norway.
|
|
#
|
|
# The file 'test_kw_function.py' is part of ERT - Ensemble based Reservoir Tool.
|
|
#
|
|
# ERT is free software: you can redistribute it and/or modify
|
|
# it under the terms of the GNU General Public License as published by
|
|
# the Free Software Foundation, either version 3 of the License, or
|
|
# (at your option) any later version.
|
|
#
|
|
# ERT is distributed in the hope that it will be useful, but WITHOUT ANY
|
|
# WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
|
# FITNESS FOR A PARTICULAR PURPOSE.
|
|
#
|
|
# See the GNU General Public License at <http://www.gnu.org/licenses/gpl.html>
|
|
# for more details.
|
|
import os
|
|
import random
|
|
from ecl.ecl import EclKW, EclDataType, EclGrid , Ecl3DKW
|
|
from ecl.util import IntVector
|
|
from ecl.test import ExtendedTestCase
|
|
|
|
class KWFunctionTest(ExtendedTestCase):
|
|
|
|
def test_region_filter(self):
|
|
nx = 10
|
|
ny = 10
|
|
nz = 1
|
|
actnum = IntVector( initial_size = nx*ny*nz , default_value = 1 )
|
|
actnum[nx*ny - 1] = 0
|
|
|
|
grid = EclGrid.createRectangular( (nx,ny,nz) , (1,1,1) , actnum = actnum)
|
|
self.assertEqual( grid.getNumActive() , nx*ny*nz - 1 )
|
|
|
|
kw = Ecl3DKW.create( "REGIONS" , grid , EclDataType.ECL_INT , global_active = True )
|
|
kw.assign( 0 )
|
|
kw[0:nx*ny/2] = 1
|
|
kw[5,2,0] = 0
|
|
kw[0,9,0] = 2
|
|
|
|
kw.fixUninitialized( grid )
|
|
|
|
# Not assigned because they are in contact with a '2'; these
|
|
# two are problem cells.
|
|
self.assertEqual( kw[0,ny - 2,0] , 0)
|
|
self.assertEqual( kw[1,ny - 1,0] , 0)
|
|
|
|
# Not assigned because it is inactive
|
|
self.assertEqual( kw[nx - 1,ny - 1,0] , 0)
|
|
|
|
self.assertEqual( kw[5,2,0] , 1 )
|
|
for j in range(5,10):
|
|
self.assertEqual( kw[5,j,0] , 1 )
|
|
|
|
for i in range(10):
|
|
self.assertEqual( kw[i,7,0] , 1 )
|
|
|