C NCLFORTSTART
      SUBROUTINE get_ij( lat_array,long_array,
     &           lat,longitude,i_point,j_point,ny,nx )
      IMPLICIT NONE
      INTEGER nx , ny
      REAL    lat_array(nx,ny), long_array(nx,ny)
      REAL    lat, longitude
      integer i_point,j_point
C NCLEND
      REAL    longd, latd
      integer i,j
      real ir,jr
      real dist_min, dist
      real w00, w01, w02, w10, w11, w12, w20, w21, w22
      dist_min = 1.e+20
      do j=1,ny
      do i=1,nx
        latd = (lat_array(i,j)-lat)**2
        longd = amin1((long_array(i,j)-longitude)**2,
     &                (long_array(i,j)+longitude)**2 )
        dist = sqrt(latd + longd)
        if(dist_min .gt. dist) then
          dist_min = dist
          ir = float(i)
          jr = float(j)
        end if
      enddo
      enddo
      i_point = nint(ir)
      j_point = nint(jr)

! we will just return the nearest point at present
      return
      end