In [34]:

    

import netCDF4
import matplotlib.pyplot as plt
import matplotlib.tri as Tri
import numpy as np
%matplotlib inline


    

In [35]:

    

url='http://geoport.whoi.edu/thredds/dodsC/usgs/vault0/models/tides/fvcom/spectral_tides.nc'


    

In [36]:

    

nc = netCDF4.Dataset(url)
ncv = nc.variables


    

In [38]:

    

# read node locations
lat = ncv['lat'][:]
lon = ncv['lon'][:]
# read connectivity array
nv = ncv['nv'][:].T - 1


    

In [39]:

    

# create a triangulation object, specifying the triangle connectivity array
tri = Tri.Triangulation(lon,lat, triangles=nv)


    

In [40]:

    

ncv['h'].shape


    

    
Out[40]:





(40992,)

In [41]:

    

# plot depth using tricontourf
h = nc.variables['h'][:]
fig = plt.figure(figsize=(12,8))
ax = fig.add_subplot(111,aspect=1.0/np.cos(lat.mean() * np.pi / 180.0))
plt.tricontourf(tri,-h,levels=range(-300,10,10))
plt.colorbar();