docs: def _get_lat_lon_indices(lsm_lat_array, lsm_lon_array, lat, lon):
      """
      Determines the index in the array (1D or 2D) where the
      lat/lon point is
      """
      if lsm_lat_array.ndim == 2 and lsm_lon_array.ndim == 2:
          lsm_lat_indices_from_lat, lsm_lon_indices_from_lat = \
              np.where((lsm_lat_array == lat))
          lsm_lat_indices_from_lon, lsm_lon_indices_from_lon = \
              np.where((lsm_lon_array == lon))
  
          index_lsm_grid_lat = np.intersect1d(lsm_lat_indices_from_lat,
                                              lsm_lat_indices_from_lon)[0]
          index_lsm_grid_lon = np.intersect1d(lsm_lon_indices_from_lat,
                                              lsm_lon_indices_from_lon)[0]
  
      elif lsm_lat_array.ndim == 1 and lsm_lon_array.ndim == 1:
          index_lsm_grid_lon = np.where(lsm_lon_array == lon)[0][0]
          index_lsm_grid_lat = np.where(lsm_lat_array == lat)[0][0]
      else:
          raise IndexError("Lat/Lon lists have invalid dimensions. "
                           "Only 1D or 2D arrays allowed ...")
  
      return index_lsm_grid_lat, index_lsm_grid_lon