diff --git a/src/phase_mechanical_plastic_phenopowerlaw.f90 b/src/phase_mechanical_plastic_phenopowerlaw.f90
index 985a93c44..fe2fc2ef3 100644
--- a/src/phase_mechanical_plastic_phenopowerlaw.f90
+++ b/src/phase_mechanical_plastic_phenopowerlaw.f90
@@ -545,7 +545,7 @@ associate(prm => param(ph), stt => state(ph), dlt => deltastate(ph))
         neighbor_e = geom(ph)%IPneighborhood(1,n,en)
         
         if (any(dNeq(phase_O_0(ph)%data(en)%asQuaternion(), &
-        phase_O_0(ph)%data(neighbor_e)%asQuaternion()))) then
+        phase_O_0(ph)%data(neighbor_e)%asQuaternion()))) then                                             !> Selection of grain boundary elements
           Ability_Nucleation: if(stt%f_twin(twin_var,en)>(stt%fmc_twin(twin_var,en)+prm%checkstep)) then
     
             stt%fmc_twin(twin_var,en) = stt%fmc_twin(twin_var,en)+prm%checkstep