处理3a12数据，怎么提取不出来经纬度？

南宫蓝天

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
fil=addfile(/3A12.19980101.7.HDF","r")
pr=fil->surfacePrecipitation(499:619,79:159)
lon=fil->surfacePrecipitation&nlon(499:619)
lat=fil->surfacePrecipitation&nlat(79:159)
wks=gsn_open_wks("x11","dengzhixian")
res=True

res@mpLimitMode       = "Corners"
res@mpLeftCornerLonF  = lon(0)
res@mpRightCornerLonF = lon(dimsizes(lon)-1)
res@mpLeftCornerLatF  = lat(0)
res@mpRightCornerLatF = lat(dimsizes(lat)-1)

res@mpOutlineBoundarySets       = "National"         
res@mpGridAndLimbOn             = True     
res@pmTickMarkDisplayMode  = "Always"         

res@sfXArray         = lon
res@sfYArray         = lat

plot = gsn_contour_map(wks,pr,res)

end

fatal:(nlon) is no  a coordinate variable for file (3A12.19980101.7)
fatal:["Execute.c":8575]:Execute: Error occurred at or near line 7 in file Precipitation.ncl

talkd

那就是没有咯

南宫蓝天

talkd 发表于 2015-5-10 11:00
那就是没有咯

dimensions里有nlon，变量里，surfaceprecipitation（nlon,nlat）

talkd

要在Coordinates里有才能读取出来