diff --git a/src/make_matrix.f90 b/src/make_matrix.f90
index bb5a16c533b59129429eb1396e52c77ebaf33bdf..ec2406c4a395fdd4fa04989a9693d5b96c4ef9c6 100644
--- a/src/make_matrix.f90
+++ b/src/make_matrix.f90
@@ -111,7 +111,7 @@ double precision weight
 !(((((((((((((((( declaration of the subroutine internal variables )))))))))))))
 !------------------------------------------------------------------------------|
 logical is_plastic_temp,flag_vrm_pb
-integer i,k,k1,k2,k3,ii,jj,ij,iint,nint,nb,err,ic
+integer i,j,k,k1,k2,k3,ii,jj,ij,iint,nint,nb,err,ic
 integer iproc,nproc,ierr,mpe
 double precision eps,yield,fixt,yieldnow
 double precision r,s,t,w,volume,vol
diff --git a/src/pressure_cut.f90 b/src/pressure_cut.f90
index d4e23c14c940d179e91bec6180cbfbf5ea3216ba..6d132ab03896d1b266ef19e705d33031e1b2ee56 100644
--- a/src/pressure_cut.f90
+++ b/src/pressure_cut.f90
@@ -471,7 +471,7 @@ double precision,intent(inout) :: pressure
 !(((((((((((((((( declaration of the subroutine internal variables )))))))))))))
 !------------------------------------------------------------------------------|
 
-integer :: err,i,k,ii,jj,kk,levelp,jcut
+integer :: err,i,j,k,ii,jj,kk,levelp,jcut
 double precision,allocatable :: lsfp(:,:)
 double precision :: r(params%mpe),s(params%mpe),t(params%mpe),h(params%mpe)
 double precision :: r0p,s0p,t0p,rstp