
<div dir="ltr"><div class="gmail_default" style="font-size:small">Michele,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I think the error might be coming from trying to panel the &quot;id&quot; created from calling gsn_add_annotation, rather than the actual plot. Try this:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8000001907349px">Â  gsn_panel(wks,isoplot,(/1,1/</span><span style="font-size:12.8000001907349px">),False)</span><br></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8000001907349px"><br></span></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8000001907349px">You might also look into maximize_output</span></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8000001907349px"><br></span></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8000001907349px">--Mary</span></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8000001907349px"><br></span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jul 8, 2015 at 5:58 AM, Michele Petrini <span dir="ltr">&lt;<a href="mailto:mpetrini139@yahoo.it" target="_blank">mpetrini139@yahoo.it</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">


  <div 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&#39;t a resource in this object<br>


      warning:NhlGetValues:Error retrieving pmOverlaySequenceIds<br>


      fatal:[&quot;Execute.c&quot;:8578]:Execute: Error occurred at or near line


      1927 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl<br>


      <br>


      fatal:[&quot;Execute.c&quot;:8578]:Execute: Error occurred at or near line


      6847 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl<br>


      <br>


      fatal:[&quot;Execute.c&quot;:8578]:Execute: Error occurred at or near line


      7546 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl<br>


      <br>


      fatal:[&quot;Execute.c&quot;: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Â Â Â Â Â Â Â Â Â  = &quot;Variation in bedrock elevation


    (30-0 kyrs BP) (m)&quot;<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Â Â Â  = &quot;Always&quot;<br>


    Â  resL@lgPerimOnÂ Â Â Â Â Â Â Â Â Â Â Â Â  = False<br>


    Â  resL@xyExplicitLegendLabels = &quot;Topography variation (m)&quot;<br>


    Â  resL@pmLegendSideÂ Â Â Â Â Â Â Â Â Â  = &quot;Bottom&quot;Â Â Â  <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Â Â Â Â Â Â Â Â Â Â Â Â Â Â  =


    &quot;Explicit&quot;Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 


    ; 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Â Â Â Â Â Â Â Â Â Â Â Â  = (/&quot;A&quot;,&quot;B&quot;/)Â Â Â Â Â Â Â Â Â Â Â Â Â Â  <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Â Â Â Â Â Â Â Â Â  = &quot;Variation in ice thickness (30-0


    kyrs BP) (m)&quot;<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Â Â Â  = &quot;Always&quot;<br>


    Â Â  resR@lgPerimOnÂ Â Â Â Â Â Â Â Â Â Â Â Â  = False<br>


    Â Â  resR@xyExplicitLegendLabels = &quot;Ice thickness variation (m)&quot;<br>


    Â Â  resR@pmLegendSideÂ Â Â Â Â Â Â Â Â Â  = &quot;Bottom&quot;<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Â Â Â Â Â Â Â Â Â Â Â  = &quot;Corners&quot;Â Â Â Â Â Â Â Â Â Â Â  ; 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Â Â Â Â Â Â Â Â Â Â  = &quot;LambertEqualArea&quot;<br>


    Â  res1@mpCenterLatFÂ Â Â Â Â Â Â Â Â Â  = 90<br>


    Â  res1@mpCenterLonFÂ Â Â Â Â Â Â Â Â Â  = 0<br>


    Â  res1@mpOutlineOnÂ Â Â Â Â Â Â Â Â Â Â  = True<br>


    Â  res1@tfDoNDCOverlayÂ Â Â Â Â Â Â Â  = True<br>


    Â  <br>


    Â  res1@cnFillOnÂ Â Â Â Â Â Â Â Â Â Â Â Â Â  = True<br>


    Â  res1@cnLevelSelectionModeÂ Â  = &quot;ManualLevels&quot;<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Â Â Â Â Â Â  = &quot;black&quot;<br>


    Â  res2@gsLineThicknessFÂ  = 3.0<br>


    Â  res2@gsLineLabelString = &quot;A B&quot;<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Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â  = &quot;CenterLeft&quot;<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 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 href="mailto:mpetrini139@yahoo.it" target="_blank">mpetrini139@yahoo.it</a>


Skype: michele.petrins


Mobile: <a href="tel:%2B39%203398367372" value="+393398367372" target="_blank">+39 3398367372</a> 


</pre>


  </div>


<br>_______________________________________________<br>


ncl-talk mailing list<br>


<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>


List instructions, subscriber options, unsubscribe:<br>


<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>


<br></blockquote></div><br></div>