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

#file  = 'thermal-moist_PTempAll_rank000000.nc'
#var   = 'PTempAll'

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

# added
IntegPeriod  = 4320000.0e0
tn = 20
dt = IntegPeriod / tn
##IntegPeriod  = 432000.0e0
##tn = 2
##dt = IntegPeriod / tn

##for t in [3600,4200]
#for t in timearray
##for t in timearray[0..2]
##for t in timearray[6*24*5..timearray.size-1]
for timestep in 1..tn
#timestep = 20

#  p t
  p timestep

##  system('ruby draw_core.rb ' + file0 + ' ' + var0 + ' ' + file + ' ' + var + ' ' + t.to_s )
##  system('ruby draw_core3.rb ' + t.to_s )
##  system('ruby draw_core4.rb ' + t.to_s )
#  system('ruby draw_core4_mod.rb ' + t.to_s )

#  system('ruby draw_xt.rb ' + timestep.to_s + ' ' + tn.to_s )
#  system('ruby draw_tz.rb ' + timestep.to_s )
#  system('ruby draw_tz_exnermeanx.rb ' )
#  system('ruby draw_xt_sfcheatflux.rb ' + timestep.to_s )
#  system('ruby draw_xt_sfcheatflux.rb ' )
#  system('ruby draw_xt_velx_at500m.rb ' + timestep.to_s )
#  system('ruby draw_xt_temp_at500m.rb ' + timestep.to_s )
#  system('ruby draw_xt_watervapormixingratio_at500m.rb ' + timestep.to_s )
#  system('ruby draw_xt_watervapormixingratio_at1000m.rb ' + timestep.to_s )
#  system('ruby draw_xt_velz_at500m.rb ' + timestep.to_s )
#  system('ruby draw_tz_temp.rb ' )
#  system('ruby draw_tz_watervapormixingratio.rb ' )
  system('ruby draw_xt_rerativehumidity_at500m.rb ' + timestep.to_s )

  ##system("pstoimg -flip r270 -type gif dcl.ps")
#  system("convert -rotate 90 dcl.ps dcl_H2O-l-Rain_FallFluxAtLB_xt" + (timestep-1).to_s + ".gif" )
  system("convert -rotate 90 dcl.ps dcl.gif" )

##  降水量
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-l-Rain_FallFluxAtLB_xt%0#{3}d.gif", timestep )
##  海面からの蒸発量
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-g_SfcMassFlux_xt%0#{3}d.gif", timestep )
##  水平平均した水平速度
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_VelX_tz.gif" )
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_VelX_tz%0#{3}d.gif", timestep )
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_VelX_tz2.gif" )
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_VelX_tz-contour%0#{3}d.gif", timestep )
#  水平平均したエクスナー関数
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner_tz.gif" )
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner_tz2.gif" )
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_Exner_tz-contour.gif" )
##  水平平均した温度
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_Temp_tz.gif")
##  水平平均した水蒸気混合比
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_tz.gif" )
##  海面からの熱フラックス
#  com = sprintf( "mv dcl.gif figs/dcl_SfcHeatFlux_xt%0#{3}d.gif", timestep )
##  高度 500 m あたり (450 m) での水平風
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_VelX_xt-at500m%0#{3}d.gif", timestep )
##  高度 500 m あたり (450 m) での温度
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_Temp_xt-at500m%0#{3}d.gif", timestep )
##  高度 500 m あたり (450 m) での水蒸気混合比
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_xt-at500m%0#{3}d.gif", timestep )
##  高度 1000 m あたり (1050 m) での水蒸気混合比
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-gAll_xt-at1000m%0#{3}d.gif", timestep )
##  高度 5000 m あたり (4950 m) での鉛直流
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_VelZ_xt-at5000m%0#{3}d.gif", timestep )
##  高度 500 m あたり (450 m) での鉛直流
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_VelZ_xt-at500m%0#{3}d.gif", timestep )
##  高度 500 m あたり (450 m) での相対湿度
  com = sprintf( "mv dcl.gif zzz-figs/dcl_RevHumid_xt-at500m%0#{3}d.gif", timestep )
  system(com)

end