I   irecord, nx_in, ny_in, mythid)

      do i=-1,nx_in+2

       arrayin(0,j)       = arrayin(nx_in,j)

       arrayin(-1,j)      = arrayin(nx_in-1,j)

       arrayin(nx_in+1,j) = arrayin(1,j)

       arrayin(nx_in+2,j) = arrayin(2,j)

      do i=-1,nx_in+2

         do i=1,nx_in

         NorthValue = NorthValue / nx_in

         do i=-1,nx_in+2

         do i=-1,nx_in+2

     &          xG(i,j,bi,bj) .ge. x_in(nx_in+1)   .or.

              print*,'nx_in,ny_in'    ,nx_in        ,ny_in

              print*,'x_in(0,nx_in+1)',x_in(0)      ,x_in(nx_in+1)

     I   nx_in, ny_in, method, mythid)

      integer       filePrec, irecord, nx_in, ny_in

      _RL      arrayin(-1:nx_in+2 ,      -1:ny_in+2)

      _RL      x_in   (-1:nx_in+2), y_in(-1:ny_in+2)

     I                irecord, nx_in, ny_in, myThid )

        length_of_rec=MDS_RECLEN( filePrec, nx_in*ny_in, myThid )

          READ(ioUnit,rec=irecord) (buffer_r4(i),i=1,nx_in*ny_in)

          CALL MDS_BYTESWAPR4(nx_in*ny_in,buffer_r4)

          READ(ioUnit,rec=irecord) (buffer_r8(i),i=1,nx_in*ny_in)

          CALL MDS_BYTESWAPR8(nx_in*ny_in,buffer_r8)

           CALL MPI_BCAST(buffer_r4,nx_in*ny_in,MPI_REAL,

           CALL MPI_BCAST(buffer_r8,nx_in*ny_in,MPI_DOUBLE_PRECISION,

          DO i=1,nx_in

          DO i=1,nx_in

          DO i=1,nx_in

          DO i=1,nx_in

          ijs = (j-1)*nx_in

          DO i=1,nx_in

          ijs = (j-1)*nx_in

          DO i=1,nx_in

      INTEGER       filePrec, irecord, nx_in, ny_in

       _RL          arrayin( -1:nx_in+2 , -1:ny_in+2 )

      Real*8, target ::  buffer_r8(nx_in,ny_in)

      Real*4, target ::  buffer_r4(nx_in,ny_in)

      Real*8   buffer_r8(nx_in,ny_in)

      Real*4   buffer_r4(nx_in,ny_in)

      IF ( nx_in*ny_in .GT. exf_interp_bufferSize ) THEN