require "numru/ggraph"
include NumRu

gphys1 = GPhys::IO.open("./src/nu_010.nc","u").cut("lat"=>-0..90).cut("t"=>1000).cut("height"=>8)
gphys2 = GPhys::IO.open("./src/nu_050.nc","u").cut("lat"=>-0..90).cut("t"=>1000).cut("height"=>8)
gphys3 = GPhys::IO.open("./src/nu_250.nc","u").cut("lat"=>-0..90).cut("t"=>1000).cut("height"=>8)
gphys4 = GPhys::IO.open("./src/theory.nc","u_H_theo").cut("lat"=>0..90)

DCL.gropn(1)

DCL.uzfact(0.6)

GGraph.set_fig("itr"=>1,"viewport"=>[0.1,0.7,0.2,0.6])
GGraph.set_axes("ytitle"=>"zonal velocity (m/s)","yunits"=>"","xtitle"=>"Latitude (deg)","xunits"=>"")

rmiss = DCL.glpget('rmiss')

GGraph.line(gphys1,true,"exchange"=>false,"annotate"=>false,"title"=>"","index"=>20,"min"=>0,"max"=>90)
GGraph.line(gphys2,false,"exchange"=>false,"annotate"=>false,"index"=>40)
GGraph.line(gphys3,false,"exchange"=>false,"annotate"=>false,"index"=>30)
GGraph.line(gphys4,false,"exchange"=>false,"annotate"=>false,"index"=>12)

DCL.grcls

