OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #1064 from GitPaean/fixing_thp_initialization

Browse Source 
when no THP control/limit, thp value will be set to zero.
This commit is contained in:
Atgeirr Flø Rasmussen 2017-03-10 10:23:15 +01:00 committed by GitHub
commit 81ebba2df8
1 changed files with 16 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
opm/autodiff/WellStateFullyImplicitBlackoil.hpp
View File

@ -171,6 +171,22 @@ namespace Opm
currentControls()[ newIndex ] = 0; currentControls()[ newIndex ] = 0;
} }
} }
// If in the new step, there is no THP related target/limit anymore, its thp value should be
// set to zero.
const WellControls* ctrl = wells->ctrls[w];
const int nwc = well_controls_get_num(ctrl);
int ctrl_index = 0;
for (; ctrl_index < nwc; ++ctrl_index) {
if (well_controls_iget_type(ctrl, ctrl_index) == THP) {
break;
}
}
// not finding any thp related control/limits
if (ctrl_index == nwc) {
thp()[w] = 0.;
}
} }
} }
} }