diff --git a/src/constitutive.f90 b/src/constitutive.f90
index 697275367..7ddceb098 100644
--- a/src/constitutive.f90
+++ b/src/constitutive.f90
@@ -432,11 +432,11 @@ subroutine constitutive_hooke_SandItsTangents(S, dS_dFe, dS_dFi, &
  
   E = 0.5_pReal*(matmul(transpose(Fe),Fe)-math_I3)                                                  !< Green-Lagrange strain in unloaded configuration
   S = math_mul3333xx33(C,matmul(matmul(transpose(Fi),E),Fi))                                        !< 2PK stress in lattice configuration in work conjugate with GL strain pulled back to lattice configuration
- 
-  forall (i=1:3, j=1:3)
+
+  do i =1, 3;do j=1,3
     dS_dFe(i,j,1:3,1:3) = matmul(Fe,matmul(matmul(Fi,C(i,j,1:3,1:3)),transpose(Fi)))                !< dS_ij/dFe_kl = C_ijmn * Fi_lm * Fi_on * Fe_ko
     dS_dFi(i,j,1:3,1:3) = 2.0_pReal*matmul(matmul(E,Fi),C(i,j,1:3,1:3))                             !< dS_ij/dFi_kl = C_ijln * E_km * Fe_mn
-  end forall
+  enddo; enddo
 
 end subroutine constitutive_hooke_SandItsTangents