# -*- coding: euc-jp -*-
require "numru/ggraph"
include NumRu

#files = "thermal-moist_ExnerThermExp.nc"
files = "zzz-thermal-moist_ExnerThermExp.nc"
gphys1 = GPhys::IO.open(files, 'ExnerThermExp')


#-- 圧力方程式における熱膨張項の寄与の水平平均値を計算 --#
gphysdisp0 = gphys1.mean( 'x' ).cut( 't'=>2592000 )

print "\n Contribution of Thermal Expansion Term in Pressure Eq. (30 day) :\n"
p gphysdisp0.val