& mpi_myXGlobalLo, mpi_myYGlobalLo
INTEGER mpi_myYGlobalLo(nPx*nPy)
bjG = (mpi_myYGlobalLo(np1)-1)/sNy+bj
jG = mpi_myYGlobalLo(npe+1)-1+(bi-1)*sNy+j
bjG = (mpi_myYGlobalLo(np)-1)/sNy+bj
mpi_myYGlobalLo(np) = itemp(2)