=begin

=exner-cont.rb

$B%(%/%9%J!<4X?t$N(B x-t $B%3%s%?!<?^$rIA$/(B Ruby $B%9%/%j%W%H(B

$B3+H/HG(B : $B3+H/L\I8$O%+%i!<%P!<$r:n@.$9$k$3$H(B

=end

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

include NumRu

filename = ARGV[0]
varname = "Exner"
gphys0 = GPhys::NetCDF_IO.open(filename,varname)
gphys0 = gphys0[0..-1,0,0..-1]

if ARGV.index("-ps") then
  DCL::gropn(2)
elsif
  DCL::gropn(4)
end

xmin = 0.0
xmax = 10000.0
ymin = 0.0
ymax = 20.0

valmin = 0.998
valmax = 1.008

vxmin = 0.2
vxmax = 0.8
vymin = 0.2
vymax = 0.8

x_title = "x"
x_unit = "m"
y_title = "t"
y_unit = "s"

DCL::sglset("LCNTL", false )
DCL::udlset("LMSG", false )
DCL::uglset("LMSG", false )
DCL::gllset("LMISS", true )

# $B<jF0$G%+%i!<%H!<%s$r;XDj(B
#
#cltstep = 10
#delta = (valmax - valmin) / cltstep
#for k in 1..10
#  tlev1 = valmin + (k-1)*delta
#  tlev2 = tlev1 + delta
#  ipat = (13+8*k)*1000 + 999
#  DCL::uestlv(tlev1,tlev2,ipat)
#end

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(gphys0.coord(0).val)
DCL::uwsgya(gphys0.coord(1).val)
DCL::ueitlv
DCL::uegtla(valmin, valmax, 0 )
DCL::uetone(gphys0.val)
DCL::udgcla(valmin, valmax, 0 )
DCL::udcntz(gphys0.val)
DCL::ussttl(x_title, x_unit, y_title, y_unit)
DCL::usdaxs
DCL::sglset("LCLIP", false )
i=0
DCL::sgtxzr(vxmax+0.01,vymax-0.03*i,"z=5050 m",0.02,0,-1,1)
DCL::uzrset("ROFFXT", 0.06)
title = "Exner(1)"
DCL::uxsttl("t", title, 0 )


#--- $B%+%i!<%P!<(B 

MyTool::colorbar(valmin,valmax)


DCL::grcls

