Use a 'P' sign that is somehow related to the equation we're solving

The former definition was inherited from the hybrid solver in MRST and
not really applicable to compressible TPFA.
This commit is contained in:
Bård Skaflestad 2010-11-05 16:58:19 +01:00
parent 8780df0a55
commit 223681b8c4

View File

@ -58,7 +58,7 @@ compr_accum_term(size_t nc,
size_t c;
for (c = 0; c < nc; c++) {
P[c] = -totcompr[c] * porevol[c] / dt;
P[c] = totcompr[c] * porevol[c] / dt;
}
}
@ -76,6 +76,6 @@ compr_src_add_press_accum(size_t nc,
size_t c;
for (c = 0; c < nc; c++) {
src[c] -= P[c] * p0[c];
src[c] = P[c] * p0[c];
}
}