require "numru/ggraph"
include NumRu

time=1000
max=1045.0

gphys = GPhys::IO.open("./src/test.nc","msf").cut("lat"=>0..60).cut("t"=>time)

DCL.gropn(1)
DCL.uzfact(0.6)
GGraph.set_axes("xlabelint"=>5,"ylabelint"=>1,"xtitle"=>"latitude (deg)","ytitle"=>"height (km)","xunits"=>"","yunits"=>"")
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