In [1]:

    

import netCDF4


    

In [2]:

    

url='http://testbedapps.sura.org/thredds/dodsC/inundation/uf/slosh/ike/2Dvrwoww'
nc = netCDF4.Dataset(url)
nc.variables.keys()


    

    
Out[2]:





[u'time', u'lon', u'lat', u'z', u'eta', u'etamax', u'u', u'v']

In [3]:

    

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)

In [4]:

    

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)

In [5]:

    

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)

In [8]:

    

figure(figsize(10,6))
itime=50
pcolor(lon[:,:],lat[:,:],eta[itime,:,:])
colorbar()


    

    
Out[8]:





<matplotlib.colorbar.Colorbar instance at 0x4c6e368>






    

In [ ]: