require "numru/ggraph"
include NumRu

a=6.4*(10**6)
cint=0.2   # 等値線間隔 phi0=0,2 のときは 0.1, phi0=6 のときは 0.2

gphys = GPhys::IO.open("./src/phi0_06.nc","msf").cut("t"=>1000)
gphys = gphys*a/((10**10))

DCL.gropn(1)

DCL.uzfact(0.6)
GGraph.set_fig("itr"=>1,"viewport"=>[0.1,0.7,0.2,0.6])

rmiss = DCL.glpget('rmiss')

GGraph.tone(gphys,true,"exchange"=>false,"title"=>"","annotate"=>false,"interval"=>cint)
GGraph.color_bar
GGraph.contour(gphys,false,"exchange"=>false,"interval"=>cint)
DCL.grcls