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

file = 'BS1998_PTempRad.nc'
dens0 = GPhys::IO.open("./BS1998_restart.nc",   'DensBZ'  ).val
#file = 'gphys.nc'

PTRad   = GPhys::IO.open(file,   'PTempRad'  )
ExnerBZ   = GPhys::IO.open("BS1998_restart.nc",   'ExnerBZ'  ).val
denscut = dens0[1,0,5..104]
exnercut = ExnerBZ[1,0,5..104]
#PTRad = PTRad * denscut *891.0 * exnercut

DCL.gropn(2)
GGraph.line( PTRad.cut('x'=>50.0,'y'=>100.0,'t'=>0), true, 'exchange'=>true)
DCL.grcls
