Large scale 3D data read and write using MPI

Question

I am trying to combine two 3D matrices (A and B) into one 3D matrix (C). Due to the nature of my problem, it\'s suitable to separate the total job into px * py