[ncl-talk] Defined polyline is missing on map

[Beata Szabo]

Dear All,

I intend to denote eight subregions in a European map by gsn_add_polyline
function. I tried to carry out it by the following way:


  ;-- create points for box
;************************************************

   ybi = (/ 59.0, 59.0, 50.0,  50.0,  59.0/)
   xbi = (/-10.0,  2.0,  2.0, -10.0, -10.0/)

   yip = (/ 44.0, 44.0, 36.0,  36.0,  44.0/)
   xip = (/-10.0,  3.0,  3.0, -10.0, -10.0/)

   yfr = (/ 50.0, 50.0, 44.0,  44.0,  50.0/)
   xfr = (/ -5.0,  5.0,  5.0,  -5.0,  -5.0/)

   yme = (/ 55.0, 55.0, 48.0,  48.0,  55.0/)
   xme = (/  2.0, 16.0, 16.0,   2.0,   2.0/)

   ysc = (/ 70.0, 70.0, 55.0,  55.0,  70.0/)
   xsc = (/  5.0, 30.0, 30.0,   5.0,   5.0/)

   yal = (/ 48.0, 48.0, 44.0,  44.0,  48.0/)
   xal = (/  5.0, 15.0, 15.0,   5.0,   5.0/)

   ymd = (/ 44.0, 44.0, 36.0,  36.0,  44.0/)
   xmd = (/  3.0, 25.0, 25.0,   3.0,   3.0/)

   xea = (/ 55.0, 55.0, 44.0,  44.0,  55.0/)
   yea = (/ 16.0, 30.0, 30.0,  16.0,  16.0/)

; add the boxes
;************************************************
  resp                  = True                      ; polyline mods desired
  resp at gsLineColor      = "black"                     ; color of lines
  resp at gsLineThicknessF = 5.0                       ; thickness of lines


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

 dumbi = new(4,graphic)
  dumip = new(4,graphic)
  dumfr = new(4,graphic)
  dumme = new(4,graphic)
  dumsc = new(4,graphic)
  dumal = new(4,graphic)
  dummd = new(4,graphic)
  dumea = new(4,graphic)

; draw each line separately. Each line must contain two points.
  do i = 0 , 3
    dumbi(i)=gsn_add_polyline(wks,plot,xbi(i:i+1),ybi(i:i+1),resp)
    dumip(i)=gsn_add_polyline(wks,plot,xip(i:i+1),yip(i:i+1),resp)
    dumfr(i)=gsn_add_polyline(wks,plot,xfr(i:i+1),yfr(i:i+1),resp)
    dumme(i)=gsn_add_polyline(wks,plot,xme(i:i+1),yme(i:i+1),resp)
    dumsc(i)=gsn_add_polyline(wks,plot,xsc(i:i+1),ysc(i:i+1),resp)
    dumal(i)=gsn_add_polyline(wks,plot,xal(i:i+1),yal(i:i+1),resp)
    dummd(i)=gsn_add_polyline(wks,plot,xmd(i:i+1),ymd(i:i+1),resp)
    dumea(i)=gsn_add_polyline(wks,plot,xea(i:i+1),yea(i:i+1),resp)

  end do

On the resulted map the first seven boxes are presented but the last one
(xea, yea) is missing (the map is attached). I do not understand why. Can
someone inform me what I did wrong and how I should solve this issue?

Thank you for your help in advance!
Beata

Dr. Beata Szabo-Takacs
Global Change Research Institute CAS
Bělidla 986/4a
60300 Brno
Czech Republic
-------------- next part --------------
A non-text attachment was scrubbed...
Name: MAP.000002.png
Type: image/png
Size: 216140 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160926/fb9f7d0a/attachment.png