docs: def get_mesh(sitecol, complete=True):
      """
      :returns:
          a lon-lat or lon-lat-depth array depending if the site collection
          is at sea level or not
      """
      sc = sitecol.complete if complete else sitecol
      if sc.at_sea_level():
          mesh = numpy.zeros(len(sc), [('lon', F64), ('lat', F64)])
          mesh['lon'] = sc.lons
          mesh['lat'] = sc.lats
      else:
          mesh = numpy.zeros(len(sc), [('lon', F64), ('lat', F64),
                                       ('depth', F64)])
          mesh['lon'] = sc.lons
          mesh['lat'] = sc.lats
          mesh['depth'] = sc.depths
      return mesh