求助！新建数组计算日变压出现错误

nightelfxu

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/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"
begin
f = addfile("I:/05-data_daily/slp.1967.nc","r")
; 读取变量
    pressure = f->slp
    ; 获取维度信息
    time = f->time
    lat = f->lat
    lon = f->lon

    ; 获取维度大小
    n_time = dimsizes(time)
    n_lon = dimsizes(lon)
    n_lat = dimsizes(lat)
    ; 创建一个数组来存储逐日变压数据
    daily_pressure_diff = new((/n_time-1, n_lat, n_lon/), typeof(pressure))
    ; 计算逐日变压数据
    do i = 1, n_time-1
        daily_pressure_diff(i,:,:) = pressure(i,:,:) - pressure(i-1,:,:)
    end do
    ; 打印逐日变压数据的维度信息
    printVarSummary(daily_pressure_diff)
end



fatal:Subscript out of range, error in subscript #0
fatal:["Execute.c":8640]:Execute: Error occurred at or near line 26 in file /cygdrive/i/03-2019NCL/JB/testslp.ncl

nightelfxu

错误说我下标超出范围，没有弄明白哪里错了，求助大神！

安达充

daily_pressure_diff(i-1,:,:) = pressure(i,:,:) - pressure(i-1,:,:)

nightelfxu

安达充 发表于 2023-8-10 11:43
daily_pressure_diff(i-1,:,:) = pressure(i,:,:) - pressure(i-1,:,:)

感谢 就是这个问题！