# -*- coding: euc-jp -*-
require "numru/ggraph"
include NumRu

timestep  = ARGV[0].to_i

x1 = 0
x2 = 256000
#x2 = 500 * 1000 - 500 * 10
z1 = 0
z2 = 30000
t2= ARGV[0].to_i

IntegPeriod  = 4320000.0e0
tn = 25
dt = IntegPeriod / tn

files = 'thermal-moist_H2O-g_Cond.nc'
gphys  = GPhys::IO.open(files, 'H2O-g_Cond')


#DCL.gropn(1)
DCL.gropn(2)
#DCL.sldiv('y',1,2)
DCL.sgpset('lfull',true)
DCL.sgpset('lcntl', false)
DCL.uzfact(0.7)
#DCL.uzfact(0.3)
#GGraph.set_fig( 'viewport'=>[0.25,0.7,0.15,0.6] )
# アスペクト比をそれっぽくする
GGraph.set_fig( 'viewport'=>[0.1,0.85,0.15,0.6] ) 
#GGraph.set_fig( 'viewport'=>[0.25,0.8,0.10,0.275] )

gphysdisp0 = gphys
gphysdisp0 = gphysdisp0.mean( 'x' )
gphysdisp0 = gphysdisp0.mean( 'y' )
gphysdisp0 = gphysdisp0.cut( 't'=>0..t2 )
#gphysdisp0 = gphysdisp0.cut( 't'=>dt*(timestep-1)..dt*timestep )
gphysdisp  = gphysdisp0

#GGraph.contour( gphysdisp, true, 'transpose'=>true )
#GGraph.tone( gphysdisp, true, 'transpose'=>true )
GGraph.tone( gphysdisp, true, 'transpose'=>true, 
#            'lev'=>[0,0.5e-3,1.0e-3,1.5e-3,2.0e-3,2.5e-3,3.0e-3,3.5e-3,4.0e-3,4.5e-3,5.0e-3,5.5e-3,6.0e-3,6.5e-3,7.0e-3],
#             # レベル＆パターンを陽に指定
#             'pat'=>[10999,12999,15999,20999,25999,30999,35999,40999,55999,65999,70999,75999,80999,90999,95999]
#             'lev'=>[-2.5e-8,-2.0e-8,-1.5e-8,-1.0e-8,-0.5e-8,-1.0e-10,0,1.0e-10,0.5e-8,1.0e-8,1.5e-8,2.0e-8,2.5e-8 ],
             'lev'=>[   -3.0e-8,-2.5e-8,-2.0e-8,-1.5e-8,-1.0e-8,-0.5e-8,-0.25e-8,-5.0e-10,5.0e-10,0.5e-8,1.0e-8,2.0e-8],
             'pat'=>[10999,  15999,   20999,  30999,  40999, 50999,  65999,    70999,    1,    75999, 80999, 85999, 90999]
#             'pat'=>[10999,20999,   25999,  30999,35999,   40999,45999,1,    50999,67999,70999,    80999,85999,90999]
)
GGraph.color_bar
#GGraph.color_bar( 'vlength'=>0.07, 'inffact'=>0.5 )
#GGraph.color_bar( 'vlength'=>0.12, 'inffact'=>1.0 )
DCL.grcls



### めも ###
###---------------------------------
###
### ゼロしろ
##
#            'pat'=>[1,10999,15999,20999,      30999,35999,40999,            55999,65999,70999,75999,80999,90999,95999] # 14 色
#            'pat'=>[1,10999,15999,20999,      30999,      40999,      50999,      60999,70999,75999,80999,90999]       # 12 色
# (降水量, 蒸発量)
#            'pat'=>[1,10999,      20999,      30999,      40999,            55999,65999,70999,75999,80999,90999]       # 11 色
###
### non ゼロしろ
##
# (エクスナー関数)
#            'pat'=>[  10999,15999,20999,25999,30999,35999,40999,45999,50999,      67999,70999,80999,85999,90999]
# (水平風, 鉛直風まま)
#            'pat'=>[  10999,15999,20999,25999,30999,      40999,            55999,      70999,75999,80999,90999,95999] # 12 色
# (地表面熱フラックス)
#            'pat'=>[  10999,15999,20999,25999,30999,      40999,            55999,65999,70999,      80999,90999]       # 11 色
#            'pat'=>[  10999,15999,20999,      30999,      40999,            55999,65999,70999,75999,80999,90999]       # 11 色
#            'pat'=>[  10999,      20999,      30999,      40999,            55999,65999,70999,75999,80999,90999]       # 10 色
###
### 試作中
##
# (水蒸気混合比まま)
#            'pat'=>[  10999,15999,20999,      30999,      40999,            55999,65999,70999,75999,80999,90999]       # 11 色
###
###
# (水蒸気混合比たまに)
#            'pat'=>[  10999,15999,20999,      30999,      40999,                  65999,70999,75999,80999,90999]       # 10 色
#            'pat'=>[  10999,15999,20999,      30999,      40999,                  65999,70999,75999,80999,90999]       # 9 色
###
###
### でふぉ
##
# (温度)
#            'pat'=>[  10999,12999,15999,20999,25999,30999,35999,40999,55999,65999,70999,75999,80999,90999,95999] # 15 色