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

file0 = 'odaka1998_restart.nc'
file1 = 'odaka1998_1_PTemp.nc'
file2 = 'odaka1998_2_PTemp.nc'
file3 = 'odaka1998_3_PTemp.nc'
file4 = 'odaka1998_4_PTemp.nc'
file5 = 'odaka1998_5_PTemp.nc'
file6 = 'odaka1998_6_PTemp.nc'
file7 = 'odaka1998_7_PTemp.nc'
file8 = 'odaka1998_8_PTemp.nc'
file9 = 'odaka1998_9_PTemp.nc'
file10 = 'odaka1998_10_PTemp.nc'

varname0 = 'PTempBZ'
varname1 = 'PTemp'
ltime = 172800.0

var0  = GPhys::IO.open(file0, varname0)
#var1  = GPhys::IO.open([file1,file2,file3,file4,file5,file6,file7,file8,file9,file10], varname1)
#time  = GPhys::IO.open([file1,file2,file3,file4,file5,file6,file7,file8,file9,file10], 't').val

var1  = GPhys::IO.open(file1, varname1)
time  = GPhys::IO.open(file1, 't').val

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

VarAll = var1.cut(true,true,true,0.0) + var0[5..504,0..-1,5..104]
VarMean = VarAll.average('x').average('y')

GGraph.line( VarMean, 
             true, 'exchange'=>true ,
             'index'=>2, 'type'=>1, 'label'=>'00',
             'title'=>'Mean Potemtial Temp.',
             'max'=>280.0, 'min'=>220.0)    

for i in 1..6 do

time = 3600.0*i*2
#time = 86400.0*i*2

VarAll = var1.cut(true,true,true,time) + var0[5..504,0..-1,5..104]
VarMean = VarAll.mean(0,1)

str = (i*2).to_s

GGraph.line( VarMean, 
               false, 'exchange'=>true ,
               'index'=>2, 'type'=>1, 'label'=>str,
               'max'=>280.0, 'min'=>220.0)    
end

# <Close DCL>
DCL.grcls