diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidFields.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidFields.H
index 337b88acb2..1be5924d42 100644
--- a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidFields.H
+++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidFields.H
@@ -214,9 +214,6 @@ forAll(fluidRegions, i)
         )
     );
 
-    // Force p_rgh to be consistent with p
-    p_rghFluid[i] = thermoFluid[i].p() - rhoFluid[i]*ghFluid[i] - pRefFluid[i];
-
     fluidRegions[i].setFluxRequired(p_rghFluid[i].name());
 
     initialMassFluid[i] = fvc::domainIntegrate(rhoFluid[i]).value();
diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidPressureControls.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidPressureControls.H
index e2e7e46fde..5c58bf8aef 100644
--- a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidPressureControls.H
+++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidPressureControls.H
@@ -12,4 +12,17 @@ forAll(fluidRegions, i)
             false
         )
     );
+
+    hydrostaticInitialisation
+    (
+        p_rghFluid[i],
+        thermoFluid[i].p(),
+        rhoFluid[i],
+        UFluid[i],
+        ghFluid[i],
+        ghfFluid[i],
+        pRefFluid[i],
+        thermoFluid[i],
+        pimples.pimple(i).dict()
+    );
 }