&                    .AND. exch2_jHi(ns,is).EQ.(tNy+1) ) THEN

            np(1) = ns

     &                    .AND. exch2_jLo(ns,is).EQ. 0

     &                    .AND. exch2_jHi(ns,is).EQ. 0 ) THEN

            np(2) = ns

     &                    .AND. exch2_iLo(ns,is).EQ.(tNx+1)

     &                    .AND. exch2_iHi(ns,is).EQ.(tNx+1) ) THEN

            np(3) = ns

     &                    .AND. exch2_iLo(ns,is).EQ. 0

     &                    .AND. exch2_iHi(ns,is).EQ. 0 ) THEN

            np(4) = ns

         DO ns=1,exch2_nNeighbours(is)

     &     ' ns:',ns,' it=',exch2_neighbourId(ns,is),

     &     ', iLo,iHi=', exch2_iLo(ns,is), exch2_iHi(ns,is),

     &     ', jLo,jHi=', exch2_jLo(ns,is), exch2_jHi(ns,is)

      INTEGER i, j, js, jp, jt, ii, is, it, ns, nt, k, tx, ty

         DO ns=1,exch2_nNeighbours(is)

     &                    .AND. exch2_jLo(ns,is).EQ.(tNy+1)

               ns = MIN(nbNeighb,W2_maxNeighbours)

               exch2_neighbourId(ns,is) = it

               tile_edge2edge(ns,is) = 10*i + ii

               exch2_pij(1,ns,is) = 1

               exch2_pij(2,ns,is) = 0

               exch2_pij(3,ns,is) = 0

               exch2_pij(4,ns,is) = 1

               exch2_oi(ns,is) = 0

               exch2_oj(ns,is) = 0

               exch2_iLo(ns,is) = ii1 - ddi - exch2_tBasex(is)

               exch2_iHi(ns,is) = ii2 + ddi - exch2_tBasex(is)

               exch2_jLo(ns,is) = jj1 - ddj - exch2_tBasey(is)

               exch2_jHi(ns,is) = jj2 + ddj - exch2_tBasey(is)

               ns = MIN(nbNeighb,W2_maxNeighbours)

               exch2_neighbourId(ns,is)  = it

               tile_edge2edge(ns,is) = 10*i + ii

                 exch2_pij(k,ns,is) = facet_pij(k,i,js)

               exch2_oi(ns,is) = facet_oi(i,js)

               exch2_oj(ns,is) = facet_oj(i,js)

               exch2_iLo(ns,is) = isbd1 - ddi - exch2_tBasex(is)

               exch2_iHi(ns,is) = isbd2 + ddi - exch2_tBasex(is)

               exch2_jLo(ns,is) = jsbd1 - ddj - exch2_tBasey(is)

               exch2_jHi(ns,is) = jsbd2 + ddj - exch2_tBasey(is)

         DO ns=1,MIN(nbNeighb,W2_maxNeighbours)

     &     ' ns:',ns,' it=',exch2_neighbourId(ns,is),

     &     ', iLo,iHi=', exch2_iLo(ns,is), exch2_iHi(ns,is),

     &     ', jLo,jHi=', exch2_jLo(ns,is), exch2_jHi(ns,is)

       DO ns=1,exch2_nNeighbours(is)

        i  = tile_edge2edge(ns,is)/10

        ii = MOD(tile_edge2edge(ns,is),10)

          exch2_neighbourDir(ns,is) = i

        it = exch2_neighbourId(ns,is)

           IF ( exch2_opposingSend(ns,is).EQ.0 ) THEN

            exch2_opposingSend(ns,is) = nt

            nn = exch2_opposingSend(ns,is)

     &      ns,' (',tile_edge2edge(ns,is),' ) has multiple connection'

        IF ( exch2_opposingSend(ns,is).EQ.0 ) THEN

     &      ns,' (',tile_edge2edge(ns,is),' ) no connection from',it

       DO ns=1,exch2_nNeighbours(is)

        it = exch2_neighbourId(ns,is)

        nt = exch2_opposingSend(ns,is)

        IF ( ii.NE.is .OR. nn.NE.ns ) THEN

     &      ns,' (',tile_edge2edge(ns,is),' ) connected'

      INTEGER is, js, ns, it, jt, nt, tx, ty

       DO ns=1,W2_maxNeighbours

        tile_edge2edge(ns,is) = 0

        exch2_neighbourDir(ns,is) = 0

      DO ns=1,nbStrips

       DO j=jMinUpd(ns),jMaxUpd(ns)

      DO j=jMinUpd(ns),jMaxUpd(ns)

      DO j=jMinUpd(ns),jMaxUpd(ns)

      DO j=jMinUpd(ns),jMaxUpd(ns)

      INTEGER ns, nbStrips

      DO ns=1,nbStrips

        DO i=iMinUpd(ns),iMaxUpd(ns)

       DO i=iMinUpd(ns),iMaxUpd(ns)

       DO i=iMinUpd(ns),iMaxUpd(ns)

       DO i=iMinUpd(ns),iMaxUpd(ns)

      INTEGER ns, nbStrips

      integer ndv, did, ns,ne, n1,n2, indg, indv

        ns = IFNBLNK(name)

        CALL MNC_DIM_INIT_ALL(fname,name(ns:ne),dlen,'N', myThid)

          IF (((ne-ns) .EQ. (n2-n1)) 

     &         .AND. (mnc_d_names(j)(ns:ne) .EQ. name(ns:ne))) THEN