#
# Test program for GGraph, implemented rectangular curvelinear coordinate. 
# (2004/01/25 S. Takehiro)
#
#  Drawing initial temperature field for dynamo benchmark test. 
#
#require "numru/ggraph"                                                      
require "./ggraph"
include NumRu

gphys = GPhys::NetCDF_IO.open('case1_initial.nc', 'temp')


# 子午面断面
gp = gphys.cut('lon'=>0.0,'t'=>0.0 )
DCL.gropn(1)
GGraph.set_fig('itr'=>5 )
GGraph.tone( gp, true, {'transpose'=>true} )
GGraph.contour( gp, false, {'transpose'=>true} )

# 赤道面断面
gp = gphys.cut('lat'=>0.0,'t'=>0.0 )
GGraph.set_fig('itr'=>5 )
GGraph.tone( gp, true, {'transpose'=>true} )
GGraph.contour( gp, false, {'transpose'=>true} )

# 赤道面断面(ちょっと大きめ)
GGraph.set_fig('itr'=>5 )
GGraph.set_fig('similar'=>[0.25,-0.1,0.05] )
GGraph.tone( gp, true, {'transpose'=>true} )
GGraph.contour( gp, false, {'transpose'=>true} )

# 正射投影
gp = gphys.cut('rad'=>1.0,'t'=>0.0 )
GGraph.set_fig('itr'=>30 )
GGraph.set_fig('map_axis'=>[0,30,20] )
GGraph.tone( gp, true )
GGraph.contour( gp, false )

DCL.grcls