# -*- coding: euc-jp -*-
#
# (dcpam 用) ある年における惑星表面の物理量の東西平均時間変化を描画する ctl スクリプト
#
# usage: draw_ctl_ts.rb を実行する.
#
# * 2012/01/05 drawzm_clim_1panel_ctl.rb を基に作成
#

require "numru/ggraph"
include NumRu

#data = ARGV[1].to_s

#if data == "dcpam"  then

   dir = '..'
   vname = ARGV[0]
   exp = ARGV[1]

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

   iyrs = 10
   iyre = 10

=begin
elsif data == "NCEP" then

   dir = '/home2/yukai-home2/NCEP'
   vname = ARGV[0]

   p vname

   gphys = GPhys::NetCDF_IO.open( dir + '/' +  'ulwrftop.mon.mean' + '.nc', 'ulwrf')
   tm = gphys.shape[0] # number of elements for time      dimension

   iyrs = 1988
   iyre = 2007
   
end   
=end

p vname, exp
p iyrs, iyre

figdir = 'figs_zm_ts'

if FileTest::directory?(figdir) == false then  

   Dir::mkdir( figdir )
   Dir::mkdir( figdir+'/eps' )

end

    system('ruby drawzm_ts.rb ' + dir + ' ' + vname + ' ' + iyrs.to_s + ' ' + iyre.to_s + ' ' + exp )

# png ファイルの作成

    com = sprintf( "convert -rotate 90 dcl.ps " +vname+ "_zm_ts.png" )
    system(com)

# eps ファイルの作成

    com = sprintf( "convert -rotate 90 dcl.ps " +vname+ "_zm_ts.eps" )
    system(com)

    com = sprintf( "mv *.png "+figdir )
    system(com)
    com = sprintf( "mv *.eps  "+figdir+"/eps" )
    system(com)