Name fields to avoid duplicate registration

applications/solvers/compressible/rhoCentralFoam/centralCourantNo.H

@@ -22,16 +22,13 @@ License
     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
 
 Global
-    compressibleCourantNo
+    centralCourantNo
 
 Description
     Calculates the mean and maximum wave speed based Courant Numbers.
 
 \*---------------------------------------------------------------------------*/
 
-scalar CoNum = 0.0;
-scalar meanCoNum = 0.0;
-
 if (mesh.nInternalFaces())
 {
     surfaceScalarField amaxSfbyDelta

applications/solvers/compressible/rhoCentralFoam/createFields.H

@@ -8,7 +8,6 @@ psiThermo& thermo = pThermo();
 
 volScalarField& p = thermo.p();
 volScalarField& e = thermo.he();
-const volScalarField& T = thermo.T();
 const volScalarField& psi = thermo.psi();
 const volScalarField& mu = thermo.mu();
 

applications/solvers/compressible/rhoCentralFoam/rhoCentralFoam.C

@@ -53,6 +53,10 @@ int main(int argc, char *argv[])
 
     dimensionedScalar v_zero("v_zero", dimVolume/dimTime, 0.0);
 
+    // Courant numbers used to adjust the time-step
+    scalar CoNum = 0.0;
+    scalar meanCoNum = 0.0;
+
     Info<< "\nStarting time loop\n" << endl;
 
     while (runTime.run())
@@ -160,7 +164,7 @@ int main(int argc, char *argv[])
         // estimated by the central scheme
         amaxSf = max(mag(aphiv_pos), mag(aphiv_neg));
 
-        #include "compressibleCourantNo.H"
+        #include "centralCourantNo.H"
         #include "readTimeControls.H"
         #include "setDeltaT.H"
 
@@ -199,8 +203,6 @@ int main(int argc, char *argv[])
         U.correctBoundaryConditions();
         rhoU.boundaryField() = rho.boundaryField()*U.boundaryField();
 
-        volScalarField rhoBydt(rho/runTime.deltaT());
-
         if (!inviscid)
         {
             solve