def METHOD_NAME(filename):
    # Open the file using the NetCDF4 library
    nc = Dataset(filename)
    # analyze_dataset(nc)
    latlon = nc.variables['geospatial_lat_lon_extent']
    # print('geospatial_lat_lon_extent:', latlon)
    # Original example uses Brightness Temperature values, 'CMI'
    # But that's not in our dataset.
    whichvar = 'Rad'
    data = nc.variables[whichvar][:]
    prettyname = nc.variables[whichvar].long_name
    # Show it without basemap:
    # plt.imshow(data, cmap='Greys')
    # plt.show()
    # Create the basemap reference for the Satellite Projection
    bmap = Basemap(projection='geos',
                   # Extents
                   llcrnrlon=latlon.geospatial_westbound_longitude,
                   llcrnrlat=latlon.geospatial_southbound_latitude,
                   urcrnrlon=latlon.geospatial_eastbound_longitude,
                   urcrnrlat=latlon.geospatial_northbound_latitude,
                   # Should these be nadir, or center?
                   lon_0=latlon.geospatial_lon_nadir,
                   lat_0=latlon.geospatial_lat_nadir,
                   satellite_height=35786023.0, ellps='GRS80')
    # Plot GOES-16 Channel using 170 and 378 as the temperature thresholds
    bmap.imshow(data, origin='upper', cmap='Greys')
    # Draw the coastlines, countries, parallels and meridians
    bmap.drawcoastlines(linewidth=0.3, linestyle='solid', color='blue')
    bmap.drawcountries(linewidth=0.3, linestyle='solid', color='blue')
    # These don't seem to do anything
    bmap.drawparallels(np.arange(-90.0, 90.0, 10.0),
                       linewidth=0.1, color='yellow')
    bmap.drawmeridians(np.arange(0.0, 360.0, 10.0),
                       linewidth=0.1, color='yellow')
    # Insert the legend
    bmap.colorbar(location='bottom', label=prettyname)
    # Show the plot
    plt.show()