#!/usr/bin/ruby
# coding: utf-8

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

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

#使用するデータの指定
#緯度：0度(赤道),高度：100hPa-1hPa,東西風速は経度方向に平均した
gp = gp.cut('lat'=>0).mean('lon').cut('level'=>100..1)
#時間は2012年から2020年を切り出した
gp = gp.cut('time'=>DateTime.parse("2012-01-01 00:00:0.0")..DateTime.parse("2020-12-31 00:00:0.0"))

#図をpdfファイルとして保存
DCL.gropn(2)

DCL.sgpset('lfull',true)  #描画範囲を最大に設定
DCL.uzfact(0.6)           #フォントサイズの設定
DCL.sgpset('lclip',true)  #軸範囲を超えた領域を描画しない
DCL.glpset('lmiss',true)  #欠損値処理

#描画画面の設定(横軸は線形軸,縦軸は対数軸に指定．図の描画範囲の設定)
GGraph.set_fig('itr'=>2,'viewport'=>[0.1,0.85,0.2,0.55])

#グラフを色づけ(トーン)と等値線で描画, 縦軸と横軸の入れ替え
GGraph.tone_and_contour(gp,true,"exchange"=>true)
GGraph.color_bar          #カラーバーを表示


#画面を閉じる.
DCL.grcls