require "numru/ggraph"
include NumRu

time_start=360
time_end=720

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

mu_max = GPhys::IO.open("./src/c_5d-4/mu_minmax.nc","mu_max").cut("t"=>time_start..time_end)
mu_min = GPhys::IO.open("./src/c_5d-4/mu_minmax.nc","mu_min").cut("t"=>time_start..time_end)

data1 = mu_max.val
data2 = mu_min.val
data3 = NArray.float(361).fill(0.0) # 加熱中心のサイン緯度 \mu_0
data4 = NArray.float(361).fill(0.0)  # 縦軸=0の線

time = mu_max.axis('t').pos.val
time = time-360
t = VArray.new(time,{"long_name"=>"day","units"=>"day"},"day")
axis_t = Axis.new.set_pos(t)

for i in 0..360 do
data3[i] = 0.2*sin((2.0*PI/360.0)*time[i])
end

data1 = VArray.new(data1,{"long_name"=>"position","units"=>"1"},"mu")
data2 = VArray.new(data2,{"long_name"=>"position","units"=>"1"},"mu")
data3 = VArray.new(data3,{"long_name"=>"position","units"=>"1"},"mu")
data4 = VArray.new(data4,{"long_name"=>"position","units"=>"1"},"mu")

gphys1 = GPhys.new(Grid.new(axis_t),data1)
gphys2 = GPhys.new(Grid.new(axis_t),data2)
gphys3 = GPhys.new(Grid.new(axis_t),data3)
gphys4 = GPhys.new(Grid.new(axis_t),data4)

DCL.gropn(1)
DCL.uzfact(0.6)
GGraph.set_fig("itr"=>1,"viewport"=>[0.1,0.6,0.2,0.6])
GGraph.set_axes("xtitle"=>"day","xunits"=>"","yunits"=>"","ytitle"=>"Position of Psi_(Max,Min)")
rmiss = DCL.glpget('rmiss')

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

DCL.grcls