#!/usr/bin/env ruby =begin = mkfig.rb -- make figure for zonal mean for 25 means relative:25MEANS.1979-2003.NCEP/ 以下の netCDF ファイルから 25 年平均値を取り出しプロット. 現在 mkfig_25years.rb に機能をマージしたため, 未使用. 100hPa*.nc: shum 300hPa*.nc: strm uwnd t dse dse_bar dse_dash lhe lhe_bar lhe_dash GAUSSIAN*.nc: nswrs nlwrs nswrf ulsrf lhtfl shtfl lheat 各図の物理量は (a)オイラー平均質量流線関数, (b)東西風, (c)温度, (d)比湿, (e)放射収支(黒線:大気上端正味短波フラックス, 赤線:大気上端上向き長波フラックス, 緑線:地表面正味短波フラックス, 青線:地表面正味長波フラックス), (f)熱収支(黒線:地表面正味潜熱フラックス, 赤線:地表面正味顕熱フラックス, 緑線:降水), (g)南北乾燥静的エネルギー輸送量 (黒線:合計, 赤線:帯状平均成分, 緑線:帯状擾乱成分), (h)南北潜熱輸送量 (黒線:合計, 赤線:帯状平均成分, 緑線:帯状擾乱成分), (i)EP フラックスとその発散 (矢羽: EP フラックス, トーン:発散), (j)残差循環の質量流線関数 を一枚にまとめた図. =end require "getopts" # for option_parse require "numru/ggraph" require "libgphys-n" include NumRu include Misc::EMath ## enable to the option min, max module NumRu module GGraph module_function def vector_for_dcchart(fx, fy, newframe=true, options=nil) if ! defined?(@@vector_for_dcchart_options) @@vector_for_dcchart_options = Misc::KeywordOptAutoHelp.new( ['newfig', true, 'if false, do not cleared before figure setting.'], ['title', nil, 'Title of the figure(if nil, internally determined)'], ['annotate', true, 'if false, do not put texts on the right margin even when newframe==true'], ['transpose', false, 'if true, exchange x and y axes'], ['xintv', 1, '(Effective only if flow_vect) interval of data sampling in x'], ['yintv', 1, '(Effective only if flow_vect) interval of data sampling in y'], ['xfact', 1, 'scale factor for x'], ['yfact', 1, 'scale factor for y'], ['itr', 1, 'log or linear'], ['unit_vect', false, 'Show the unit vector'], ['max_unit_vect', false, '(Effective only if flow_vect && unit_vect) If true, use the maximum arrows to scale the unit vector; otherwise, normalize in V coordinate.'] ) end opts = @@vector_for_dcchart_options.interpret(options) fx = fx.first2D.copy fy = fy.first2D.copy sh = fx.shape if sh != fy.shape raise ArgumentError, "shapes of fx and fy do not agree with each other" end fx = fx.transpose(1,0) if opts['transpose'] fy = fy.transpose(1,0) if opts['transpose'] if ((xi=opts['xintv']) >= 2) idx = NArray.int(sh[0]/xi).indgen!*xi # [0,xi,2*xi,..] fx = fx[idx, true] fy = fy[idx, true] end if ((yi=opts['yintv']) >= 2) idx = NArray.int(sh[1]/yi).indgen!*yi # [0,yi,2*yi,..] fx = fx[true, idx] fy = fy[true, idx] end xax = fx.coord(0) yax = fy.coord(1) if newframe nextfig = @@next_fig.dup if ( @@next_fig != nil ) # backup next_fig fig(xax, yax, {'itr'=>itr}) axes(xax, yax) if opts['title'] ttl = opts['title'] else fxnm = fx.data.get_att('long_name') || fx.name fynm = fy.data.get_att('long_name') || fy.name ttl = '('+fxnm+','+fynm+')' end title( ttl ) annotate(fx.lost_axes) if opts['annotate'] @@next_fig = nextfig if ( @@next_fig != nil ) end fig(xax, yax, {'new_frame'=>false, 'itr'=>1, 'yreverse'=>false}) \ if (opts['newfig']) DCL.uwsgxa(xax.val) DCL.uwsgya(yax.val) xfact = opts["xfact"]; yfact = opts["yfact"] before2=DCLExt.ug_set_params( {'LNRMAL'=>false, 'LMSG'=>false, 'XFACT1'=>xfact, 'YFACT1'=>yfact} ) # before=DCLExt.ug_set_params({'lunit'=>true}) if opts['unit_vect'] DCL.ugvect(fx.val[true, -1..0], fy.val[true, -1..0]) DCLExt.unit_vect(xfact, yfact, nil, nil, "vxuoff"=>0.01, "vyuoff"=>-0.1) if opts['unit_vect'] # DCLExt.ug_set_params(before) if opts['unit_vect'] nil end end end def GGraph::line(gphys, newframe=true, options=nil) if newframe!=true && newframe!=false raise ArgumentError, "2nd arg (newframe) must be true or false" end if ! defined?(@@line_options) @@line_options = Misc::KeywordOptAutoHelp.new( ['title', nil, 'Title of the figure(if nil, internally determined)'], ['annotate', true, 'if false, do not put texts on the right margin even when newframe==true'], ['exchange', false, 'whether to exchange x and y axes'], ['index', 1, 'line/mark index'], ['type', 1, 'line type'], ['label', nil, 'if a String is given, it is shown as the label'], ['max', nil, 'maximam value for draw line'], ['min', nil, 'minimam value for draw line'] ) end opts = @@line_options.interpret(options) gp = gphys.first1D if !opts['exchange'] x = gp.coord(0) y = gp.data xax = x if opts['min'] ymin = opts['min'].to_f else ymin = gp.data.val.min end if opts['max'] ymax = opts['max'].to_f else ymax = gp.data.val.max end yax = VArray.new(NArray[ymin, ymax], gp.data, gp.data.name) else y = gp.coord(0) x = gp.data yax = y if opts['min'] xmin = opts['min'] else xmin = gp.data.val.min end if opts['max'] xmax = opts['max'] if opts['max'] else xmax = gp.data.val.max end xax = VArray.new(NArray[xmin, xmax], gp.data, gp.data.name) end if newframe fig(xax, yax) axes(xax, yax) title( (opts['title'] || gp.data.get_att('long_name')) ) annotate(gp.lost_axes) if opts['annotate'] end if opts['label'] lcharbk = DCL.sgpget('lchar') DCL.sgpset('lchar',true) DCL.sgsplc(opts['label']) end DCL.uulinz(x.val, y.val, opts['type'], opts['index']) DCL.sgpset('lchar',lcharbk) if opts['label'] nil end ##----------------------- # ラインメソッド再定義? def plot_line_main(gphys_array, line_opts={}) line_index_ary = Array.new; name_ary = Array.new default_index = 12 line_hash = { "index"=> default_index }.update(line_opts) line_index_ary.push(default_index) name_ary.push(gphys_array[0].data.get_att("long_name")) GGraph.line( gphys_array[0], true, line_hash ) gphys_array.size.times{ |num| unless num == 0 line_hash = { "index"=> (num*10 + default_index)} line_index_ary.push(num*10 +default_index) name_ary.push(gphys_array[num].data.get_att("long_name")) GGraph.line( gphys_array[num], false, line_hash ) end } return name_ary, line_index_ary end ##----------------------- # ラインインデックスの種類を表示 def __show_line_index(str_ary,line_index_ary, x=0.15, y=0.1) charsize = 0.85 * DCL.uzpget('rsizec1') dvx = 0.01 dvy = charsize*1.6 raise TypeError,"Array expected" if ! str_ary.is_a?(Array) vxmin,vxmax,vymin,vymax = DCL.sgqvpt vx = x vy = y + charsize/2 str_ary.size.times{|num| DCL::sgtxzv(vx,vy,"--- #{str_ary[num]}", charsize, 0, -1, line_index_ary[num]) vy -= dvy if num == 6 vx = 0.30 vy = 0.12 - charsize/2 end } nil end ##----------------------- ## 25 年平均のオブジェクトを作成 def make_mean_gphys(phys, src, season, year, var=phys) pathary = [] y1 = year[0]; y2 = year[1] path = "../../../../#{phys}.#{src}/" season.each do |m| y1.upto(y2) do |y| y = y-1 if m == "12" fn = path + "#{phys}.#{y}.#{src}/#{phys}_#{y}-#{m}_#{src}.nc" pathary << fn end end # gp = mean_gphys_date_weight(pathary, phys.downcase) gp = mean_gphys(pathary, var.downcase) return gp end ## GGraph の margin_info を参考に, 余白の上下につけるようにした. def margin_info(program=nil, data_source=nil, char_height=nil, date=nil, xl=0.0, xr=0.0, yb=nil, yt=nil) program = File.expand_path($0) if !program if date program = program + ' ' + Date.today.to_s else program = program if date.nil? program = program + ' ' + Date.today.to_s end end if !data_source data_source = Dir.pwd data_source = '' if data_source == program.sub(/\/[^\/]+$/,'') end sz = 1.0/( program.length + data_source.length ) char_height = [0.008, sz].min if !char_height yb = 2.0 * char_height if !yb yt = 2.0 * char_height if !yt DCL.slmgn(xl, xr, yb, yt) # 上につける DCL.slsttl(program, 't', -1.0, 1.0, char_height, 1) DCL.slsttl(data_source, 't', 1.0, 1.0, char_height, 2) # 下につける DCL.slsttl(program, 'b', -1.0, 1.0, char_height, 3) DCL.slsttl(data_source, 'b', 1.0, 1.0, char_height, 4) end ################################################################## unless getopts( "", "range:", "kikan:" # "annual", "winter", "spring", "summer", "fall" を指定 ) print "#{$0}:illegal options. please exec this program with -h/--help. \n" exit 1 end ## define Constants Infty = 1.0/0.0 L = UNumeric.new(2.5008e6, 'J.Kg-1') Rad = UNumeric.new(6.37E6, 'm') G = UNumeric.new(9.81 , 'm.s-2') # make gphys objects ## set kikan winter = "DJF"; spring = "MAM" summer = "JJA"; fall = "SON" annual = "ANNUAL" case $OPT_kikan when "annual" kikan = annual when "winter" kikan = winter when "spring" kikan = spring when "summer" kikan = summer when "fall" kikan = fall else raise "#{$OPT_kikan} is invalid value." end ##### 気候値取得 ## 質量流線関数 strm_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/100hPa_PHYS_1979-2003_#{kikan}_NCEP.nc", "strm") ## 東西風速 uwnd_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/100hPa_PHYS_1979-2003_#{kikan}_NCEP.nc", "uwnd") ## 温度場 temp_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/100hPa_PHYS_1979-2003_#{kikan}_NCEP.nc", "t") ## 比湿 shum_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/300hPa_PHYS_1979-2003_#{kikan}_NCEP.nc", "shum") ## 放射 nswrf_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/GAUSSIAN_PHYS_1979-2003_#{kikan}_NCEP.nc", "uswrf") ulwrf_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/GAUSSIAN_PHYS_1979-2003_#{kikan}_NCEP.nc", "ulwrf") nlwrs_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/GAUSSIAN_PHYS_1979-2003_#{kikan}_NCEP.nc", "nlwrs") nswrs_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/GAUSSIAN_PHYS_1979-2003_#{kikan}_NCEP.nc", "nswrs") ## 熱 #### 潜熱フラックス lhtfl_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/GAUSSIAN_PHYS_1979-2003_#{kikan}_NCEP.nc", "lhtfl") #### 顕熱フラックス shtfl_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/GAUSSIAN_PHYS_1979-2003_#{kikan}_NCEP.nc", "shtfl") #### 降水(潜熱に換算したもの) lheat_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/GAUSSIAN_PHYS_1979-2003_#{kikan}_NCEP.nc", "prate") ## 乾燥静的エネルギー dse_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/100hPa_PHYS_1979-2003_#{kikan}_NCEP.nc", "dsefl") dse_bar_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/100hPa_PHYS_1979-2003_#{kikan}_NCEP.nc", "dsefl_bar") dse_dash_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/100hPa_PHYS_1979-2003_#{kikan}_NCEP.nc", "dsefl_dash") ## 潜熱輸送量 lhe_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/100hPa_PHYS_1979-2003_#{kikan}_NCEP.nc", "lhefl") lhe_bar_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/100hPa_PHYS_1979-2003_#{kikan}_NCEP.nc", "lhefl_bar") lhe_dash_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/100hPa_PHYS_1979-2003_#{kikan}_NCEP.nc", "lhefl_dash") ## EP フラックス epflx_phi_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/100hPa_PHYS_1979-2003_#{kikan}_NCEP.nc", "epflx_phi") epflx_p_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/100hPa_PHYS_1979-2003_#{kikan}_NCEP.nc", "epflx_p") epflx_div_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/100hPa_PHYS_1979-2003_#{kikan}_NCEP.nc", "epflx_div") ## 残差循環 strm_rmean_climat = GPhys::IO.open("../25MEANS.1979-2003.NCEP/100hPa_PHYS_1979-2003_#{kikan}_NCEP.nc", "strm_rmean") ##### ------------ 各年 ## set kikan year = $OPT_range.split(/:/).collect{|y| y.to_i } winter = ["12", "01", "02"]; spring = ["03", "04", "05"] summer = ["06", "07", "08"]; fall = ["09", "10", "11"] annual = winter + spring + summer + fall case $OPT_kikan when "annual" kikan = annual when "winter" kikan = winter when "spring" kikan = spring when "summer" kikan = summer when "fall" kikan = fall else raise "#{$OPT_kikan} is invalid value." end ## 質量流線関数 strm = make_mean_gphys('STRM', 'NCEP', kikan, year).cut( 'level'=>1000..100 ) ## 東西風速 uwnd = make_mean_gphys('UWND', 'NCEP', kikan, year).mean('lon').cut( 'level'=>1000..100 ) ## 温度場 temp = make_mean_gphys('T', 'NCEP', kikan, year).mean('lon').cut( 'level'=>1000..100 ) ## 比湿 shum = make_mean_gphys('SHUM', 'NCEP', kikan, year).mean('lon').cut( 'level'=>1000..100 ) ## 放射 dswrf = make_mean_gphys('DSWRF', 'NCEP', kikan, year).mean('lon') uswrf = make_mean_gphys('USWRF', 'NCEP', kikan, year).mean('lon') nswrf = uswrf - dswrf # 正味の短波放射フラックス nswrf.data.set_att('long_name', 'net short wave radiation flux at nominal top of atmosphere') ulwrf = make_mean_gphys('ULWRF', 'NCEP', kikan, year).mean('lon') ulwrf.data.set_att('long_name', 'upward long wave radiation flux at nominal top of atmosphere') nlwrs = make_mean_gphys('NLWRS', 'NCEP', kikan, year).mean('lon') # 正味の長波放射フラックス nlwrs.data.set_att('long_name', 'net long wave radiation flux at surface') nswrs = make_mean_gphys('NSWRS', 'NCEP', kikan, year).mean('lon') # 地表面正味の短波放射フラックス nswrs.data.set_att('long_name', 'net short wave radiation flux at surface') ## 熱 #### 潜熱フラックス lhtfl = make_mean_gphys('LHTFL', 'NCEP', kikan, year).mean('lon') lhtfl.data.set_att('long_name', 'net latent heat flux at surface') #### 顕熱フラックス shtfl = make_mean_gphys('SHTFL', 'NCEP', kikan, year).mean('lon') shtfl.data.set_att('long_name', 'net sensible heat flux at surface') #### 降水(潜熱に換算したもの) prate = make_mean_gphys('PRATE', 'NCEP', kikan, year).mean('lon') lheat = prate*L lheat.data.set_att('long_name', 'latent heat with precipitation') ## 乾燥静的エネルギー dsefl = make_mean_gphys('DSEFL', 'NCEP', kikan, year, 'dsefl') dsefl_bar = make_mean_gphys('DSEFL', 'NCEP', kikan, year, 'dsefl_bar') dsefl_dash = make_mean_gphys('DSEFL', 'NCEP', kikan, year, 'dsefl_dash') ## 潜熱輸送量 lhefl = make_mean_gphys('LHEFL', 'NCEP', kikan, year, 'lhefl') lhefl_bar = make_mean_gphys('LHEFL', 'NCEP', kikan, year, 'lhefl_bar') lhefl_dash = make_mean_gphys('LHEFL', 'NCEP', kikan, year, 'lhefl_dash') # get cos lat = lhefl.coord('lat') grid = Grid.new(lhefl.axis('lat')) cos_lat = GPhys.new(grid, cos(lat*2*PI/360.0)) dse = dsefl.integrate('level') /G * Rad * cos_lat * 2 * PI * 100 dse_bar = dsefl_bar.integrate('level') /G * Rad * cos_lat * 2 * PI * 100 dse_dash = dsefl_dash.integrate('level') /G * Rad * cos_lat * 2 * PI * 100 dse.data.set_att('units', 'W') dse_bar.data.set_att('units', 'W') dse_dash.data.set_att('units', 'W') lhe = lhefl.integrate('level') /G * Rad * cos_lat * 2 * PI * 100 lhe_bar = lhefl_bar.integrate('level') /G * Rad * cos_lat * 2 * PI * 100 lhe_dash = lhefl_dash.integrate('level') /G * Rad * cos_lat * 2 * PI * 100 lhe.data.set_att('units', 'W') lhe.data.set_att('long_name', 'latent heat energy transport (total)') lhe_bar.data.set_att('units', 'W') lhe_dash.data.set_att('units', 'W') ## EP フラックス epflx_phi = make_mean_gphys('EPFLX', 'NCEP', kikan, year, 'epflx_phi') epflx_p = make_mean_gphys('EPFLX', 'NCEP', kikan, year, 'epflx_p') epflx_div = make_mean_gphys('EPFLX', 'NCEP', kikan, year, 'epflx_div') ## 残差循環 strm_rmean = make_mean_gphys('STRM_RMEAN', 'NCEP', kikan, year) ##### ------------ 差分 ## 質量流線関数 strm = strm - strm_climat ## 東西風速 uwnd = uwnd - uwnd_climat ## 温度場 temp = temp - temp_climat ## 比湿 shum = shum - shum_climat ## 放射 nswrf = nswrf - nswrf_climat ulwrf = ulwrf - ulwrf_climat nlwrs = nlwrs - nlwrs_climat nswrs = nswrs - nswrs_climat radiation = [ nswrf, ulwrf, nswrs, nlwrs ] ## 熱 #### 潜熱フラックス lhtfl = lhtfl - lhtfl_climat #### 顕熱フラックス shtfl = shtfl - shtfl_climat #### 降水(潜熱に換 lheat = lheat - lheat_climat heat = [ lhtfl, shtfl, lheat ] ## 乾燥静的エネルギ dse = dse - dse_climat dse_bar = dse_bar - dse_bar_climat dse_dash = dse_dash - dse_dash_climat dry_ene = [ dse, dse_bar, dse_dash ] ## 潜熱輸送量 lhe = lhe - lhe_climat lhe_bar = lhe_bar - lhe_bar_climat lhe_dash = lhe_dash - lhe_dash_climat lat_ene = [ lhe, lhe_bar, lhe_dash ] ## EP フラックス epflx_phi = epflx_phi - epflx_phi_climat epflx_p = epflx_p - epflx_p_climat epflx_div = epflx_div - epflx_div_climat ## 残差循環 strm_rmean = strm_rmean - strm_rmean_climat ############################################################################# # マージンに載せる年月情報 mgn_str = { 'annual' => 'ANNUAL', 'winter' => 'DJF', 'spring' => 'MAM', 'summer' => 'JJA', 'fall' => 'SON' } # set User Path for dcldatabase DCL.glcset('DUPATH','/home/daktu32/.dcldir/') DCL.swpset('lsep', true) # ページ毎別々のファイルに落す DCL.gropn(-2) margin_info($0, mgn_str[$OPT_kikan]+" MEAN (#{year.uniq})") DCL.sgpset('lcntl', false) # DCL.uzfact(0.75) # DCL.sgpset('lcorner',false) # コーナーを取っ払う DCL.sgpset('lfprop',true) # # DCL.udpset('lmsg',false) # DCL.uscset('cyspos', 'B' ) # y 軸の単位の位置を下方へ DCL.sgpset('lfull',true) # フルスクリーン DCL.sldiv('y',2,5) ############################################################################## ##-------------------------------------- ## HEAT BALANCE (seasonal) # ビューポート設定 # vpt00 = NArray[0.2, 0.8, 0.14, 0.44] # vpt01 = NArray[0.2, 0.8, 0.14, 0.44] vpt00 = NArray[0.12, 0.88, 0.14, 0.5] vpt01 = NArray[0.11, 0.87, 0.14, 0.5] vpt10 = vpt00 vpt11 = vpt01 vpt20 = vpt00 vpt21 = vpt01 vpt30 = vpt00 vpt31 = vpt01 vpt40 = vpt00 vpt41 = vpt01 itr = 2 # dcl に欠損値情報を教える before = DCLExt.gl_set_params('lmiss'=>true,'rmiss'=>temp.data.get_att("missing_value")[0]) # 質量流線関数(左上の絵) level = [-1.0e12, 0, +1.0e12]; pattern = [35999, 35999, 70999, 70999] GGraph.next_linear_contour_options('int'=>2.5e9) ## draw GGraph.set_fig('viewport'=>vpt00, 'itr'=>itr) GGraph.set_axes('xunits'=>'', 'yunits'=>'', 'xtitle'=>'', 'ytitle'=>'') DCL.uzpset('labelxb',false) GGraph.tone( strm, true , 'annot'=>false, 'titl'=>'', 'pat'=>pattern, 'lev'=>level ) GGraph.contour( strm, false , 'annot'=>false, 'titl'=>'' ) GGraph::color_bar('constwidth'=>true) DCL.uzpset('pad1',0.2) ; DCL.uxsttl('t',"(a) Mass Stream Function(#{strm.data.units})", -1) ; DCL.uzpset('pad1',0.7) # 東西風速(右上の絵) ## set param GGraph.next_linear_contour_options('min'=>-10,'int'=>0.5, 'max'=>10) level = [-5, -2.5, 0, 2.5, 5] pattern = [30999, 35999, 40999, 70999, 75999, 85999] ## draw GGraph.set_fig('viewport'=>vpt01, 'new_frame'=>true, 'itr'=>itr) GGraph.set_axes('ytitle'=>'') GGraph.tone( uwnd, true , 'annot'=>false, 'titl'=>'', 'pat'=>pattern, 'lev'=>level ) GGraph.contour( uwnd, false , 'annot'=>false, 'titl'=>'' ) GGraph::color_bar('constwidth'=>true) DCL.uzpset('pad1',0.2) ; DCL.uxsttl('t',"(b) Zonal Wind(#{uwnd.data.units})", -1) ; DCL.uzpset('pad1',0.7) # 温度場(左下の絵) ## set param GGraph.next_linear_contour_options( 'min'=>-5, 'max'=>5, 'int'=>0.5) level = [-5, -2.5, 0, 2.5, 5] pattern = NArray[25999, 30999, 40999, 70999, 75999, 85999] ## draw GGraph.set_fig('viewport'=>vpt10, 'new_frame'=>true, 'itr'=>itr) GGraph.set_axes('ytitle'=>'','xtitle'=>nil) GGraph.tone( temp, true , 'annot'=>false, 'titl'=>'', 'pat'=>pattern, 'lev'=>level ) GGraph.contour( temp, false , 'annot'=>false, 'titl'=>'' ) GGraph::color_bar('constwidth'=>true) DCL.uzpset('pad1',0.2) ; DCL.uxsttl('t',"(c) Temperature(#{temp.data.units})", -1) ; DCL.uzpset('pad1',0.7) # 比湿(右下の絵) ## set param GGraph.next_linear_contour_options('int'=>5.0e-5) level = NArray[0.006, 0.008, 0.01, 0.012, 0.014, 0.016]*0.1 - 0.001 pattern = NArray[25999, 30999, 35999, 40999, 70999, 75999, 85999] ## draw GGraph.set_fig('viewport'=>vpt11, 'new_frame'=>true, 'itr'=>itr) GGraph.set_axes('ytitle'=>'') GGraph.tone( shum, true , 'annot'=>false, 'titl'=>'', 'pat'=>pattern, 'lev'=>level ) GGraph.contour( shum, false , 'annot'=>false, 'titl'=>'' ) GGraph::color_bar('constwidth'=>true) DCL.uzpset('pad1',0.2) ; DCL.uxsttl('t',"(d) Specific Humidity(#{shum.data.units})", -1) ; DCL.uzpset('pad1',0.7) itr = 1 # 放射(左下の絵) ## draw GGraph.set_fig('viewport'=>vpt20, 'new_frame'=>true, 'itr'=>itr) GGraph.set_axes('ytitle'=>'','xtitle'=>nil) names, idxs = plot_line_main(radiation, 'annot'=>false, 'titl'=>'', 'min'=>-20, 'max'=>20) __show_line_index(names,idxs, x=0.22) DCL.uzpset('pad1',0.2) ; DCL.uxsttl('t',"(e) Radiation Budget(#{radiation[0].data.units})", -1) ; DCL.uzpset('pad1',0.7) # 熱(右下の絵) GGraph.set_fig('viewport'=>vpt21, 'new_frame'=>true, 'itr'=>itr) GGraph.set_axes('ytitle'=>'') names, idxs = plot_line_main(heat, 'annot'=>false, 'titl'=>'', 'min'=>-20, 'max'=>20) __show_line_index(names,idxs) DCL.uzpset('pad1',0.2) ; DCL.uxsttl('t',"(f) Heat Budget(#{heat[0].data.units})", -1) ; DCL.uzpset('pad1',0.7) # 乾燥静的エネルギーフラックス(左下の絵) GGraph.set_fig('viewport'=>vpt30, 'new_frame'=>true, 'itr'=>itr) GGraph.set_axes('ytitle'=>'','xtitle'=>nil) names, idxs = plot_line_main(dry_ene, 'annot'=>false, 'titl'=>'', 'min'=>-10.0e14, 'max'=>+10.0e14) __show_line_index(names,idxs, x=0.3) DCL.uzpset('pad1',0.2) ; DCL.uxsttl('t','(g) Dry Static Energy Transport(W)',-1) ; DCL.uzpset('pad1',0.7) # 南北エネルギーフラックス(左下の絵) GGraph.set_fig('viewport'=>vpt31, 'new_frame'=>true, 'itr'=>itr) GGraph.set_axes('ytitle'=>'') names, idxs = plot_line_main(lat_ene, 'annot'=>false, 'titl'=>'', 'min'=>-8.0e14, 'max'=>+8.0e14) __show_line_index(names,idxs) DCL.uzpset('pad1',0.2) ; DCL.uxsttl('t','(h) Latent Heat Transport(W)',-1) ; DCL.uzpset('pad1',0.7) itr = 2 # EP フラックス(左下の絵) ## settings level = NArray[-5.0e1, -1.0e1, 0, 1.0e1, 5.0e1]*1.0e-6 pattern = NArray[30999, 35999, 40999, 55999, 75999, 85999] vect_ratio = NArray[1.0, 100.0]*8e-12 vect_fact = 50.0 ## draw GGraph.set_fig('viewport'=>vpt40, 'new_frame'=>true, 'itr'=>itr) GGraph.set_axes('ytitle'=>'','xtitle'=>nil) DCL.uzpset('labelxb',true) GGraph.tone( epflx_div, true , 'annot'=>false, 'titl'=>'', 'pat'=>pattern, 'lev'=>level ) GGraph::vector_for_dcchart(epflx_phi, - epflx_p, false, 'xfact'=>vect_ratio[0]*vect_fact, 'yfact'=>vect_ratio[-1]*vect_fact, 'unit'=>true, 'annot'=>false, 'itr'=>2, 'xintv'=>2 ) GGraph::color_bar('constwidth'=>true) DCL.uzpset('pad1',0.2) ; DCL.uxsttl('t','(i) EP Flux(kg.s-2), Div F(m.s-2)',-1) ; DCL.uzpset('pad1',0.7) # 残差循環(左下の絵) ## settings level = [-1.0e12, 0, +1.0e12]; pattern = [35999, 35999, 70999, 70999] GGraph.next_linear_contour_options('int'=>2.5e9) ## draw GGraph.set_fig('viewport'=>vpt41, 'new_frame'=>true, 'itr'=>itr) GGraph.set_axes('ytitle'=>'') GGraph.tone( strm_rmean, true , 'annot'=>false, 'titl'=>'', 'pat'=>pattern, 'lev'=>level ) GGraph::color_bar('constwidth'=>true) GGraph.contour( strm_rmean, false , 'annot'=>false, 'titl'=>'' ) DCL.uzpset('pad1',0.2) ; DCL.uxsttl('t','(j) Residual Mass Stream Function(kg/s)',-1) ; DCL.uzpset('pad1',0.7) DCL.grcls #######################################################