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

z2 = 48000
zkm = 48

# スイッチ
for val in ["DthetaDt", "DqDt", "All"]

  ## 時間積分した熱膨張項への寄与 (dthetadt のみ, エクスナー関数)
  if val == "DthetaDt"

    p "### exner-thermalexpansionterm dthetadt meanx ###"
    system('ruby makenetCDF_xzt_thermalexpansionterm-dthetadt-plussfcmassflux.rb ' + z2.to_s )

    system('ruby draw_z_exner-thermalexpansiontermdthetadtmeanx.rb' + ' no' + ' 1' )
    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_ExnerThermExpDthetaDt_z_1-5_z%0#{2}dkm.gif", zkm )
    system(com)

    system('ruby draw_z_exner-thermalexpansiontermdthetadtmeanx.rb' + ' yes' + ' 1' )
    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_ExnerThermExpDthetaDt_z_z=0-12km_1-5_z%0#{2}dkm.gif", zkm )
    system(com)

    ### All ###

    if z2 == 30000

    system('ruby draw_z_exner-thermalexpansiontermdthetadtmeanx.rb' + ' no' + ' all' )
    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_ExnerThermExpDthetaDt_z_1-5_all.gif" )
    system(com)

    system('ruby draw_z_exner-thermalexpansiontermdthetadtmeanx.rb' + ' yes' + ' all' )
    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_ExnerThermExpDthetaDt_z_z=0-12km_1-5_all.gif" )
    system(com)

    end
    p "### exner-thermalexpansionterm dthetadt meanx ###"

  ## 時間積分した熱膨張項への寄与 (dqdt のみ, エクスナー関数)
  elsif val == "DqDt"

    p "### exner-thermalexpansionterm dqdt meanx ###"
    system('ruby makenetCDF_xzt_thermalexpansionterm-dqdt-plussfcmassflux.rb ' + z2.to_s )

    system('ruby draw_z_exner-thermalexpansiontermdqdtmeanx.rb' + ' no' + ' 1' )
    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_ExnerThermExpDqDt_z_1-5_z%0#{2}dkm.gif", zkm )
    system(com)

    system('ruby draw_z_exner-thermalexpansiontermdqdtmeanx.rb' + ' yes' + ' 1' )
    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_ExnerThermExpDqDt_z_z=0-12km_1-5_z%0#{2}dkm.gif", zkm )
    system(com)

    ### All ###

    if z2 == 30000

      system('ruby draw_z_exner-thermalexpansiontermdqdtmeanx.rb' + ' no' + ' all' )
      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_ExnerThermExpDqDt_z_1-5_all.gif" )
      system(com)

      system('ruby draw_z_exner-thermalexpansiontermdqdtmeanx.rb' + ' yes' + ' all' )
      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_ExnerThermExpDqDt_z_z=0-12km_1-5_all.gif" )
      system(com)

    end
    p "### exner-thermalexpansionterm dqdt meanx ###"

  ## 時間積分した熱膨張項への寄与 (エクスナー関数)
  elsif val == "All"

    p "### exner-thermalexpansionterm meanx ###"
    system('ruby draw_z_exner-thermalexpansiontermmeanx.rb' + ' no' + ' 1' )
    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_ExnerThermExp_z_1-5_z%0#{2}dkm.gif", zkm )
    system(com)


    system('ruby draw_z_exner-thermalexpansiontermmeanx.rb' + ' yes' + ' 1' )
    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_ExnerThermExp_z_z=0-12km_1-5_z%0#{2}dkm.gif", zkm )
    system(com)

    ### All ###

    if z2 == 30000

      system('ruby draw_z_exner-thermalexpansiontermmeanx.rb' + ' no' + ' all' )
      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_ExnerThermExp_z_1-5_all.gif" )
      system(com)

      system('ruby draw_z_exner-thermalexpansiontermmeanx.rb' + ' yes' + ' all' )
      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_ExnerThermExp_z_z=0-12km_1-5_all.gif" )
      system(com)

    end
    p "### exner-thermalexpansionterm meanx ###"

  end

end