Polygon use on map, multiple boxes

Dear NCL, I tried using the polygon example
http://www.ncl.ucar.edu/Applications/Scripts/polyg_4.ncl
 to place several boxes over a map. It worked for the first box, yet
when I try to place two more boxes on the map, I get these errors
regarding elements. Could I have some advice on how to get past these
errors? I looked in the ncl archives and I don't see any question like
this. Code and terminal history are below.

Thank you in advance
Erik

;================================================;
; maponly_2.ncl
;================================================;
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/wrf/WRFUserARW.ncl"
; ================================================;
begin
  wks = gsn_open_wks("ps","maponly") ; open a ps file
  res = True
  res_at_pmTickMarkDisplayMode = "Always" ; turn on built-in tickmar
  res_at_gsnDraw = False ; don't draw yet
  res_at_gsnFrame = False ; don't advance yet
  res_at_mpFillOn = False ; turn off gray fill
  res_at_mpOutlineBoundarySets = "National" ; turn on country boundaries
  res_at_mpGeophysicalLineColor = "Navy" ; color of cont. outlines
  res_at_mpGeophysicalLineThicknessF = 1.5 ; thickness of outlines
; res_at_gsnMaximize = True
  res_at_mpGridSpacingF = 0.50
  res_at_mpMaxLatF = 35 ; choose subregion
  res_at_mpMinLatF = -20
  res_at_mpMaxLonF = 35
  res_at_mpMinLonF = -35
  plot = gsn_csm_map_ce(wks,res) ; draw map

  ;************************************************
  ; create points for box1
  ;************************************************
   ypts = (/ 18.0, 18.0, 5.0, 5.0, 18.0/)
   xpts = (/-30.0, -20.0,-20.0, -30.0,-30.0/)
  ;************************************************
  ; add the box
  ;************************************************
  res1 = True ; polyline mods desired
  res1_at_gsLineColor = "red" ; color of lines
  res1_at_gsLineThicknessF = 2.0 ; thickness of lines
  res1_at_gsLineLabelString= "L1" ; adds a line label string

  ; create array of dummy graphic variables. This is required, b/c each line
  ; must be associated with a unique dummy variable.

  dum = new(4,graphic)
  ; draw each line separately. Each line must contain two points.
  do i = 0 , 3
    dum(i)=gsn_add_polyline(wks,plot,xpts(i:i+1),ypts(i:i+1),res1)
  end do
  ;************************************************
  ; label the box with additional text
  ;************************************************
  tres = True
  tres_at_txFontHeightF = 0.02
  gsn_text(wks,plot,"L1",-25.0,12.0,tres)

  draw(plot)
  frame(wks)
  delete(xpts)
  delete(ypts)
  delete(dum)
  ;************************************************
  ; create points for box2
  ;************************************************
   ypts = (/ 18.0, 18.0, 5.0, 5.0, 18.0/)
   xpts = (/-20.0, -10.0,-10.0, -20.0,-20.0/)
  ;************************************************
  ; add the box
  ;************************************************
  res2 = True ; polyline mods desired
  res2_at_gsLineColor = "red" ; color of lines
  res2_at_gsLineThicknessF = 2.0 ; thickness of lines
  res2_at_gsLineLabelString= "L2" ; adds a line label string

  ; create array of dummy graphic variables. This is required, b/c each line
  ; must be associated with a unique dummy variable.

  dum = new(4,graphic)
  ; draw each line separately. Each line must contain two points.
  do i = 0 , 3
    dum(i)=gsn_add_polyline(wks,plot,xpts(i:i+1),ypts(i:i+1),res2)
  end do
  ;************************************************
  ; label the box with additional text
  ;************************************************
  tres2 = True
  tres2_at_txFontHeightF = 0.02
  gsn_text(wks,plot,"L2",-15.0,12.0,tres2)

  draw(plot)
  frame(wks)
   delete(xpts)
  delete(ypts)
  delete(dum)

   ;create points for box3
  ;************************************************
   ypts = (/ 18.0, 18.0, 5.0, 5.0, 18.0/)
   xpts = (/-10.0, 10.0,10.0, -10.0,-10.0/)
  ;************************************************
  ; add the box
  ;************************************************
  res3 = True ; polyline mods desired
  res3_at_gsLineColor = "red" ; color of lines
  res3_at_gsLineThicknessF = 2.0 ; thickness of lines
  res3_at_gsLineLabelString= "L3" ; adds a line label string

  ; array of dummy graphic variables. This is required, b/c each line
  ; be associated with a unique dummy variable.

  dum = new(4,graphic)
  ; draw each line separately. Each line must contain two points.
  do i = 0 , 3
    dum(i)=gsn_add_polyline(wks,plot,xpts(i:i+1),ypts(i:i+1),res3)
  end do
  ;************************************************
  ; label the box with additional text
  ;************************************************
  tres2 = True
  tres2_at_txFontHeightF = 0.02
  gsn_text(wks,plot,"L2",-15.0,12.0,tres2)

  draw(plot)
  frame(wks)

    delete(xpts)
  delete(ypts)
  delete(dum)

end

"makemap.ncl" 129 lines, 5085 characters

ENOBLE_at_ATHENA: /usr/people/enoble> ncl makemap.ncl
 Copyright (C) 1995-2007 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 4.3.1
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.
warning:TransformPostDraw: tfPolyDrawList element 0 is invalid
warning:TransformPostDraw: tfPolyDrawList element 1 is invalid
warning:TransformPostDraw: tfPolyDrawList element 2 is invalid
warning:TransformPostDraw: tfPolyDrawList element 3 is invalid
warning:TransformPostDraw: tfPolyDrawList element 0 is invalid
warning:TransformPostDraw: tfPolyDrawList element 1 is invalid
warning:TransformPostDraw: tfPolyDrawList element 2 is invalid
warning:TransformPostDraw: tfPolyDrawList element 3 is invalid
warning:TransformPostDraw: tfPolyDrawList element 4 is invalid
warning:TransformPostDraw: tfPolyDrawList element 5 is invalid
warning:TransformPostDraw: tfPolyDrawList element 6 is invalid
warning:TransformPostDraw: tfPolyDrawList element 7 is invalid
ENOBLE_at_ATHENA: /usr/people/enoble>
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jul 21 2008 - 13:13:32 MDT