diff --git a/ebos/eclbaseaquifermodel.hh b/ebos/eclbaseaquifermodel.hh
index 66234324c..7b92faba5 100644
--- a/ebos/eclbaseaquifermodel.hh
+++ b/ebos/eclbaseaquifermodel.hh
@@ -57,6 +57,13 @@ public:
         : simulator_(simulator)
     {}
 
+    /*!
+     * \brief Called once the problem has been fully initialized and the initial
+     *        condition has been applied.
+     */
+    void initialSolutionApplied()
+    { }
+
     /*!
      * \brief This method is called when a new episode (report step) starts.
      */
diff --git a/ebos/eclproblem.hh b/ebos/eclproblem.hh
index 555fe7f8e..5124f566d 100644
--- a/ebos/eclproblem.hh
+++ b/ebos/eclproblem.hh
@@ -1298,7 +1298,6 @@ public:
      */
     void initialSolutionApplied()
     {
-
         if (!GET_PROP_VALUE(TypeTag, DisableWells)) {
             // initialize the wells. Note that this needs to be done after initializing the
             // intrinsic permeabilities and the after applying the initial solution because
@@ -1315,6 +1314,8 @@ public:
         this->simulator().vanguard().releaseEquilGrid();
 
         updateCompositionChangeLimits_();
+
+        aquiferModel_.initialSolutionApplied();
     }
 
     /*!