# -*- 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

file1 = 'odaka1998_1_VelZ.nc'
file2 = 'odaka1998_2_VelZ.nc'
file3 = 'odaka1998_3_VelZ.nc'
file4 = 'odaka1998_4_VelZ.nc'
file5 = 'odaka1998_5_VelZ.nc'
file6 = 'odaka1998_6_VelZ.nc'
file7 = 'odaka1998_7_VelZ.nc'
file8 = 'odaka1998_8_VelZ.nc'
file9 = 'odaka1998_9_VelZ.nc'
file10 = 'odaka1998_10_VelZ.nc'

varname = 'VelZ'

var  = GPhys::IO.open([file1,file2,file3,file4,file5,file6,file7,file8,file9,file10], varname)
time  = GPhys::IO.open([file1,file2,file3,file4,file5,file6,file7,file8,file9,file10], 't').val

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(4)

str = (0).to_s

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

for i in 1..5 do

time = 86400.0*i*2

str = (i*2).to_s

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

end

# <Close DCL>
DCL.grcls