require "numru/ggraph"
include NumRu

ncfn = "zzz/momfluxspe_deepconv-\?.nc"
ncfn = "out/a3_h100m_momfluxspe_deepconv.nc"
ncfn = "out_additional/a3_h100m_momfluxspe_deepconv.nc"
ncfn = "out_additional/a3_h100m_L080km_momfluxspe_deepconv.nc"
ncfn = "out_additional/a3_h100m_L160km_momfluxspe_deepconv.nc"
#ncfn = "out_additional/a3_h100m_H020km_momfluxspe_deepconv.nc"
#ncfn = "zzz/a3_h100m_momfluxspe_deepconv-\?.nc"
#ncfn = "momfluxspe_deepconv.nc"

ncfn_anal = "out/a3_h100m_momfluxspe_analytic.nc"

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]

vname = "Mflux"
url = ncfn + "@" + vname
gp = GPhys::IO.open_gturl( url )

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

zs=0
ze=5e3
ze=na_z[-1]
ts=0
te=na_time[-1]
#te=1200


GGraph.set_fig 'itr'=> 1, 'viewport'=>[0.1,0.85,0.1,0.65], 'window'=>[-55e-4,0,zs,ze]
flag_first = true
for i in 1..5
  time = ((te-ts)/5)*i
  gpout = gp.cut('z'=>zs..ze,'t'=>time)
  GGraph.line gpout, flag_first, 'index'=>i*10, 'exchange'=>true
  flag_first = false
end


gpout = gp.cut('z'=>zs..ze,'t'=>ts..te)
GGraph.set_fig 'itr'=> 1, 'viewport'=>[0.1,0.85,0.1,0.65], 'window'=>[ts,te,zs,ze]
GGraph.tone gpout, true    , 'transpose'=>true, 'nlev'=>50, 'min'=>-6e-3, 'max'=>0.5e-3 #, 'lev'=>levels, 'pat'=>patterns
GGraph.contour gpout, false, 'transpose'=>true, 'lev'=>[0]#, 'min'=>-2.25, 'max'=>0.75 #, 'lev'=>levels
GGraph.color_bar


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]
#te=1200

gpout = gp.cut('wn'=>wns..wne,'z'=>zs..ze,'t'=>te)
gpout = gpout*2
GGraph.set_fig 'itr'=> 1, 'viewport'=>[0.1,0.85,0.1,0.65], 'window'=>[wns,wne,zs,ze]
GGraph.tone gpout, true    , 'min'=>-4.5, 'max'=>1.5, 'nlev'=>50 #, 'lev'=>levels, 'pat'=>patterns
GGraph.contour gpout, false, 'lev'=>[0]
GGraph.color_bar

gpout = gp.cut('wn'=>wns..wne,'z'=>510,'t'=>ts..te)
gpout = gpout*2
GGraph.set_fig 'itr'=> 1, 'viewport'=>[0.1,0.85,0.1,0.65], 'window'=>[wns,wne,ts,te]
GGraph.tone gpout, true    , 'min'=>-4.5, 'max'=>1.5, 'nlev'=>50 #, 'lev'=>levels, 'pat'=>patterns
GGraph.contour gpout, false, 'lev'=>[0]
GGraph.color_bar

gpout = gp.cut('wn'=>wns..wne,'z'=>510,'t'=>te)
gpout = gpout*2
GGraph.set_fig 'itr'=> 1, 'viewport'=>[0.1,0.85,0.1,0.65], 'window'=>[wns,wne,-5,1]
GGraph.line gpout, true
ncfn = ncfn_anal
vname = "MfluxSpe"
url = ncfn + "@" + vname
gp = GPhys::IO.open_gturl( url )
GGraph.line gp, false, 'index'=>20



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

DCL.grcls