# -*- coding: euc-jp -*- require "numru/ggraph" include NumRu # added (07/19) timestep = ARGV[0].to_i # 引数で渡そうと思ったがうまくいかない... IntegPeriod = 4320000.0e0 tn = 20 #tn = 50 dt = IntegPeriod / tn #IntegPeriod = 432000.0e0 #tn = 2 #dt = IntegPeriod / tn x1 = 0 x2 = 512000 #x2 = 500 * 1000 - 500 * 10 z1 = 0 z2 = 30000 files = 'thermal-moist_RevHumid.nc' gphys = GPhys::IO.open(files, 'RevHumid') #DCL.gropn(1) DCL.gropn(2) #DCL.sldiv('y',1,3) 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] ) gphysdisp0 = gphys gphysdisp0 = gphysdisp0.mean( 'y' ) gphysdisp0 = gphysdisp0.mean( 'x' ) #gphysdisp0 = gphysdisp0.cut( 'z'=>450 ) gphysdisp = gphysdisp0 gphysdisp.long_name = 'Rerative Humidity' zax = gphysdisp.coord(0) # <- gphysdisp は mean や cut して 2 次元になった tax = gphysdisp.coord(1) # <- gphysdisp は mean や cut して 2 次元になった gphysdisp = gphysdisp0 #GGraph.fig( xax, tax, 'new_frame'=>true ) GGraph.fig( tax, zax, 'new_frame'=>true ) #GGraph.contour( gphysdisp, true ) #GGraph.tone( gphysdisp, true ) #GGraph.tone( gphysdisp, false ) #GGraph.tone( gphysdisp, true, 'exchange'=>true, GGraph.tone( gphysdisp, false, 'exchange'=>true, # # レベル＆パターンを陽に指定 # 相対湿度 全時刻統一 ver. 'lev'=>[ 0, 0.1e-2,0.10,0.20, 0.30, 0.40, 0.50, 0.60, 0.70, 0.80, 0.90, 1.0 ], '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 色 ) #GGraph.axes( xax, tax ) GGraph.axes( tax, zax ) GGraph.title( gphysdisp.long_name ) GGraph.color_bar 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 色 ###---------------------------------