# -*- coding: euc-jp -*- require "numru/ggraph" include NumRu #timestep = ARGV[0].to_i x1 = 0 #x2 = 80000 x2 = 256000 #x2 = 512000 #x2 = 500 * 1000 - 500 * 10 z1 = 0 #z2 = 20000 z2 = 48000 #t2 = 10800 t2 = 86400 * 30 t2= ARGV[0].to_i files = "thermal-moist_VelZ.nc" gphys = GPhys::IO.open(files, 'VelZ') #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.1,0.85,0.2,0.45] ) # 前の残り #GGraph.set_fig( 'viewport'=>[0.25,0.8,0.10,0.275] ) gphysdisp0 = gphys gphysdisp0 = gphysdisp0.mean( 'x' ) gphysdisp0 = gphysdisp0.mean( 'y' ) gphysdisp = gphysdisp0.cut('t'=>0..t2) 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], ## # レベル＆パターンを陽に指定 ### デフォルト？ ## 'lev'=>[ -35.0,-30.0,-25.0,-20.0,-15.0,-10.0,-5.0, 0, 5.0, 10.0, 15.0, 20.0, 25.0, 30.0], 'lev'=>[ -60.0,-45.0,-30.0,-15.0, 0, 15.0, 30.0, 45.0, 60.0, 75.0, 90.0, 105.0], 'pat'=>[10999,20999,30999,35999,40999,45999,50999,67999,70999,80999,85999,90999] ### サーマル用 ## 'lev'=>[ -2.4e-3,-1.2e-3,0,1.2e-3,2.4e-3,3.6e-3,4.8e-3,6.0e-3,7.2e-3], # 'lev'=>[ 0, 0.02, 0.04, 0.06, 0.08, 0.10, 0.12, 0.14, 0.16, 0.18, 0.20], # 'pat'=>[10999,20999,30999,35999,40999,45999,50999,67999,70999,80999,85999,90999] # 'pat'=>[10999,20999,35999,55999,65999,70999,75999,80999,85999,90999] ## 'pat'=>[10999,12999,15999,20999,25999,30999,35999,40999,55999,65999,70999,75999,80999,90999,95999] ## 'pat'=>[10999,15999,20999,25999,30999,35999,40999,1,45999,50999,67999,70999,80999,85999,90999] ) #GGraph.color_bar( 'vlength'=>0.07, 'inffact'=>0.5 ) #GGraph.color_bar( 'vlength'=>0.12, 'inffact'=>1.0 ) GGraph.color_bar #gphysdisp = gphysdisp0.cut('t'=>dt*1..dt*2) ##GGraph.contour( 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], ## # レベル＆パターンを陽に指定 # 'lev'=>[ -35.0,-30.0,-25.0,-20.0,-15.0,-10.0,-5.0, 0, 5.0, 10.0, 15.0, 20.0, 25.0, 30.0], # 'pat'=>[10999,12999,15999,20999,25999,30999,35999,40999,55999,65999,70999,75999,80999,90999,95999] ## 'pat'=>[10999,15999,20999,25999,30999,35999,40999,1,45999,50999,67999,70999,80999,85999,90999] #) ##GGraph.color_bar( 'vlength'=>0.07, 'inffact'=>0.5 ) ##GGraph.color_bar( 'vlength'=>0.12, 'inffact'=>1.0 ) #GGraph.color_bar #gphysdisp = gphysdisp0.cut('t'=>dt*2..dt*3) #GGraph.contour( gphysdisp, true, 'transpose'=>true ) # #gphysdisp = gphysdisp0.cut('t'=>dt*3..dt*4) #GGraph.contour( gphysdisp, true, 'transpose'=>true ) ### めも ### ###--------------------------------- ### ### ゼロしろ ## # '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 色 # level がこれのとき # -> 'lev'=>[ -60.0,-45.0,-30.0,-15.0, 0, 15.0, 30.0, 45.0, 60.0, 75.0, 90.0, 105.0], # 'pat'=>[ 10999,15999,20999, 30999, 40999, 55999,65999,70999,75999,80999,85999,90999,95999] # 13 色 # (地表面熱フラックス) # '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 色 DCL.grcls