OilfieldToolsNet/opm-common
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Don't Assume Increasing WLIST Position In Restart Files

Browse Source 
The well list positions for individual wells (i.e,, IWLS) need not,
and typically are not, strictly increasing in the restart files.
Therefore we must only 'resize()' the well name vectors if the
'well_order' is strictly larger than the current size.  Otherwise,
we lose information about which wells are in which well lists.
This commit is contained in:
Bård Skaflestad 2022-06-15 11:29:29 +02:00
parent 5495266115
commit 042830e860
1 changed files with 10 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
src/opm/io/eclipse/rst/state.cpp
View File

@ -380,13 +380,17 @@ void RstState::add_wlist(const std::vector<std::string>& zwls,
const auto& well_name = this->wells[well_index].name;
for (auto wlist_index = 0*this->header.max_wlist; wlist_index < this->header.max_wlist; wlist_index++) {
int well_order = iwls[iwls_offset + wlist_index];
if (well_order != 0) {
const auto& wlist_name = zwls[zwls_offset + wlist_index];
auto& wlist = this->wlists[wlist_name];
wlist.resize( well_order );
wlist[well_order - 1] = well_name;
const auto well_order = iwls[iwls_offset + wlist_index];
if (well_order < 1) {
continue;
}
auto& wlist = this->wlists[zwls[zwls_offset + wlist_index]];
if (wlist.size() < static_cast<std::vector<std::string>::size_type>(well_order)) {
wlist.resize(well_order);
}
wlist[well_order - 1] = well_name;
}
}
}