require "numru/ggraph"
include NumRu

ivar_target = 2
vnames = ["VelX", "VelZ", "PTemp"]
minmaxs = [[9,11], [-1,1], [-1,1]]
datafn0 = "../test_a3_H020km/lee-wave_test_a3_"
datafn0 = "../flux_ver2/output_data/a3_H020km/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 (na_time.size-1)-1..(na_time.size-1)
for itime in 0..(na_time.size-1)
  print itime, "\n"

  ivar = ivar_target

  gp = gps[ivar]
  gp = gps[ivar].cut('t'=>na_time[itime])

  GGraph.set_fig 'itr'=> 1, 'viewport'=>[0.1,0.9,0.15,0.65]

  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.grcls