diff --git a/src/IO.f90 b/src/IO.f90
index b6dbf5664..9ea35503b 100644
--- a/src/IO.f90
+++ b/src/IO.f90
@@ -604,7 +604,7 @@ subroutine IO_warning(warning_ID,el,ip,g,ext_msg)
       msg = 'read only the first document'
     case default
       msg = 'unknown warning number'
-    end select
+  end select
 
   !$OMP CRITICAL (write2out)
   write(IO_STDERR,'(/,a)')                ' ┌'//IO_DIVIDER//'┐'
@@ -658,7 +658,7 @@ subroutine selfTest
   if(any([1,1,1]     /= IO_stringPos('a')))         error stop 'IO_stringPos'
   if(any([2,2,3,5,5] /= IO_stringPos(' aa b')))     error stop 'IO_stringPos'
 
-  str=' 1.0 xxx'
+  str = ' 1.0 xxx'
   chunkPos = IO_stringPos(str)
   if(dNeq(1.0_pReal,IO_floatValue(str,chunkPos,1))) error stop 'IO_floatValue'
 
diff --git a/src/YAML_types.f90 b/src/YAML_types.f90
index 3d83831e6..2a6bd64f9 100644
--- a/src/YAML_types.f90
+++ b/src/YAML_types.f90
@@ -266,7 +266,7 @@ subroutine selfTest
         if(any(dNeq(l2%get_as1dFloat(1),[2.0_pReal,3.0_pReal]))) error stop 'byIndex_as1dFloat'
         call l2%append(l3)
         x = l2%as2dFloat()
-        if(x(2,1)/= 4.0_pReal)                                   error stop 'byKey_as2dFloat'
+        if(dNeq(x(2,1),4.0_pReal))                               error stop 'byKey_as2dFloat'
         if(any(dNeq(pack(l2%as2dFloat(),.true.),&
                [2.0_pReal,4.0_pReal,3.0_pReal,5.0_pReal])))      error stop 'byKey_as2dFloat'
         n => l2
diff --git a/src/homogenization.f90 b/src/homogenization.f90
index 9be7c0abf..e5d3d034e 100644
--- a/src/homogenization.f90
+++ b/src/homogenization.f90
@@ -39,8 +39,6 @@ module homogenization
     thermal_type                                                                                    !< thermal transport model
   integer(kind(DAMAGE_none_ID)),              dimension(:),   allocatable :: &
     damage_type                                                                                     !< nonlocal damage model
-  integer(kind(HOMOGENIZATION_undefined_ID)), dimension(:),   allocatable :: &
-    homogenization_type                                                                             !< type of each homogenization
 
   type, private :: tNumerics_damage
     real(pReal) :: &
diff --git a/src/homogenization_mechanical.f90 b/src/homogenization_mechanical.f90
index 3b75a66f0..c19695d3d 100644
--- a/src/homogenization_mechanical.f90
+++ b/src/homogenization_mechanical.f90
@@ -71,6 +71,9 @@ submodule(homogenization) mechanical
 
   end interface
 
+  integer(kind(HOMOGENIZATION_undefined_ID)), dimension(:),   allocatable :: &
+    homogenization_type                                                                             !< type of each homogenization
+
 contains
 
 !--------------------------------------------------------------------------------------------------
diff --git a/src/phase.f90 b/src/phase.f90
index 81e5ab250..c95e828b4 100644
--- a/src/phase.f90
+++ b/src/phase.f90
@@ -332,8 +332,7 @@ contains
 subroutine phase_init
 
   integer :: &
-    ph, &                                                                                            !< counter in phase loop
-    so                                                                                               !< counter in source loop
+    ph
   class (tNode), pointer :: &
     debug_constitutive, &
     materials, &
@@ -475,7 +474,6 @@ subroutine crystallite_init()
     co, &                                                                                           !< counter in integration point component loop
     ip, &                                                                                           !< counter in integration point loop
     el, &                                                                                           !< counter in element loop
-    so, &
     cMax, &                                                                                         !< maximum number of  integration point components
     iMax, &                                                                                         !< maximum number of integration points
     eMax                                                                                            !< maximum number of elements