require "numru/ggraph"
require "narray"
include NumRu

# 縦軸の値の範囲 
min=-3.0
max=3.0
#phi0   00   02   06    08
#min  -3.0 -4.5 -8.0 -10.0
#max   3.0  4.5  8.0  10.0

gphys1 = GPhys::IO.open("./src/phi0_00.nc","u").cut("t"=>1000).cut("height"=>0)
data=gphys1.val

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

for i in 0..191 do
mu[i] = sin(Lat[i]*(Math::PI/180.0))
end

Mu = VArray.new(mu,{"long_name"=>"sin Latitude","units"=>""},"sinLat")
data = VArray.new(data,{"long_name"=>"zonal wind","units"=>"m/s"},"u")
line = VArray.new(line,{"long_name"=>"zonal wind","units"=>"m/s"},"u")
axis_mu = Axis.new.set_pos(Mu)
gphys1 = GPhys.new( Grid.new(axis_mu),data)
gphys2 = GPhys.new( Grid.new(axis_mu),line)

DCL.gropn(1)

DCL.uzfact(0.6)
GGraph.set_axes("xlabelint"=>0.2,"xtitle"=>"sin(Lat)")
GGraph.set_fig("itr"=>1,"viewport"=>[0.1,0.5,0.2,0.6])

rmiss = DCL.glpget('rmiss')

GGraph.line(gphys1,true,"exchange"=>false,"title"=>"","annotate"=>false,"index"=>20,"min"=>min,"max"=>max)
GGraph.line(gphys2,false,"exchange"=>false,"index"=>10)

DCL.grcls