require "numru/ggraph"
include NumRu

time=1000

gphys = GPhys::IO.open("./src/nu_005.nc","msf").cut("t"=>800..1000).cut("lat"=>0..60)
gphys = gphys.mean("t")
#gphys = GPhys::IO.open("./src/nu_010.nc","msf").cut("t"=>time).cut("lat"=>0..60)

max=gphys.max("lat","height")

DCL.gropn(1)

DCL.uzfact(0.6)
GGraph.set_axes("xlabelint"=>5,"ylabelint"=>1)
GGraph.set_fig("itr"=>1,"viewport"=>[0.1,0.7,0.2,0.6])
rmiss = DCL.glpget('rmiss')

GGraph.tone(gphys,true,"exchange"=>false,"title"=>"","annotate"=>false,"levels"=>[rmiss,max*0,max*0.1,max*0.2,max*0.3,max*0.4,max*0.5,max*0.6,max*0.7,max*0.8,max*0.9,max,rmiss])
GGraph.color_bar
GGraph.contour(gphys,false,"exchange"=>false,"title"=>"","levels"=>[rmiss,max*0,max*0.1,max*0.2,max*0.3,max*0.4,max*0.5,max*0.6,max*0.7,max*0.8,max*0.9,max,rmiss])

DCL.grcls

