require "numru/ggraph"
include NumRu

yranges=[
  [-50.0,0.5],
  [-7.0,0.5],
  [-7.0,0.5],
  [-7.0,0.5]
]
ncfns = [
  "out/a1_h100m_momfluxspe_deepconv.nc",
  "out/a2_h100m_momfluxspe_deepconv.nc",
#  "out/a3_h100m_momfluxspe_deepconv.nc",
#  "out_additional/a3_h100m_momfluxspe_deepconv.nc",
#  "zzz/a3_h100m_momfluxspe_deepconv-\?.nc",
#  "out_additional/a3_h100m_H020km_momfluxspe_deepconv.nc",
#  "out_additional/a3_h100m_H030km-\?_momfluxspe_deepconv.nc",
  "out_additional/a3_h100m_L080km_H020km_SpongeVelX10Zt10kmTau1800s_momfluxspe_deepconv.nc",
  "out/a4_h100m_momfluxspe_deepconv.nc",
]
anal_ncfns = [
  "out/a1_h100m_momfluxspe_analytic.nc",
  "out/a2_h100m_momfluxspe_analytic.nc",
  "out/a3_h100m_momfluxspe_analytic.nc",
  "out/a4_h100m_momfluxspe_analytic.nc",
]


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

DCL.sldiv('y',2,2)
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]

timee = 18000
timee = 9000
dtime = 60
for it in 0..(timee/dtime-1)
  time = dtime*it

  vname = "MfluxSpe"
  for ifile in 0..(ncfns.size-1)
    wns=0
    wne=5e-3
    z=510
    y1=yranges[ifile][0]
    y2=yranges[ifile][1]

    GGraph.set_fig 'itr'=> 1, 'viewport'=>[0.1,0.9,0.1,0.65], 'window'=>[wns,wne,y1,y2]

    ncfn = ncfns[ifile]
    url = ncfn + "@" + vname
    gp = GPhys::IO.open_gturl( url )
    gpout = gp.cut('wn'=>wns..wne,'z'=>z,'t'=>time)
    gpout = gpout*2
    GGraph.line gpout, true, 'index'=>20

    ncfn = anal_ncfns[ifile]
    url = ncfn + "@" + vname
    gp = GPhys::IO.open_gturl( url )
    gpout = gp.cut('wn'=>wns..wne)
    GGraph.line gpout, false, 'index'=>40


#    label = 'z=' + z.to_s + ', time=' + time.to_s
#    rsize = DCL.sgqtxs()
#    DCL.sgstxs( rsize*0.25 )
#    DCL.sgtxv( 0.85, 0.2, label )
#    DCL.sgstxs( rsize )
  end
end


DCL.grcls