# -*- 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_000000000-000150000"

filehead = "BS1998_"
cuttime = 100000..150000
#file0 = 'BS1998_restart_rank000000.nc'
file_base = 'restart_long-DensBZ.nc'
ptemp_file = "../#{filedir}/#{filehead}PTemp_rank000000.nc"
ptemprad_file = "../#{filedir}/#{filehead}PTempRad_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"

comax =0.000010
comin = -0.00010
densBZ  = GPhys::IO.open(file_base, 'DensBZ')
ptemp  = GPhys::IO.open(ptemp_file, 'PTemp')
ptemprad  = GPhys::IO.open(ptemprad_file, 'PTempRad')

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

rhoQ = densBZ * ptemprad
rhoQ_XYZmean = rhoQ.mean('x').mean('y').mean('z')

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

# <Open DCL>
DCL.gropn(4)
#var0 = var0[5..104,0..-1,5..104]
p rhoQ_XYZmean.cut('t'=>100000).val.to_f
p rhoQ_XYZmean.cut('t'=>150000).val.to_f
#VarAll = var0.cut(true,true,true,0.0) + var_base

GGraph.line( rhoQ_XYZmean.cut('t'=>cuttime), 
             true, 'exchange'=>false ,
             'index'=>2, 'type'=>1, #'label'=>'10',
             'title'=>'Mean Potemtial Temp.',
             'max'=>comax, 'min'=>comin)    

DCL.grcls