#!/usr/bin/ruby
# coding: utf-8
#使用するらいぶらりの読込
require "numru/ggraph"
include NumRu

#NetCDF
#格納
gp = GPhys::IO.open("rhum.2019.nc","rhum")

#cutmethodを使うところ

#画面を開く
#  引数の1が画面への描写、2がファイル出力（デフォルトでpdfファイル化される）
DCL.swpset("ifl",1)
DCL.gropn(2)

#各種設定
DCL.sgpset('lfull',true) # 描画範囲を最大に設定 (ウィンドウのギリギリま で使用)
DCL.uzfact(0.6)          # フォントサイズ (元の値に対する比率)
DCL.sgpset('lclip',true) # 軸範囲を超えた領域を描画しない
DCL.glpset('lmiss',true) # 欠損値処理
GGraph.set_axes('xlabelint'=>30) # x 軸にラベルを書く間隔
GGraph.set_axes('ylabelint'=>30) # y 軸にラベルを書く間隔

#画面の準備
#1縦横線形。、2たてが対数、3横が対数、4が縦横対数
GGraph.set_fig('itr'=>10)

GGraph.set_map('coast_world'=>true,'grid'=>true)

#forを使うことで連続させる
for k in 1..12
for i in 5..10   

  #time_label = "2019-" + i.to_s + " " + j.to_s + k.to_s + " 00:00"

  #6時間毎のデータの表示    
    time_label = "2019-" + k.to_s + "-" + i.to_s

  #繰返し回数
  print k," : ",time_label,"￥n"

#湾岸線を書いて下さい
#GGraph.set_map('coast_world'=>true,'grid'=>true)

#色付描写
#第一引数はGPhysおぶじぇくと
#第二引数はtrue faulseで決める
GGraph.tone(gp.cut('time'=>DateTime.parse(time_label)),true,'annotate'=>true,'map_axes'=>true,'nlev'=>30)
#カラーバーの描画
GGraph.color_bar

#等値線を描く
GGraph.contour(gp,false)

#ベクトルの描画

end
end


#おしまい
DCL.grcls