# -*- 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_000050000-000100000"
filedir = "time_000180000-000280000"
filehead = "BS1998_PTemp"
#restart_filehead = "test_restart"
#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"
#restart_file0 = "../#{restart_filehead}_rank000000.nc"
#restart_file1 = "../#{restart_filehead}_rank000001.nc"
#restart_file2 = "../#{restart_filehead}_rank000002.nc"
#restart_file3 = "../#{restart_filehead}_rank000003.nc"
#restart_file4 = "../#{restart_filehead}_rank000004.nc"
#restart_file5 = "../#{restart_filehead}_rank000005.nc"

varname_conv = 'PTemp'
varunit_conv = 'K'
varname0 = 'PTempBZ'
varname1 = 'PTemp'
#cuttime  = 
comax = 250.50
comin = 240.50
var_base  = GPhys::IO.open(file_base, varname0)

var0  = GPhys::IO.open([file0,file1,file2,file3,file4,file5], varname1)
#var_base  = GPhys::IO.open([restart_file0,restart_file1,restart_file2,restart_file3,restart_file4,restart_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) + var_base
VarAll = var0.cut(true,true,true,180000) + var_base
#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'=>'0',
             'title'=>'PTemp',
             'max'=>comax, 'min'=>comin)    
#=begin
for i in 1..4 do

time = 3600.0*i*6 + 180000

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

str = (i*6).to_s

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