diff --git a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/30CO6H2_Hematite/DEM/in.liggghts_init b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/30CO6H2_Hematite/DEM/in.liggghts_init
index 30fee690..d0d48418 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/30CO6H2_Hematite/DEM/in.liggghts_init
+++ b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/30CO6H2_Hematite/DEM/in.liggghts_init
@@ -24,7 +24,7 @@ fix         m3 all property/global coefficientRestitution peratomtypepair 1 0.3
 fix         m4 all property/global coefficientFriction peratomtypepair 1 0.5
 
 # pair style
-pair_style  gran model hertz tangential history # Hertzian without cohesion
+pair_style  gran model hertz tangential history
 pair_coeff  * *
 
 # timestep, gravity
@@ -48,11 +48,13 @@ compute         rke all erotate/sphere
 thermo_style    custom step atoms ke c_rke vol 
 thermo          1000
 thermo_modify   lost ignore norm no
-compute_modify  thermo_temp dynamic yes
+compute_modify  thermo_temp dynamic no extra 0
 
 # insert the first particles so that dump is not empty
 run     1
-dump    dmp all custom 100 post/dump.liggghts_init id type x y z vx vy vz fx fy fz radius mass 
+dump    dmp all custom 100 ../DEM/post/dump.liggghts_init id type x y z vx vy vz fx fy fz radius mass 
+
+run     10 upto
+
+write_restart ../DEM/post/restart/liggghts.restart
 
-run        10 upto 
-write_restart   post/restart/liggghts.restart
diff --git a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/30CO6H2_Hematite/DEM/in.liggghts_run b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/30CO6H2_Hematite/DEM/in.liggghts_run
index 7a6832d6..406515c8 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/30CO6H2_Hematite/DEM/in.liggghts_run
+++ b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/30CO6H2_Hematite/DEM/in.liggghts_run
@@ -26,7 +26,7 @@ fix         m3 all property/global coefficientRestitution peratomtypepair 1 0.3
 fix         m4 all property/global coefficientFriction peratomtypepair 1 0.5
 
 # pair style
-pair_style  gran model hertz tangential history # Hertzian without cohesion
+pair_style  gran model hertz tangential history
 pair_coeff  * *
 
 # timestep, gravity
@@ -92,6 +92,7 @@ variable time equal time
 variable m1   equal mass[1]
 variable rp   equal radius[1]
 variable rho1 equal mass[1]/((4/3)*PI*radius[1]*radius[1]*radius[1])
+
 fix printmass all print ${WI} "${time} ${m1} ${rho1} ${rp}" file mass_rho_rad.txt title "#time mass rho rad_par"
 
 # compute kinetic energy of particles to see if they are moving changing
@@ -207,8 +208,9 @@ compute         rke all erotate/sphere
 thermo_style    custom step atoms ke c_rke vol
 thermo          10
 thermo_modify   lost ignore norm no
-compute_modify  thermo_temp dynamic yes
+compute_modify  thermo_temp dynamic no extra 0
 
 dump    dmp all custom 100 ../DEM/post/dump*.liggghts_run id type x y z vx vy vz fx fy fz radius mass f_fracRed[1] f_fracRed[2] f_fracRed[3]
 
 run     1
+
diff --git a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/33CO3H2_Hematite/DEM/in.liggghts_init b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/33CO3H2_Hematite/DEM/in.liggghts_init
index 161b494d..6590c6e2 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/33CO3H2_Hematite/DEM/in.liggghts_init
+++ b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/33CO3H2_Hematite/DEM/in.liggghts_init
@@ -48,12 +48,13 @@ compute         rke all erotate/sphere
 thermo_style    custom step atoms ke c_rke vol
 thermo          1000
 thermo_modify   lost ignore norm no
-compute_modify  thermo_temp dynamic yes
+compute_modify  thermo_temp dynamic no extra 0
 
 # insert the first particles so that dump is not empty
 run     1
-dump    dmp all custom 100 post/dump.liggghts_init id type x y z vx vy vz fx fy fz radius mass 
+dump    dmp all custom 100 ../DEM/post/dump.liggghts_init id type x y z vx vy vz fx fy fz radius mass
 
 run     10 upto 
 
-write_restart   ../DEM/post/restart/liggghts.restart
+write_restart ../DEM/post/restart/liggghts.restart
+
diff --git a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/33CO3H2_Hematite/DEM/in.liggghts_run b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/33CO3H2_Hematite/DEM/in.liggghts_run
index 9b486956..2a0221ad 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/33CO3H2_Hematite/DEM/in.liggghts_run
+++ b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/33CO3H2_Hematite/DEM/in.liggghts_run
@@ -208,7 +208,9 @@ compute         rke all erotate/sphere
 thermo_style    custom step atoms ke c_rke vol
 thermo          10
 thermo_modify   lost ignore norm no
-compute_modify  thermo_temp dynamic yes
+compute_modify  thermo_temp dynamic no extra 0
 
 dump        dmp all custom 100 ../DEM/post/dump*.liggghts_run id type x y z vx vy vz fx fy fz radius mass f_fracRed[1] f_fracRed[2] f_fracRed[3]
+
 run         1
+
diff --git a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/GOD0k1_Hematite/DEM/in.liggghts_init b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/GOD0k1_Hematite/DEM/in.liggghts_init
index c56b3306..972be3fc 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/GOD0k1_Hematite/DEM/in.liggghts_init
+++ b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/GOD0k1_Hematite/DEM/in.liggghts_init
@@ -48,11 +48,13 @@ compute         rke all erotate/sphere
 thermo_style    custom step atoms ke c_rke vol
 thermo          1000
 thermo_modify   lost ignore norm no
-compute_modify  thermo_temp dynamic yes
+compute_modify  thermo_temp dynamic no extra 0
 
 # insert the first particles so that dump is not empty
 run     1
 dump    dmp all custom 100 ../DEM/post/dump.liggghts_init id type x y z vx vy vz fx fy fz radius mass
 
-run     10 upto 
-write_restart   ../DEM/post/restart/liggghts.restart
+run     10 upto
+
+write_restart ../DEM/post/restart/liggghts.restart
+
diff --git a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/GOD0k1_Hematite/DEM/in.liggghts_run b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/GOD0k1_Hematite/DEM/in.liggghts_run
index a3e52fc5..7d30e22e 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/GOD0k1_Hematite/DEM/in.liggghts_run
+++ b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/GOD0k1_Hematite/DEM/in.liggghts_run
@@ -196,8 +196,9 @@ compute         rke all erotate/sphere
 thermo_style    custom step atoms ke c_rke vol
 thermo          10
 thermo_modify   lost ignore norm no
-compute_modify  thermo_temp dynamic yes
+compute_modify  thermo_temp dynamic no extra 0
 
 dump    dmp all custom 100 ../DEM/post/dump*.liggghts_run id type x y z vx vy vz fx fy fz radius mass f_fracRed[1] f_fracRed[2] f_fracRed[3]
 
 run     1
+
diff --git a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/ISO4695/DEM/in.liggghts_init b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/ISO4695/DEM/in.liggghts_init
index 60a70170..972be3fc 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/ISO4695/DEM/in.liggghts_init
+++ b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/ISO4695/DEM/in.liggghts_init
@@ -48,7 +48,7 @@ compute         rke all erotate/sphere
 thermo_style    custom step atoms ke c_rke vol
 thermo          1000
 thermo_modify   lost ignore norm no
-compute_modify  thermo_temp dynamic yes
+compute_modify  thermo_temp dynamic no extra 0
 
 # insert the first particles so that dump is not empty
 run     1
@@ -56,4 +56,5 @@ dump    dmp all custom 100 ../DEM/post/dump.liggghts_init id type x y z vx vy vz
 
 run     10 upto
 
-write_restart   ../DEM/post/restart/liggghts.restart
+write_restart ../DEM/post/restart/liggghts.restart
+
diff --git a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/ISO4695/DEM/in.liggghts_run b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/ISO4695/DEM/in.liggghts_run
index 00466c6f..7860971b 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/ISO4695/DEM/in.liggghts_run
+++ b/tutorials/cfdemSolverRhoPimpleChem/SingleParticleCases/ISO4695/DEM/in.liggghts_run
@@ -196,8 +196,9 @@ compute         rke all erotate/sphere
 thermo_style    custom step atoms ke c_rke vol
 thermo          10
 thermo_modify   lost ignore norm no
-compute_modify  thermo_temp dynamic yes
+compute_modify  thermo_temp dynamic no extra 0
 
 dump    dmp all custom 100 ../DEM/post/dump*.liggghts_run id type x y z vx vy vz fx fy fz radius mass f_fracRed[1] f_fracRed[2] f_fracRed[3]
 
 run     1
+