Temporary arrays on large grids -> seg-fault

In certain calculations when the mesh-grid is extensively big a seg-fault will occur.

By tracking down it happens here:

        do iPart = 1,nParts
          partSize = sizeSum(iPart) - sizeSum(iPart-1)
          trsfBuff(sizeSum(iPart-1)+1:sizeSum(iPart)) = &
            reshape( srcData(srcComBox(1,1,iPart):srcComBox(2,1,iPart), &
                             srcComBox(1,2,iPart):srcComBox(2,2,iPart), &
                             srcComBox(1,3,iPart):srcComBox(2,3,iPart),:), &
                     (/partSize/) )
        end do ! iPart

(line 2617 of mesh3d.F90)

It is because the reshape creates a temporary array upon copy.

If all processors hit this point at the same time (unlikely, but may occur), the mesh grid will be duplicated which may turn out to be too much memory.

I think that all array assignments should be made into explicit loops to circumvent possible array-temporaries.