require "numru/ggraph"
include NumRu

vnames = ["VelX", "VelZ", "PTemp"]
minmaxs = [[9,11], [-1,1], [-1,1]]
datafn0 = "../test_a3/lee-wave_test_a3_"
datafn0 = "../flux_ver2/output_data/a3/lee-wave_test_a3_"

vname_momflux = "Mflux"
ncfn_momflux = "a3_h100m_H020km_momfluxspe_deepconv.nc"

gps = []
for ivar in 0..(vnames.size-1)
  vname = vnames[ivar]
#  ncfn = datafn0 + vname + "_rank00000\?.nc"
  ncfn = datafn0 + vname + "_zcoordinate.nc"

  url = ncfn + "@" + vname
  gp = GPhys::IO.open_gturl( url )
  gps << gp
end
na_time = gps[0].coord('t').val

vname = vname_momflux
ncfn = ncfn_momflux
url = ncfn + "@" + vname
gp_momflux = GPhys::IO.open_gturl( url )


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]

#for itime in 0..(na_time.size-1)
for itime in (na_time.size-1)-1..(na_time.size-1)

#for itime in 0..1
  for ivar in 0..(vnames.size-1)
    gp = gps[ivar]
    gp = gps[ivar].cut('t'=>na_time[itime])

    GGraph.set_fig 'itr'=> 1, 'viewport'=>[0.15,0.9,0.2,0.65]
#    DCL.sgscmn(4)  # blue-cyan-white-yellow-red

    print minmaxs[ivar], "\n"
    p minmaxs[ivar][0]
    p minmaxs[ivar][1]
    print minmaxs[ivar][0], ', ', minmaxs[ivar][1]#, ', ', ((minmaxs[0][ivar]+minmaxs[1][ivar])/2), "\n"
    GGraph.tone    gp, true , 'min'=>minmaxs[ivar][0], 'max'=>minmaxs[ivar][1], 'nlev'=>40
    #, 'lev'=>levels, 'pat'=>patterns
    GGraph.contour gp, false, 'lev'=>[(minmaxs[ivar][0]+minmaxs[ivar][1])/2] #'min'=>minmaxs[ivar][0], 'max'=>minmaxs[ivar][1], 'nlev'=>40
    #, 'lev'=>levels
    GGraph.color_bar( 'landscape'=>true, 'tickintv'=>0 )
  end
#  DCL.grfrm

  gp = gp_momflux
  GGraph.set_fig 'itr'=> 1, 'viewport'=>[0.15,0.9,0.2,0.65]
#  DCL.sgscmn(1)
  GGraph.tone    gp_momflux, true #, 'transpose'=>true, 'min'=>-6e-3, 'max'=>0.5e-4#, 'lev'=>levels, 'pat'=>patterns
  GGraph.contour gp_momflux, false#, 'transpose'=>true, 'min'=>-6e-3, 'max'=>0.5e-4#, 'lev'=>[0] #, 'lev'=>levels
  GGraph.color_bar( 'landscape'=>true, 'tickintv'=>0 )

end

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

DCL.grcls