[WIP] painfully trying to plot isofill with gengrid

[doutriaux1]

@danlipsa @sankhesh @aashish24 I need some help here.

I'm trying to get isofill (then isoline) to plot a generic grid w/o vertices info. In "theory" it shouldbe doable since we use only pointdata for isofill.

So far I create the unstructured grid and it appears I set the data on it, but the code complains there's no input data.

Here's my script and I'm attaching the grids as well

from __future__ import print_function
import vcs
import cdms2
import os
import numpy
import vtk

data = cdms2.open("/Users/doutriaux1/PeterC/ne30_TS.nc")("TS")
gridFile = cdms2.open("/Users/doutriaux1/PeterC/ne30np4_latlon.091226.nc")
lats = gridFile("lat").filled()
lons = gridFile("lon").filled()
ncols = len(lats)
latAxis = cdms2.auxcoord.TransientAuxAxis1D(lats)
lonAxis = cdms2.auxcoord.TransientAuxAxis1D(lons)
grid = cdms2.gengrid.TransientGenericGrid(latAxis,lonAxis) 
data.setAxis(-1,grid.getAxisList()[0])
data.setGrid(grid)
import vcs
x=vcs.init()
gm = vcs.createisofill()
x.plot(data,gm)
x.png("isof")

grid.zip