JUst a reminder: this writeDomain call will do an MPI_Allreduce internally to detect where each participant should store its information in the global target buffer. If you do not want this and you already know the size/offset of each rank, use writeDomain(currentStep, globalSize, globalOffset, ...) interface instead.