require "numru/ggraph"
include NumRu

ncfn = "out/a3_h100m_x200km_momfluxspe_deepconv.nc"
ncfn = "zzz/momfluxspe_deepconv-\?.nc"
ncfn = "out_additional/a3_h100m_L080km_H020km_SpongeVelX10Zt10kmTau1800s_momfluxspe_deepconv.nc"
vname = "MfluxSpe"
url = ncfn + "@" + vname
gp = GPhys::IO.open_gturl( url )

na_time = gp.coord('t').val

wns=0
wne=5e-3
zs=0
ze=5e3
ts=0
te=na_time[-1]


iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
DCL.gropn(iws)

#DCL.sldiv('y',1,4)
DCL.sgpset('isub', 96)     # control character of subscription: '_' --> '`'
DCL.sgpset('lfull',true)
DCL.uzfact(0.6)            # font size
DCL.sgpset('lclip', true)  # clip graph
DCL.glpset('lmiss',true)   # missing value

#DCL.sgscmn(4)  # blue-cyan-white-yellow-red
#DCL.sgscmn(5)  # gray-scale
#DCL.sgscmn(14) # blue-white-red

#levels = [110,120,130,140,150,160,170,180,190,200,210,220,230,240]
#patterns = [10999,15999,20999,25999,30999,35999,40999,50999,60999,65999,70999,75999,80999,90999,95999]

for it in 0..(na_time.size-1)
  gpout = gp.cut('wn'=>wns..wne,'z'=>zs..ze,'t'=>na_time[it])
  GGraph.set_fig 'itr'=> 1, 'viewport'=>[0.1,0.9,0.1,0.65], 'window'=>[wns,wne,zs,ze]
  GGraph.tone gpout, true    , 'min'=>-4, 'max'=>1.5, 'nlev'=>50 #, 'lev'=>levels, 'pat'=>patterns
  GGraph.contour gpout, false, 'lev'=>[0]
  GGraph.color_bar
end



#GGraph.line gphys, true, 'exchange'=>false

DCL.grcls