diff --git a/code/constitutive_nonlocal.f90 b/code/constitutive_nonlocal.f90
index 31c07292c..7afc99a24 100644
--- a/code/constitutive_nonlocal.f90
+++ b/code/constitutive_nonlocal.f90
@@ -860,7 +860,7 @@ do i = 1,maxNinstance
                                                 + 3.0_pReal*constitutive_nonlocal_Cslip_66(4,4,i) )                                          ! (C11iso-C12iso)/2 with C11iso=(3*C11+2*C12+4*C44)/5 and C12iso=(C11+4*C12-2*C44)/5
   constitutive_nonlocal_nu(i) =   ( constitutive_nonlocal_Cslip_66(1,1,i) + 4.0_pReal*constitutive_nonlocal_Cslip_66(1,2,i) &
                                     - 2.0_pReal*constitutive_nonlocal_Cslip_66(1,2,i) ) &
-                                / ( 4.0_pReal*constitutive_nonlocal_Cslip_66(1,1,i) + 6.0_pReal*constitutive_nonlocal_Cslip_66(1,2,i) &
+                     / ( 4.0_pReal*constitutive_nonlocal_Cslip_66(1,1,i) + 6.0_pReal*constitutive_nonlocal_Cslip_66(1,2,i) &
                                     + 2.0_pReal*constitutive_nonlocal_Cslip_66(4,4,i) )    
   constitutive_nonlocal_Cslip_66(1:6,1:6,i) = math_Mandel3333to66(math_Voigt66to3333(constitutive_nonlocal_Cslip_66(1:6,1:6,i)))
   constitutive_nonlocal_Cslip_3333(1:3,1:3,1:3,1:3,i) = math_Voigt66to3333(constitutive_nonlocal_Cslip_66(1:6,1:6,i))
diff --git a/code/lattice.f90 b/code/lattice.f90
index 5d32d9e28..4a9a63d59 100644
--- a/code/lattice.f90
+++ b/code/lattice.f90
@@ -752,19 +752,15 @@ pure function lattice_symmetrizeC66(structName,C66)
  select case(structName(1:3))
    case ('iso')
      forall(k=1_pInt:3_pInt)
-       forall(j=1_pInt:3_pInt) 
-         lattice_symmetrizeC66(k,j) = C66(1,2)
-       end forall
+       forall(j=1_pInt:3_pInt) lattice_symmetrizeC66(k,j) = C66(1,2)
        lattice_symmetrizeC66(k,k) = C66(1,1)
        lattice_symmetrizeC66(k+3,k+3) = 0.5_pReal*(C66(1,1)-C66(1,2))
      end forall
    case ('fcc','bcc')
      forall(k=1_pInt:3_pInt)
-       forall(j=1_pInt:3_pInt)
-         lattice_symmetrizeC66(k,j) =   C66(1,2)
-         lattice_symmetrizeC66(k,k) =     C66(1,1)
-         lattice_symmetrizeC66(k+3_pInt,k+3_pInt) = C66(4,4)
-       end forall
+       forall(j=1_pInt:3_pInt) lattice_symmetrizeC66(k,j) =   C66(1,2)
+       lattice_symmetrizeC66(k,k) =     C66(1,1)
+       lattice_symmetrizeC66(k+3_pInt,k+3_pInt) = C66(4,4)
      end forall    
    case ('hex')
      lattice_symmetrizeC66(1,1) = C66(1,1)