#!/usr/bin/env ruby
=begin
=draw_mass.rb

  * 領域全体の質量の時間変化を描く
  * 初期の値に対する比を描いている

==USAGE

 draw_mass.rb file

 * ARGUMENTS
  * file: 出力 NetCDF ファイル名

=end

require "numru/ggraph"
require "numru/dcl"

include NumRu

file = ["./nc/arare_diag-odaka1998-dx100-20051130.nc",
        "./nc/arare_diag-odaka1998-dx100-20051130-2.nc" ]

#file = ARGV[0]
Inittime = 0.0

# open file and initialize GPhys object

gp_mass = GPhys::IO.open(file, "MassTotal")

# disturb from value at initial time

gp_mass_ratio = ( gp_mass ) / gp_mass.cut("t"=>Inittime)

DCL.gropn(4)
DCL.sgpset("lcntl", false)
DCL.uzfact(0.7)
GGraph.line(gp_mass_ratio, true, 
              "index"=>13, 
              "title"=>"Mass Change from t=0")
DCL.grcls