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

set thp() value to zero when no THP constraints exist

Browse Source 
in updateWellState()
This commit is contained in:
Kai Bao
2017-03-09 15:51:03 +01:00
parent 7a7dcdf3ba
commit cb85319a72
1 changed files with 7 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
opm/autodiff/StandardWellsDense_impl.hpp
View File

@@ -1426,7 +1426,8 @@ namespace Opm {
const WellControls* wc = wells().ctrls[w]; const WellControls* wc = wells().ctrls[w];
const int nwc = well_controls_get_num(wc); const int nwc = well_controls_get_num(wc);
// Looping over all controls until we find a THP constraint // Looping over all controls until we find a THP constraint
for (int ctrl_index = 0; ctrl_index < nwc; ++ctrl_index) { int ctrl_index = 0;
for ( ; ctrl_index < nwc; ++ctrl_index) {
if (well_controls_iget_type(wc, ctrl_index) == THP) { if (well_controls_iget_type(wc, ctrl_index) == THP) {
// the current control // the current control
const int current = well_state.currentControls()[w]; const int current = well_state.currentControls()[w];
@@ -1479,6 +1480,11 @@ namespace Opm {
break; break;
} }
} // end of for loop for seaching THP constraints } // end of for loop for seaching THP constraints
// no THP constraint found
if (ctrl_index == nwc) { // not finding a THP contstraints
well_state.thp()[w] = 0.0;
}
} // end of for (int w = 0; w < nw; ++w) } // end of for (int w = 0; w < nw; ++w)
} }