#!/usr/bin/ruby
#使用するライブラリの読み込み
require "numru/ggraph"
include NumRu

#NetCDFファイル"dlwrf.sfc.gauss.day.ltm.1991-2020.nc"から変数"dlwrf"を読み、GPhysオブジェクトgp1に格納
gp1 = GPhys::IO.open( "dlwrf.sfc.gauss.day.ltm.1991-2020.nc", "dlwrf" )
#NetCDFファイル"ulwrf.ntat.gauss.day.ltm.1991-2020.nc"から変数"ulwrf"を読み、GPhysオブジェクトgp2に格納
gp2 = GPhys::IO.open( "ulwrf.ntat.gauss.day.ltm.1991-2020.nc", "ulwrf" )

#"dlwrf"の設定
#経度(lon)軸に対して平均する
gp1 = gp1.mean('lon')
#時間(time)に対して平均する
gp1 = gp1.mean('time')

#"ulwrf"の設定
#緯度(lon)軸に対して平均する
gp2 = gp2.mean('lon')
#時間(time)に対して平均する
gp2 = gp2.mean('time')

# 画面を開く (open)
#   引数の 1 は画面への描画を表す
#          2 はファイルへの出力を表す
#            (デフォルトでは出力は pdf 形式でファイル名は dcl.pdf)
DCL.gropn(2)

# 描画に関わる各種の指定
DCL.sgpset('lfull',true) # 描画範囲を最大に設定 (ウィンドウのギリギリま で使用)
DCL.uzfact(0.3)          # フォントサイズ (元の値に対する比率)
DCL.sgpset('lclip',true) # 軸範囲を超えた領域を描画しない
DCL.glpset('lmiss',true) # 欠損値処理

# 描画画面を準備
#   itr の 1 は横軸, 縦軸ともに線形を表す
#          2 は横軸が線形軸, 縦軸が対数軸を表す
#          3 は横軸が対数軸, 縦軸が線形軸を表す
#          4 は横軸, 縦軸ともに対数軸を表す
GGraph.set_fig( 'itr'=> 1, 'viewport'=>[0.1,0.9,0.2,0.55], 'window'=>[-90,90,100,500] )

#折れ線グラフを描画
GGraph.line( gp1, true )
GGraph.line( gp2, false, 'index'=>20 )


# 画面を閉じる (close)
DCL.grcls