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

# NetCDF ファイル "air.2m.gauss.2021.nc" から変数 "air" を読み, GPhys オブジェクト gpa に格納
gpa = GPhys::IO.open( "air.2m.gauss.2021.nc", "air" )
gpa = gpa.cut('time'=>DateTime.parse("2021-01-01 00:00:0.0")..DateTime.parse("2021-12-31 00:00:0.0"))
gpa = gpa.cut('lon'=>135)
gpa = gpa.cut('lat'=>35.2375)

# NetCDF ファイル "csdsf.sfc.gauss.2021.nc" から変数 "csdsf" を読み, GPhys オブジェクト gpc に格納
gpc = GPhys::IO.open( "csdsf.sfc.gauss.2021.nc", "csdsf" )
gpc = gpc.cut('time'=>DateTime.parse("2021-01-01 00:00:0.0")..DateTime.parse("2021-12-31 00:00:0.0"))
gpc = gpc.cut('lon'=>135)
gpc = gpc.cut('lat'=>35.2375)

# 画面を開く (open)
#   引数の 2 はファイルへの出力を表す
DCL.swcset('fname','flux1')
DCL.swpset("ifl",1)
DCL.gropn(2)

# 描画画面を準備
#   itr の 1 は横軸, 縦軸ともに線形を表す
GGraph.set_fig( 'itr'=> 1, 'window'=>[230,320,0,440] )

# 散布図の点を描画
#   第一引数は描画する点の x 成分の GPhys オブジェクト
#   第二引数は描画する点の y 成分の GPhys オブジェクト
GGraph.scatter( gpa, gpc, true)

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