f = cdms2.open(sys.prefix+"/sample_data/clt.nc")
# Get the variable
#s = f("TREFHT",latitude=(90,0))
s=f("clt",latitude=(90,30))
# Initialize a canvas
x = vcs.init()
# Create and initilize an isofill
iso = x.createisofill()
# Draw the extension arrows on either side of the color legend
iso.ext_1 = "y"
iso.ext_2 = "y"
iso.datawc_y1=90
iso.datawc_y2=30
# Set the fill colors to the ones we retrieved earlier
p=x.createprojection()
p.type=-3
iso.projection=p
# Plot the cfsr variable from lat 25 to lat 75.6, and from lon 25 to lon 85 using the isofill defined above
x.plot(s, iso)
currently this produces
where it should really produce
