ncl等值线图加点求助!

对方的凤飞飞

  ;-----------------------------------------------------------------
;
;  DKRZ - NCL Example:  NCL_map_settings.ncl
;
;  30.10.14 kmf
;-----------------------------------------------------------------
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin
  setvalues NhlGetWorkspaceObjectId()
  "wsMaximumSize": 600000000
  end setvalues
  f1 =  addfile("ETOPO2v2g_f4.nc", "r")
  var1  =  f1->z(:,:)
  lon = f1->x(:)
  lon@units = "degrees_east"
  lat = f1->y(:)
  lat@units = "degrees_north"
  var1!0 = "lat" ;第一维的名称
  var1&lat = lat ;第一维的值，直接从文件中的变量u中提取
  var1!1 = "lon" ;第二维的名称
  var1&lon = lon  ;第二维的值
  wks = gsn_open_wks("png","plot3")
  gsn_define_colormap(wks,"WhiteBlueGreenYellowRed");WhiteBlueGreenYellowRed

  res                       =  True
  res@gsnMaximize           =  True            ;-- maximize graphics output
  res@gsnLeftString=" "
  res@gsnRightString=" "
  res@tmXBMode     = "Explicit"
  res@gsnDraw        = False      ; Will draw after overlay
  res@gsnFrame       = False
  res@lbLabelStride         =  1               ;-- every other label
  res@lbBoxMinorExtentF     =  0.15   
  res@gsnAddCyclic=False         ;-- decrease the height of the labelbar
  res@pmLabelBarOrthogonalPosF = -0.07         ;-- move the labelbar upward
  res@gsnPolar="SH"
  res@cnFillOn              =  True            ;-- turn on contour fill
  res@cnLinesOn             =  False           ;-- turn off contour lines
  res@cnLineLabelsOn        =  False           ;-- turn off line labels
  res@cnLevelSelectionMode  = "ManualLevels"   ;-- set contour levels manually
  res@cnMinLevelValF        =  0      ;-- minimum contour level
  res@cnMaxLevelValF        =4500         ;-- maximum contour level
  res@cnLevelSpacingF       = 200              ;-- contour level spacing
  res@mpCenterLonF          = 0              ;-- center at lon=10
  res@mpMaxLatF          = -60   
  res@mpDataBaseVersion     = "RANGS"      ;-- map resolution "MediumRes"  
  ;res@mpGridAndLimbOn       =  True            ;-- plot grid lines
  res@mpGridLineColor       = "grey30"         ;-- set grid line color
  res@mpPerimOn             =  False           ;-- don't draw the box around the plot
  res@mpLandFillColor             = "transparent"
  res@mpOceanFillColor            = "white"
  res@mpFillDrawOrder             = "PostDraw"
  res@tiMainString          = " "   ;-- title string
  res@tiMainFontHeightF     =  0.02            ;-- title font size

plot = gsn_csm_contour_map_polar(wks, var1, res)

  res1                       =  True
  res1@gsnDraw        = False      ; Will draw after overlay
  res1@gsnFrame       = False
  res1@gsnPolar="SH"
  res1@mpMaxLatF          = -60   
  res1@mpDataBaseVersion     = "RANGS"     
  map = gsn_csm_map(wks,res1)
  lons = (/76.37/)
  lats=(/-69.37/)
;添加text和点
  mkres                       =  True
  mkres@gsMarkerColor = "Black"
  mkres@gsMarkerSizeF = 15.        ; Increase marker sizes.

  dum = gsn_add_polymarker(wks,map,lons,lats,mkres)
  txres               = True
  txres@txFontHeightF = 0.015
  txres@txJust        = "CenterLeft"
  
  text_id = gsn_add_text(wks,map,"Zhongshan",lons,lats,txres)
  overlay(plot,map)
  draw(plot)
  frame(wks)
全文代码如下,就是想在等值线图上添加一个站点和text,是不是不能overlay呀,一直报错fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member,我改成不overlay也不行,那个地方出了问题呢,思考了一个下午!请大神告知!

AXD

加站点及文本确实不需要overlay，map那部分可以在plot部分实现，把map的删了，需要的属性在plot里面加就行，然后把overlay(plot,map)也删了。

对方的凤飞飞

大神,我本来是这样写的  plot = gsn_csm_contour_map_polar(wks, var1, res)
  lats  = (/-79.37/)
  lons  = (/76.37/)
  mkres                       =  True
  mkres@gsMarkerColor = "Blue"
  mkres@gsMarkerSizeF = 15.        ; Increase marker sizes.
  dum = gsn_add_polymarker(wks,plot,lons,lats,mkres)
  txres               = True
  txres@txFontHeightF = 0.02
  txres@txJust        = "CenterLeft"
  text_id = gsn_add_text(wks,plot,"Zhongshan",lons,lats,txres)
但是结果不显示点呀

对方的凤飞飞

AXD 发表于 2020-9-2 19:10
加站点及文本确实不需要overlay，map那部分可以在plot部分实现，把map的删了，需要的属性在plot里面加就行 ...

大神,我本来是这样写的  plot = gsn_csm_contour_map_polar(wks, var1, res)
  lats  = (/-79.37/)
  lons  = (/76.37/)
  mkres                       =  True
  mkres@gsMarkerColor = "Blue"
  mkres@gsMarkerSizeF = 15.        ; Increase marker sizes.
  dum = gsn_add_polymarker(wks,plot,lons,lats,mkres)
  txres               = True
  txres@txFontHeightF = 0.02
  txres@txJust        = "CenterLeft"
  text_id = gsn_add_text(wks,plot,"Zhongshan",lons,lats,txres)
但是结果不显示点呀

AXD

对方的凤飞飞 发表于 2020-9-2 20:28
大神,我本来是这样写的  plot = gsn_csm_contour_map_polar(wks, var1, res)
  lats  = (/-79.37/)
  l ...

我没你的数据，简单画了个没数据的，可以画，代码如下：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin

  wks  = gsn_open_wks("png","maponly")       ; send graphics to PNG file
  res                       =  True
  res@gsnDraw        = False      ; Will draw after overlay
  res@gsnFrame       = False
  res@gsnPolar="SH"

                               
登录/注册后可看大图

  plot = gsn_csm_map_polar(wks,res)        ; draw global map

  lons = (/76.37/)
  lats=(/-69.37/)
;添加text和点
  mkres                       =  True
  mkres@gsMarkerColor = "red"
  mkres@gsMarkerSizeF = 60.        ; Increase marker sizes.
  mkres@gsMarkerIndex = 1

  dum = gsn_add_polymarker(wks,plot,lons,lats,mkres)

  txres               = True
  txres@txFontHeightF = 0.015
  txres@txJust        = "CenterLeft"
  text_id = gsn_add_text(wks,plot,"Zhongshan",lons,lats,txres)

  draw(plot)
  frame(wks)

  end

对方的凤飞飞

AXD 发表于 2020-9-3 07:44
我没你的数据，简单画了个没数据的，可以画，代码如下：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/g ...

你好,我也运行了你这个代码,但是我一旦用上高程数据就不行,能否告诉我你的邮箱,我发给你数据,其实数据是公开的高程数据,但是有点大,大几百兆吧.

对方的凤飞飞

AXD 发表于 2020-9-3 07:44
我没你的数据，简单画了个没数据的，可以画，代码如下：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/g ...

谢谢搞定了.