ncl这样读取dat文件正确吗

隽今、

已经尝试了好几遍不同的读法，但绘出的降水空间分布图还是不符合实际？有人能帮我看看吗，谢谢。
dset J:\Datasets\Obs\GPCP\GPCP_PEN\GPCP_PEN_merged.dat
title variable
undef -9.99e+33
xdef 144 linear 1.25 2.5
ydef 72 linear -88.75 2.5
zdef 1 levels 1000
tdef 3285 linear 00Z03Jan1979 7200mn
vars 1
gpcp=>gpcp 0 t,z,y,x variable
endvars

fili            = "~/bylw/GPCP_PEN_merged.dat"

    setfileoption("bin","ReadByteOrder","LittleEndian")

    gpcp            = fbindirread(fili,0,(/3285,1,72,144/),"float")

    gpcp@_FillValue = -9.99e+33

    gpcp            = where(isnan_ieee(gpcp),gpcp@_FillValue,gpcp)

    gpcp!0          = "time"

    gpcp!1          = "lev"

    gpcp!2          = "lat"

    gpcp!3          = "lon"

    gpcp&time       = ispan(0,3284,1)

    gpcp&lat        = fspan(-88.75,88.75,72)

    gpcp&lon        = fspan(1.25,358.75,144)

    test    = a(::-1,:)

    printVarSummary(test)

    printMinMax(test,True)

    test!0  ="lat"

    test!1  ="lon"

    test&lat=gpcp&lat

    test&lon=gpcp&lon

    test&lat@units="degrees_north"

    test&lon@units="degrees_east"

    wks  = gsn_open_wks("png","/mnt/d/ncl/aaa/op/test.png")   

    gsn_define_colormap(wks,"BlueRed")

    ;;场设置

    res = True

    res@cnFillOn       = True

    res@gsnAddCyclic   = True

    ;res@trYReverse     = True

    res@cnMaxLevelValF = 11.

    res@cnMinLevelValF = 0.

    res@cnLevelSpacingF= 1

    plot1 = gsn_csm_contour_map(wks,test,res)

隽今、

才发现漏了一行代码：a       = dim_avg_n_Wrap(gpcp(:,0,:,:),0)
好痛苦，我感觉读取方式是正确的，为什么这个降水空间分布一点也不符合现实，失去了所有力气……