ax = plt.subplot(111)

plt.title(title)

m = Basemap(resolution='l',
            projection='laea',
            llcrnrlat=y1,
            urcrnrlat=y2,
            llcrnrlon=x1,
            urcrnrlon=x2,
            lon_0=lon,
            lat_0=lat,
            lat_ts=(y1 + y2) / 2)
m.drawcountries(linewidth=0.5)
m.drawcoastlines(linewidth=0.5)
m.drawmeridians(np.arange(Round(x1, False), Round(x2, False), step),
                labels=[0, 0, 0, 1],
                color='black',
                linewidth=0.2)
m.drawparallels(np.arange(Round(y1, False),
                          Round(y2) + step, step),
                labels=[1, 0, 0, 0],
                color='black',
                linewidth=0.2)
m.shadedrelief()

# Loop through country list and add label for each
missing = []
for c in COUNTRIES.split(','):
    try:
        plotCountry(m, ax, c)
plt.title(title)

m = Basemap(resolution='i',
            projection='laea',
            llcrnrlat=y1,
            urcrnrlat=y2,
            llcrnrlon=x1,
            urcrnrlon=x2,
            lon_0=lon,
            lat_0=lat,
            lat_ts=(y1 + y2) / 2)
m.drawcountries(linewidth=0.5)
m.drawcoastlines(linewidth=0.5)
m.drawmeridians(np.arange(
    Round(x1, False, step) - step,
    Round(x2, True, step) + step, step),
                labels=[0, 0, 0, 1],
                color='black',
                linewidth=0.2)
m.drawparallels(np.arange(
    Round(y1, False, step) - step,
    Round(y2) + step, step),
                labels=[1, 0, 0, 0],
                color='black',
                linewidth=0.2)
m.shadedrelief()

# Loop through country list and add label for each
missing = []
for c in COUNTRIES.split(','):
plt.title(title)

m = Basemap(resolution='i',
            projection='laea',
            llcrnrlat=y1,
            urcrnrlat=y2,
            llcrnrlon=x1,
            urcrnrlon=x2,
            lon_0=lon,
            lat_0=lat,
            lat_ts=(y1 + y2) / 2)
m.drawcountries(linewidth=0.5)
m.drawcoastlines(linewidth=0.5)
m.drawmeridians(np.arange(
    Round(x1, False, step) - step,
    Round(x2, False, step) + 2 * step, step),
                labels=[0, 0, 0, 1],
                color='black',
                linewidth=0.2)
m.drawparallels(np.arange(
    Round(y1, False, step) - step,
    Round(y2) + step, step),
                labels=[1, 0, 0, 0],
                color='black',
                linewidth=0.2)
m.shadedrelief()

# Loop through country list and add label for each
missing = []
for c in COUNTRIES.split(','):
ax = plt.subplot(111)

plt.title(title)

m = Basemap(resolution='l',
            projection='laea',
            llcrnrlat=y1,
            urcrnrlat=y2,
            llcrnrlon=x1,
            urcrnrlon=x2,
            lon_0=lon,
            lat_0=lat,
            lat_ts=(y1 + y2) / 2)
m.drawcountries(linewidth=0.5)
m.drawcoastlines(linewidth=0.5)
m.drawmeridians(np.arange(Round(x1, False), Round(x2, False), 10),
                labels=[0, 0, 0, 1],
                color='black',
                linewidth=0.2)
m.drawparallels(np.arange(Round(y1, False),
                          Round(y2) + 10, 10),
                labels=[1, 0, 0, 0],
                color='black',
                linewidth=0.2)
m.shadedrelief()

# Loop through country list and add label for each
missing = []
for c in COUNTRIES.split(','):
    try:
        plotCountry(m, ax, c)