极端降水阈值求助

caojinqiu

请问极端降水阈值有没有NCL的程序或者python的程序。95%百分位。

caojinqiu

http://www.ncl.ucar.edu/Applicat ... percentile_plot.ncl
undef("percent_to_value")
function percent_to_value( \
  i_data           : numeric, \
  i_percentiles

: numeric  \
)
local None
begin

  retVal = new( dimsizes(i_percentiles), float )

  data1d = ndtooned( i_data )
  notMissing = data1d( ind(.not.ismissing(data1d) ) )

  qsort(notMissing)

  do p = 0, dimsizes(i_percentiles)-1
;   pInd = round( i_percentiles(p) * .01 * dimsizes(notMissing) + 0.5, 3 ) -1
;   pInd = where( pInd.ge.dimsizes(notMissing), dimsizes(notMissing)-1, pInd )
    floatInd = i_percentiles(p) * .01 * dimsizes(notMissing) - 0.5
    floorInd = toint( floor(floatInd) )
    floorInd = where( floorInd.lt.0, 0, floorInd )
    ceilInd = toint( ceil(floatInd) )
    ceilInd = where( ceilInd.ge.dimsizes(notMissing), \
        dimsizes(notMissing)-1, ceilInd )
;   print(pInd + " " + dimsizes(notMissing))
    if( ceilInd.eq.floorInd ) then
      retVal(p) = notMissing(floorInd)
    else
      retVal(p) = notMissing(floorInd) * ( ceilInd - floatInd ) \
          + notMissing(ceilInd) * ( floatInd - floorInd )
    end if
  end do

  return(retVal)
end ; percent_to_value

caojinqiu

杰西卡·芬 发表于 2020-9-13 11:17
代码没看懂，楼主可以解释一下吗？

我直接调用的，（1）load "/data1/***/***/***/box_percentile_plot.ncl"；把ncl官网上的ncl脚本放在自己的路径下；（2）prec(j,i,k)=percent_to_value(pr(365*j:(365*j+364),i,k),95)，然后做循环，，，你自己去ncl
官网看吧

caojinqiu

已找到

caojinqiu

如果不是这个函数，请大家指正

言深深

这个程序并不难，计算出来序列的均值和方差就ok了，通过查表获得95%分为上面的值（1.96）就可以编写了
abs(x(i)-x平均)>1.96*x方差，即可

zhangyichi123

感谢楼主！

浅色记忆

感谢楼主！

杰西卡·芬

代码没看懂，楼主可以解释一下吗？

YQY优优

caojinqiu 发表于 2020-09-15 22:18
我直接调用的，（1）load "/data1/***/***/***/box_percentile_plot.ncl"；把ncl官网上的ncl脚本放在自己的路径下；（2）prec(j,i,k)=percent_to_value(pr(365*j:(365*j+364),i,k),95)，然后做循环，，，你自己去ncl
官网看吧

你好啊，我运行的时候总是官网的脚本报错，Subscripts out of range 。真是不知道咋回事了

来自苹果APP客户端