Since this is gridded data, using CSV as output and then using interpolate is a lot of extra work. ERDDAP supports many gridded file format output as well, and since you want to access this data later, NetCDF is a great option. You can do this by just using the extension ‘.nc’ on the ERDDAP request. Here’s an example IPython Notebook: http://nbviewer.ipython.org/gist/rsignell-usgs/a4f418c3f588ca32e116