#! /usr/local/bin/ruby

ncfile = "../2d-dg2004-run1.nc"
tstart = 0.0
tend = 1500.0
dt = 1.0
n = ((tend-tstart)/dt).to_i + 1

(1..n-1).each do |i|
  time = i*dt
#  print "proccessing #{time}\n"
  fn = sprintf("%06d", i)
  system "gpview --nocont --srange -59:59 --wsn 2 #{ncfile}@vor,t=#{time} > /dev/null"
  system "convert -density 144 -geometry 50% -rotate 90 dcl.ps pngs/t#{fn}.png"
end
system "rm dcl.ps"