Merge pull request #2105 from GitPaean/aqunum_activating_cells

AQUNUM can activate inactive cells
This commit is contained in:
Joakim Hove
2020-11-13 17:23:14 +01:00
committed by GitHub
3 changed files with 93 additions and 2 deletions

View File

@@ -2148,6 +2148,69 @@ BOOST_AUTO_TEST_CASE(TESTCP_ACTNUM_UPDATE) {
}
BOOST_AUTO_TEST_CASE(TESTCP_ACTNUM_AQUNUM) {
const std::string deckData = R"(
RUNSPEC
DIMENS
3 2 1 /
GRID
COORD
2000.0000 2000.0000 2000.0000 1999.9476 2000.0000 2002.9995
2049.9924 2000.0000 2000.8726 2049.9400 2000.0000 2003.8722
2099.9848 2000.0000 2001.7452 2099.9324 2000.0000 2004.7448
2149.9772 2000.0000 2002.6179 2149.9248 2000.0000 2005.6174
2000.0000 2050.0000 2000.0000 1999.9476 2050.0000 2002.9995
2049.9924 2050.0000 2000.8726 2049.9400 2050.0000 2003.8722
2099.9848 2050.0000 2001.7452 2099.9324 2050.0000 2004.7448
2149.9772 2050.0000 2002.6179 2149.9248 2050.0000 2005.6174
2000.0000 2100.0000 2000.0000 1999.9476 2100.0000 2002.9995
2049.9924 2100.0000 2000.8726 2049.9400 2100.0000 2003.8722
2099.9848 2100.0000 2001.7452 2099.9324 2100.0000 2004.7448
2149.9772 2100.0000 2002.6179 2149.9248 2100.0000 2005.6174 /
ZCORN
2000.0000 2000.8726 2000.8726 2001.7452 2001.7452 2002.6179
2000.0000 2000.8726 2000.8726 2001.7452 2001.7452 2002.6179
2000.0000 2000.8726 2000.8726 2001.7452 2001.7452 2002.6179
2000.0000 2000.8726 2000.8726 2001.7452 2001.7452 2002.6179
2002.9995 2003.8722 2003.8722 2004.7448 2004.7448 2005.6174
2002.9995 2003.8722 2003.8722 2004.7448 2004.7448 2005.6174
2002.9995 2003.8722 2003.8722 2004.7448 2004.7448 2005.6174
2002.9995 2003.8722 2003.8722 2004.7448 2004.7448 2005.6174 /
ACTNUM
0 0 1 1 0 1 /
PORO
6*0.15 /
AQUNUM
1 1 1 1 1000000.0 10000 0.25 400 2585.00 285.00 1 1 /
1 2 1 1 1000000.0 10000 0.25 400 2585.00 285.00 1 1 /
/
AQUCON
1 3 3 2 2 1 1 'J+' 1.00 1 /
/
EDIT)";
Opm::Parser parser;
const auto deck = parser.parseString( deckData) ;
Opm::EclipseGrid grid( deck);
const std::vector<int>& grid_actnum = grid.getACTNUM();
const std::vector<int> desired_actnum = {1, 1, 1, 1, 0, 1};
BOOST_CHECK( grid_actnum.size() == desired_actnum.size() );
for (size_t n=0; n< grid_actnum.size(); n++) {
BOOST_CHECK_EQUAL( grid_actnum[n], desired_actnum[n] );
}
Opm::EclipseState es(deck);
const auto& grid_actnum2 = es.getInputGrid().getACTNUM();
BOOST_CHECK( desired_actnum.size() == grid_actnum2.size());
for (size_t n=0; n< grid_actnum.size(); n++) {
BOOST_CHECK_EQUAL( grid_actnum2[n], desired_actnum[n] );
}
}
BOOST_AUTO_TEST_CASE(TEST_altGridConstructors) {
const char* deckData =