require "numru/ggraph"
include NumRu

file  = 'EPT.nc'
var   = 'EPT'

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

### netcdf の時間でまわす

for t in timearray[0..90] #時間割る描画間隔

### t を出す
  p t

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


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

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

end