# -*- coding: euc-jp -*-
require "numru/ggraph"
include NumRu

#file0 = 'thermal-moist_restart_rank000000.nc'
#var0  = 'PTempBZ'
#file  = 'thermal-moist_PTemp_rank000000.nc'
#var   = 'PTemp'
var   = 'PTemp_H2O-g'

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

dt = 3600      # (1 hour)          # 絵を出力する時間刻み幅
tt = 4320000   # (50 days)         # 全積分時間
#tt = 36000   # (10 hours)          # 全積分時間
tnmin = 0                          # 繰り返し回数の最小値 (固定)
tnmax = tt/dt  # (1200)            # 繰り返し回数の最大値

####--- ファイル名の通し番号を入れておく配列を作る ---###
inum=[0]
for tn in tnmin..tnmax 
  inum[tn] = tn*dt
end


###--- 繰り返しはじめ ---###
##for t in [3600, 4200]
#for t in timearray
##for t in timearray[6*24*3+1-1..timearray.size-1]

#for tn in tnmin..tnmax 
for tn in 1011..tnmax 

#  p t
  p tn

#  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_core2_sub.rb ' + t.to_s )
  system('ruby draw_core2_sub.rb '  + tn.to_s + ' ' + dt.to_s + ' ')

  ###--- ファイルの rename & png を gif に変換 ---###
    ###--- ファイル名を変数に入れておく ---###
#  filename = printf("dcl_" + var + "_t%0#{7}d.gif\n", inum[tn] )
#  p filename

#  com = sprintf( "convert dcl_001.png dcl_t%0#{7}d.gif", inum[tn] )
#  com = sprintf( "convert dcl_001.png" + filename)
  com = sprintf( "convert dcl_001.png dcl_" + var + "_t%0#{7}d.gif\n", inum[tn] )
  system(com)
end

###--- アニメーションを作る ---###
#filename0 = printf("dcl_" + var + "_anim.gif")
com = sprintf("gifsicle --delay=20 --loopcount=1 dcl_" + var + "_t*.gif > dcl_" + var + "_anim.gif\n")
system(com)