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)