<div dir="ltr"><div><div><div><div><div>Hi Barry,<br><br></div>I think you can try the following changes (note ncities need be defined before these lines):<br> <br>  markid = new((/4,ncities/),graphic)<br>  textid = new((/4,ncities/),graphic)<br><br></div>and in the loop:<br><br>    do n = 0,3<br>      markid(n,:) = gsn_add_polymarker(wks,plot(n)<wbr>,lon,lat,mkres)<br>      textid(n,:) = gsn_add_text(wks,plot(n),<wbr>cities,lon,lat,txres)<br>    end do<br><br><br></div>A simplified version of your script that works is attached (just maps without data):<br><br></div>HTH<br></div>Rashed<br><div><div><br></div></div></div>