【求助】干旱指数SPEI计算过程中出现无效值，求解决方案（附个人改写的IDL代码）

水禹刀

http://bbs.06climate.com/forum.php?mod=viewthread&tid=14971提供的spei.exe文件在计算指数时对于部分站点的数据会得到无效值（不同时间尺度时都会发生），本人对C不熟又不会写makefie，所以只能把C源码改写成IDL。

以甘肃某站点（52323）为例（IDL代码、该站1961-2010年月值均温和降水见附件），对于1963年3月的一个月尺度SPEI输出结果-NaN,同时IDL显示：Program caused arithmetic error: Floating illegal operand

经检查发现问题出现在log-logistic概率计算过程中，3月份该站的概率分布参数γ，α，β依次为-2.7648521；4.3525177 ；3.2208551，概率分布函数F(x)的计算式中的[α/(x-γ)]^β对应1963年3月份的值就导致了（-2.603）^3.22这样的计算问题，查看该站50年的P-PET序列后发现-NaN问题值对应的都是序列中的极端值，似乎是概率函数拟合的问题，导致对站点序列中的极端值计算指数出错，但是本人在之前算SPI的时候从没遇过这种情况，请问该如何解决，谢谢~

SPEItest.rar (7.75 KB, 下载次数: 52)

2013-7-28 00:47 上传

点击文件名下载附件
阅读权限: 10

weirunchu

我也正在学习！

lxdcdut

围观解决方案

Hibiscus

感谢分享！

刘布袋子

求解决方案

jane1124

求fortran版的

maggiezy

我最近也在算，总是算到B（贝塔值）时出错，大神求解

concise

请问有无FORTRAN 程序的？

xjgis

这个很不错，好好学习中！！！！

xiaosashapi

求解决方案~