require "numru/ggraph"
include NumRu

time=510

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

gphys = GPhys::IO.open("./src/c_5d-1/msf_mu.nc","msf").cut("t"=>time)
max=gphys.max("mu","height")
min=gphys.min("mu","height")

DCL.gropn(1)
DCL.uzfact(0.6)
GGraph.set_axes("xlabelint"=>0.1,"xunits"=>"","yunits"=>"","ytitle"=>"z/H","xtitle"=>"sin(Lat)")
GGraph.set_fig("itr"=>1,"viewport"=>[0.1,0.7,0.2,0.7])
rmiss = DCL.glpget('rmiss')

GGraph.tone(gphys,true,"exchange"=>false,"title"=>"","annotate"=>false,"levels"=>[rmiss,min,min*0.9,min*0.8,min*0.7,min*0.6,min*0.5,min*0.4,min*0.3,min*0.2,min*0.1,0.0,max*0.2,max*0.4,max*0.6,max*0.8,max*0.99,max,rmiss])
GGraph.color_bar
GGraph.contour(gphys,false,"exchange"=>false,"title"=>"","annotate"=>false,"levels"=>[rmiss,min,min*0.9,min*0.8,min*0.7,min*0.6,min*0.5,min*0.4,min*0.3,min*0.2,min*0.1,0.0,max*0.2,max*0.4,max*0.6,max*0.8,max*0.99,max,rmiss])

DCL.grcls