ncl画grd数据的问题

liyuan3970

用ncl画一个grd数据数据是96*73*65的ctl是*Dec-Feb frequency of heat wave
dset hwf.grd
undef -999.0
xdef 96 linear -180 3.75
ydef 73 linear -90 2.5
zdef 1 linear 1 1
tdef 65 linear jan1950 1yr
vars 1
ap 0 0 rain
endvars
我画一个气候场写ncl的脚本是
;*******************导入库函数******************************
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/wrf/WRFUserARW.ncl"
;************************************************
begin
path = "hwf.grd"
dims = (/65,96,73/)
f=fbindirread(path,0,dims,"float")
a=new((/65,96,73/),"float")
do i=0,64
do j=0,95
do k=0,72
a(i,j,k)=where(f(i,j,k).eq.-999.0,0,f(i,j,k))
end do
end do
end do   
a!0="time"
a!1="lon"
a!2="lat"
a&time=ispan(0,64,1)
a&lon=fspan(0,356.25,96)
a&lat=fspan(-87.5,90,73)
a&lon@unit="degrees_east"
a&lat@unit="degrees_north"
b=dim_avg_n(a,0)
b!0="lon"
b!1="lat"
b&lon=fspan(0,356.25,96)
b&lat=fspan(-87.5,90,73)
b&lon@unit="degrees_east"
b&lat@unit="degrees_north"


printVarSummary(b)
wks = gsn_open_wks("png","avg")
res=True
;;;;;;;;;;;限定地图区域绘图;;;;;;地图和数据要分别设置，互相独立;;;;;;;;;;;;;
res@gsnAddCyclic        = False
res@mpLimitMode       = "LatLon"
res@mpMinLatF         = 0      
res@mpMaxLatF         = 90
res@mpMinLonF         = 60
res@mpMaxLonF         = 130
;;;;;;;;;;;;;只显示中国区域的图;;;;;;;;;;;;;;;;;;;;;;
res@mpDataBaseVersion="Ncarg4_1"                ;中等分辨率
res@mpDataSetName="Earth..4"                        ;第4版地图，有中国边界数据
res@mpOutlineOn            = True
res@mpOutlineSpecifiers=(/"China:states","Taiwan"/)        ;中国边界有问题，缺藏南、台湾
res@mpOutlineBoundarySets ="NoBoundaries"

;设置等值线
res@cnLevelSelectionMode="ManualLevels"
res@cnSmoothingOn=True
;res@
res@cnMinLevelValF=0
res@cnMaxLevelValF=12
res@cnLevelSpacingF=0.5
plot = gsn_csm_contour_map(wks,b,res)      
end


提示的错误是
   check_for_y_lat_coord: Warning: Data either does not contain a valid latitude coordinate array or doesn't contain one at all.
(0)    A valid latitude coordinate array should have a 'units' 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 a valid longitude coordinate array or doesn't contain one at all.
(0)    A valid longitude coordinate array should have a 'units' attribute equal to one of the following values:
(0)        'degrees_east' 'degrees-east' 'degree_east' 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'

第一个循环是去缺省值的，感觉ncl不擅长处理grd数据这样出的图等值线非常的密很难看求大神指导

liyuan3970

没有设置坐标变量导致的问题，现在想想这个太简单了

handsoman22

看我幸运不？！呵呵呵～

南信陶局

感谢您的分享

朱丽飞

具体是怎么解决的？麻烦您能解释一下嘛？