diff --git a/opm/core/pressure/IncompTpfa.cpp b/opm/core/pressure/IncompTpfa.cpp index 23e82d92..9417b7eb 100644 --- a/opm/core/pressure/IncompTpfa.cpp +++ b/opm/core/pressure/IncompTpfa.cpp @@ -121,7 +121,7 @@ namespace Opm pressure.resize(grid_.number_of_cells); faceflux.resize(grid_.number_of_faces); - ifs_tpfa_solution soln = { NULL, NULL }; + ifs_tpfa_solution soln = { NULL, NULL, NULL, NULL }; soln.cell_press = &pressure[0]; soln.face_flux = &faceflux[0]; @@ -204,7 +204,7 @@ namespace Opm pressure.resize(grid_.number_of_cells); faceflux.resize(grid_.number_of_faces); - ifs_tpfa_solution soln = { NULL, NULL }; + ifs_tpfa_solution soln = { NULL, NULL, NULL, NULL }; soln.cell_press = &pressure[0]; soln.face_flux = &faceflux[0]; diff --git a/opm/core/pressure/tpfa/ifs_tpfa.h b/opm/core/pressure/tpfa/ifs_tpfa.h index fbe5a4f6..f2ce3851 100644 --- a/opm/core/pressure/tpfa/ifs_tpfa.h +++ b/opm/core/pressure/tpfa/ifs_tpfa.h @@ -42,6 +42,9 @@ struct ifs_tpfa_data { struct ifs_tpfa_solution { double *cell_press; double *face_flux ; + + double *well_press; /* BHP */ + double *well_flux ; /* Perforation (total) fluxes */ }; struct ifs_tpfa_forces {