求取90百分位数

一顾留白

求取三维数组的90百分位数出错，出错行标红了，错误显示：fatal:Dimension sizes of left hand side and right hand side of assignment do not match

do i =0,90,1
    do j=0,180,1  
       h1 = ndtooned(h(:,i,j))
       h2 = h1(ind(.not.ismissing(h1)))
        m = dimsizes(h2)
        qsort(h2)

        a=toint(floor(0.9*m))
        b=toint(ceil(0.9*m))

        h3(i,j)=(h2(a)+h2(b))/2
       print(h3(i,j))

    end do
end do   

我的思路是每个格点得到一个时间序列，选出非缺测值排序，然后找到90百分位的数值存储，最后得到lat×lon的数组h3，麻烦大佬们看一下是什么原因，或者大家有什么别的方法求取90百分位数的吗？  

一顾留白

问题解决啦，下面是正确的脚本

      if (.not.all(ismissing(h)))
        h1 = ndtooned(h)
        h2 := h1(ind(.not.ismissing(h1)))
       end if
   
      m = dimsizes(h2)
      qsort(h2)
      
      if ((toint(floor(0.9*m))) .eq. (0.9*m))then
        h3 =  h2(toint(0.9*m))
      else
        a=toint(floor(0.9*m))
        b=toint(ceil(0.9*m))
        h3=(h2(a)+h2(b))/2  
       end if
      print(h3)

曾璨_stella

还没学这个

Adiligirl

h2 := h1(ind(.not.ismissing(h1)))

一顾留白

谢谢啦，这样改了原来的错误没有了，但是出现了新的错误
fatal:Subscript out of range, error in subscript #0
fatal:An error occurred reading h1
很奇怪，不知道是什么原理

chen--

遇到同样的问题，请问您解决了吗

Ethan87

检查一下数组的范围，应该是提取部分数组数据时出了问题。

想个名字

{:eb502:}{:eb502:}

duoci01

您好！我的数据数组是[年64*日365*纬度145*经度288]四维构成的，现在老板让我做第90个百分位求取极端气温，不知道怎么弄，麻烦大佬指点

一顾留白

duoci01 发表于 2023-3-6 23:09
您好！我的数据数组是[年64*日365*纬度145*经度288]四维构成的，现在老板让我做第90个百分位求取极端气温， ...

我还不太清楚怎么上传文件，麻烦你留个邮箱我把脚本发你~