require "numru/ggraph"
include NumRu

PI = acos(-1.0)
a=6.4*(10**6)
H=15000
Ro=0.0283
Omega=(2.0*PI)/(8.64*(10**4))
U=2.0*Omega*a*Ro

gphys = GPhys::IO.open("./src/steady/msf_ave.nc","msf")
gphys=gphys/(H*U)

DCL.gropn(1)
DCL.uzfact(0.6)
GGraph.set_axes("xlabelint"=>0.1,"xunits"=>"","yunits"=>"","ytitle"=>"z/H","xtitle"=>"sin(Lat)")
GGraph.set_fig("itr"=>1,"viewport"=>[0.1,0.7,0.2,0.7])
rmiss = DCL.glpget('rmiss')

GGraph.tone(gphys,true,"exchange"=>false,"title"=>"","annotate"=>false)
GGraph.color_bar
GGraph.contour(gphys,false,"exchange"=>false,"title"=>"")

DCL.grcls