NCL画出的图没有显示经纬度度数

Geo-Meteo巡天er

现在刚开始学习NCL，但是在画气压图时，坐标轴上总是不显示经纬度的度数（如图），下面是程序，请教大家看看怎么回事？

begin
;************************************************
; read in netCDF file
;************************************************
  a = addfile("slp.nc","r")
;************************************************
; read in zonal winds
;************************************************
;  slp = short2flt(a->slp)
  t = a->slp(0,:,:)                                ; read July zonal winds
;************************************************
; create plot

  res                   = True

  res@mpLimitMode = "LatLon"
  res@mpMinLatF = 10
  res@mpMaxLatF = 60
  res@mpMinLonF = 70
  res@mpMaxLonF = 140

  res@mpProjection      = "Mercator"        ; choose projection
  res@mpGridAndLimbOn   = True              ; turn on lat/lon lines
  res@mpPerimOn         = True             ; turn off box around plot
  res@mpGridLatSpacingF = 10               ; spacing for lat lines
  res@mpGridLonSpacingF = 10               ; spacing for lon lines
  res@mpFillOn          = False

  res@gsnLeftString  = " "
  res@gsnRightString = " "

  res@cnFillOn          = True              ; color plot desired
  res@cnLineLabelsOn    = False             ; turn off contour lines

;  res@pmTickMarkDisplayMode = "Always"
  res@tmXBOn     = True
  res@tmXBValues = (/70,80,90,100,110,120,130,140/)
  res@tmXBLabels = (/"70~S~o~N~E","80~S~o~N~E","90~S~o~N~E","100~S~o~N~E","110~S~o~N~E","120~S~o~N~E","130~S~o~N~E","140~S~o~N~E"/)
  res@tmYROn     = True
  res@tmYLValues = (/15,25,35,45,55/)
  res@tmYLLabels = (/"15~S~o~N~N","25~S~o~N~N","35~S~o~N~N","45~S~o~N~N","55~S~o~N~N"/)

  res@vpXF            = 0.1                 ; make plot bigger
  res@vpYF            = 0.9         
  res@vpWidthF        = 0.8
  res@vpHeightF       = 0.8

;色标位置颜色设定
  res@lbOrientation       = "horizontal"
  res@lbLabelFontColor    = "black"
  res@lbLabelFontHeightF  = 0.02           ; label bar font height
  res@pmLabelBarOrthogonalPosF = 0.05
  res@lbBoxLinesOn        = True
  res@lbBoxLineThicknessF = 0.3
  res@pmLabelBarHeightF   = 0.06

;色标名称设定
  res@lbTitleOn          = True
  res@lbTitleString      = "~F25~ hPa"
  res@lbTitleFontColor   = "black"
  res@lbTitlePosition    = "Bottom"
  res@lbTitleFontHeightF = 0.018

;主标题
  res@tiMainString       = "~F25~ Sea Level Pressure"         ; add a title
  res@tiMainFontHeightF  = 0.03                               ; font height

  wks = gsn_open_wks("png","SLP")          ; send graphics to PNG file
  gsn_define_colormap(wks,"rainbow")

  contour = gsn_csm_contour_map(wks,t,res)  ; create the plot

;  draw(wks)
;  frame(contour)

end

Geo-Meteo巡天er

已自行解决该问题。

晴天小雅wyy

Geo-Meteo巡天er 发表于 2019-9-4 11:00
已自行解决该问题。

请问楼主怎么解决的？

lc2016

怎么解决的呢

chasing

请问楼主怎么解决的啊