require "numru/ggraph"
include NumRu


itr = 3
x1 = -7; x2 = 0; y1 = -32 ; y2=-24
DCL.gropn(1)
DCL.sgpset('lcntl', false) ; DCL.uzfact(0.7)
#GGraph.set_fig 'itr'=>itr, 'viewport'=>[svx1,svx2,svy1,svy2], 'window'=>[x1,x2,y1,y2]
#GGraph.set_fig 'itr'=>itr, 'window'=>[x1,x2,y1,y2]
GGraph.set_fig 'itr'=>itr

#gpview --overplot=5 --overplot_color --var=BinAbsCoef,MolNum=1,Temp=300,Press=1e5,VMRMol2=100e-6,VMRMol3=0,BinWaveNum=6050e2 kdactbl-2017-05-24-2.nc ../prog08_calc_rte_check_tbl-x2/kdactbl-2017-05-24-2-04-noO3-x2.nc ../prog08_calc_rte_check_tbl-x4/kdactbl-2017-05-24-2-02-noO3-x4.nc ../prog08_calc_rte_check_tbl-x8/kdactbl-tmp.nc
ncfn = '../prog08_calc_rte_check/kdactbl-2017-05-24-2.nc'
#ncfn = '../prog08_calc_rte_check_tbl-x2/kdactbl-2017-05-24-2-04-noO3-x2.nc'
#ncfn = '../prog08_calc_rte_check_tbl-x2/kdactbl-2017-05-24-2-04-noO3-x2.nc'
gp = GPhys::IO.open(ncfn, 'BinAbsCoef')
#gp = gp.cut('MolNum'=>1).cut('Temp'=>300).cut('Press'=>1e5).cut('VMRMol2'=>100e-6).cut('VMRMol3'=>0).cut('BinWaveNum'=>6050e2).cut('VMRMol1'=>0.9e-7..1)
#gp = gp.cut('MolNum'=>1).cut('Press'=>1e5).cut('VMRMol2'=>100e-6).cut('VMRMol3'=>0).cut('BinWaveNum'=>6050e2).cut('VMRMol1'=>1e-3)
gp = gp.cut('MolNum'=>1).cut('Temp'=>300).cut('Press'=>0.4e-4..1e6).cut('VMRMol2'=>100e-6).cut('VMRMol3'=>0).cut('BinWaveNum'=>6050e2).cut('VMRMol1'=>1e-3)
gp = gp.log10
#vmr = gp.coord('VMRMol1').val
#vmr = NMath::log10(vmr)
#va = VArray.new( vmr, { "long_name"=>'VMRMol1', "units"=>'1' }, "VMRMol1" )
#gp.axis('VMRMol1').set_pos(va)
GGraph.line( gp )

ncfn = '../prog08_calc_rte_check_tbl-x2/kdactbl-2017-05-24-2-04-noO3-x2.nc'
gp = GPhys::IO.open(ncfn, 'BinAbsCoef')
#gp = gp.cut('MolNum'=>1).cut('Temp'=>300).cut('Press'=>1e5).cut('VMRMol2'=>100e-6).cut('VMRMol3'=>0).cut('BinWaveNum'=>6050e2).cut('VMRMol1'=>0.9e-7..1)
#gp = gp.cut('MolNum'=>1).cut('Press'=>1e5).cut('VMRMol2'=>100e-6).cut('VMRMol3'=>0).cut('BinWaveNum'=>6050e2).cut('VMRMol1'=>1e-3)
gp = gp.cut('MolNum'=>1).cut('Temp'=>300).cut('Press'=>0.4e-4..1e6).cut('VMRMol2'=>100e-6).cut('VMRMol3'=>0).cut('BinWaveNum'=>6050e2).cut('VMRMol1'=>1e-3)
gp = gp.log10
GGraph.line( gp, false, 'index'=>20 )

ncfn = '../prog08_calc_rte_check_tbl-x4/kdactbl-2017-05-24-2-02-noO3-x4.nc'
gp = GPhys::IO.open(ncfn, 'BinAbsCoef')
#gp = gp.cut('MolNum'=>1).cut('Temp'=>300).cut('Press'=>1e5).cut('VMRMol2'=>100e-6).cut('VMRMol3'=>0).cut('BinWaveNum'=>6050e2).cut('VMRMol1'=>0.9e-7..1)
#gp = gp.cut('MolNum'=>1).cut('Press'=>1e5).cut('VMRMol2'=>100e-6).cut('VMRMol3'=>0).cut('BinWaveNum'=>6050e2).cut('VMRMol1'=>1e-3)
gp = gp.cut('MolNum'=>1).cut('Temp'=>300).cut('Press'=>0.4e-4..1e6).cut('VMRMol2'=>100e-6).cut('VMRMol3'=>0).cut('BinWaveNum'=>6050e2).cut('VMRMol1'=>1e-3)
gp = gp.log10
GGraph.line( gp, false, 'index'=>30 )

DCL.grcls