=begin
 表題  shallow_zd_case5.rb

 使い方  ruby shallow_zd_case5.rb

 履歴  2004/06/16 小高正嗣  gave より作成 
=end

require "numru/gphys"
require "numru/dcl"
require "colorbar.rb"

include NumRu

filename = "shallow_zd_case5.nc"

varname = "zeta"
Zeta0 = GPhys::NetCDF_IO.open(filename,varname)

varname = "h"
H0 = GPhys::NetCDF_IO.open(filename,varname)

iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
DCL::gropn(-iws)
DCL::sldiv('Y', 2, 4)

xmin = 0.0
xmax = 354.375
ymin = -85.76058197
ymax = 85.76058197

vxmin = 0.1
vxmax = 0.9
vymin = 0.15
vymax = 0.55

x_title = "longitude"
x_unit = "deg."
y_title = "latitude"
y_unit = "deg."

# カラーバーの設定 ###################################################
$cbar_conf = {
  "vx0"=>0.05,            # カラーバーの左下角の x 座標
  "vy0"=>0.05,            # カラーバーの左下角の y 座標
  "vxlength"=>0.25,      # カラーバーの 横の長さ
  "vylength"=>0.02,      # カラーバーの 縦の長さ
  "tick1"  => 1,         # 目盛(大)を付ける間隔
  "tick2"  => 2,        # 目盛(小)を付ける間隔
  "eqlev"  => false,     # 
  "nobound"=> false      # 
  }
#########################################################################

DCL::sglset("LCORNER", false )
DCL::sglset("LCNTL", false )
DCL::udlset("LMSG", false )
DCL::uglset("LMSG", false )
DCL::gllset("LMISS", true )
DCL::uzrset('RSIZEC1', 0.021)
DCL::uzrset('RSIZEC2', 0.028)
DCL::uzrset('RSIZEL1', 0.014)
DCL::uzrset('RSIZEL2', 0.021)
DCL::sglset('LFULL', true)

for i in [0, 5, 10, 15]

Zeta = Zeta0[0..-1,0..-1,i]
H = H0[0..-1,0..-1,i]

DCL::grfrm
DCL::grswnd(xmin,xmax,ymin,ymax)
DCL::grsvpt(vxmin,vxmax,vymin,vymax)
DCL::grstrn(1)
DCL::grstrf

DCL::sglset("LCLIP", true )
DCL::uwsgxa(Zeta.coord(0).val)
DCL::uwsgya(Zeta.coord(1).val)
DCL::ueitlv
DCL::uegtla(-3.520313248e-05, 3.669281487e-05, 0 )
DCL::uetone(Zeta.val)
DCL::udgcla(-3.520313248e-05, 3.669281487e-05, 0 )
DCL::udcntz(Zeta.val)
DCL::ussttl(x_title, x_unit, y_title, y_unit)
DCL::usdaxs
DCL::sglset("LCLIP", false )
DCL::uzrset("ROFFXT", 0.06)
title = "Zeta(1/s)"
DCL::uxsttl("t", title, 0 )
DCL::Util::color_bar($cbar_conf)

DCL::grfrm
DCL::grswnd(xmin,xmax,ymin,ymax)
DCL::grsvpt(vxmin,vxmax,vymin,vymax)
DCL::grstrn(1)
DCL::grstrf

DCL::sglset("LCLIP", true )
DCL::uwsgxa(H.coord(0).val)
DCL::uwsgya(H.coord(1).val)
DCL::ueitlv
DCL::uegtla(4996.775762, 5957.739839, 0 )
DCL::uetone(H.val)
DCL::udgcla(4996.775762, 5957.739839, 0 )
DCL::udcntz(H.val)
DCL::ussttl(x_title, x_unit, y_title, y_unit)
DCL::usdaxs
DCL::sglset("LCLIP", false )
DCL::uzrset("ROFFXT", 0.06)
title = "h(m)"
DCL::uxsttl("t", title, 0 )
DCL::Util::color_bar($cbar_conf)

end

DCL::grcls