怎么把点和图画在一张上面

tianmeng

画了一张温度图，想点几个站点位置在上面，结果就变成了一张图画得是站点，一张图画得是温度平面图，怎么把他们结合在一起？

;   Example script to produce plots for a WRF real-data run,
;   with the ARW coordinate dynamics option.

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
;load "./WRFUserARW.ncl"

begin
;
; The WRF ARW input file.
; This needs to have a ".nc" appended, so just do it.
  a = addfile("./wrfout.nc","r")
; We generate plots, but what kind do we prefer?
  type = "png"
; type = "pdf"
; type = "ps"
; type = "ncgm"
  wks = gsn_open_wks(type,"L")
  gsn_define_colormap(wks,"default")

; Set some basic resources
  res = True
  res@MainTitle                   = "REAL-TIME WRF"
  pltres = True
  mpres = True


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

; What times and how many time steps are in the data set?

    it=0

; First get the variables we will need

     L = wrf_user_getvar(a,"T2",it)
; Plotting options for T
      opts = res
      opts@cnFillOn = True
      opts@ContourParameters = (/ 0., 24., 1./)
      contour_tc = wrf_contour(a,wks,L,opts)
      delete(opts)
   ;Plotting Provincial boundary
     mpres@mpDataBaseVersion="MediumRes"
     mpres@mpDataSetName="Earth..4"
     mpres@mpGeophysicalLineColor = "Black"
     mpres@mpNationalLineColor    = "Black"
     mpres@mpGridLineColor        = "Black"
     mpres@mpLimbLineColor        = "Black"
mpres@mpPerimLineColor       = "Black"
     mpres@mpUSStateLineThicknessF     = 2
     mpres@mpNationalLineThicknessF    = 2
     mpres@mpUSStateLineColor="Black"
     mpres@mpOutlineSpecifiers="China:states"

    ; MAKE PLOTS
      plot = wrf_map_overlays(a,wks,(/contour_tc/),pltres,mpres)

    ; plot
      lat =(/38.93/)
      lon = (/100.43/)
      mres = True
mres@gsMarkerSizeF = 10.
      mres@gsMarkerIndex = 16
      mres@gsMarkerColor = "black"
      dum1 =gsn_add_polymarker(wks,plot,lon,lat,mres)
      draw(wks);
      frame(wks)


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
end

niuda

overlay............

dataiyang

增加页面设置的属性呀！
mpres@gsnDraw=False
mpres@gsnFrame=False

janweis

由于wrf_map_overlays为NCL_WRF的功能，
原函数为wrf_overlays (nc_file, wks, (/graphics/), pltres)
在设置中需要加上pltres@FramePlot＝False
而后在后面加上frame(wks)