diff --git a/src/constitutive.f90 b/src/constitutive.f90
index 4df97b240..23ae3f07b 100644
--- a/src/constitutive.f90
+++ b/src/constitutive.f90
@@ -476,19 +476,11 @@ subroutine constitutive_LpAndItsTangents(Lp, dLp_dS, dLp_dFi, &
 
  end select plasticityType
 
-#if defined(__INTEL_COMPILER) || defined(__PGI)
- forall(i = 1_pInt:3_pInt, j = 1_pInt:3_pInt)
-#else
- do concurrent(i = 1_pInt:3_pInt, j = 1_pInt:3_pInt)
-#endif
+ do i=1,3; do j=1,3
    dLp_dFi(i,j,1:3,1:3) = matmul(matmul(Fi,S),transpose(dLp_dMp(i,j,1:3,1:3))) + &
                           matmul(matmul(Fi,dLp_dMp(i,j,1:3,1:3)),S)
-   dLp_dS(i,j,1:3,1:3)  = matmul(matmul(transpose(Fi),Fi),dLp_dMp(i,j,1:3,1:3))       ! ToDo: @PS: why not:   dLp_dMp:(FiT Fi)
-#if defined(__INTEL_COMPILER) || defined(__PGI)
- end forall
-#else
- enddo
-#endif
+   dLp_dS(i,j,1:3,1:3)  = matmul(matmul(transpose(Fi),Fi),dLp_dMp(i,j,1:3,1:3))                     ! ToDo: @PS: why not:   dLp_dMp:(FiT Fi)
+ enddo; enddo
 
 end subroutine constitutive_LpAndItsTangents