想问下cdo怎么去气候态啊

wulalalala

在cdo里面我想对一个1979-2010年的数据做去除气候态的处理，想问下有没有直接的命令能做呢，我翻了半天手册也没找到
我用cdo  ydaysub  infile1979-2010.nc  -ydayavg   infile1979-2010.nc   outfile 会发生报错：Fields have different size (vfarsub)
想问下要怎么办呀

FDU.EDU

请问楼主解决了吗，遇到了同样的问题

wulalalala

FDU.EDU 发表于 2021-11-23 09:18
请问楼主解决了吗，遇到了同样的问题

没有欸，自己重新写了一个ncl搞的

沃德玛雅柯

分享一个计算月平均场相对于某个气候态的异常，算日数据的话应该换成对应的day的计算就行：首先用selyear和cat命令得到基准期所有的数据（如1850-1900年，每年12个月），然后利用ymonmean命令计算气候态时期各个月的多年平均值，并用ymonsub命令计算原数据与气候态的差值：
# =======1.to select the climatology period data===========
year=2024
month=10
cd /global_st/best/ST+LT/
cdo -f nc selyear,1850 best_monthly_1850-`echo ${year}``echo ${month}`_time.nc best_monthly_1850-1900.nc

for iyear in `seq 1851 1900`
   do
            cdo -f nc selyear,`echo ${iyear}` best_monthly_1850-`echo ${year}``echo ${month}`_time.nc temp.nc
            cdo -f nc cat temp.nc best_monthly_1850-1900.nc
            echo "------------------------------------------"
            echo "    Monthly data of ${iyear}.nc file over !  "
            echo "------------------------------------------"
            rm -f temp.nc
   done  
# =======2.to calcualte the climatology mean and substract it to gain the climotology of data===========
cdo ymonsub best_monthly_1850-`echo ${year}``echo ${month}_time`.nc -ymonmean best_monthly_1850-1900.nc best_monthly_1850-`echo ${year}``echo ${month}`_clim-1850-1900.nc