NCL

大气CYX

NCL读取WRF文件，然后提取风的变量，然后怎么计算风速？我写的程序老是报错，注释掉的也试过，不行
sst        = new((/dimsizes(wrffile),nlat,mlon/),float)
do n = 0, dimsizes(wrffile)-1
print("Working On" + wrffile(n))
ff = addfile(wrffile(n),"r")
U = ff->U(8:56,:,:,:)
V = ff->V(8:56,:,:,:)
W = ff->W(8:56,:,:,:)
speed = U*U+V*V+W*W
;(/(/V/)*(/V/)+(/U/)*(/U/)+(/W/)*(/W/)/)
;speed = wind_speed(U,V)
sst(n,:,:) = sst(sqrt(speed),:,:)