图做出来无经纬度，想请教一下原因

Tirpitz

使用全球地形起伏模型(ETOPO2)的数据想要绘制地形等高线，但是使用gsn_csm_contour_map最后得到的图上没有经纬度，请问是什么原因呢？
原始文件信息、部分代码以及样图如下。

原始文件信息：

1. #ETOPO2v2g_f4.nc
   
2. dimensions:
   
3.         x = 10801 ;
   
4.         y = 5401 ;
   
5. variables:
   
6.         float x(x) ;
   
7.                 x:long_name = "x" ;
   
8.                 x:actual_range = -180., 180. ;
   
9.         float y(y) ;
   
10.                 y:long_name = "y" ;
    
11.                 y:actual_range = -90., 90. ;
    
12.         float z(y, x) ;
    
13.                 z:long_name = "z" ;
    
14.                 z:_FillValue = NaNf ;
    
15.                 z:actual_range = -10722.f, 8046.f ;
    
16. 
    
17. // global attributes:
    
18.                 :Conventions = "COARDS" ;
    
19.                 :title = "z" ;
    
20.                 :history = "grdreformat ETOPO2v2g_GMT_int.grd ETOPO2v2g_GMT_flt2.grd" ;
    
21.                 :node_offset = 0 ;
    
22. }

复制代码

绘图代码：

1. begin
   
2. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
3. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   
4. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   
5. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
   
6. load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"
   
7. load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/ut_string.ncl"
   
8. 
   
9. ;选定范围               
   
10.                 f2 = addfile("ETOPO2v2g_f4.nc","r")
    
11.                 z = f2->z((2700+34*30):(2700+44*30),(5400+87*30):(5400+108*30))
    
12.                 y = f2->lat((2700+34*30):(2700+44*30))
    
13.                 x = f2 ->lon((5400+87*30):(5400+108*30))
    
14. ;为二维数组z添加属性               
    
15.                 z!0 = "lat"
    
16.                 z!1 = "lon"
    
17.                 x@units = "degrees_east"
    
18.                 y@units = "degrees_north"
    
19.                 z&lat = y
    
20.                 z&lon = x
    
21.                 z&lat@units = "degrees_north"
    
22.                 z&lon@units = "degrees_east"
    
23. 
    
24.                 wks = gsn_open_wks("png","2018_uv10")
    
25. 
    
26. ;set basic resource
    
27.         res                = True
    
28.         res@gsnDraw        = False
    
29.         res@gsnFrame       = False
    
30.         res@gsnMaximize    = True
    
31. 
    
32.                 res@gsnLeftString  = ""
    
33.                 res@gsnRightString = ""
    
34.                 res@vpWidthF       = 0.6           ;窗口大小
    
35.                 res@vpHeightF      = 0.4
    
36.         res@pmLabelBarOrthogonalPosF = 0.05
    
37. 
    
38. ;------Topography
    
39.                 ret = res
    
40.                 ret@cnFillOn = True
    
41.                 ret@cnLinesOn = False
    
42.                 ret@cnFillPalette = "MPL_Greys"
    
43. 
    
44.                 ret@mpDataBaseVersion                        = "MediumRes"
    
45.         ret@mpOutlineOn                                        = True                       
    
46.                 ret@mpDataSetName                                = "Earth..4"
    
47.                 ret@mpOutlineSpecifiers                        = (/"China:states","Gansu"/)
    
48.         ret@mpGeophysicalLineThicknessF        = 3
    
49.         ret@mpNationalLineThicknessF        = 3
    
50.                 ret@mpFillDrawOrder                                = "PostDraw"
    
51.                 ret@mpFillOn                                        = False
    
52.                 ret@mpProjection                = "mercator"      
    
53.                 ret@mpLimitMode                 = "Corners"
    
54.                 ret@mpFillAreaSpecifiers                = (/"water","land" /)
    
55.                 ret@mpSpecifiedFillColors                = (/"white","white"/)
    
56.                 ret@mpSpecifiedFillColors                = (/100,0/)
    
57.                 ret@mpLeftCornerLatF            = min(y)
    
58.                 ret@mpLeftCornerLonF            = min(x)
    
59.                 ret@mpRightCornerLatF           = max(y)      
    
60.                 ret@mpRightCornerLonF           = max(x)
    
61. 
    
62.                 ret@lbLabelPosition = "left"
    
63.                 ret@lbOrientation = "Vertical"
    
64.                 ret@pmLabelBarOrthogonalPosF =-1.25
    
65.                
    
66.                 ret@gsnAddCyclic = False                ;未覆盖全地球，值设置为False
    
67. 
    
68.                 ret@tmXBTickSpacingF=4.            ;坐标轴间距
    
69.                 ret@tmYLTickSpacingF=2.
    
70. 
    
71.                 top = gsn_csm_contour_map(wks,z,ret)
    
72. 
    
73. draw(top)
    
74.         frame(wks)
    
75. 
    
76. end
    

复制代码