diff --git a/solvers_post/SLFMFoam/Mixturefraction.H b/solvers_post/SLFMFoam/Mixturefraction.H
index 79d73d1..f4c7a20 100644
--- a/solvers_post/SLFMFoam/Mixturefraction.H
+++ b/solvers_post/SLFMFoam/Mixturefraction.H
@@ -3,7 +3,7 @@ fvScalarMatrix mfEqn
 (
 
       fvm::div(phi, mf)
-    - fvm::laplacian(1.47*turbulence->mut(), mf) //let 1/Sc = 1.47
+    - fvm::laplacian((1/Sc)*turbulence->mut(), mf) //let 1/Sc = 1.47
 );
 
 
diff --git a/solvers_post/SLFMFoam/MixturefractionVar.H b/solvers_post/SLFMFoam/MixturefractionVar.H
index 2db84b5..b05b582 100644
--- a/solvers_post/SLFMFoam/MixturefractionVar.H
+++ b/solvers_post/SLFMFoam/MixturefractionVar.H
@@ -7,9 +7,9 @@ fvScalarMatrix mfVarEqn
 (
 
      fvm::div(phi, mfVar)
-    - fvm::laplacian(1.47*turbulence->mut(), mfVar) //let 1/Sc = 1.47
+    - fvm::laplacian((1/Sc)*turbulence->mut(), mfVar) 
     ==
-    2 * ( 1.47 * turbulence->mut() ) * (Gradmf & Gradmf)
+    2 * ( (1/Sc) * turbulence->mut() ) * (Gradmf & Gradmf)
     - 2 * rho * SDR
 );
 
diff --git a/solvers_post/SLFMFoam/SLFMFoam.C b/solvers_post/SLFMFoam/SLFMFoam.C
index 70e2244..d434ede 100644
--- a/solvers_post/SLFMFoam/SLFMFoam.C
+++ b/solvers_post/SLFMFoam/SLFMFoam.C
@@ -25,10 +25,14 @@ Application
     SLFMFoam
 
 Description
-    Steady state solver for turbulent combustion.
+    Steady Laminar Flamelet Model(SLFM) solver for turbulent combustion. SLFM 
+    assumes a turbulent flame composed of thin stretched laminar flamelets.
+    Each flamelet library is generated by external program and imported by 
+    the solver.
     
 References
-    N. Peters, Turbulent Combustion, Cambridge University Press; 2000.
+    A.Y. Klimenko, R.W. Bilger, Progress in Energy and Combustion Science 25 (1999) 595-687
+    N. Peters, Turbulent Combustion, Cambridge University Press (2000)
 
 Contact
     POSTECH combustion lab.
@@ -55,7 +59,7 @@ int main(int argc, char *argv[])
     #include "createControl.H"
     #include "createFields.H"
     #include "createFieldRefs.H"
-    #include "BetaPDF.H" //make detailed integration space and evaluate beta-pdf integration
+    #include "BetaPDF.H" //make detailed integration space and perform beta-pdf integration
 
     #include "createFvOptions.H"
     #include "initContinuityErrs.H"
@@ -81,8 +85,8 @@ int main(int argc, char *argv[])
         // --- Pressure-velocity SIMPLE corrector loop
         {
             #include "UEqn.H"
-            #include "Mixturefraction.H" //mixture fraction
-            #include "MixturefractionVar.H" //mixture fraction variance 
+            #include "Mixturefraction.H" //mean mixture fraction
+            #include "MixturefractionVar.H" //mean mixture fraction variance 
             #include "setSDR.H" //calculate scalar dissipation rate
             #include "updateT_RHO.H" //update temperature and density
             #include "pEqn.H"
diff --git a/solvers_post/SLFMFoam/createFields.H b/solvers_post/SLFMFoam/createFields.H
index e5b53c9..0b611ba 100644
--- a/solvers_post/SLFMFoam/createFields.H
+++ b/solvers_post/SLFMFoam/createFields.H
@@ -308,7 +308,9 @@ scalar stoiMF(readScalar(SLFMdict.lookup("stoiMF")));
 scalar stretchFac(readScalar(SLFMdict.lookup("stretchFactor")));
 scalar lowerN(readScalar(SLFMdict.lookup("lowerN")));
 scalar upperN(readScalar(SLFMdict.lookup("upperN")));
+scalar Sc(readScalar(SLFMdict.lookup("Sc")));
 scalar etamax(lowerN+upperN);
+word outletName(SLFMdict.lookup("outletName"));
 
 scalarField etaValue(etamax+1,0);
 scalarField etaIndex(etamax+1,0);
diff --git a/solvers_post/SLFMFoam/readSLFMlib.H b/solvers_post/SLFMFoam/readSLFMlib.H
index dc26a22..f58e724 100644
--- a/solvers_post/SLFMFoam/readSLFMlib.H
+++ b/solvers_post/SLFMFoam/readSLFMlib.H
@@ -30,12 +30,12 @@ scalarFieldArray2d W_SLFM //reaction rate [1/sec]
         scalarField(etamax+1,0)
     )
 );
-scalarFieldArray1d rho_SLFM //density [g/cm3 ?]
+scalarFieldArray1d rho_SLFM //density [g/cm3]
 (
     NstList.size(),
     scalarField(etamax+1, 0)
 );
-scalarFieldArray1d h_SLFM //enthalpy [?]
+scalarFieldArray1d h_SLFM //enthalpy [erg/gm]
 (
     NstList.size(),
     scalarField(etamax+1, 0)
diff --git a/solvers_post/SLFMFoam/updateT_RHO.H b/solvers_post/SLFMFoam/updateT_RHO.H
index cf6e310..3099a32 100644
--- a/solvers_post/SLFMFoam/updateT_RHO.H
+++ b/solvers_post/SLFMFoam/updateT_RHO.H
@@ -63,7 +63,7 @@ const fvPatchList& patches = mesh.boundary();
 forAll(patches, patchI)
 {
     const fvPatch& curPatch = patches[patchI];
-    if(curPatch.name() == "outlet")
+    if(curPatch.name() == outletName)
     {
         forAll(curPatch, faceI)
         {