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

prettification

Browse Source
This commit is contained in:
Kjetil Olsen Lye 2012-04-13 13:21:33 +02:00
parent 5b21e74999
commit 47dc72812b
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -420,8 +420,8 @@ namespace Opm
double cell_depth = grid.cell_centroids[3*cell+2]; double cell_depth = grid.cell_centroids[3*cell+2];
double saturation_sum = 0.0; double saturation_sum = 0.0;
for(size_t i = 0; i < densities.size(); i++) { for(size_t p = 0; p < np; p++) {
saturation_sum += saturations[densities.size()*cell + i]; saturation_sum += saturations[np*cell + p];
} }
if(saturation_sum == 0) { if(saturation_sum == 0) {
saturation_sum = 1.0; saturation_sum = 1.0;
@ -429,7 +429,7 @@ namespace Opm
double density = 0.0; double density = 0.0;
for(size_t p = 0; p < np; p++) { 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[np*cell + p] * densities[p]; density += saturations[np*cell + p] * densities[p] / saturation_sum;
} }
// Is the sign correct? // Is the sign correct?