SciTools/iris — pull request #422 by pelson (2013-03-26 13:46:26+00:00)
Title: Working OpenDAP functionality.

With this change, the following works:
```
import numpy
import matplotlib.pyplot as plt

import iris
import iris.quickplot as qplt

bathy = iris.load_cube(['http://geoport.whoi.edu/thredds/dodsC/bathy/gom15'])

# from http://colorbrewer2.org/index.php?type=sequential&scheme=Greens&n=9
earth_colors = [(247, 252, 245),
(229, 245, 224), (199, 233, 192), (161, 217, 155),
(116, 196, 118), (65, 171, 93), (35, 139, 69),
(0, 109, 44), (0, 68, 27)]

# from http://colorbrewer2.org/index.php?type=sequential&scheme=Blues&n=7
sea_colors = [(239, 243, 255), (198, 219, 239), (158, 202, 225), (107, 174, 214), (66, 146, 198), (33, 113, 181), (8, 69, 148)]

colors = numpy.array(sea_colors[:1:-1] + earth_colors[2:], dtype=numpy.float32)
colors /= 256

levels = [-4000, -2500, -400, -145, -10, 0, 10, 145, 400, 800, 1200, 1600, 2000]
qplt.contourf(bathy, levels, colors=colors, extend='both')

plt.show()
```

![topo_dap](https://f.cloud.github.com/assets/810663/303068/58d80746-961b-11e2-85d9-d7f3cc3fa16a.png)