# -*- coding: euc-jp -*-
# 
# NCEP/NCAR 再解析データの 1988-2007 年における東西平均した各緯度における
# OLR, OSR の値を出力する
#
# usage: 出力したい変数に合わせてコメントアウトを追加/削除してから,
#        このスクリプトを実行する
#
# 2012/01/10 drawgm.rb を元に作成
#


require "numru/ggraph"
include NumRu


dir    = '/home2/yukai-home2/NCEP'
vname1 = 'sst_T62'
vname2 = 'ulwrf'
vname3 = 'uswrf'
vname4 = 'dswrf'

file2 = 'ulwrftop.mon.mean'
file3 = 'uswrftop.mon.mean'
file4 = 'dswrftop.mon.mean'


#
# NECP の データは 1988-2007 年度の各月平均を使用する
#
ts2 = 17417760
te2 = 17592240

#
# 各 nc ファイルの読み込み
#
gphys2 = GPhys::NetCDF_IO.open(dir+'/'+file2+".nc", vname2)
gphys3 = GPhys::NetCDF_IO.open(dir+'/'+file3+".nc", vname3)
gphys4 = GPhys::NetCDF_IO.open(dir+'/'+file4+".nc", vname4)

gphys5 = gphys4 - gphys3 # 射出太陽放射
gphys5.name      = 'outgoing shortwave'
gphys5.long_name = 'outgoing shortwave'

#
# 変数 time の定義
#
time = gphys2.coord('time').val

lat  = gphys2.coord('lat').val

for i in 0..(lat.size-1)

  #p i

  gphystmp2 = ( gphys2.cut('lat'=>lat[i]).mean('lon').cut('time'=>ts2..te2).mean('time') )
  gphystmp5 = ( gphys5.cut('lat'=>lat[i]).mean('lon').cut('time'=>ts2..te2).mean('time') )

#  p 'lat' + ': ' + lat[i].to_s + '  OLR' + ': ' + gphystmp2.to_s
  p gphystmp2.to_s

#  p 'lat' + ': ' + lat[i].to_s + '  OSR' + ': ' + gphystmp5.to_s
  p gphystmp5.to_s


end