diff --git a/code/IO.f90 b/code/IO.f90
index 86961a785..59632cd0a 100644
--- a/code/IO.f90
+++ b/code/IO.f90
@@ -262,7 +262,6 @@ endsubroutine
    IO_hybridIA(1,i) = deltas(1)*(mod(bin/(steps(3)*steps(2)),steps(1))+center)  ! phi1
    IO_hybridIA(2,i) = deltas(2)*(mod(bin/ steps(3)          ,steps(2))+center)  ! Phi
    IO_hybridIA(3,i) = deltas(3)*(mod(bin                    ,steps(3))+center)  ! phi2
-!    write(6,'(i6,x,a3,i5,x,a5,x,3(f8.5,x))') i,'bin',bin,'euler',IO_hybridIA(:,i)
    binSet(j) = binSet(i)
  enddo
  close(999)
diff --git a/code/crystallite.f90 b/code/crystallite.f90
index 560515477..3973c3b19 100644
--- a/code/crystallite.f90
+++ b/code/crystallite.f90
@@ -86,8 +86,7 @@ subroutine crystallite_init(Temperature)
  use debug, only:            debug_info, &
                              debug_reset
  use math, only:             math_I3, &
-                             math_EulerToR, &
-                             math_RToEuler
+                             math_EulerToR
  use FEsolving, only:        FEsolving_execElem, &
                              FEsolving_execIP
  use mesh, only:             mesh_element, &
@@ -1044,7 +1043,7 @@ subroutine crystallite_stressAndItsTangent(updateJaco)
               k = kl(1,mod(cycleCounter/iJacoStiffness,9)+1,g,i,e)
               l = kl(2,mod(cycleCounter/iJacoStiffness,9)+1,g,i,e)      
               crystallite_dPdF(:,:,k,l,g,i,e) = (crystallite_P(:,:,g,i,e) - storedP(:,:,g,i,e))/pert_Fg                             ! ... use tangent dP_ij/dFg_kl
-            elseif (.not. storedConvergenceFlag(g,i,e)) then                                                                        ! if crystallite didn�t converge before...
+            elseif (.not. storedConvergenceFlag(g,i,e)) then                                                                        ! if crystallite didnÕt converge before...
               crystallite_dPdF(:,:,:,:,g,i,e) = crystallite_fallbackdPdF(:,:,:,:,g,i,e)                                             ! ... use (elastic) fallback
             endif
       enddo; enddo; enddo
@@ -1550,8 +1549,6 @@ use prec, only:                       pInt, &
                                       pReal
 use math, only:                       math_pDecomposition, &
                                       math_RtoQuaternion, &
-                                      math_RtoEuler, &
-                                      math_QuaternionToEuler, &
                                       math_QuaternionDisorientation, &
                                       inDeg, &
                                       math_qConj
@@ -1560,7 +1557,6 @@ use FEsolving, only:                  FEsolving_execElem, &
 use IO, only:                         IO_warning
 use material, only:                   material_phase, &
                                       homogenization_Ngrains, &
-                                      material_EulerAngles, &
                                       phase_constitution
 use mesh, only:                       mesh_element, &
                                       mesh_ipNeighborhood, &
@@ -1592,7 +1588,6 @@ logical error
 
 
 !$OMP PARALLEL DO
-  n = 0_pInt
   do e = FEsolving_execElem(1),FEsolving_execElem(2)
     do i = FEsolving_execIP(1,e),FEsolving_execIP(2,e)
       do g = 1,homogenization_Ngrains(mesh_element(3,e))
diff --git a/code/material.f90 b/code/material.f90
index a42207c3c..d841001f0 100644
--- a/code/material.f90
+++ b/code/material.f90
@@ -667,9 +667,6 @@ subroutine material_populateGrains()
            constituentGrain = constituentGrain + myNorientations
 
          endif
-!          do j = 1,myNorientations
-!            write(6,'(i6,x,a5,x,3(f8.5,x))') j,'euler',orientationOfGrain(:,grain+j)
-!          enddo
 ! ----------------------------------------------------------------------------
          symExtension = texture_symmetry(textureID) - 1_pInt
          if (symExtension > 0_pInt) then                                      ! sample symmetry
@@ -690,7 +687,6 @@ subroutine material_populateGrains()
        enddo  ! constituent
 
 ! ----------------------------------------------------------------------------
-
        if (.not. microstructure_elemhomo(micro)) then                           ! unless element homogeneous, reshuffle grains
          do i=1,myNgrains-1                                                     ! walk thru grains
            call random_number(rnd)
@@ -724,7 +720,6 @@ subroutine material_populateGrains()
                material_phase(g,i,e) = phaseOfGrain(grain+(i-1)*dGrains+g)
                material_EulerAngles(:,g,i,e) = orientationOfGrain(:,grain+(i-1)*dGrains+g)
              end forall
-
              grain = grain + FE_Nips(mesh_element(2,e)) * dGrains               ! wind forward by Nips*NgrainsPerIP
            endif
            ! write (6,*) e