# -*- coding: euc-jp -*-

# 圧力方程式における熱膨張項の寄与の値を計算するスクリプト
# (x 平均)

require "numru/ggraph"
include NumRu

#files = "thermal-moist_ExnerThermExp.nc"
files = "thermal-moist_Exner-ThermExpTerm-DthetaDt.nc"
gphys1 = GPhys::IO.open(files, 'ExnerThermExpTerm')


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

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