NCL修改横纵坐标标签求助

lyricskeko

我这是某卫星反演后的纬向月均风速资料，我现在想绘制纬向和高度的剖面图，但是数据里横纵坐标的标签是错的，想把横坐标标签从0-40改成50°S-50°N，每25度一间隔，纵坐标高度改为18-50km，每5km一间隔，请问如何修改代码
begin
a = addfile("/cygdrive/d/wind_datemean.nc", "r")

; 提取变量
height = a->height
u = a->zonal_wind   ; 风速数据
latitude = a->latitude
u_transposed = transpose(u)  ; 转置 u 数据
; 绘制填色图
wks = gsn_open_wks("png", "h_lat_saber")
res = True
res@cnFillOn = True
;res@gsnLeftString = "Latitude vs Height"
res@tiYAxisString = "Height (km)"        ; 在纵轴旁添加 "Height (km)"
res@cnFillPalette = "BlWhRe"
;res@gsnRightString = "Data"
res@trYMinF = 0
res@trYMaxF = 32
res@cnLevelSelectionMode = "ManualLevels"       ; manually select levels
res@cnLevelSpacingF      = 5.0                  ; contour spacing
res@cnMinLevelValF       = -40.                 ; min level
res@cnMaxLevelValF       =  40.                 ; max level
res@cnLineLabelsOn = True

; 绘制图形
plot = gsn_csm_contour(wks, u_transposed, res)

; 关闭绘图
end