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

include NumRu

filename = "2dMars_dx200dz100_VelZ.nc"
varname = "VelZ"
gphys1 = GPhys::IO.open(filename,varname)
gphys1 = gphys1[0..99,0..99,12]


DCL::gropn(-2)

xmin = 100.0
xmax = 19900.0
zmin = 50.0
zmax = 9950.0

vxmin = 0.1
vxmax = 0.9
vzmin = 0.3
vzmax = 0.7

x_title = ""
x_unit = "m"
z_title = ""
z_unit = "m"

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

DCL::grfrm
DCL::grswnd(xmin,xmax,zmin,zmax)
DCL::grsvpt(vxmin,vxmax,vzmin,vzmax)
DCL::grstrn(1)
DCL::grstrf

DCL::sglset("LCLIP", true )
DCL::uwsgxa(gphys1.coord(0).val)
DCL::uwsgya(gphys1.coord(1).val)
DCL::ueitlv
DCL::uegtla(-10.0, 15.0, 0.5 )
DCL::udiclv
DCL::udgcla(-10.0, 15.0, 5.0 )
DCL::uetone(gphys1.val)
DCL::udcntz(gphys1.val)
DCL::ussttl(x_title, x_unit, z_title, z_unit)
DCL::usdaxs
DCL::sglset("LCLIP", false )
DCL::uzrset("ROFFXT", 0.06)

title = "Vertical Velocity (2D)"
DCL::uxsttl("t", title, 0 )
GGraph.color_bar

DCL::grcls