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

Opm::computeWDP(): Don't shadow (well) loop counter in phase loop.

Browse Source
This commit is contained in:
Bård Skaflestad 2012-04-12 22:43:33 +02:00
parent 8ab03f2288
commit e66a2572c6
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
opm/core/utility/miscUtilities.cpp
View File

@ -408,6 +408,8 @@ namespace Opm
void computeWDP(const Wells& wells, const UnstructuredGrid& grid, const std::vector<double>& saturations, void computeWDP(const Wells& wells, const UnstructuredGrid& grid, const std::vector<double>& saturations,
const std::vector<double>& densities, std::vector<double>& wdp) const std::vector<double>& densities, std::vector<double>& wdp)
{ {
const size_t np = densities.size();
// Simple for now: // Simple for now:
for(int i = 0; i < wells.number_of_wells; i++) { for(int i = 0; i < wells.number_of_wells; i++) {
double depth_ref = wells.depth_ref[i]; double depth_ref = wells.depth_ref[i];
@ -418,9 +420,9 @@ namespace Opm
double cell_depth = grid.cell_centroids[3*cell+2]; double cell_depth = grid.cell_centroids[3*cell+2];
double density = 0.0; double density = 0.0;
for(size_t i = 0; i < densities.size(); i++) { for(size_t p = 0; p < np; p++) {
// Is this a smart way of doing it? // Is this a smart way of doing it?
density += saturations[densities.size()*cell+i]*densities[i]; density += saturations[np*cell + p] * densities[p];
} }
// Is the sign correct? // Is the sign correct?