N America cartopy blank using Natural earth for fill in

namap.py

"""
Features
--------

A demonstration of some of the built-in Natural Earth features found
in cartopy.

modified from caropy examples/features.py
"""
__tags__ = ['Lines and polygons']

import cartopy.crs as ccrs
import cartopy.feature as cfeature
import matplotlib.pyplot as plt


def main():
    fig = plt.figure()
    ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())
    ax.set_extent([-166, -105, 71, 35], crs=ccrs.PlateCarree())

    states_provinces = cfeature.NaturalEarthFeature(
        category='cultural',
        name='admin_1_states_provinces_lines',
        scale='50m',
        facecolor='none')

    roads = cfeature.NaturalEarthFeature(
        category='cultural',
        name='roads',
        scale='10m',
        facecolor='none')

    cstline = cfeature.NaturalEarthFeature(
        category="physical",
        name="land",
        scale="10m"
    )

    # ax.add_feature(cfeature.LAND)
    ax.add_feature(cfeature.OCEAN)
    # ax.add_feature(cfeature.COASTLINE)
    ax.add_feature(cstline, edgecolor="black", facecolor="tan")
    ax.add_feature(cfeature.BORDERS, linestyle=':')
    ax.add_feature(states_provinces, edgecolor="black")
    ax.add_feature(roads, edgecolor="gray")
    # ax.add_feature(cfeature.LAKES, alpha=0.5)
    # ax.add_feature(cfeature.RIVERS)

    plt.show()


if __name__ == '__main__':
    main()