In [1]: import vcs
x=
In [2]: x=vcs.init()
In [3]: f=cdms2.open("/Users/doutriaux1/uvcmetrics-test-data/obs/NCEP_09_climo.nc")
In [4]: a=f("Z3")
In [5]: p=x.getprojection("lambert")
In [6]: b=x.createboxfill()
In [7]: b.projection=p
In [8]: x.plot(a(latitude=(20,60),longitude=(-160,-120)),b)