diff --git a/applications/solvers/dfHighSpeedFoam/createFieldsSave.H b/applications/solvers/dfHighSpeedFoam/createFieldsSave.H
index c3fd5cd4..37723319 100644
--- a/applications/solvers/dfHighSpeedFoam/createFieldsSave.H
+++ b/applications/solvers/dfHighSpeedFoam/createFieldsSave.H
@@ -41,6 +41,6 @@ else if(ddtSchemes == "RK3SSP")
 {
     rkcoe1[0]=1.0; rkcoe2[0]=0.0; rkcoe3[0]=1.0;
     rkcoe1[1]=0.75; rkcoe2[1]=0.25; rkcoe3[1]=0.25;
-    rkcoe1[2]=0.33; rkcoe2[2]=0.66; rkcoe3[2]=0.66;
+    rkcoe1[2]=1.0/3.0; rkcoe2[2]=2.0/3.0; rkcoe3[2]=2.0/3.0;
     rk=3;
 }
diff --git a/applications/solvers/dfHighSpeedFoam/rhoYEqn.H b/applications/solvers/dfHighSpeedFoam/rhoYEqn.H
index fac11306..ba634503 100644
--- a/applications/solvers/dfHighSpeedFoam/rhoYEqn.H
+++ b/applications/solvers/dfHighSpeedFoam/rhoYEqn.H
@@ -110,7 +110,7 @@ time_monitor_corrDiff += double(end - start) / double(CLOCKS_PER_SEC);
                     rhoYi_rhs = 
                         (
                             turbName == "laminar"
-                            ?  (fvc::laplacian(DEff(), Yi) - fvc::div(phiUc,Yi,"div(phic,Yi)")) 
+                            ?  (fvc::laplacian(DEff(), Yi) - fvc::div(phiUc,Yi,"div(phi,Yi_h)")) 
                             :  fvc::laplacian(DEff(), Yi)
                         );