Must subset the mesh.  Closes #24.

2400x2400 seems to be too much for pcolormesh, subsetting to 1200x1200
seems to work, though.
Merge branch 'issue24' into devel