diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/Allrun.sh b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/Allrun.sh
index 20153244..73759c80 100755
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/Allrun.sh
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/Allrun.sh
@@ -27,5 +27,5 @@ else
 fi
 
 #- run parallel CFD-DEM in new terminal
-gnome-terminal -e "bash $casePath/parCFDDEMrun.sh" 
-#bash $casePath/parCFDDEMrun.sh
+#gnome-terminal -e "bash $casePath/parCFDDEMrun.sh" 
+bash $casePath/parCFDDEMrun.sh
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/combustionProperties b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/combustionProperties
index c81ec351..1f0c2fbf 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/combustionProperties
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/combustionProperties
@@ -15,21 +15,17 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-//combustionModel  noCombustion<psiThermoCombustion>;//laminar<psiChemistryCombustion>;////noCombustion<psiThermoCombustion>;//PaSR<psiChemistryCombustion>;
 combustionModel  PaSR<rhoChemistryCombustion>;
 
 active  false; //true;
 
-
 noCombustionCoeffs
-{
-    
-}
+{}
+
 PaSRCoeffs
 {
     Cmix                1.0;
     turbulentReaction   off;
 }
 
-
 // ************************************************************************* //
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/couplingProperties b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/couplingProperties
index 230f7919..f1da1954 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/couplingProperties
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/couplingProperties
@@ -102,16 +102,7 @@ reactionHeatProps
 speciesProps
 {
 	ChemistryFile 		"$casePath/CFD/constant/foam.inp";
-	tempFieldName 		"T";
-	partTempName  		"partTemp";
-	densityFieldName 	"rho";
-	partRhoName		"partRho";
-        voidfractionFieldName   "voidfraction";
-        totalMoleFieldName      "molarConc";
-        partMoleName            "partMolarConc";
-        pressureFieldName       "p";
-        partPName               "partP";
-	Nevery			10;
+	Nevery			1;
 }
 
 SyamlalThermCondProps
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/foam.inp b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/foam.inp
index e14dbf33..c504c918 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/foam.inp
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/foam.inp
@@ -7,6 +7,5 @@ species
 
 reactions
 {
-
 }
 
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/controlDict b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/controlDict
index 72bf9029..9fffc4b1 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/controlDict
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/controlDict
@@ -29,7 +29,7 @@ deltaT          0.0001;
 
 writeControl    timeStep;
 
-writeInterval   10000;
+writeInterval   1000;
 
 purgeWrite      0;
 
@@ -54,7 +54,6 @@ maxDeltaT       0.1;
 
 // ************************************************************************* //
  libs (  
-   //     "libsimpleFunctionObjects.so"
         "libfieldFunctionObjects.so"
      );
 
@@ -63,9 +62,7 @@ functions
     probes1
     {
 	type probes;
-
         functionObjectLibs ("libsampling.so");
-
 	#include "probesDict";
     }
   
@@ -110,189 +107,6 @@ functions
     {
         type	rhoReactionThermoMoleFractions;
     }
-
-/*    Average-O2
-    {
-        type            fieldAverage;
-        functionObjectLibs ( "libfieldFunctionObjects.so" );
-        enabled         true;
-        outputControl   outputTime; //timeStep; //
-        //outputInterval  0.50;
-        log             true;
-        restartOnRestart    false;
-        restartOnOutput     false;
-        
-        
-        fields
-        (
-            O2
-            {
-                mean        on;
-                prime2Mean  on;
-                base        time;
-            }
-        );
-    }
-    
-   species_O2 // Name also used to identify output folder 
-    { 
-        type            cellSource; 
-        functionObjectLibs ("libfieldValueFunctionObjects.so"); 
-        enabled         true; 
-        outputControl   timeStep;//outputTime; 
-        outputInterval  1,
-        log             true;    
-        valueOutput     false;   //true; 
-        source          all;   // Type of cell source or all for the full domain
-        sourceName      O2; 
-        operation       weightedAverage; // volAverage; //volIntegrate; //weightedAverage; //none //sum;
-        fields 
-        ( 
-            O2
-        ); 
-    } 
-    Average-CO2
-    {
-        type            fieldAverage;
-        functionObjectLibs ( "libfieldFunctionObjects.so" );
-        enabled         true;
-        outputControl   outputTime; //timeStep; //
-        //outputInterval  0.50;
-        log             true;
-        restartOnRestart    false;
-        restartOnOutput     false;
-        
-        
-        fields
-        (
-           CO2
-            {
-                mean        on;
-                prime2Mean  on;
-                base        time;
-            }
-        );
-    }
-    
-   species_CO2 // Name also used to identify output folder 
-    { 
-        type            cellSource; 
-        functionObjectLibs ("libfieldValueFunctionObjects.so"); 
-        enabled         true; 
-        outputControl   timeStep;//outputTime; 
-        outputInterval  1,
-        log             true;    
-        valueOutput     false;   //true; 
-        source          all;   // Type of cell source or all for the full domain
-        sourceName      CO2; 
-        operation       weightedAverage; // volAverage; //volIntegrate; //weightedAverage; //none //sum;
-        fields 
-        ( 
-            CO2
-        ); 
-    }
-*/ 
 } 
 
 // ************************************************************************* //
- 
- /* // IF swak4Foam is installed following can be used;
-   volumetricFlowSurfaceField
-   {
-        type patchExpression;
-        outputControlMode timeStep;
-        outputInterval    1;
-        verbose true;
-        accumulations ( sum  );
-        patches ( outlet );
-        expression "CO2";
-    }
-    
-    massFlowSurfaceField
-    {
-        $volumetricFlowSurfaceField;
-        patches ( outlet );
-        variables ("Area=0.25;"
-                    
-        );
-        expression "CO2*rho*Area*U";
-    }
-
-    CO2_Mass_in_Volume
-    {
-        type                swakExpression;
-        valueType           internalField; //patch;
-        //aliases{alpha alpha.water;}
-        //patchName           outlet;
-        verbose             true;
-        expression          "rho*voidfraction*CO2*vol()";
-        accumulations       ( sum );
-        outputControlMode   timeStep; //outputTime;
-        outputInterval      1;
-    }
-    
-    OutletMassFlow-CO2
-    {
-        type patchExpression;
-        accumulations (
-            sum
-        );
-        patches (
-            outlet
-        );
-        expression "CO2";
-        verbose true;
-        autoInterpolate true;
-    }
-    
-    CO2_Concentration_in_Volume
-    {
-        type                swakExpression;
-        valueType           internalField; //patch;
-        //aliases{alpha alpha.water;}
-        //patchName           outlet;
-        verbose             true;
-        expression          "CO2";
-        accumulations       ( sum );
-        outputControlMode   timeStep; //outputTime;
-        outputInterval      1;
-    }
-    
-    O2_Mass_in_Volume
-    {
-        type                swakExpression;
-        valueType           internalField; //patch;
-        //aliases{alpha alpha.water;}
-        //patchName           outlet;
-        verbose             true;
-        expression          "rho*voidfraction*O2*vol()";
-        accumulations       ( sum );
-        outputControlMode   timeStep; //outputTime;
-        outputInterval      1;
-    }
-   
-    O2_Concentration_in_Volume
-    {
-        type                swakExpression;
-        valueType           internalField; //patch;
-        //aliases{alpha alpha.water;}
-        //patchName           outlet;
-        verbose             true;
-        expression          "O2";
-        accumulations       ( sum );
-        outputControlMode   timeStep; //outputTime;
-        outputInterval      1;
-    }
-    
-    OutletMassFlow-O2
-    {
-        type patchExpression;
-        accumulations (
-            sum
-        );
-        patches (
-            outlet
-        );
-        expression "O2";
-        verbose true;
-    } */
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/probesDict b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/probesDict
index d48c2713..12a28583 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/probesDict
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/probesDict
@@ -34,10 +34,9 @@ writeInterval  10;
 // Locations to be probed.
 probeLocations
 (
-	(7.5 2.5 2.5)
-        
-
-
+	(7.0 2.5 2.5)
+        (8.0 2.5 2.5)
+        (1.0 2.5 2.5)
 );
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //