diff --git a/src/mesh_marc.f90 b/src/mesh_marc.f90
index 9f6c1bcba..0d082a0f1 100644
--- a/src/mesh_marc.f90
+++ b/src/mesh_marc.f90
@@ -82,9 +82,8 @@ subroutine mesh_init(ip,el)
   if (debug_e < 1 .or. debug_e > nElems)    call IO_error(602,ext_msg='element')
   if (debug_i < 1 .or. debug_i > elem%nIPs) call IO_error(602,ext_msg='IP')
  
-  FEsolving_execElem = [ 1,nElems ]                                                                 ! parallel loop bounds set to comprise all DAMASK elements
-  allocate(FEsolving_execIP(2,nElems), source=1)                                                    ! parallel loop bounds set to comprise from first IP...
-  FEsolving_execIP(2,:) = elem%nIPs
+  FEsolving_execElem =        [1,nElems]
+  FEsolving_execIP   = spread([1,elem%nIPs],2,nElems)
  
   allocate(calcMode(elem%nIPs,nElems),source=.false.)                                               ! pretend to have collected what first call is asking (F = I)
   calcMode(ip,mesh_FEasCP('elem',el)) = .true.                                                      ! first ip,el needs to be already pingponged to "calc"