Python：SPEI指数的计算

1.准备工作

将nc文件修改为想要格式：参考sample_precip.nc和sample_tavg.nc
(1）根据ERA5小时数据求月总降水和月平均气温
在python中

precip_sum = precip.resample(time="1M").sum()
tavg = tavg.resample(time="1M").mean()

(2) 在nco中修改nc文件的格式

ncatted -a units,precip,c,c,millimeter in.nc   #给precip变量添加单位属性为“millimeter”
ncatted -a units,tavg,c,c,K in.nc    #给tavg变量添加单位属性为“K”
ncatted -O units,precip,o,c,"millimeter" in.nc   #修改precip的单位为millimeter
ncatted -a ,precip,d,, in,nc    #删除precip变量的所有属性
ncrename -v(-d) old_name,new_name .nc   #重命名变量（属性）
##修改维度在python中实现
precip_monthly = precip_monthly.rename({'longitude':'lon','latitude':'lat'})
ncpdq -v precip -a lat,lon,time in.nc out.nc   #将precip变量的纬度顺序更改为lat,lon,time

2.在Anaconda prompt终端运行

首先激活软件包环境
conda activate indices_env
（1）求pet

process_climate_indices --index pet