Fixed bug in gauss seidel gravity solver.

This commit is contained in:
Xavier Raynaud
2012-05-22 09:51:32 +02:00
parent 48adf98b5a
commit 8e0936fadd

View File

@@ -998,7 +998,7 @@ namespace Opm
{
nbcell[0] = -1;
gf[0] = 0.0;
if (cell > 0) {
if (pos > 0) {
nbcell[0] = cells[pos - 1];
gf[0] = -gravflux[pos - 1];
}