code: def circle(m, centerlon, centerlat, radius, *args, **kwargs):
      
      glon1 = centerlon
      glat1 = centerlat
      X = []
      Y = []
      for azimuth in range(0, 360):
          glon2, glat2, baz = _gccalc(glon1, glat1, azimuth, radius)
          X.append(glon2)
          Y.append(glat2)
      X.append(X[0])
      Y.append(Y[0])
  
      proj_x, proj_y = m(X,Y)
      return zip(proj_x, proj_y)
docs: Return lon, lat tuples of a "circle" which matches the chosen Basemap projection
      Takes the following arguments:
      m = basemap instance
      centerlon = originating lon
      centrelat = originating lat
      radius = radius