为什么会读数据出错！！

xiaotuanyuan

我想画200hpa风场的多年平均，下载的资料是月平均资料，1948年1月到至今
Variable: U
Type: float
Total Size: 582575040 bytes
            145643760 values
Number of Dimensions: 4
Dimensions and sizes:        [time | 815] x [level | 17] x [lat | 73] x [lon | 144]
Coordinates:
            time: [1297320..1891944]
            level: [1000..10]
            lat: [90..-90]
            lon: [ 0..357.5]
上面是资料，时间维数那里是有815个的，是对的但是写脚本的时候写成：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/cd_string.ncl"
begin
in = addfile("DATA/uwnd.mon.mean.nc","r")
f = addfile("DATA/vwnd.mon.mean.nc","r")
u = in->uwnd
v = f->vwnd
wks = gsn_open_wks("pdf","uvmain")
res = True
res@vcRefMagnitudeF = 10.0
res@vcRefLengthF = 0.045
res@vcMinDistanceF = 0.017
do i=1,815
u(i,9,:,:)=u(i-1,9,:,:)+u(i,9,:,:)
v(i,9,:,:)=v(i-1,9,:,:)+v(i,9,:,:)
i=i+1
end do
u(i,9,:,:)=u(i,9,:,:)/815
v(i,9,:,:)=v(i,9,:,:)/815
plot = gsn_csm_vector_map_ce(wks,u(i,9,:,:),v(i,9,:,:),res)
end
写成这个循环，出现了
fatal:Subscript out of range, error in subscript #0
fatal:An error occurred reading u
fatal:["Execute.c":8575]:Execute: Error occurred at or near line 16 in file test1

为什么读u的时候是错的，如果我把815改成12 ，就是只读前12个月的，就可以，这是哪里错了，还有算多年平均的风场这种，ncl应该怎么实现，我觉得我这个办法很蠢，希望大家多多指教  谢谢大家

尽头的尽头

ncl的第一维是0，你循环改成0，814就可以了，而且ncl有自己的函数计算平均

气象小新

xiaotuanyuan

尽头的尽头 发表于 2015-12-22 08:35
ncl的第一维是0，你循环改成0，814就可以了，而且ncl有自己的函数计算平均

我的是i-1,所以其实起始维数是0，有计算函数啊，我就说这样很蠢  我去看看函数