monthly anomaly

以下のサンプルはpythonからcdoを呼び出して、月毎平均値の適切な期間の気候値と、さらに偏差を求めています。偏差の月別のファイルも出力しているのは、この研究での後の処理を容易にするためです。

#--------------------

import os

import numpy as np

file_mth='../org_plev3d.nc/hgt.mon.mean.nc'

file_clm='hgt.mcl.nc'

file_anm_all='hgt.man_all.nc'

str_cmd='cdo ymonmean -selyear,1980/2010 %s %s'%(file_mth, file_clm)

print(str_cmd)

os.system(str_cmd)

str_cmd='cdo ymonsub -selvar,hgt %s -selvar,hgt %s %s'%(file_mth,file_clm,file_anm_all)

print(str_cmd)

os.system(str_cmd)

for mth in np.arange(1,13):

file_anm_1mo='hgt.man_%2.2i.nc'%mth

str_cmd='cdo selmon,%i %s %s'%(mth,file_anm_all,file_anm_1mo)

print(str_cmd)

os.system(str_cmd)