OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-01-03 21:16:54 -06:00
Code Issues Packages Projects Releases Wiki Activity

Guard against 'bc' being NULL.

Browse Source 
This is in preparation of introducing compressible boundary conditions
as represented by the 'compr_bc' module.
This commit is contained in:
Bård Skaflestad 2011-10-25 19:34:32 +02:00
parent d1026c82d5
commit 51111176db
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/cfs_tpfa_residual.c
View File

@ -612,7 +612,8 @@ compute_fpress(grid_t *G,
c1 = G->face_cells[2*f + 0];
c2 = G->face_cells[2*f + 1];
if (((c1 < 0) || (c2 < 0)) && (bc->type[f] == PRESSURE)) {
if (((c1 < 0) || (c2 < 0)) &&
(bc != NULL) && (bc->type[f] == PRESSURE)) {
fpress[f] = bc->bcval[f];
}
}
@ -638,7 +639,8 @@ compute_flux(grid_t *G,
c1 = G->face_cells[2*f + 0];
c2 = G->face_cells[2*f + 1];
if (((c1 < 0) || (c2 < 0)) && (bc->type[f] == FLUX)) {
if (((c1 < 0) || (c2 < 0)) &&
(bc != NULL) && (bc->type[f] == FLUX)) {
fflux[f] = bc->bcval[f];
continue;
}
@ -652,7 +654,7 @@ compute_flux(grid_t *G,
if ((c1 >= 0) && (c2 >= 0)) {
dp = cpress[c1] - cpress[c2];
} else if (bc->type[f] == PRESSURE) {
} else if ((bc != NULL) && (bc->type[f] == PRESSURE)) {
if (c1 < 0) {
dp = bc->bcval[f] - cpress[c2];
/* g = -g; */