import vcs
import cdms2
import sys
import os
import MV2

f=cdms2.open(os.path.join(vcs.prefix,"sample_data","clt.nc"))
s=f("clt",time=slice(0,1),squeeze=1)
s=MV2.masked_less(s,65.)
x=vcs.init()
gm=x.createisofill()
#gm=x.createboxfill()
gm.missing = 252
#gm.levels = [65,70,75,80,85,90,95,100]
x.plot(s,gm)
x.png(gm.g_name)
x.interact()



