require "numru/ggraph"
include NumRu

level=2500.0
time=43200.0

InputFile1='3dMars_dx200dz200_02_VorZ.nc'
InputFile2='3dMars_dx200dz200_02_VelX.nc'
InputFile3='3dMars_dx200dz200_02_VelY.nc'
InputFile4='3dMars_dx200dz200_02_VelZ.nc'
InputFile5='3dMars_dx200dz200_02_PotTempDev.nc'


# データの読み込み
#
VorZ = GPhys::NetCDF_IO.open(InputFile1, 'VorZ' )
VelX = GPhys::NetCDF_IO.open(InputFile2, 'VelX' )
VelY = GPhys::NetCDF_IO.open(InputFile3, 'VelY' )
VelZ = GPhys::NetCDF_IO.open(InputFile4, 'VelZ' )
PotTemp = GPhys::NetCDF_IO.open(InputFile5, 'PotTemp' )


DCL.gropn(1)
DCL.sgpset('lcntl', false) 

GGraph.tone( VorZ.cut('z'=>level, 't'=>time), true,
             'max'=>0.02, 'min'=>-0.02,
             'interval'=>0.001,	
             'title'=>'Vorticity and Wind Vector' )
GGraph.color_bar

GGraph.vector( VelX.cut('z'=>level, 't'=>time), 
               VelY.cut('z'=>level, 't'=>time), false,
               'xintv'=>5, 'yintv'=>5, 
               'unit_vect'=>true )

GGraph.tone( PotTemp.cut('z'=>level, 't'=>time), true,
             'max'=>2.0, 'min'=>-1.0,
             'interval'=>0.1,	
             'title'=>'Potential Temp. and Vertical Velocity' )
GGraph.color_bar

GGraph.contour( VelZ.cut('z'=>level, 't'=>time), false,
               'interval'=>2.5 )




DCL.grcls