# 表題  make_tempdev.rb
#
# 履歴  2007/10/14 小高正嗣
#
# 内容  温位の水平平均偏差データを作成する
#
require "numru/gphys"
include NumRu

# データの読み込み
#
PotTempBasic = GPhys::NetCDF_IO.open('2dMars_dx200dz100_BasicZ.nc', 'PotTempBasicZ')
PotTemp = GPhys::NetCDF_IO.open('2dMars_dx200dz100_PotTemp.nc', 'PotTemp')

# 水平平均と偏差の計算
#
PotTempAll  = PotTemp + PotTempBasic
PotTempMean = PotTempAll.mean(0)
PotTempDev  = PotTempAll - PotTempMean

outfile = NetCDF.create('2dMars_dx200dz100_PotTempDev.nc')
GPhys::NetCDF_IO.write( outfile, PotTempDev )
outfile.close

outfile = NetCDF.create('2dMars_dx200dz100_PotTempMean.nc')
GPhys::NetCDF_IO.write( outfile, PotTempMean )
outfile.close