def add_boundaries(ax, sector, lw=1, color='k'):
latn = sector['latn']
latx = sector['latx']
lonn = sector['lonn']
lonx = sector['lonx']
m = Basemap(projection='cyl',
lat_0=35,
lon_0=-130,
resolution='c',
area_thresh=0.1,
llcrnrlon=lonn,
llcrnrlat=latn,
urcrnrlon=lonx,
urcrnrlat=latx)
coastline = m.coastpolygons
for line in coastline:
ax.plot(line[0], line[1], linewidth=lw, color=color)
cnts, types = m._readboundarydata('countries')
countries = LineCollection(cnts, color=color, linewidth=lw)
ax.add_collection(countries)
