NCL where 求助！！！！！

暮夜冬雨

a=addfile("wrfout_d01_2014-08-24_90","r")

aa= a->QVAPOR

ba= a->QCLOUD

ca=a->QRAIN

da=a->QICE

ea=a->QSNOW

fa=a->QGRAUP

s=(aa+ba+ca+da+ea+fa)*1000

wa=a->W

do i=0,79

w=new((/25,80,100,100/),float)

    w(:,i,:,:)=0.50*(wa(0:24,i,0:99,0:99)+wa(0:24,i+1,0:99,0:99))

end do

s_1=dimsizes(w)

w_1=dimsizes(w)

sum_1=dimsizes(w)

s_1 = where (s.gt.0.01,s_1*2,s_1)

;w_1 = where (w.gt.1.00,1,0)

我是想通过where计算出s_1，w_1大于阈值的个数，其中s，w，s_1，w_1都是4维25*80*100*100的数组，请问各位大神，我要怎么修改啊

暮夜冬雨

s_1=dimsizes(w)
w_1=dimsizes(w)
sum_1=dimsizes(w)
s_1 = where (s.gt.0.01,1,s_1@_FillValue)
;w_1 = where (w.gt.1.00,1,0)

暮夜冬雨

是这样子的，上面是我胡乱改动测试的

御风者

你用的函数不太对症，有更简单的。用ind函数。
; 已知s_1，w_1都是4维数组，想取各自超过某一阈值(设为1000)的个数。
; 这里以s_1这个数组为例吧，代码如下：

temp_s_1 = ndtooned(s_1)
indices = ind(temp_s_1.gt.1000.)
n_gt1000 =  dimsizes(indices)
三行代码搞定。

暮夜冬雨

多谢帮助，已解决！！感谢