require "numru/ggraph"
include NumRu

# gpview OLR.nc@OLR --mean=lon,lat SLR.nc@SLR --mean=lon,lat Rain.nc@Rain --mean=lon,lat EVap.nc@EVap --mean=lon,lat Sens.nc@Sens --mean=lon,lat --overplot=5 --range=0:500

varname    = 'OLR'
filename   = varname + '.nc'
gphys      = GPhys::IO.open(filename, varname)
lat_weight = GPhys::IO.open(filename, "lat_weight")

gphys = gphys.mean('lon')
#gphys = (gphys * lat_weight.to_a).sum('lat') / 2.0

summation = gphys


varname    = 'OSR'
filename   = varname + '.nc'
gphys      = GPhys::IO.open(filename, varname)
lat_weight = GPhys::IO.open(filename, "lat_weight")

gphys = gphys.mean('lon')
#gphys = (gphys * lat_weight.to_a).sum('lat') / 2.0

summation = summation + gphys


varname    = 'SLR'
filename   = varname + '.nc'
gphys      = GPhys::IO.open(filename, varname)
lat_weight = GPhys::IO.open(filename, "lat_weight")

gphys = gphys.mean('lon')
#gphys = (gphys * lat_weight.to_a).sum('lat') / 2.0

summation = summation - gphys


varname    = 'SSR'
filename   = varname + '.nc'
gphys      = GPhys::IO.open(filename, varname)
lat_weight = GPhys::IO.open(filename, "lat_weight")

gphys = gphys.mean('lon')
#gphys = (gphys * lat_weight.to_a).sum('lat') / 2.0

summation = summation - gphys


varname    = 'Rain'
filename   = varname + '.nc'
gphys      = GPhys::IO.open(filename, varname)
lat_weight = GPhys::IO.open(filename, "lat_weight")

gphys = gphys.mean('lon')
#gphys = (gphys * lat_weight.to_a).sum('lat') / 2.0

summation = summation - gphys


varname    = 'EVap'
filename   = varname + '.nc'
gphys      = GPhys::IO.open(filename, varname)
lat_weight = GPhys::IO.open(filename, "lat_weight")

gphys = gphys.mean('lon')
#gphys = (gphys * lat_weight.to_a).sum('lat') / 2.0

#summation = summation - gphys


varname    = 'Sens'
filename   = varname + '.nc'
gphys      = GPhys::IO.open(filename, varname)
lat_weight = GPhys::IO.open(filename, "lat_weight")

gphys = gphys.mean('lon')
#gphys = (gphys * lat_weight.to_a).sum('lat') / 2.0

summation = summation - gphys








summation = (summation * lat_weight.to_a).sum('lat') / 2.0
#summation = summation.mean('lat') / 2.0

summation = summation.cut('time'=>1000..2000)
#summation = summation.mean('time')

p summation.to_a

#__END__


DCL.gropn(1)

GGraph.set_fig('viewport'=>[0.15,0.80,0.15,0.6])

DCL.sgpset('lfull',true)     # 全画面表示
DCL.sgpset('lcntl', false) ; DCL.uzfact(0.6)

#GGraph.set_linear_contour_options('interval'=>4)

GGraph.line( summation )
#GGraph.contour(summation)

#GGraph.line(lat_weight.sum)


DCL.grcls