In [1]:
import netCDF4
url='http://geoport.whoi.edu/thredds/dodsC/usgs/data2/rsignell/models/maurer.ncml'
nc=netCDF4.Dataset(url).variables

In [2]:
nc.keys()


Out[2]:
[u'time', u'lev', u'lat', u'lon', u'temperature', u'salinity']

In [3]:
t=nc['temperature'][1,0,::10,::10]
imshow(t)
gca().invert_yaxis()
title('temperature');



In [4]:
s=nc['salinity'][1,0,::10,::10]
imshow(s,vmin=32.,vmax=36.)
gca().invert_yaxis()
title('salinity');



In [4]: