import netCDF4
url='http://testbedapps.sura.org/thredds/dodsC/inundation/uf/slosh/ike/2Dvrwoww'
nc = netCDF4.Dataset(url)
nc.variables.keys()
[u'time', u'lon', u'lat', u'z', u'eta', u'etamax', u'u', u'v']
eta=nc.variables['eta']
print eta
<type 'netCDF4.Variable'> float64 eta(u'time', u'number_of_j_points', u'number_of_i_points') units: m cell_methods: time: point coordinates: time lat lon _FillValue: -1e+20 valid_min: -2.43829203 valid_max: 5.38797045 standard_name: water_surface_height_above_reference_datum long_name: Water Surface Height Above Reference Datum VerticalDatum: urn:ogc:def:datum:epsg::5102 unlimited dimensions = (u'time',) current size = (90, 191, 242)
lon=nc.variables['lon']
print lon
<type 'netCDF4.Variable'> float64 lon(u'number_of_j_points', u'number_of_i_points') standard_name: longitude long_name: Longitude units: degrees_east valid_min: -96.444458 valid_max: -92.8651733 unlimited dimensions = () current size = (191, 242)
lat=nc.variables['lat']
print lat
<type 'netCDF4.Variable'> float64 lat(u'number_of_j_points', u'number_of_i_points') standard_name: latitude long_name: Latitude units: degrees_north valid_min: 28.1099072 valid_max: 30.2245941 unlimited dimensions = () current size = (191, 242)
figure(figsize(10,6))
itime=50
pcolor(lon[:,:],lat[:,:],eta[itime,:,:])
colorbar()
<matplotlib.colorbar.Colorbar instance at 0x4c6e368>