#!/usr/bin/ruby require "numru/ggraph" include NumRu gp = GPhys::IO.open("air.mon.mean.nc","air") #月平均気温のデータ gp = gp.cut('lon'=>135,'lat'=>35) #東経135度、北緯35度を選択 DCL.swpset("ifl",1) #pngで出力 DCL.gropn(2) GGraph.set_fig('itr'=>1) #両軸線形グラフ yr = 1980 #年の設定 id = 10 #線の色の設定 #基準とする1980年の温度分布 gp80 = gp.cut('time'=>DateTime.parse(yr.to_s + "-08-01 00:00:0.0")) #1980年のデータで規格化された1980年の温度分布(つまり1) gpy = gp.cut('time'=>DateTime.parse(yr.to_s + "-08-01 00:00:0.0")) / gp80 GGraph.line(gpy,true,"exchange"=>true,"min"=>0.5,"max"=>2) for i in 1..4 #1990~2020 yr = yr + 10 #先のデータの十年後 id = id + 10 #線の色の変更:黒→赤→緑→青→黄 #1980年で規格化された各年の温度分布 gpy = gp.cut('time'=>DateTime.parse(yr.to_s + "-08-01 00:00:0.0")) / gp80 GGraph.line(gpy,false,"exchange"=>true,"index"=>id,"min"=>0.5,"max"=>2) end DCL.grcls