title: indices of 2D lat lon data
body: <p>I am trying to find the equivalent (if there exists one) of an NCL function that returns the indices of two-dimensional latitude/longitude arrays closest to a user-specified latitude/longitude coordinate pair.</p>
  
  <p>This is the link to the NCL function that I am hoping there is an equivalent to in python.  I'm suspecting at this point that there is not, so any tips on how to get indices from lat/lon coordinates is appreciated</p>
  
  <p><a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/getind_latlon2d.shtml" rel="nofollow noreferrer">https://www.ncl.ucar.edu/Document/Functions/Contributed/getind_latlon2d.shtml</a></p>
  
  <p>Right now , I have my coordinate values saved into an .nc file and are read by:</p>
  
  <pre><code>coords='coords.nc'
  fh = Dataset(coords, mode='r')
  lons = fh.variables['g5_lon_1'][:,:]
  lats = fh.variables['g5_lat_0'][:,:]
  rot = fh.variables['g5_rot_2'][:,:]
  fh.close()
  </code></pre>