关于南亚高压东伸脊点的计算

Jyy1108

根据摸鱼大佬提供的西太副高西伸脊点的python程序，改编计算南亚高压东伸脊点的ncl程序，但结果与文献不一致。以下是参考的python脚本和我改编的ncl脚本，有大佬计算过类似的程序吗，求指教。
python：
wpsh_wp = np.zeros((42))
for i in range(z1_array.shape[0]):
    for j in range(z1_array.shape[2]):
        if z1_array[i,:,j].max()>=5880:
            wpsh_wp = lon[j]
            break
    if wpsh_wp==0:
        for j in range(z1_array.shape[2]):
            if z1_array[i,:,j].max()>=5870:
                wpsh_wp = lon[j]
                break

ncl：
    data = addfile("F:\data\NCEP1\hgt.mon.mean.nc", "r")
    ;print(data)

    it_s = 196101
    it_e = 201512
    time = data->time
    YYYY = cd_calendar(time, -1)
    rec_s = ind(it_s.eq.YYYY)
    rec_e = ind(it_e.eq.YYYY)

    hgt0 = data->hgt(rec_s:rec_e,{200},{0:90},:)
    ;printVarSummary(hgt0)

    hgt = month_to_season(hgt0, "JJA")

    ; 提取北半球0°～180°~170°W范围内的数据
    north_hem_data_1 = hgt(:,:,{0:190})

    ; 提取北半球30°W ～ 0°范围内的数据
    north_hem_data_2 = hgt(:,:,{330:360})

    ; 合并三个范围的数据
    north_hem_data = array_append_record(north_hem_data_1({lon|:},{time|:},{lat|:}), north_hem_data_2({lon|:},{time|:},{lat|:}), 0)

    ; 打印结果
    ;printVarSummary(north_hem_data)
    lon = north_hem_data&lon
    ;print(lon)

    hgt_200 = north_hem_data({time|:},{lat|:},{lon|:})
    ;printVarSummary(hgt_200)

wpsh = new((/55/),float,0)
;print(wpsh)

do ny = 0,54
    do iy = 0,88
        if(max(hgt_200(ny,:,iy)).ge.12520) then
            wpsh(ny) = lon(iy); write branch
        end if
    end do
end do
print(wpsh)


参考文献：
东( 西) 伸脊点指数: 北半球 30°W ～ 0° ～180° ～ 170° W 范 围 内，1 250 dagpm 等 值 线 最 东( 西) 端所在经度。（夏季南亚高压年代际变化及其
对长江中下游降水的影响. 张莹等，2019）