# -*- 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_000000000-000172000"
filedir = "time_000000000-000080000"
filehead = "BS1998_PTemp"
#file0 = 'BS1998_restart_rank000000.nc'
file_base = 'restart_long-PTempBZ.nc'
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"

varname_conv = 'PTemp'
varunit_conv = 'K'
varname0 = 'PTempBZ'
varname1 = 'PTemp'
#cuttime  = 
#comax = 249.0
#comin = 247.0
comax = 0.0
comin = -1.0
var_base  = GPhys::IO.open(file_base, varname0)
var0  = GPhys::IO.open([file0,file1,file2,file3,file4,file5], varname1)

#var0_cut = var0.cut('t'=>cuttime)



#time  = GPhys::IO.open(file1, 't').val
#x_co  = GPhys::IO.open(file1, 'x').val

# <Open DCL>
DCL.gropn(2)
#var0 = var0[5..104,0..-1,5..104]

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

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

for i in 1..4 do

time = 20000.0*i

VarAll = var0.cut(true,true,true,time)
VarMean = VarAll.mean(0,1)

str = (i).to_s

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

# <Close DCL>
DCL.grcls