1. 气温单位从K转变到℃：cdo addc,-273.15 input.nc output_K2C.nc

addc - 加上一个常数；input.nc为日平均温度；output_K2C.nc为输出文件；-b F32是为了转为浮点型，否则

若input.nc中变量不是浮点型，可能会报numeric conversion not representable的错误。

2. 获取指示日平均温度是否大于0的mask：cdo -gtc,0 output_K2C.nc mask.nc

gtc - 大于常数；mask.nc为表示数据是否大于阈值（0）由0,1组成的索引文件

相似的算子有：

eq     Equal

ne     Not equal

le     Less equal

lt     Less than

ge     Greater equal

gt     Greater than

eqc     Equal constant

nec     Not equal constant

lec     Less equal constant

ltc     Less than constant

gec     Greater equal constant

gtc     Greater than constant

3. 根据mask.nc选取数据：cdo ifthen mask.nc output_K2C.nc output_gtc0.nc

ifthen表示如果 mask中变量的值不等于0，输出文件output_gtc0.nc中的变量值等于output_K2C.nc中的值，如果等于0的话，对应位置的值则为缺省值（NaN）

4.  计算每个月的融化指数：cdo monsum output_ge0.nc output_ge0_monTDD.nc

参考文献：

https://wap.sciencenet.cn/blog-1081898-1275862.html?mobile=1

https://www.shizhanxia.com/5964.html

https://cloud.tencent.com/developer/article/1618310