31   DO i=1, N
 32 ! These are the neighbors in two directions:
 33      p(i)=MODULO(i-1+1,N)+1
 34      m(i)=MODULO(i-1-1,N)+1
 35      WRITE(*,*) m(i), i, p(i)
 36   END DO