#!/usr/bin/env ruby

require "numru/ggraph"
include NumRu

varC   = "MixRtC"
gphysC = GPhys::IO.open( ARGV[0], varC ).sum('x').sum('z')

varR   = "MixRtR"
gphysR = GPhys::IO.open( ARGV[0], varR ).sum('x').sum('z') 

varV  = "MixRtV"
gphysV = GPhys::IO.open( ARGV[0], varV ).sum('x').sum('z') 

varDV  = "DelMixRtV"
gphysDV = GPhys::IO.open( ARGV[0], varDV ).sum('x').sum('z') 

varP   = "PotTemp"
gphysP = GPhys::IO.open( ARGV[0], varP ).sum('x').sum('z') 

varRV   = "RhoV"
gphysRV = GPhys::IO.open( ARGV[0], varRV ).sum('x').sum('z') 

gphys = gphysDV + gphysC + gphysR

DCL::swpset('IHEIGHT', 600 )
DCL::swpset('IWIDTH',  600 )

DCL.gropn( 2 )
DCL.sgpset('lfull',true) 
#DCL.sglset('lclip',true)
DCL.sgpset('lcntl', false) 

DCL.uzfact(0.7)
#GGraph.line( gphysRV.cut('t'=>0..3000), true, 'index'=>43)
GGraph.line( gphysRV, true, 'index'=>43, 'max'=>13.82, 'min'=>13.65 )
GGraph.line( gphysDV, true, 'index'=>43, 'max'=>3.0, 'min'=>-2.0)
GGraph.line( gphysC, false,'index'=>33)
GGraph.line( gphysR, false,'index'=>23)
GGraph.line( gphys,  true, 'index'=>43)
GGraph.line( gphysP, true, 'index'=>43)

DCL.grcls