关于ncl坐标设定问题

李乐乐

我想给一组数据添加经纬度属性，画图的时候报错，想请问各位是哪里出错了，部分程序如下


f = asciiread("/home/yxli/slp.txt",(/72,37,12,155/),"float")

;printVarSummary(f)
slp2 = f1->Press
slp1 = f
;printVarSummary(slp2)
slp = new((/72,37,12,163/),float）


slp!0 = "lon"
slp!1 = "lat"
slp!2 = "month"
slp!3 = "year"
lon_p = ispan(-180,175,5)
lat_p = ispan(-90,90,5)
month_p = ispan(1,12,1)
year_p = ispan(1850,2004,1)
lon_p@units = "degrees_east"
lat_p@units = "degrees_north"
slp&lon = lon_p
slp&lat = lat_p
;slp&lon@units = "degrees_east"
;slp&lat@units = "degrees_north"


报错为
(0) check_for_y_lat_coord: Warning: Data either does not contai
n(0) a valid latitude coordinate array or doesn't contain one at
all.(0) A valid latitude coordinate array should have a 'units'
(0) attribute equal to one of the following values:
(0)     'degrees_north' 'degrees-north' 'degree_north' 'degrees
north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'(0) check_for_lon_coord: Warning: Data either does not contain
(0) a valid longitude coordinate array or doesn't contain one a
t all.(0) A valid longitude coordinate array should have a 'units'
(0) attribute equal to one of the following values:
(0)     'degrees_east' 'degrees-east' 'degree_east' 'degrees ea
st' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'