diff --git a/src/thermophysicalModels/radiation/radiationModel/radiationModel/radiationModel.C b/src/thermophysicalModels/radiation/radiationModel/radiationModel/radiationModel.C
index cfbb2b614552e0ae7055dc2738f39705982f9bfb..64a2cfc3efd25568864d33bf6ab6186886d76128 100644
--- a/src/thermophysicalModels/radiation/radiationModel/radiationModel/radiationModel.C
+++ b/src/thermophysicalModels/radiation/radiationModel/radiationModel/radiationModel.C
@@ -143,14 +143,14 @@ Foam::tmp<Foam::fvScalarMatrix> Foam::radiation::radiationModel::Sh
 ) const
 {
     volScalarField& h = thermo.h();
-    const volScalarField cp = thermo.Cp();
+    const volScalarField Cp = thermo.Cp();
     const volScalarField T3 = pow3(T_);
 
     return
     (
         Ru()
-      - fvm::Sp(4.0*Rp()*T3/cp, h)
-      - Rp()*T3*(T_ - 4.0*h/cp)
+      - fvm::Sp(4.0*Rp()*T3/Cp, h)
+      - Rp()*T3*(T_ - 4.0*h/Cp)
     );
 }
 
@@ -161,14 +161,14 @@ Foam::tmp<Foam::fvScalarMatrix> Foam::radiation::radiationModel::Shs
 ) const
 {
     volScalarField& hs = thermo.hs();
-    const volScalarField cp = thermo.Cp();
+    const volScalarField Cp = thermo.Cp();
     const volScalarField T3 = pow3(T_);
 
     return
     (
         Ru()
-      - fvm::Sp(4.0*Rp()*T3/cp, hs)
-      - Rp()*T3*(T_ - 4.0*hs/cp)
+      - fvm::Sp(4.0*Rp()*T3/Cp, hs)
+      - Rp()*T3*(T_ - 4.0*hs/Cp)
     );
 }