Updating rho with thermo rho at the end of the pEq.H, independently of

simpleRho or limiting p.

applications/solvers/compressible/rhoPimpleFoam/pEqn.H

@@ -98,19 +98,11 @@ K = 0.5*magSqr(U);
 if (pressureControl.limit(p))
 {
     p.correctBoundaryConditions();
-    thermo.correctRho(psi*p - psip0, rhoMin, rhoMax);
-    rho = thermo.rho();
-}
-else if (pimple.SIMPLErho())
-{
-    thermo.correctRho(psi*p - psip0, rhoMin, rhoMax);
-    rho = thermo.rho();
-}
-else
-{
-    thermo.correctRho(psi*p - psip0, rhoMin, rhoMax) ;
 }
 
+thermo.correctRho(psi*p - psip0, rhoMin, rhoMax) ;
+rho = thermo.rho();
+
 if (thermo.dpdt())
 {
     dpdt = fvc::ddt(p);