diff --git a/ebos/equil/initstateequil.hh b/ebos/equil/initstateequil.hh
index c7cf0d32d..3c597a7e0 100644
--- a/ebos/equil/initstateequil.hh
+++ b/ebos/equil/initstateequil.hh
@@ -1940,7 +1940,7 @@ private:
             saturations = psat.deriveSaturations(pos, eqreg, ptable);
             if (aquifer.hasCell(global_index)) {
                 saturations = {0.0, 0.0, 1.0};
-                const auto& aqu_cell = aquifer.getCell(global_index);
+                const auto &aqu_cell = aquifer.getCell(global_index);
                 std::ostringstream ss;
                 ss << "FOR AQUIFER CELL AT { " << aqu_cell.I + 1 << " " << aqu_cell.J + 1 << " "
                    << aqu_cell.J + 1 << " } OF NUMERICAL AQUIFER " << aqu_cell.aquifer_id << " , "
@@ -1948,6 +1948,14 @@ private:
                 OpmLog::info(ss.str());
             }
             pressures   = psat.correctedPhasePressures();
+            if (aquifer.hasCell(global_index)) {
+                const auto &aqu_cell = aquifer.getCell(global_index);
+                // TODO: NOT totally sure what we should do here to empoly the pressure specified by AQUNUM
+                if (aqu_cell.init_pressure > 0.) {
+                    const double pres = aqu_cell.init_pressure;
+                    pressures = {pres, pres, pres};
+                }
+            }
 
             const auto temp = this->temperature_[cell];