# Title: Ruby script drawing contour map for deepconv/arare5 output data 
#
# History: 2011/09/27 (Masatsugu Odaka)
#
require "numru/ggraph"
include NumRu

file0 = 'odaka1998_restart.nc'
file1 = 'odaka1998_10_PTemp.nc'


varname0 = 'PTempBZ'
varname1 = 'PTemp'
ltime = 864000.0

var0  = GPhys::IO.open(file0, varname0)
var1  = GPhys::IO.open(file1, varname1)
time  = GPhys::IO.open(file1, 't').val

GGraph.set_linear_contour_options( 'min'=>-5.0, 'max'=>5.0, 'interval'=>1.0 )  
GGraph.set_linear_tone_options( 'min'=>-5.0, 'max'=>5.0, 'interval'=>0.1 )  


# <Open DCL>
DCL.swlset( 'ldump', true )
DCL.gropn(4)
#for i in 0...time.length do
# <Drawing data by GGraph>

  VarAll = var1.cut(true,true,true,ltime) + var0[5..404,0..-1,5..54]
  VarDev = VarAll.eddy(0)

  GGraph.set_fig('viewport'=>[0.1,0.8,0.5,0.7])
  DCL.udpset('lmsg', false)
  GGraph.set_axes('xunit'=>'m', 'yunit'=>'m', 'xtitle'=>'')
  GGraph.tone( VarDev.cut(0..40000.0,100.0,true) )
  GGraph.contour( VarDev.cut(0..40000.0, 100.0, true), false )

  GGraph.set_fig('viewport'=>[0.1,0.8,0.2,0.4], 'new_frame'=>false)
  DCL.udpset('lmsg', true)
  GGraph.set_axes('xunit'=>'m', 'yunit'=>'m', 'xtitle'=>'')
  GGraph.tone( VarDev.cut(40000.0..80000.0,100.0,true), true, 'title'=>'', 'annot'=> false)
  GGraph.contour( VarDev.cut(40000.0..80000.0, 100.0, true), false, 'title'=>'', 'annot'=>false)

  GGraph.color_bar( 'tickintv'=>5 )
#end

# <Close DCL>
DCL.grcls