# -*- coding: utf-8 -*-
require "numru/ggraph"
include NumRu

file  = 'thermal-moist_PTempAll.nc'
var   = 'PTempAll'
#file  = 'thermal-moist_PressAll-20121219-Default_z48km.nc'
#var   = 'PressAll'

gphys  = GPhys::IO.open( file , var  )
timearray = gphys.coord(3).val

for t in timearray

  timestep = t % ( 3600 * 3 )
  if timestep == 0 then
  p t

  system('ruby draw_xz_rainwatermixingratio.rb ' + t.to_s )
  system('ruby psfat.rb'  + ' dcl.ps > dclnew.ps' )
  system("convert -rotate 90 dclnew.ps dcl.gif" )
#  com = sprintf( "mv dcl.gif zzz-figs/xz-rainwatermixingratio/dcl_H2O-l-RainAll_xz_%0#{8}d.gif", t )
#  com = sprintf( "mv dcl.gif zzz-figs/zzzNew/dcl_H2O-l-RainAll_xz_%0#{8}d.gif", t )
#  com = sprintf( "mv dcl.gif zzz-figs/zzzNew2/xz-rain/dcl_H2O-l-RainAll_xz_1-2-3_z48km_%0#{8}d.gif", t )
  com = sprintf( "mv dcl.gif zzz-figs/zzzNew2/xz-rain/dcl_H2O-l-RainAll_xz_1-4-3_z48km_%0#{8}d.gif", t )
  system(com)

  end
end


### z = 0-20km ###

for t in timearray

  timestep = t % ( 3600 * 3 )
  if timestep == 0 then
  p t

  system('ruby draw_xz_rainwatermixingratio_z=0-20km.rb ' + t.to_s )
  system('ruby psfat.rb'  + ' dcl.ps > dclnew.ps' )
  system("convert -rotate 90 dclnew.ps dcl.gif" )
#  com = sprintf( "mv dcl.gif zzz-figs/xz-rainwatermixingratio_z=0-15km/dcl_H2O-l-RainAll_xz_%0#{8}d.gif", t )
#  com = sprintf( "mv dcl.gif zzz-figs/zzzNew/dcl_H2O-l-RainAll_xz_z=0-20km_%0#{8}d.gif", t )
#  com = sprintf( "mv dcl.gif zzz-figs/zzzNew2/xz-rain_z=0-20km/dcl_H2O-l-RainAll_xz_z=0-20km_1-2-3_z48km_%0#{8}d.gif", t )
  com = sprintf( "mv dcl.gif zzz-figs/zzzNew2/xz-rain_z=0-20km/dcl_H2O-l-RainAll_xz_z=0-20km_1-4-3_z48km_%0#{8}d.gif", t )
  system(com)

  end
end