# -*- 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"
cuttime = 180000..280000
#file0 = 'BS1998_restart_rank000000.nc'
#file_base = '../test_restart_rank000000.nc'

file_base = 'restart_long-DensBZ.nc'
ptemp_file = "../#{filedir}/#{filehead}_rank000000.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"

comax =0.0
comin = -0.030
#densBZ_a  = GPhys::IO.open(file_base, 'DensBZ')
densBZ  = GPhys::IO.open(file_base, 'DensBZ')
ptemp  = GPhys::IO.open([file0,file1,file2,file3,file4,file5], 'PTemp')
#ptemp  = GPhys::IO.open(ptemp_file, 'PTemp')
#densBZ = densBZ_a[5..204,0..-1,5..104]
#var0_cut = var0.cut('t'=>cuttime)

rho_theta = ptemp * densBZ

rho_theta_XYZmean = rho_theta.mean('x').mean('y').mean('z')

#p rho_theta_XYZmean.cut('t'=>40000).val.to_f
#p rho_theta_XYZmean.cut('t'=>80000).val.to_f

#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
#deriv = (rho_theta_XYZmean.cut('t'=>63000).val.to_f - rho_theta_XYZmean.cut('t'=>66000).val.to_f) / 3.0e3
#p deriv
GGraph.line( rho_theta_XYZmean.cut('t'=>cuttime), 
             true, 'exchange'=>false ,
             'index'=>2, 'type'=>1, #'label'=>'10',
             'title'=>'rho theta',
             'max'=>comax, 'min'=>comin)    

DCL.grcls