require "numru/gphys"
require "numru/ggraph"
require "numru/dcl"

include NumRu


filename = "3dMars_dx200dz100_02_VorZ.nc"
varname = "VorZ"
gphys0 = GPhys::IO.open(filename,varname)
gphys0 = gphys0[0..99,0..99,0,6]


DCL::gropn(2)

xmin = 100.0
xmax = 19900.0
ymin = 100.0
ymax = 19900.0

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

x_title = "x"
x_unit = "m"
y_title = "y"
y_unit = "m"

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

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(-0.05, 0.05, 0.002 )
#DCL::udiclv
#DCL::udgcla(-10.0, 10.0, 5.0 )
DCL::uetone(gphys0.val)
#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=50 m",0.02,0,-1,1)
#i=1
#DCL::sgtxzr(vxmax+0.01,vymax-0.03*i,"t=43200 s",0.02,0,-1,1)
###
DCL::uzrset("ROFFXT", 0.06)

title = "Vorticity"
DCL::uxsttl("t", title, 0 )
GGraph.color_bar

DCL::grcls