The storage is contiguous but the domain decomposition and processor assignment are controlled by the library. And it is also in fortran 95. I will try to make it work with the WithArray methods and I'll let you know how it goes :) Thanks for the quick reply!