# -*- coding: euc-jp -*-
# Title: Ruby script drawing contour map for deepconv/arare5 output data 
#
# History: 2011/09/27 (Masatsugu Odaka)
#
require "numru/ggraph"
include NumRu

filedir = "time_000180000-000280000"
filehead = "BS1998_VelZ"
file0 = "../#{filedir}/#{filehead}_rank000000.nc"
file1 = "../#{filedir}/#{filehead}_rank000001.nc"
file2 = "../#{filedir}/#{filehead}_rank000002.nc"
file3 = "../#{filedir}/#{filehead}_rank000003.nc"
file4 = "../#{filedir}/#{filehead}_rank000004.nc"
file5 = "../#{filedir}/#{filehead}_rank000005.nc"
varname1 = 'VelZ'
var  = GPhys::IO.open([file0,file1,file2,file3,file4,file5], varname1)
time  = GPhys::IO.open(file1, 't').val
cuttime = 180000
zero_time = 180000

var_square = var**2
var_square_root = var_square**0.5
var_square_root_mean = var_square_root.average('x').average('y')

# <Open DCL>
DCL.swlset( 'ldump', true )
DCL.gropn(2)

str = (0).to_s

GGraph.line( var_square_root_mean.cut(true,zero_time), 
             true, 'exchange'=>true ,
             'index'=>2, 'type'=>1, 'label'=>str,
             'max'=>5.0, 'min'=>0.0)    

for i in 1..4 do

time = 3600.0*i*6 + zero_time

str = (i*6).to_s

GGraph.line( var_square_root_mean.cut(true,time), 
             false, 'exchange'=>true ,
             'index'=>2, 'type'=>1, 'label'=>str,
             'max'=>5.0, 'min'=>0.0)    

end

# <Close DCL>
DCL.grcls