require "numru/ggraph"
include NumRu

gphys1 = GPhys::IO.open("./src/phi0_06.nc","u").cut("t"=>1000).cut("height"=>15)

Lat_mat = gphys1.axis('lat').pos.val
line = NArray.float(192).fill(0.0)

Lat_va = VArray.new(Lat_mat,{"long_name"=>"Latitude","units"=>"degree"},"lat")
line = VArray.new(line,{"long_name"=>"zonal wind","units"=>"m/s"},"u")
Lat_new = Axis.new.set_pos(Lat_va)
gphys5 = GPhys.new( Grid.new(Lat_new),line)

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"=>-10,"max"=>105)
GGraph.line(gphys5,false,"exchange"=>false,"annotate"=>false,"index"=>10)

DCL.grcls

