require "numru/ggraph" include NumRu ### ARGV[0] は引数をもってくるときつかう time = ARGV[0].to_i #time = 1200 ### 範囲指定 x1 = 10000 x2 = 30000 z1 = 0 z2 = 20000 #ymid =10000 files = "MSE.nc" gphys = GPhys::IO.open(files, 'MSE') DCL.gropn(2) #DCL.sldiv('y',1,2) # 複数枚のときにいる( 'y方向に'1x2') で分割 DCL.sgpset('lfull',true) DCL.sgpset('lcntl', false) #DCL.sgpset('lclip', true) DCL.uzfact(0.6) GGraph.set_fig( 'viewport'=>[0.3, 0.6, 0.2, 0.5] ) #xmin,xmax,ymin,ymax #-- axis じゃないとだめ. coord では to_gphys が通らない xax = gphys.axis(0) zax = gphys.axis(2) #p xax #p zax xax = xax.to_gphys.cut('x'=>x1..x2) zax = zax.to_gphys.cut('z'=>z1..z2) gphys = gphys.cut( 'x'=>x1..x2 ) gphys = gphys.cut( 'z'=>z1..z2 ) ### 時間でカット gphysdisp0 = gphys.cut( 't'=>time ) gphysdisp0 = gphysdisp0.mean( 'y' ) gphysdisp = gphysdisp0 GGraph.fig( xax, zax, 'new_frame'=>true ) ## new_frame fales の場合 'viewport'=>[0.25,0.9,vpy1,vpy2], 'window'=>[x1,x2,z1,z2] ) GGraph.tone( gphysdisp, false, # 'lev'=>[260,280,300,320,340,360,380,400,420,440,460], # 'lev'=>[-2.5,-2.0,-1.5,-1.0,-0.5,0,0.5,1.0,1.5,2.0,2.5], 'lev'=>[ 3.28e5, 3.32e5, 3.36e5, 3.4e5, 3.44e5, 3.48e5, 3.52e5, 3.56e5, 3.6e5, 3.65e5, 3.7e5, 3.8e5, 3.9e5], # レベル&パターンを陽に指定 'pat'=>[ 15999,20999,30999,35999,45999,65999,70999,75999,77999,80999,85999,90999,95999] ) GGraph.axes( xax, zax ) GGraph.title( gphysdisp.long_name ) GGraph.color_bar ( 'vlength'=>0.24, 'inffact'=>0.5 ) GGraph.annotate ( [ time.to_s + ' s' ] ) DCL.grcls