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

Adv   = GPhys::IO.open('odaka1998_10_PTempAdv.nc',   'PTempAdv'  )
NDiff = GPhys::IO.open('odaka1998_10_PTempNDiff.nc', 'PTempNDiff')
Turb  = GPhys::IO.open('odaka1998_10_PTempTurb.nc',  'PTempTurb' )
Disp  = GPhys::IO.open('odaka1998_10_PTempDisp.nc',  'PTempDisp' )
Rad   = GPhys::IO.open('odaka1998_10_PTempRad.nc',   'PTempRad'  )
Sfc   = GPhys::IO.open('odaka1998_10_PTempSfc.nc',   'PTempSfc'  )

Turb  = Turb + Sfc

# x-y mean

Adv_XYmean  = Adv.average('x').average('y').convert_units('K/day')
Turb_XYmean = Turb.average('x').average('y').convert_units('K/day')
Rad_XYmean  = Rad.average('x').average('y').convert_units('K/day')

time  = GPhys::IO.open('odaka1998_10_PTempAdv.nc', 't').val

# time mean

#Adv_Tmean = Adv_XYmean[0..-1, 67..72].average('t')
#Turb_Tmean= Turb_XYmean[0..-1, 67..72].average('t')
#Rad_Tmean = Rad_XYmean[0..-1, 67..72].average('t')


# <Open DCL>
DCL.swlset( 'ldump', true )
DCL.gropn(4)

#for i in 0...time.length do
# <Drawing data by GGraph>
#  GGraph.line( var.cut(5500,100,true,time[i]) )

  GGraph.line( Adv_XYmean[0..-1,-1], \
               true, 'exchange'=>true , \
               'index'=>22, 'type'=>1, 'label'=>'adv', \
               'title'=>'Potential Temp. Tendency', \
               'max'=>100.0, 'min'=>-100.0)    # °Üή

  GGraph.line( Turb_XYmean[0..-1,-1], \
               false, 'exchange'=>true , \
                'index'=>32, 'type'=>1, 'label'=>'turb', \
               'max'=>100.0, 'min'=>-100.0)    # Íðή³È»¶

  GGraph.line( Rad_XYmean[0..-1,-1], \
               false, 'exchange'=>true , \
               'index'=>42, 'type'=>1, 'label'=>'rad', \
               'max'=>100.0, 'min'=>-100.0)    # Êü¼Í²ÃÇ®

#end
# <Close DCL>
DCL.grcls