require "numru/ggraph"
include NumRu

#vname1 = 'SurfTemp'
#vname2 = 'Rain'
#vname3 = 'SoilMoist'
#vname4 = 'SurfSnow'

gphys = GPhys::NetCDF_IO.open("data/Ps.nc", "Ps")
tm = gphys.shape[2] # number of elements for time      dimension

num = 1 * 10
#num = 4 * 10

for t in 0..tm/num-num/num
  p 0.to_s + ' ' + (t*num).to_s + ' '+ (tm-1).to_s

  system('ruby drawzm.rb ' + (t*num).to_s + ' ' + ((t+1)*num-1).to_s )

  system("pstoimg -flip r270 dcl.ps")
  com = sprintf( "mv dcl.png figs_zm/dcl%0#{5}d.png", t*num )
  system(com)
end