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

z2 = 48000
zkm = 48
#t2 = 86400 * 8
#t2 = 86400 * 50
t2 = 86400 * 20

# スイッチ
#for val in ["Exner", "H2O-gAll", "Sekisan"]
for val in ["Exner", "Press"]
#for val in ["Sekisan"]
#for val in ["H2O-gAll"]

##### t 図 #####

  p val

  ##  地表面でのエクスナー関数
  if val == "Exner"

#    system('ruby draw_t_exner_atSfc.rb' )
    system('ruby draw_t_exner_atSfc.rb ' + t2.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/zzzNew/dcl_Exner_t_atSfc_1-2_z48km.gif" )
#    com = sprintf( "mv dcl.gif zzz-figs/zzzNew2/dcl_Exner_t_atSfc_1-2_z48km_20day.gif" )
    com = sprintf( "mv dcl.gif zzz-figs/zzzNew2/dcl_Exner_t_atSfc_1-2-3_z48km_20day.gif" )
#    com = sprintf( "mv dcl.gif zzz-figs/zzzNew/dcl_Exner_t_atSfc_hikaku.gif" )
#    com = sprintf( "mv dcl.gif zzz-figs/zzzNew2/dcl_Exner_t_atSfc_hikakau_10day.gif" )
#    com = sprintf( "mv dcl.gif zzz-figs/zzzNew2/dcl_Exner_t_atSfc_hikakau_20day.gif" )
    system(com)

  ##  地表面気圧
  elsif val == "Press"

#    system('ruby draw_t_press_atSfc.rb' )
    system('ruby draw_t_press_atSfc.rb ' + t2.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/zzzNew/dcl_Press_t_atSfc_1-2_z48km.gif" )
#    com = sprintf( "mv dcl.gif zzz-figs/zzzNew2/dcl_Press_t_atSfc_1-2_z48km_20day.gif" )
    com = sprintf( "mv dcl.gif zzz-figs/zzzNew2/dcl_Press_t_atSfc_1-2-3_z48km.gif" )
#    com = sprintf( "mv dcl.gif zzz-figs/zzzNew/dcl_Press_t_atSfc_1-4-2_z48km.gif" )
#    com = sprintf( "mv dcl.gif zzz-figs/zzzNew/dcl_Press_t_atSfc_hikaku.gif" )
#    com = sprintf( "mv dcl.gif zzz-figs/zzzNew2/dcl_Press_t_atSfc_hikaku_10day.gif" )
#    com = sprintf( "mv dcl.gif zzz-figs/zzzNew2/dcl_Press_t_atSfc_hikaku_20day.gif" )
    system(com)


  ##  地表面での水蒸気混合比
  elsif val == "H2O-gAll"

    system('ruby makenetCDF_watervapormixingratio_atSfc.rb' )
    system('ruby draw_t_watervapormixingratio_atSfc.rb' )
    system('ruby psfat.rb'  + ' dcl.ps > dclnew.ps' )
    system("convert -rotate 90 dclnew.ps dcl.gif" )
    com = sprintf( "mv dcl.gif zzz-figs/zzzNew/dcl_H2O-gAll_t_atSfc_1-5_all.gif" )
    system(com)

  ## 積算値の時間変化とか
  elsif val == "Sekisan"

  ## 水平平均水蒸気混合比の総量の時間変化
    system('ruby makenetCDF_t_TQwatervapormeanx.rb ' + z2.to_s )
    system('ruby draw_t_budgetofwatervapormeanx.rb' )
    system('ruby psfat.rb'  + ' dcl.ps > dclnew.ps' )
    system("convert -rotate 90 dclnew.ps dcl.gif" )
    com = sprintf( "mv dcl.gif zzz-figs/zzzNew/dcl_Budget-H2O-gAll_t_1-5_z%0#{2}dkm.gif", zkm )
    system(com)

  ## 水平平均降水量の積算値の時間変化
    system('ruby makenetCDF_t_sekisan-rainfallfluxatLB.rb ' + z2.to_s )
    system('ruby draw_t_sekisan-rainfallfluxatLB.rb' )
    system('ruby psfat.rb'  + ' dcl.ps > dclnew.ps' )
    system("convert -rotate 90 dclnew.ps dcl.gif" )
    com = sprintf( "mv dcl.gif zzz-figs/zzzNew/dcl_Sekisan-H2O-l-Rain_FallFluxAtLB_t_1-5_z%0#{2}dkm.gif", zkm )
    system(com)

  ## 水平平均蒸発量の積算値の時間変化
    system('ruby makenetCDF_t_sekisan-watervaporsfcmassflux.rb ' + z2.to_s )
    system('ruby draw_t_sekisan-watervaporsfcmassflux.rb' )
    system('ruby psfat.rb'  + ' dcl.ps > dclnew.ps' )
    system("convert -rotate 90 dclnew.ps dcl.gif" )
    com = sprintf( "mv dcl.gif zzz-figs/zzzNew/dcl_Sekisan-H2O-g_SfcMassFlux_t_1-5_z%0#{2}dkm.gif", zkm )
    system(com)

  end

  p "--------------------"

end



##  6 時間ごとに積算した水平平均降水量
#  system('ruby draw_t_int-rainfallfluxatLB.rb' )
#  system("convert -rotate 90 dcl.ps dcl.gif" )
#  com = sprintf( "mv dcl.gif zzz-figs/zzzNew/dcl_Int-H2O-l-Rain_FallFluxAtLB_t.gif" )
#  system(com)

##  水平平均した蒸発量
#  system('ruby draw_t_watervapormixingratio_sfcmassflux.rb' )
#  system('ruby draw_t_watervapormixingratio_sfcmassflux.rb ' + timestep.to_s )
#  system("convert -rotate 90 dcl.ps dcl.gif" )
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-g_SfcMassFlux_t.gif" )
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-g_SfcMassFlux_t_1-5_x256km_z48km_%0#{8}d.gif", timestep*dt )
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-g_SfcMassFlux_t_1-5_x256km_z48km_%0#{8}d.gif", timestep*dt )
#  system(com)

##  水平平均した降水量
#  system('ruby draw_t_rainfallfluxatLB.rb' )
#  system('ruby draw_t_rainfallfluxatLB.rb ' + timestep.to_s )
#  system("convert -rotate 90 dcl.ps dcl.gif" )
##  com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-l-Rain_FallFluxAtLB_t.gif" )
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-l-Rain_FallFluxAtLB_t_1-5_x256km_z48km_%0#{8}d.gif", timestep*dt )
#  com = sprintf( "mv dcl.gif zzz-figs/dcl_H2O-l-Rain_FallFluxAtLB_t_1-5_x256km_z48km_%0#{8}d.gif", timestep*dt )
#  system(com)

#end