【zz】ncl画散点图（with markers of different colors and sizes）

http://www.ncl.ucar.edu/Applications/scatter.shtml

;----------------------------------------------------------------------
; scatter_7.ncl
;
; Concepts illustrated:
;   - Drawing a scatter plot with markers of different colors and sizes
;   - Annotating a plot with labelbars, text, and markers
;   - Creating a color map using named colors
;   - Drawing polymarkers and text in NDC space
;----------------------------------------------------------------------
; The original version of this script was contributed by Larry
; McDaniel of IMAGe/NCAR. It had several lines of data processing
; calls before the plotting section.
;
; This script is a pared down version of the original one, with
; the data processing calls removed. The values are read from
; two simple ASCII files.
;----------------------------------------------------------------------
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"

begin
;----------------------------------------------------------------------
; Data is hard-coded here. It was originally calculated by running
; a compute-intensive script.
;----------------------------------------------------------------------
  avdataT = (/1.965796,1.777351,1.915377,1.818644,1.766584,1.56958, \
              2.101717,1.859297,1.662189/)
  avdataP = (/-8.902548,-9.490757,-9.976177,-1.000492,2.248538,\
              -2.698187,9.177564,3.478979,4.059988/)

;----------------------------------------------------------------------
; Open a PS file and define the color table
;----------------------------------------------------------------------
  colors  = (/"white","black","dodgerblue2","firebrick","forestgreen",\
              "dodgerblue2","darkorange2","forestgreen","deeppink2",\
              "darkorange2","purple"/)

  wks = gsn_open_wks ("pdf","scatter")
  gsn_define_colormap(wks,colors)

;----------------------------------------------------------------------
; Set some plot resources
;----------------------------------------------------------------------
  res                    = True

  res@gsnDraw            = False       ; Don't draw plot
  res@gsnFrame           = False       ; Don't advance frame

;
; Don't use gsnMaximize. Instead, we will control
; exactly where plot goes.
;
  res@vpXF              = 0.2
  res@vpYF              = 0.85
  res@vpWidthF          = 0.5
  res@vpHeightF         = 0.5

;---Some calculations to get min/max for Y axis.
  miny = min(avdataP)
  maxy = max(avdataP)
  minx = min(avdataT)
  maxx = max(avdataT)

  amaxy = maxy - abs(miny)
  amaxy = amaxy / 2.
  amaxy = amaxy + 25.
 
  aminy = maxy - abs(miny)
  aminy = aminy / 2.
  aminy = aminy - 25.

;---Resources for X/Y axes.
  res@trYMinF            =  aminy               ; min value on y-axis
  res@trYMaxF            =  amaxy               ; max value on y-axis
  res@trXMinF            =  0.                  ; min value on x-axis
  res@trXMaxF            =  5.                  ; max value on x-axis

;---XY marker resources
  res@xyMarkLineMode     = "Markers"
  res@xyMarkerSizeF      = 0.03
  res@xyMarkerThicknessF = 3.
  res@xyMonoMarkerColor  = False          ; Allow different colors for markers
  res@xyMarkerColors     = colors(2:)     ; It's okay to list more than
                                          ; you need here
  res@xyMarkers          = (/4,4,4,7,7,7,3,6,6/)

;---Tickmark resources
  res@tmXBMode         = "Manual"
  res@tmYLMode         = "Manual"
  res@tmXBTickSpacingF = 1.
  res@tmYLTickSpacingF = 5.

;---Title resources
  res@tiMainString       = "30 av Annual % Change P vs Delta T~C~" + \
                           "               for the WetSouth"
  res@tiXAxisString      = "Delta T"
  res@tiYAxisString      = "% Change P"
  res@tiMainFontHeightF  = 0.03

;----------------------------------------------------------------------
; Create the scatter plot, but don't draw it yet.
;----------------------------------------------------------------------
  plot = gsn_csm_xy (wks,transpose((/avdataT,avdataT/)), \
                         transpose((/avdataP,avdataP/)),res)

;----------------------------------------------------------------------
; Draw time stamp first, then the plot.
;----------------------------------------------------------------------
  infoTimeStamp(wks,0.005,"")
  draw(plot)

;----------------------------------------------------------------------
; Draw some individual labelbars.
;----------------------------------------------------------------------
  labels1  = (/"CRCM","ECP2","HRM3","MM5I","RCM3","WRFG"/)
  colors1  = (/"dodgerblue2","purple","deeppink2","firebrick",\
               "darkorange2","forestgreen"/)
  lbres                    = True          ; labelbar only resources
  lbres@vpWidthF           = 0.1           ; labelbar width
  lbres@vpHeightF          = 0.1           ; labelbar height
  lbres@lbBoxMajorExtentF  = 0.15          ; puts space between color boxes
  lbres@lbMonoFillPattern  = True          ; Solid fill pattern
  lbres@lbLabelFontHeightF = 0.015         ; font height. default is small
  lbres@lbLabelJust        = "CenterLeft"  ; left justify labels
  lbres@lbPerimOn          = False
;
; Each labelbar has just one label. This allows you to more
; easily control where the label goes.
;
  xpos = (/0.05,  0.20,   0.35,   0.50,  0.65,   0.8/)
  do i=0,5
     lbres@lbFillColors     = colors1(i)
     lbres@lbLabelFontColor = colors1(i)
     gsn_labelbar_ndc(wks,1,labels1(i),xpos(i),0.2,lbres)
  end do

;----------------------------------------------------------------------
; Draw some markers and text.
;----------------------------------------------------------------------
  labels2  = (/"CCSM","CGCM3","GFDL","HADCM3","NCEP"/)
  markers1 = (/     4,      7,     6,       3,     2/)
  xpos2    = (/  0.15,   0.30,  0.45,    0.60,  0.75/)
  xpos3    = (/  0.18,   0.33,  0.48,    0.63,  0.78/)

  mkres               = True         ; Marker resources
  txres               = True         ; Text resources
  txres@txFontHeightF = 0.015
  txres@txJust        = "CenterLeft"
  do i = 0,4
     mkres@gsMarkerThicknessF = 3.5
     mkres@gsMarkerSizeF      = 0.02
     mkres@gsMarkerIndex      = markers1(i)

     gsn_polymarker_ndc(wks,           xpos2(i),0.23,mkres)
     gsn_text_ndc      (wks,labels2(i),xpos3(i),0.23,txres)
  end do

  frame(wks)    ; Now advance the frame.
end