通过Python计算有效降水量

有效降水量是指能够提供给作物蒸发蒸腾，从而减少作物对灌溉水需求的雨量，可根据美国农业部推荐的方法进行逐月计算，本文将介绍其在Python环境中的实现方式。其公式如下：

\[P_{eff}=\begin{cases}P_{month}\times\left(125-0.2\times P_{month}\right)/125,P_{month}\leqslant250\mathrm{~mm}\\125+0.1\times P_{month},P_{month}>250\mathrm{~mm}\end{cases} \]

式中，\(P_{eff}\) 为月有效降水量 (mm/month)；\(P_{month}\) 为月总降水量 (mm/month)。

1 主要代码

（1）核心算法函数

# 根据月总降水量，计算月有效降水量
def calculate_peff(pre):
    peff = 0
    if pre <= 250:
        peff = pre * (125 - 0.2 * pre) / 125
    else:
        peff = 125 + 0.1 * pre
    return peff

（2）示例数据及运算结果

月降水量示例数据列表：
[36.53, 143.6, 65.43, 78.4, 186.41, 8.33, 177.23, 159.42, 42.51, 145.8, 68.02, 128.7, 73.01, 69.11, 115.83, 107.53, 132.82, 353.26, 41.37, 386.66, 58.2, 126.92, 111.73, 103.41, 257.8, 478.82, 200.01, 27.22, 71.32, 41.71]
月有效降水量运算结果（保留四位小数）：
[34.3949, 110.6065, 58.5803, 68.5655, 130.8121, 8.219, 126.9732, 118.7564, 39.6186, 111.7878, 60.6172, 102.1981, 64.4813, 61.4681, 94.3635, 89.0297, 104.5942, 160.326, 38.6316, 163.666, 52.7804, 101.1461, 91.7563, 86.3002, 150.78, 172.882, 136.0036, 26.0345, 63.1815, 38.9264]

2 参考文献

[1] 黄志刚, 王小立, 肖烨, 等. 气候变化对松嫩平原水稻灌溉需水量的影响[J]. 应用生态学报, 2015, 26(1): 260-268.