js  = exch2_myFace(is)

       IF ( js.NE.0 ) THEN

         IF ( js .NE. jp )

     &    '- facet', js, ' : X-size=', fNx, ' , Y-size=', fNy,

         jp = js

       js  = exch2_myFace(is)

       IF ( js.NE.0 ) THEN

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

        js  = exch2_myFace(is)

        IF ( js.NE.0 ) THEN

         IF ( js.NE.jp ) THEN

           jp = js

            nbTx = facet_dims(2*js-1)/tNx

          jt = INT(facet_link(i,js))

          ii = MOD( NINT(facet_link(i,js)*10.), 10 )

                 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)

               isbd1 = facet_pij(1,i,js)*itbd1

     &               + facet_pij(2,i,js)*jtbd1 + facet_oi(i,js)

               isbd2 = facet_pij(1,i,js)*itbd2

     &               + facet_pij(2,i,js)*jtbd2 + facet_oi(i,js)

               jsbd1 = facet_pij(3,i,js)*itbd1

     &               + facet_pij(4,i,js)*jtbd1 + facet_oj(i,js)

               jsbd2 = facet_pij(3,i,js)*itbd2

     &               + facet_pij(4,i,js)*jtbd2 + facet_oj(i,js)

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

       js  = exch2_myFace(is)

       IF ( js.NE.0 ) THEN

        js  = exch2_myFace(is)

          js = ny_in*.5

          do while (yG(i,j,bi,bj) .lt. y_in(js))

           js = (js - 1)*.5

          do while (yG(i,j,bi,bj) .ge. y_in(js+1))

           js = js + 1

          s_ind(i,j) = js

      integer  i, j, k, l, js, bi, bj, sp, interp_unit