#!/usr/bin/ruby
# 1陦檎岼縺ｯinterpreter縺ｮ謖�ｮ�(shebang)

# 繝ｩ繧､繝悶Λ繝ｪ縺ｮ隱ｭ縺ｿ霎ｼ縺ｿ
require "numru/ggraph"
include NumRu

# NetCDF繝輔ぃ繧､繝ｫ"air.2016.nc"縺九ｉ螟画焚"air"繧定ｪｭ縺ｿ縲；Phys繧ｪ繝悶ず繧ｧ繧ｯ繝�tmp縺ｫ譬ｼ邏�
tmp = GPhys::IO.open( "air.2016.nc", "air" )

# 遲牙悸髱｢繧�500hPa縺ｫ謖�ｮ�
tmp = tmp.cut('level'=>500)

# png 蠖｢蠑上〒譖ｸ縺榊�縺�
DCL.swpset("ifl",1)
DCL.gropn(2)

# 繧ｫ繝ｩ繝ｼ繝槭ャ繝励ｒ謖�ｮ�
DCL.sgscmn(3) 

# 謠冗判逕ｻ髱｢繧呈ｺ門ｙ (itr 縺ｮ 32 縺ｯ豁｣霍晄婿菴榊峙豕�)
# viewport縺ｧ謠冗判遽�峇繧呈欠螳�
GGraph.set_fig( 'itr'=> 32, 'viewport'=>[0.15,0.75,0.15,0.85] ) 

# 豬ｷ蟯ｸ邱壹→繧ｰ繝ｪ繝�ラ繧呈緒逕ｻ
GGraph.set_map( 'coast_world'=>true, 'grid'=>true )

# 繝ｬ繝吶Ν縺ｮ謖�ｮ�
levels = [220,230,240,250,260,270]

# 繝ｬ繝吶Ν豈弱�濶ｲ縺ｮ謖�ｮ�
patterns =  [63999,72999,81999,90999,99999]

# 謠冗判髢句ｧ区凾蛻ｻ繧貞､画焚縺ｫ譬ｼ邏�
time_label = DateTime.parse("2016-01-15 00:00:0.0")

# i縺�1縺九ｉ52縺ｾ縺ｧ郢ｰ繧願ｿ斐☆
 for i in 1..52

   # 繝医�繝ｳ縺ｧ謠冗判
   # 隨ｬ荳€蠑墓焚縺ｧ謠冗判縺吶ｋGphys繧ｪ繝悶ず繧ｧ繧ｯ繝医ｒ謖�ｮ�
   # 譁ｰ隕上↓謠冗判縺吶ｋ縺溘ａ隨ｬ莠悟ｼ墓焚縺ｯtrue縺ｨ縺吶ｋ
   # annotate 繧� true 縺ｫ險ｭ螳壹＠縺ｦ蝗ｳ縺ｮ蜿ｳ荳翫↓諠��ｱ繧呈嶌縺�
   # map_axes 繧� true 縺ｫ險ｭ螳壹＠縺ｦ霆ｸ縺ｮ繝ｩ繝吶Ν繧呈嶌縺�
   # lev, pat 縺ｧ 荳願ｨ倥〒謖�ｮ壹＠縺� levels, patterns繧偵Ξ繝吶Ν縺ｨ濶ｲ縺ｫ險ｭ螳� 
   GGraph.tone_and_contour( tmp.cut('time'=>time_label), true, 'annotate'=>false, 'map_axes'=>true, 'lev'=>levels, 'pat'=>patterns )

   # 繧ｫ繝ｩ繝ｼ繝舌�繧呈緒逕ｻ
   GGraph.color_bar

   time_display = time_label.to_s  # time_display 縺ｫ time_label縺ｮ譁�ｭ怜梛縺ｸ縺ｮ螟画鋤繧剃ｻ｣蜈･
   date = time_display.slice(0,10) # time_display 縺九ｉ 譌･莉倥ｒ蜿悶ｊ蜃ｺ縺�
   time_utc = time_display.slice(11,8) + " (UTC)" # time_display 縺九ｉ 譎ょ綾繧貞叙繧雁�縺�

   DCL::sgstxs(0.02)  # 繝�く繧ｹ繝医�繧ｵ繧､繧ｺ
   DCL::sgstxi(1)  # 繝ｩ繧､繝ｳ繧ｿ繧､繝�
   DCL::sgtxv(0.84, 0.8, 'level = 500 hPa') # 豌怜悸縺ｮ陦ｨ遉ｺ
   DCL::sgtxv(0.817, 0.77, date )  # 譌･莉倥�陦ｨ遉ｺ
   DCL::sgtxv(0.829, 0.74, time_utc ) # 譎ょ綾縺ｮ陦ｨ遉ｺ
   DCL::sgtxv(0.79, 0.5, '[K]') # 貂ｩ蠎ｦ縺ｮ蜊倅ｽ阪�陦ｨ遉ｺ

   # 騾ｲ謐礼憾豕√ｒ陦ｨ遉ｺ縺吶ｋ
   print i, "page/52pages", "\n"

   # 譎ょ綾繧抵ｼ匁凾髢� ( = 0.25譌･) 騾ｲ繧√ｋ
   time_label = time_label + 0.25
   
end

# 逕ｻ髱｢繧帝哩縺倥ｋ
DCL.grcls