require "numru/ggraph"
include NumRu

DCL.gropn(2)
#DCL.gropn(4)

gphys1 = GPhys::IO.open('../../netCDF/thermal-moist_H2O-g.nc', 'H2O-g')

###--- 引き数を置き変える ---###
tn = ARGV[0].to_i
dt = ARGV[1].to_i

####--- 定数の設定確認 ---###
print tn
printf ("\n")
print dt*tn
printf ("\n")

###--- Fix the tone level ---###
# 下限を -0.0015, 上限を 0.0012 とした
#  GGraph.set_tone_levels( 
#   'levels'=>[-0.0015,-0.0012,-0.0009,-0.0006,-0.0004,-0.0002,-0.0001,0,0.0001,0.0002,0.0004,0.0006,0.0009,0.0012],
#   'patterns'=>[13999,15999,20999,27999,39999,50999,69999,73999,76999,79999,86999,92999,98999,98999] )
# 下限を 0.0, 上限を 0.0012 とした
  GGraph.set_tone_levels( 
   'levels'=>[0,0.0001,0.0002,0.0003,0.0004,0.0005,0.0006,0.0007,0.0008,0.0009,0.001,0.0011,0.0012,0.0013],
   'patterns'=>[13999,15999,20999,27999,39999,50999,69999,73999,76999,79999,86999,92999,98999,98999] )
#GGraph.set_tone_levels( 
#  'levels'=>[-0.007,-0.006,-0.005,-0.004,-0.003,-0.002,-0.001,-0.0005,0,0.0005,0.001,0.002],
#  'patterns'=>[13999,15999,20999,27999,39999,50999,69999,73999,76999,79999,86999,92999] )

###--- Draw with tone ---###
  DCL.sgpset('lcntl', false)
#  GGraph.tone( ( ( gphys1 ).cut(true, true, tn*dt) + gphys2 )  -  ( ( gphys1 ).cut(true, true, tn*dt).mean( 'x' ) + ( gphys2 ).mean( 'x' ) ) )
#  GGraph.tone( ( gphys1 ).cut(true, true, tn*dt) + gphys2 ) ) 
  GGraph.tone( ( gphys1 ).cut(true, y=100, true, tn*dt) )
#  GGraph.tone( ( gphys1 ).cut(true, y=100, true, 100) )

###--- Draw side tone var ---### 
GGraph.color_bar("landscape"=>false,'labelintv'=>3)

DCL.grcls