require "numru/ggraph"
include NumRu

dir = '..'

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

num = 1 * 10
#num = 4 * 10

DataNumOfDay = 4
DaysOfMonth  = [31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31]

DaysOfYear = 0
DaysOfMonth.each do |i|
  DaysOfYear += i
end
p DaysOfYear

figdir = 'figs_zm_log_mon'
Dir::mkdir( figdir )

for iyr in 20..20
  for imon in 1..12
    t = DaysOfYear*(iyr-1)
    for i in 1..imon-1
      t = t + DaysOfMonth[i-1]
    end

    ts = t
    te = t + DaysOfMonth[imon-1]

    ts = ts * DataNumOfDay
    te = te * DataNumOfDay - 1

    if ts > tm then
      p 'ts is in appropriate'
    end
    if te > tm then
      p 'te is in appropriate'
    end

    p 'Year: ' + iyr.to_s + ', Month: ' + imon.to_s + ', t = ' + ts.to_s + ' to '+ te.to_s

    case imon
    when 1
      title = 'JAN.'
    when 2
      title = 'FEB.'
    when 3
      title = 'MAR.'
    when 4
      title = 'APR.'
    when 5
      title = 'MAY'
    when 6
      title = 'JUN.'
    when 7
      title = 'JUL.'
    when 8
      title = 'AUG'
    when 9
      title = 'SEP.'
    when 10
      title = 'OCT.'
    when 11
      title = 'NOV.'
    when 12
      title = 'DEC.'
    end

    system('ruby drawzm_log.rb ' + dir + ' ' + ts.to_s + ' ' + te.to_s + ' ' + title)

#    system("pstoimg -flip r270 dcl.ps")
#    com = sprintf( "mv dcl.png " + figdir + "/dcl_yr%0#{3}d_mon%0#{2}d.png", iyr, imon )
    system("pstoimg -type gif -flip r270 dcl.ps")
    com = sprintf( "mv dcl.gif " + figdir + "/dcl_yr%0#{3}d_mon%0#{2}d.gif", iyr, imon )
    system(com)
    com = sprintf( "mv dcl.ps " + figdir + "/dcl_yr%0#{3}d_mon%0#{2}d.ps", iyr, imon )
    system(com)
  end
end