MPI_COMM_mylocal=MPI_COMM_World
MPI_COMM_myglobal=MPI_COMM_World
call MPI_COMM_RANK( MPI_COMM_World, n, ierr )