In [1]:
import matplotlib.pyplot as plt
import netCDF4 as nc
%matplotlib inline
In [2]:
grid = nc.Dataset('https://salishsea.eos.ubc.ca/erddap/griddap/ubcSSnBathymetryV17-02')

By grid cell

In [3]:
plt.pcolormesh(grid.variables['bathymetry'][:])
plt.colorbar()
Out[3]:
<matplotlib.colorbar.Colorbar at 0x7fd0dcf634e0>

By lat/lon

In [4]:
plt.pcolormesh(grid.variables['longitude'][:], grid.variables['latitude'][:], 
              grid.variables['bathymetry'][:])
plt.colorbar()
Out[4]:
<matplotlib.colorbar.Colorbar at 0x7fd0dc133860>

First grid cell corresponds to lon/lat =

In [5]:
print(grid.variables['longitude'][0,0], grid.variables['latitude'][0,0])
-123.42943572998047 46.85966491699219
In [ ]: