Read Chlorophyll data from ERDDAP

In [1]:

    

import urllib
import netCDF4


    

In [2]:

    

url=('http://coastwatch.pfeg.noaa.gov/erddap/griddap/erdVHchla8day.nc?'
     'chla[(2014-05-29T00:00:00Z):1:(2014-05-29T00:00:00Z)]'
     '[(47.97916666666667):1:(34.89583333333334)][(-76.5625):1:(-63.47916666666666)]')
file='chla.nc'
urllib.urlretrieve (url, file)
nc = netCDF4.Dataset(file)
ncv = nc.variables
ncv.keys()


    

    
Out[2]:





[u'time', u'latitude', u'longitude', u'chla']

In [6]:

    

figure(figsize=(10,8))
pcolormesh(ncv['longitude'][:],ncv['latitude'][:],log(ncv['chla'][0,:,:]),vmin=-2,vmax=2);


    

In [ ]: