#!/usr/bin/env ruby

require "numru/ggraph"      # ggraph library
include NumRu

# file

file_d = "./nc/arare-odaka1998-dx200-20050918.nc"  # output netCDF file (disturb field value)
file_b = "./nc/arare-odaka1998-dx200-20050918-init.nc"  # output netCDF file (basic field value)

# definition of constant value

Sfctemp = 270.0 # ground tempareture [K]
Min = 240       # minimum value of draw range
Max = 280       # maximum value of draw range
Xmin = 0       # minimum value of x-region
Xmax = 25600   # maximum value of x-region
Zmin = 0       # minimum value of z-region
Zmax = 10000   # maximum value of z-region

# GPhys object initialization

gp_pottemp_d = GPhys::IO.open(file_d, "PotTemp")
gp_pottemp_b = GPhys::IO.open(file_b, "PotTempBasicZ")

# uniform array size and dimension

gp_pottemp_b = gp_pottemp_b.cut('x'=>Xmin..Xmax, 'z'=>Zmin..Zmax)
gp_pottemp_d = gp_pottemp_d.cut('x'=>Xmin..Xmax, 'z'=>Zmin..Zmax)

# draw

DCL.gropn(4)
DCL.sgpset('lcntl', false)
DCL.uzfact(0.7)

gp_pottemp = gp_pottemp_b + gp_pottemp_d.cut('t'=>0)
GGraph.tone( gp_pottemp, true)
GGraph.contour( gp_pottemp, false)

DCL.grcls