require "numru/ggraph"
include NumRu

time_1=500
time_2=750
time_3=1000
z=4
min=10**-(16)
max=10**2

gphys1 = GPhys::IO.open("./src/test_hordiff.nc","engspm").cut("t"=>time_1).cut("height"=>z)
gphys2 = GPhys::IO.open("./src/test_hordiff.nc","engspm").cut("t"=>time_2).cut("height"=>z)
gphys3 = GPhys::IO.open("./src/test_hordiff.nc","engspm").cut("t"=>time_3).cut("height"=>z)

DCL.gropn(1)
DCL.uzfact(0.6)
GGraph.set_axes("xunits"=>"","yunits"=>"","ytitle"=>"KE spectrum")
GGraph.set_fig("itr"=>4,"viewport"=>[0.1,0.5,0.2,0.6])
rmiss = DCL.glpget('rmiss')

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

DCL.grcls