<div dir="ltr"><div>Hi all,</div><div><br></div><div>I am struggling with overlaying a single contour plot (called &#39;tmp0&#39; below) over multiple plots in a panel. The script below results in the error &quot;fatal:NhlAddOverlay: tranform is already an annotation or overlay...&quot;</div><div><br></div><div><div>  plot(1)  = gsn_csm_contour(wks,dz_m,res1)</div><div>  plot(3)  = gsn_csm_contour(wks,dy_m,res1)</div><div>  res1@tiXAxisString          = &quot;Latitude&quot;</div><div>  plot(5)  = gsn_csm_contour(wks,dy_e,res1)</div><div>  delete(res1@tiXAxisString)</div><div><br></div><div>  res2@cnFillOn          = False</div><div>  res2@cnLinesOn         = True</div><div>  res2@cnLineThicknessF  = 2.</div><div>  res2@cnLevelSelectionMode = &quot;ManualLevels&quot;</div><div>  res2@cnLevelSpacingF = 5.e10;</div><div>  tmp0  = gsn_csm_contour(wks,psi1,res2)</div><div><br></div><div>  res2@cnLevelSpacingF = 1.e10</div><div>  tmp1  = gsn_csm_contour(wks,dpsi,res2)</div><div><br></div><div>  overlay(plot(0),tmp0)</div><div>  overlay(plot(2),tmp0) ;&lt;-fatal:NhlAddOverlay</div><div>  overlay(plot(4),tmp0) ;&lt;-fatal:NhlAddOverlay</div><div>  overlay(plot(1),tmp1)</div><div>  overlay(plot(3),tmp1) ;&lt;-fatal:NhlAddOverlay</div><div>  overlay(plot(5),tmp1) ;&lt;-fatal:NhlAddOverlay</div></div><div>...<br></div><div>  gsn_panel(wks,plot,(/3,2/),resP)<br></div><div><br></div><div>In terms of finding a solution, I would prefer to keep all the overlays contained in the graphics variable &#39;plot&#39; so that I can easily perform a panel plot. I don&#39;t think the solution provided by Mary Haley in <a href="https://www.ncl.ucar.edu/Support/talk_archives/2013/3802.html">https://www.ncl.ucar.edu/Support/talk_archives/2013/3802.html</a> would work for me, since gsn_panel draws all the plots at once.<br></div><div><br></div><div>Any help is greatly appreciated!</div><div><br></div><div>Adam</div><div><br></div><div><br></div><div><br></div></div>