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

# NetCDF ファイルから東西風と南北風の変数を読み, GPhys オブジェクトに格納
gpu = GPhys::IO.open( "uwnd.2023.nc", "uwnd" )
gpv = GPhys::IO.open( "vwnd.2023.nc", "vwnd" )

#強い高度である200mbar(hPa)のデータを切り出す (cut)
gpu = gpu.cut('level'=>200)
gpv = gpv.cut('level'=>200)

# GPhys オブジェクト gpu (gpv)内の変数 "uwnd"("vwnd") から7月のデータを切り出す (cut)
#時間(time)軸に対して平均する
gpu = gpu.cut('time'=>DateTime.parse("2023-07-01 00:00:0.0")..DateTime.parse("2023-08-01 00:00:0.0")).mean('time')
gpv = gpv.cut('time'=>DateTime.parse("2023-07-01 00:00:0.0")..DateTime.parse("2023-08-01 00:00:0.0")).mean('time')

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

# 描画画面を準備
#   itr の10 は正距円筒図法 (いわゆる経度-緯度座標) を表す
GGraph.set_fig( 'itr'=> 10 )

# 地図情報の指定. ここでは海岸線を描くように設定.
GGraph.set_map( 'coast_world'=>true, 'grid'=>true )

# ベクトルを描画
#   第一引数は描画するベクトルの x 成分の GPhys オブジェクト
#   第二引数は描画するベクトルの y 成分の GPhys オブジェクト
#   第三引数の true は新しい図を表す
GGraph.vector( gpu, gpv, true, "xintv"=>2, "yintv"=>2, "factor"=>2, "unit_vect"=>true,'map_axes'=>true,'index'=>40 )

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