有人会用NCL做循环吗？（循环算不同时间、气压层）

fengyy789

一直不太会用NCL的循环语句。。
NCEPII的数据，格式为：859（monthly,1979年-共859个时间数据）*17（气压层）*74（lat）*144(lon);
想计算每一个时间点、每一个气压层，的水汽通量散度。

用do做了两层循环： m = 0,858 (时间点)
                              n = 0,16 （气压层）
但是仿佛并不能用   qu(m,n,:,:) = q(m,n,:,:)*u(m,n,:,:)做计算。。。
求指导，感激不尽！

                               
登录/注册后可看大图

begin

a=addfile("rhum.mon.mean.nc", "r")
b=addfile("uwnd.mon.mean.nc", "r")
c=addfile("vwnd.mon.mean.nc", "r")


lat=a->lat
lon=a->lon
level = a ->level
time = a ->time

q = a->rhum
u = b->uwnd
v = c->vwnd

do m = 0,858
do n = 0,16
q_cal = q(m,n,:,:)
u_cal = u(m,n,:,:)
v_cal = v(m,n,:,:)

qu = u_cal*q_cal/9.8*1000
qv = v_cal*q_cal/9.8*1000
copy_VarCoords(u_cal,qu)
copy_VarCoords(v_cal,qv)

vapord_out=uv2dv_cfd(qu,qv,lat,lon,2)
copy_VarCoords(q_cal,vapord_out)
result(m,n,:,:) = vapord_out
end
end