diff --git a/Src/SiestaXC/mesh3d.F90 b/Src/SiestaXC/mesh3d.F90
index 148a6f9..1bd256f 100644
--- a/Src/SiestaXC/mesh3d.F90
+++ b/Src/SiestaXC/mesh3d.F90
@@ -2991,7 +2991,8 @@ subroutine setMeshDistr( distrID, nMesh, box, firstNode, nNodes, &
     ! Collect all node boxes and store them
     allocate( nodeBoxes(2,3,0:totNodes-1) )
     call gatherBoxes( box, nodeBoxes )
-    distr%box = nodeBoxes
+    distr%box(:,:,:) = nodeBoxes(:,:,:)
+    deallocate( nodeBoxes )
     goto 999  ! Exit, since no other arguments must be considered in this case
   end if ! (present(box))