diff --git a/src/grid/grid_mech_FEM.f90 b/src/grid/grid_mech_FEM.f90
index bd24ddadc..eb38e3e65 100644
--- a/src/grid/grid_mech_FEM.f90
+++ b/src/grid/grid_mech_FEM.f90
@@ -46,7 +46,7 @@ module grid_mech_FEM
 
   type(tNumerics), private :: num
   logical, private:: &
-    debug_rotation 
+    debugRotation 
 
 !--------------------------------------------------------------------------------------------------
 ! PETSc data
@@ -128,7 +128,7 @@ subroutine grid_mech_FEM_init
 !-----------------------------------------------------------------------------------------------
 ! debugging options
   debug_grid => debug_root%get('grid', defaultVal=emptyList)
-  debug_rotation = debug_grid%contains('rotation')
+  debugRotation = debug_grid%contains('rotation')
  
 !-------------------------------------------------------------------------------------------------
 ! read numerical parameter and do sanity checks
@@ -518,7 +518,7 @@ subroutine formResidual(da_local,x_local, &
   newIteration: if (totalIter <= PETScIter) then
     totalIter = totalIter + 1
     write(6,'(1x,a,3(a,i0))') trim(incInfo), ' @ Iteration ', num%itmin, '≤',totalIter+1, '≤', num%itmax
-    if (debug_rotation) &
+    if (debugRotation) &
       write(6,'(/,a,/,3(3(f12.7,1x)/))',advance='no') &
               ' deformation gradient aim (lab) =', transpose(params%rotation_BC%rotate(F_aim,active=.true.))
     write(6,'(/,a,/,3(3(f12.7,1x)/))',advance='no') &
diff --git a/src/grid/grid_mech_spectral_basic.f90 b/src/grid/grid_mech_spectral_basic.f90
index a0392e452..ad57f36c0 100644
--- a/src/grid/grid_mech_spectral_basic.f90
+++ b/src/grid/grid_mech_spectral_basic.f90
@@ -45,7 +45,7 @@ module grid_mech_spectral_basic
   type(tNumerics) :: num                                                                            ! numerics parameters. Better name?
 
   logical, private:: &
-    debug_rotation 
+    debugRotation 
 
 !--------------------------------------------------------------------------------------------------
 ! PETSc data
@@ -123,7 +123,7 @@ subroutine grid_mech_spectral_basic_init
 !-------------------------------------------------------------------------------------------------
 ! debugging options
   debug_grid => debug_root%get('grid', defaultVal=emptyList)
-  debug_rotation = debug_grid%contains('rotation')
+  debugRotation = debug_grid%contains('rotation')
   
 !-------------------------------------------------------------------------------------------------
 ! read numerical parameters and do sanity checks
@@ -478,7 +478,7 @@ subroutine formResidual(in, F, &
   newIteration: if (totalIter <= PETScIter) then
     totalIter = totalIter + 1
     write(6,'(1x,a,3(a,i0))') trim(incInfo), ' @ Iteration ', num%itmin, '≤',totalIter, '≤', num%itmax
-    if (debug_rotation) &
+    if (debugRotation) &
       write(6,'(/,a,/,3(3(f12.7,1x)/))',advance='no') &
               ' deformation gradient aim (lab) =', transpose(params%rotation_BC%rotate(F_aim,active=.true.))
     write(6,'(/,a,/,3(3(f12.7,1x)/))',advance='no') &
diff --git a/src/grid/grid_mech_spectral_polarisation.f90 b/src/grid/grid_mech_spectral_polarisation.f90
index 52f6cd06f..555d5e0e1 100644
--- a/src/grid/grid_mech_spectral_polarisation.f90
+++ b/src/grid/grid_mech_spectral_polarisation.f90
@@ -50,7 +50,7 @@ module grid_mech_spectral_polarisation
 
   type(tNumerics) :: num                                                                            ! numerics parameters. Better name?
 
-  logical, private :: debug_rotation
+  logical, private :: debugRotation
 
 !--------------------------------------------------------------------------------------------------
 ! PETSc data
@@ -133,7 +133,7 @@ subroutine grid_mech_spectral_polarisation_init
 !------------------------------------------------------------------------------------------------
 ! debugging options
   debug_grid => debug_root%get('grid',defaultVal=emptyList)
-  debug_rotation = debug_grid%contains('rotation')
+  debugRotation = debug_grid%contains('rotation')
 
 !-------------------------------------------------------------------------------------------------
 ! read numerical parameters
@@ -562,7 +562,7 @@ subroutine formResidual(in, FandF_tau, &
   newIteration: if (totalIter <= PETScIter) then
     totalIter = totalIter + 1
     write(6,'(1x,a,3(a,i0))') trim(incInfo), ' @ Iteration ', num%itmin, '≤',totalIter, '≤', num%itmax
-    if(debug_rotation) &
+    if(debugRotation) &
       write(6,'(/,a,/,3(3(f12.7,1x)/))',advance='no') &
               ' deformation gradient aim (lab) =', transpose(params%rotation_BC%rotate(F_aim,active=.true.))
     write(6,'(/,a,/,3(3(f12.7,1x)/))',advance='no') &