diff --git a/solvers_post/LagrangianCMCFoam/Flamegroup.H b/solvers_post/LagrangianCMCFoam/Flamegroup.H
index a4da813..07ca3d4 100644
--- a/solvers_post/LagrangianCMCFoam/Flamegroup.H
+++ b/solvers_post/LagrangianCMCFoam/Flamegroup.H
@@ -19,3 +19,20 @@
         F_total += F[k];
     }
     
+    forAll(rho , celli)
+    {
+        if(F_total[celli] > 0)
+        {
+            for(label k=0 ; k<=group ; k++)
+            {
+                Ffrac[k][celli] = F[k][celli]/F_total[celli];
+            }
+        }
+        else if(F_total[celli] <= 0)
+        {
+            for(label k=0 ; k<=group ; k++)
+            {
+                Ffrac[k][celli] = 0.0;
+            }
+        }
+    }
diff --git a/solvers_post/LagrangianCMCFoam/setSDRCMC.H b/solvers_post/LagrangianCMCFoam/setSDRCMC.H
index 7e11afa..2aacc41 100644
--- a/solvers_post/LagrangianCMCFoam/setSDRCMC.H
+++ b/solvers_post/LagrangianCMCFoam/setSDRCMC.H
@@ -31,20 +31,6 @@
                                 CSDRCELL[j][celli] = C2 * AMC(eta);
                         }
                 }
-		if(F_total[celli] > 0)
-               	{
-			for(label k=0 ; k<=group ; k++)
-			{
-                		Ffrac[k][celli] = F[k][celli]/F_total[celli];
-			}
-              	}
-		else if(F_total[celli] <= 0)
-		{
-                        for(label k=0 ; k<=group ; k++)
-                        {
-                                Ffrac[k][celli] = 0.0;
-                        }			
-		}
         }
 
 //get density weighted averaged value of conditional SDR for each eta and flame-group