Include WPIMULT for CTFAC output

This commit is contained in:
Lars Petter Øren Hauge 2018-10-29 09:12:50 +01:00
parent 7ebadbb055
commit d958da6542
3 changed files with 20 additions and 4 deletions

View File

@ -442,8 +442,8 @@ inline quantity srate( const fn_args& args ) {
const auto& name = well->name();
if( args.wells.count( name ) == 0 ) return zero;
const auto& well_data = args.wells.at( name );
const auto& well_data = args.wells.at( name );
const auto& segment = well_data.segments.find(segNumber);
if( segment == well_data.segments.end() ) return zero;
@ -485,7 +485,7 @@ inline quantity trans_factors ( const fn_args& args ) {
if( connection == connections.end() ) return zero;
const auto& v = connection->CF();
const auto& v = connection->CF() * connection->wellPi();
return { v, measure::transmissibility };
}

View File

@ -675,5 +675,16 @@ COMPDAT
W_4 1 1 3 3 /
/
WPIMULT
W_1 0.5 /
/
TSTEP
10 10 /
10 /
WPIMULT
W_1 0.5 /
/
TSTEP
10 /

View File

@ -1052,6 +1052,8 @@ BOOST_AUTO_TEST_CASE(BLOCK_VARIABLES) {
writer.add_timestep( 0, 0 * day, cfg.es, cfg.schedule, cfg.wells , {},{}, block_values);
writer.add_timestep( 1, 1 * day, cfg.es, cfg.schedule, cfg.wells , {},{}, block_values);
writer.add_timestep( 2, 2 * day, cfg.es, cfg.schedule, cfg.wells , {},{}, block_values);
writer.add_timestep( 3, 2 * day, cfg.es, cfg.schedule, cfg.wells , {},{}, block_values);
writer.add_timestep( 4, 2 * day, cfg.es, cfg.schedule, cfg.wells , {},{}, block_values);
writer.write();
auto res = readsum( cfg.name );
@ -1074,6 +1076,9 @@ BOOST_AUTO_TEST_CASE(BLOCK_VARIABLES) {
BOOST_CHECK_CLOSE( 2.6788413524627686 , ecl_sum_get_well_completion_var( resp, 1, "W_2", "CTFAC", 102) , 1e-5);
BOOST_CHECK_CLOSE( 2.7855057716369629 , ecl_sum_get_well_completion_var( resp, 1, "W_3", "CTFAC", 3) , 1e-5);
BOOST_CHECK_CLOSE( 50 , ecl_sum_get_well_completion_var( resp, 3, "W_1", "CTFAC", 1) , 1e-5);
BOOST_CHECK_CLOSE( 25 , ecl_sum_get_well_completion_var( resp, 4, "W_1", "CTFAC", 1) , 1e-5);
// Cell is not active
BOOST_CHECK( !ecl_sum_has_general_var( resp , "BPR:2,1,10"));
}