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)