#!/usr/bin/ruby
# 使用するライブラリの読み込み. (以下 2 行は「決まり文句」.)
require "numru/ggraph"
include NumRu

# NetCDF ファイルから東西風の変数を読み, GPhys オブジェクトに格納
gp = GPhys::IO.open( "uwnd.mon.mean.nc", "uwnd" )

# GPhys オブジェクト gp 内の変数 "uwnd" から lon=0のデータを切り出す (cut)
#Monthly mean u wind=14m/sのデータを切り出す
#latで平均
gp = gp.cut('lon'=>0).mean('lat')

# 画面を開く (open)
#          2 はファイルへの出力を表す
DCL.gropn(2)

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

# 描画画面を準備
#   itr の 1 は横軸, 縦軸ともに線形を表す
#   viewport は描画可能範囲の中での図の描画範囲を表す
#            (四隅 [x1,x2,y1,y2] の全体に対する比率)
#   window は描画する軸の範囲を表す (四隅 [x1,x2,y1,y2] の値)
GGraph.set_fig( 'itr'=> 1, 'viewport'=>[0.1,0.9,0.2,0.55], 'window'=>[DateTime.parse("2000-01-01 00:00:0.0"),DateTime.parse("2021-01-01 00:00:0.0"),0,1000]
)

# 折れ線グラフを描画
#   第一引数は描画するデータの GPhys オブジェクト
#   第二引数の true はこの時点では「決まり文句」
# 縦軸と横軸を入れ替える
GGraph.line( gp,true,"exchange"=>true )

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