# coding: utf-8
require "numru/ggraph"
include NumRu

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
na_z    = gps[0].coord('z').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)/2-1..(na_time.size-1)/2

  gp = gp_momflux
  GGraph.set_fig 'itr'=> 1, 'viewport'=>[0.1,0.9,0.15,0.65]
  GGraph.tone    gp, true , 'transpose'=>true, 'min'=>-6e-3, 'max'=>0.5e-4, 'nlev'=>50 #, 'lev'=>levels, 'pat'=>patterns
#  GGraph.contour gp, false, 'transpose'=>true, 'min'=>-6e-3, 'max'=>0.5e-4#, 'lev'=>[0] #, 'lev'=>levels

  a_x = [na_time[itime],na_time[itime]]
  a_y = [na_z[0],na_z[-1]]
  DCL::uulin(a_x, a_y)  # 線を引く

  GGraph.color_bar( 'landscape'=>true, 'tickintv'=>0 )

end

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

DCL.grcls