require "numru/ggraph"
include NumRu

file  = 'thermal1_VelZ.nc'
var   = 'VelZ'

gphys  = GPhys::IO.open( file , var  )
# PTemp の nc は y 軸なし-> t 軸は 2
timearray = gphys.coord(3).val    #gphysの軸(0):x, (1):y, (2):z, (3):t
#gphys.close

### netcdf の時間でまわす
#t=0
for t in timearray[0..90] #時間割る描画間隔

### t を出す

### よびだし
  system('ruby draw_2dthermal_VelZ.rb ' + t.to_s )

### とりあえず convert
  system("convert -rotate 90 dcl.ps dcl.gif" )

### 名前付け
###  
  com = sprintf( "mv dcl.gif figure/Vel/2dmoist40k20z-VelZ%0#{5}d.gif", t )
  system(com)

end