code: def crust_type_at(lat=None, lon=None):
      
      # Get lon into appropriate format
  
      lats = np.array(lat)
      lons = np.array(lon%360)
  
      iVals = ((90.0-lats)%180).astype(np.int)
      jVals = (lons%360.0).astype(int)
  
      # i = int((-lat+90.0)%180)
      # j = int(lon)
  
      t = _c1_crust_type_lat_lon[iVals,jVals]
  
      # t = _c1_crust_type_lat_lon[i,j]
      # des = litho.c1_region_descriptor[t]
  
      return t
  
  
  
  
  
      return t
docs: lat, lon (degrees)