<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=utf-8">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Dear helpdesk,<br>
I am trying to attach a small map as a legend of a XY plot (with
double different Y axes). The map I need to attach must contain also
a line, added with gsn_add_polyline. The two single plots (the XYY
plot and the map with the polyline) are correct (I drew them
separately first), but as I use gsn_add_annotation I get the
following error message:<br>
<br>
<br>
<font color="#000066">fatal:Invalid plot ID=282 passed to NhlGetBB<br>
warning:pmOverlaySequenceIds isn't a resource in this object<br>
warning:NhlGetValues:Error retrieving pmOverlaySequenceIds<br>
fatal:["Execute.c":8578]:Execute: Error occurred at or near line
1927 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl<br>
<br>
fatal:["Execute.c":8578]:Execute: Error occurred at or near line
6847 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl<br>
<br>
fatal:["Execute.c":8578]:Execute: Error occurred at or near line
7546 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl<br>
<br>
fatal:["Execute.c":8578]:Execute: Error occurred at or near line
647 in file Exercise1.ncl</font><br>
<br>
Do you have any idea the reason why I get such an error? Here is the
code.<br>
<br>
--------------------------------------------------------------------------------------------------------------------------------<br>
;***************************************set resources for the
isostasy plot (Left)<br>
<br>
resL = True<br>
<br>
resL@tiYAxisString = "Variation in bedrock elevation
(30-0 kyrs BP) (m)"<br>
resL@tiXAxisFontHeightF = 0.02<br>
resL@tiYAxisFontHeightF = 0.02<br>
<br>
resL@xyLineColors =
20 ; Define line
colors.<br>
resL@xyLineThicknessF =
3.0 ;
Define line thickness.<br>
resL@xyDashPattern =
0
; Make curves all solid <br>
<br>
resL@trYMaxF = -150.<br>
resL@trYMinF = -850.<br>
<br>
resL@tmXTOn = False<br>
<br>
resL@tmYLMinorOn = False<br>
resL@tmYLLabelFont = 21<br>
resL@tmYLLabelFontHeightF = 0.02<br>
<br>
resL@pmLegendDisplayMode = "Always"<br>
resL@lgPerimOn = False<br>
resL@xyExplicitLegendLabels = "Topography variation (m)"<br>
resL@pmLegendSide = "Bottom" <br>
resL@pmLegendOrthogonalPosF = 0.01<br>
resL@pmLegendWidthF = 0.4<br>
resL@pmLegendHeightF = 0.04<br>
<br>
resL@vpHeightF = 0.4<br>
resL@vpWidthF = 0.4<br>
<br>
<br>
resL@lgAutoManage = False <br>
resL@tmXBMode =
"Explicit"
; Define your own tick mark labels.<br>
;resL@tmXBLabelFont =
21
; Change font of labels.<br>
;resL@tmXBLabelFontHeightF = 0.02<br>
resL@tmXBMinorOn =
False
; No minor tick marks.<br>
resL@tmXBValues =
(/xiso(0),xiso(52)/)
; Location to put tick mark labels<br>
resL@tmXBLabels = (/"A","B"/) <br>
<br>
resL@gsnDraw = False ; do not
draw picture<br>
resL@gsnFrame = False ; do not
advance frame<br>
<br>
<br>
;******************************************set resources for ice
thickness plot (right) <br>
<br>
resR = True<br>
<br>
resR@tiYAxisString = "Variation in ice thickness (30-0
kyrs BP) (m)"<br>
resR@tiXAxisFontHeightF = 0.02<br>
resR@tiYAxisFontHeightF = 0.02<br>
<br>
resR@xyLineColors =
21
; Define line colors.<br>
resR@xyLineThicknessF =
3.0
; Define line thickness.<br>
resR@xyDashPattern =
1
; Make curves all solid<br>
<br>
resR@trYMaxF = 3100<br>
resR@trYMinF = 600<br>
<br>
resR@tmXBOn = False<br>
resR@tmXTOn = False<br>
<br>
resR@tmYRMinorOn = False<br>
resR@tmYRLabelFont = 21<br>
resR@tmYRLabelFontHeightF = 0.02<br>
<br>
resR@vpHeightF = 0.4<br>
resR@vpWidthF = 0.4<br>
<br>
<br>
resR@pmLegendDisplayMode = "Always"<br>
resR@lgPerimOn = False<br>
resR@xyExplicitLegendLabels = "Ice thickness variation (m)"<br>
resR@pmLegendSide = "Bottom"<br>
resR@pmLegendOrthogonalPosF = 0.06<br>
resR@pmLegendWidthF = 0.4<br>
resR@pmLegendHeightF = 0.04<br>
resR@lgAutoManage = False<br>
resR@gsnDraw = False ; do not
draw picture<br>
resR@gsnFrame = False ; do not
advance frame<br>
<br>
<br>
;******************************************set resources for legend
plot<br>
res1 = True<br>
res1@gsnDraw = False<br>
res1@gsnFrame = False<br>
<br>
<br>
res1@mpLimitMode = "Corners" ; choose range
of map<br>
res1@mpLeftCornerLatF = lat(41,1)<br>
res1@mpLeftCornerLonF = lon(41,1)<br>
res1@mpRightCornerLatF = lat(226,139)<br>
res1@mpRightCornerLonF = lon(226,139)<br>
<br>
res1@mpProjection = "LambertEqualArea"<br>
res1@mpCenterLatF = 90<br>
res1@mpCenterLonF = 0<br>
res1@mpOutlineOn = True<br>
res1@tfDoNDCOverlay = True<br>
<br>
res1@cnFillOn = True<br>
res1@cnLevelSelectionMode = "ManualLevels"<br>
res1@cnMinLevelValF = 0<br>
res1@cnMaxLevelValF = 5000<br>
res1@cnLevelSpacingF = 50<br>
res1@gsnDraw = False ; do not draw
picture<br>
res1@gsnFrame = False ; do not
advance frame<br>
<br>
res1@cnLinesOn = False ; contour
lines<br>
res1@cnLineLabelsOn = False ; no contour
labels<br>
res1@gsnSpreadColors = True ; use total
colormap<br>
res1@gsnSpreadColorStart = 2<br>
res1@gsnSpreadColorEnd = 27<br>
res1@cnInfoLabelOn = False ; no contour
info (right bottom)<br>
res1@lbLabelBarOn = False<br>
<br>
res1@vpHeightF = 0.2<br>
res1@vpWidthF = 0.2<br>
<br>
;****************************************************end resources<br>
<br>
isoplot =
gsn_csm_xy2(wks,xiso,brockplot1(25:77),icetplot1(25:77),resL,resR) <br>
<br>
mapannotation =
gsn_csm_contour_map(wks,Hicetfmask(45:226,1:140),res1)<br>
<br>
ypts = (/lat(145,25),lat(145,77)/)<br>
xpts = (/lon(145,25),lon(145,77)/)<br>
<br>
res2 = True<br>
res2@gsLineColor = "black"<br>
res2@gsLineThicknessF = 3.0<br>
res2@gsLineLabelString = "A B"<br>
<br>
dum = new(1,graphic)<br>
dum = gsn_add_polyline(wks,mapannotation,xpts(0:1),ypts(0:1),res2)<br>
<br>
amres = True<br>
amres@amParallelPosF = -1.15<br>
amres@amOrthogonalPosF = -0.4<br>
amres@amJust = "CenterLeft"<br>
<br>
map_final = gsn_add_annotation(isoplot,mapannotation,amres)<br>
<br>
gsn_panel(wks,map_final,(/1,1/),False)<br>
<br>
delete(resL)<br>
delete(resR)<br>
delete(res1)<br>
delete(res2) <br>
delete(amres)<br>
----------------------------------------------------------------------------------------------------------------------------------<br>
<br>
Thank you for you help,<br>
<br>
Michele<br>
<pre class="moz-signature" cols="72">--
***
Michele Petrini
Ph.D. student in Earth Science and Fluid Mechanics
Università degli studi di Trieste,
Dipartimento di Matematica e Geoscienze
Palazzina C - via Weiss 1, 34128 Trieste, Italy
Email: <a class="moz-txt-link-abbreviated" href="mailto:mpetrini139@yahoo.it">mpetrini139@yahoo.it</a>
Skype: michele.petrins
Mobile: +39 3398367372
</pre>
</body>
</html>