### 値の取得
###

var  = "VelX"
VelX = GPhys::IO.open(ARGV[0],var).cut('z'=>0..zlev)

if ( $OPT_max )
  maximum = $OPT_max.to_f
else
  maximum = 50.0
end

if ( $OPT_min )
  minimum = $OPT_min.to_f
else
  minimum = -50.0
end


###
### 図の作成
###
WindowOpen( aspect, $OPT_dump, true )

GGraph.set_axes('ytitle'=>'Height','xtitle'=>'')
GGraph.set_axes('xunits'=>'km','yunits'=>'km')

for i in 0..num
  time = time0 + i * deltime
  p time
  
  #水平速度
  GGraph.tone( VelX.cut('t'=>time), 
               true, 
               'annotate'=>infosw,
               'levels' =>[-1000,minimum],
               'pat'=>[14999]
               )
  GGraph.tone( VelX.cut('t'=>time), 
               false, 
               'annotate'=>infosw,
               'levels' =>[maximum,1000],
               'pat'=>[96999]
               )
  GGraph.tone( VelX.cut('t'=>time), 
               false, 
               'title' => 'Horizontal Velocity',
               'annotate'=>infosw,
               'nlev' => 50,
               'max' => maximum,
               'min' => minimum
               )
#  GGraph.color_bar('vlength'=>0.275)
#  GGraph.color_bar( 'landscape' => true )
  GGraph.color_bar

end

DCL.grcls


if ( $OPT_dump )

  DumpPNG( time0, deltime, num )
  
end