diff --git a/src/homogenization.f90 b/src/homogenization.f90
index 1f9fb2221..8c9ecaecb 100644
--- a/src/homogenization.f90
+++ b/src/homogenization.f90
@@ -237,7 +237,7 @@ subroutine homogenization_mechanical_response(Delta_t,FEsolving_execIP,FEsolving
     doneAndHappy
 
 
-  !$OMP PARALLEL DO PRIVATE(ce,en,ho,NiterationMPstate,converged,doneAndHappy)
+  !$OMP PARALLEL DO PRIVATE(ce,co,en,ho,NiterationMPstate,converged,doneAndHappy)
   do el = FEsolving_execElem(1),FEsolving_execElem(2)
 
     do ip = FEsolving_execIP(1),FEsolving_execIP(2)
diff --git a/src/phase.f90 b/src/phase.f90
index c4e7142f8..9943c777a 100644
--- a/src/phase.f90
+++ b/src/phase.f90
@@ -554,7 +554,7 @@ subroutine crystallite_init()
   flush(IO_STDOUT)
 
 
-  !$OMP PARALLEL DO PRIVATE(ce)
+  !$OMP PARALLEL DO PRIVATE(ce,ph,en)
   do el = 1, eMax
     do ip = 1, iMax
       ce = (el-1)*discretization_nIPs + ip