require "numru/netcdf"
require "numru/dcl"
include NumRu

upx = [ 1.0, 5.0, 10.0, 50.0, 100.0 ]
upy = [ ]

file1 = NetCDF.open("obratmos_jupiter_x1_441.0K-2.0d6Pa.nc")   
var_temp = file1.var("stab")               # 変数をオープン
temp = var_temp.get
upy.push( temp.max )
file1.close 

file1 = NetCDF.open("obratmos_jupiter_x5_425.0K-2.0d6Pa.nc")   
var_temp = file1.var("stab")               # 変数をオープン
temp = var_temp.get
upy.push( temp.max )
file1.close 

file1 = NetCDF.open("obratmos_jupiter_x10_407.0K-2.0d6Pa.nc")   
var_temp = file1.var("stab")               # 変数をオープン
temp = var_temp.get
upy.push( temp.max )
file1.close 

file1 = NetCDF.open("obratmos_jupiter_x50_736.0K-5.0d7Pa.nc")   
var_temp = file1.var("stab")               # 変数をオープン
temp = var_temp.get
upy.push( temp.max )
file1.close 

file1 = NetCDF.open("obratmos_jupiter_x100_642.0K-1.0d8Pa.nc")   
var_temp = file1.var("stab")               # 変数をオープン
temp = var_temp.get
upy.push( temp.max )
file1.close 

p upy


DCL::swpset('LDUMP', true) 
DCL::swpset('LWAIT',false)
DCL::swpset('LWAIT1',false)  
DCL::swpset('IPOSX', 50)
DCL::swpset('IPOSY',50)
DCL::swpset('IHEIGHT', 600 )
DCL::swpset('IWIDTH',  600 )

DCL::gropn(1)
DCL::grfrm
DCL::grstrn(1)                      
DCL::usspnt(upx, upy)
DCL::uspfit
DCL::grstrf                          # 確定

DCL::ussttl('abundance', 'x solar', 'stability', 's^-2')
DCL::usdaxs

DCL.uusmkt(3)
DCL.uusmks(0.03)
DCL.uusmki(3)
DCL::uumrk(upx, upy)            

DCL::uxsttl('T', 'Jupiter', 0.0)

DCL::grcls