require "numru/ggraph"
include NumRu

nu=3.5
time=510

gphys = GPhys::IO.open("./src/c_5d-3/dudz.nc","dudz").cut("t"=>time)
gphys=nu*gphys

DCL.gropn(1)
DCL.uzfact(0.6)
GGraph.set_axes("xlabelint"=>0.1,"xunits"=>"","yunits"=>"","xtitle"=>"sin (Lat)","ytitle"=>"z/H")
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,-0.032,-0.024,-0.016,-0.008,0,0.008,0.016,0.024,0.032,rmiss])
GGraph.color_bar
GGraph.contour(gphys,false,"exchange"=>false,"title"=>"","levels"=>[rmiss,-0.032,-0.024,-0.016,-0.008,0,0.008,0.016,0.024,0.032,rmiss])

DCL.grcls