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

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

var0  = GPhys::IO.open(file0, varname0)
var1  = GPhys::IO.open([file1,file2,file3,file4], varname1)
time  = GPhys::IO.open([file1,file2,file3,file4], 't').val

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

VarAll = var1.cut(true,true,true,0.0) + var0[5..104,0..-1,5..54]
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)    

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

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

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

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

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

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

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

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

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

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

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

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


# <Close DCL>
DCL.grcls