require "numru/gphys"
require "numru/ggraph"
include NumRu

file1 = '3dMars_dx200dz100_02_VelX.nc'
file2 = '3dMars_dx200dz100_02_VelY.nc'
file3 = '3dMars_dx200dz100_02_VelZ.nc'
file4 = '3dMars_dx200dz100_02_VorZ.nc'

U = GPhys::NetCDF_IO.open(file1, 'VelX')
V = GPhys::NetCDF_IO.open(file2, 'VelY')
W = GPhys::NetCDF_IO.open(file3, 'VelZ')
VorZ = GPhys::NetCDF_IO.open(file4, 'VorZ')

U1 = U[0..-1,0..-1,0,0..-1]   
V1 = V[0..-1,0..-1,0,0..-1]   
W1 = W[0..-1,0..-1,0,0..-1]   
VorZ1= VorZ[0..-1,0..-1,0,0..-1]   

Ua = (U1**2 + V1**2)**0.5
dudz = Ua/50


DCL.gropn(1)
GGraph.tone( dudz[0..-1,0..-1,6] )
#GGraph.color_bar

#GGraph.tone( VorZ1[0..-1,0..-1,6] )
#GGraph.contour( W1[0..-1,0..-1,6] )

GGraph.color_bar

DCL.grcls