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

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

=end

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

include NumRu

file = "/home/kitamo/arare3/arare_diag-odaka1998-dx200-20050917.nc" 
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