diff --git a/src/homogenization_mechanical_RGC.f90 b/src/homogenization_mechanical_RGC.f90
index 345592790..cd10aa5e2 100644
--- a/src/homogenization_mechanical_RGC.f90
+++ b/src/homogenization_mechanical_RGC.f90
@@ -778,9 +778,9 @@ pure function interfaceNormal(intFace,ho,en) result(n)
   associate (dst => dependentState(ho))
 
     n = 0.0_pReal
-    n(abs(intFace(1))) = real(intFace(1)/abs(intFace(1)),pReal)                                    ! get the normal vector w.r.t. cluster axis
+    n(abs(intFace(1))) = real(intFace(1)/abs(intFace(1)),pReal)                                     ! get the normal vector w.r.t. cluster axis
 
-    n = matmul(dst%orientation(1:3,1:3,en),n)                             ! map the normal vector into sample coordinate system (basis)
+    n = matmul(dst%orientation(1:3,1:3,en),n)                                                       ! map the normal vector into sample coordinate system (basis)
 
   end associate
 
diff --git a/src/phase_damage.f90 b/src/phase_damage.f90
index b69969898..b3f2c7c22 100644
--- a/src/phase_damage.f90
+++ b/src/phase_damage.f90
@@ -386,7 +386,7 @@ module function phase_K_phi(co,ce) result(K)
   real(pReal), dimension(3,3) :: K
   real(pReal), parameter :: l = 1.0_pReal
   
-  K = crystallite_push33ToRef(co,ce,param(material_phaseID(co,ce))%D) \
+  K = crystallite_push33ToRef(co,ce,param(material_phaseID(co,ce))%D) &
     * l**2.0_pReal
 
 end function phase_K_phi