<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body dir="auto"><div>Sebastian,</div><div><br></div><div>Try adding pltres@FramePlot = False</div><div><br></div><div>I think it is still framing the plot before you add the annotations, gsnDraw and gsnFrame are more related to the plot than the map. &nbsp;See the link below for where the information comes from. I think this is the only change you need in your script.</div><div><br></div><div><a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_map_overlays.shtml">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_map_overlays.shtml</a></div><div><br></div><div><br></div><div>Hope that helps,</div><div>-Alex</div><div><br></div><div><br><br><br></div><div><br>On Apr 29, 2015, at 7:21 PM, Sebastian Otarola-Bustos &lt;<a href="mailto:Sebastian.F.Otarola-Bustos.1@nd.edu">Sebastian.F.Otarola-Bustos.1@nd.edu</a>&gt; wrote:<br><br></div><blockquote type="cite"><div><div dir="ltr"><div><br></div><div>Hi all,&nbsp;</div><div><br></div><div>I'm trying to add the location of a city and it's name in the same plot of a contour-terrain heights drawn over a map.</div><div><br></div><div>I'm editing the &nbsp;<b style="font-style:italic">wrf_wps_ter6.ncl </b>example script. And it works fine plotting the map and the topography, but it doesn't shows me the polymarker and the name of the city.&nbsp;</div><div><br></div><div>If someone see my mistake, I would really aprecciate his help.<br><br>Best regards,<br>Sebastián.&nbsp;</div><div><br><br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT<i>/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"</i><br>begin<br>&nbsp;a=addfile("/afs/<a href="http://crc.nd.edu/user/s/sotarola/WRF/Santiago/Sim/1/geo_em.d01.nc">crc.nd.edu/user/s/sotarola/WRF/Santiago/Sim/1/geo_em.d01.nc</a>","r")<br><br>&nbsp; type = "pdf"<br>&nbsp; wks = gsn_open_wks(type,"Norte") &nbsp; &nbsp; &nbsp; &nbsp; ; Create a plot workstation<br><br>&nbsp; opts = True &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; Set some Basic Plot options<br>&nbsp; opts@MainTitle = "GEOGRID FIELDS"<br>&nbsp; opts@InitTime = False &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; Do not plot time or footers<br>&nbsp; opts@Footer = False<br>&nbsp; ter = wrf_user_getvar(a,"HGT_M",0) &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; Read the variable to memory<br>&nbsp; lat &nbsp;= wrf_user_getvar(a,"lat",-1)<br>&nbsp; lon &nbsp;= wrf_user_getvar(a,"lon",-1)<br>&nbsp; opt = True<br>&nbsp; loc1 = wrf_user_ll_to_ij(a,-70.63929,-33.35079,opt) &nbsp;; Chanaral<br>&nbsp; loc1 = loc1-1;<br>&nbsp; res = opts &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;<br>&nbsp;res@cnFillOn = True &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; Create a color fill plot<br>&nbsp; res@ContourParameters = (/ 0., 5000., 50. /) ; Set the levels<br>&nbsp;&nbsp;<br>&nbsp; contour = wrf_contour(a,wks,ter,res)<br>&nbsp; pltres = True<br>&nbsp; pltres@PanelPlot = True &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; Set plot options<br>&nbsp; mpres = True &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; Set map options<br>&nbsp; mpres@mpGeophysicalLineColor &nbsp; &nbsp; &nbsp;= "Black" &nbsp;; Overwrite basic map settings<br>&nbsp; mpres@mpGridLineColor &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = "Black"<br>&nbsp; mpres@mpLimbLineColor &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = "Black"<br>&nbsp; mpres@mpNationalLineColor &nbsp; &nbsp; &nbsp; &nbsp; = "Black"<br>&nbsp; mpres@mpPerimLineColor &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= "Black"<br>&nbsp; mpres@mpUSStateLineColor &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= "Black"<br>&nbsp; mpres@gsnDraw = False<br>&nbsp; mpres@gsnFrame =False<br>&nbsp; plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres) ; Plot field over map background<br>&nbsp;; Poly res desired<br>&nbsp; resp = True<br>&nbsp; resp@gsMarkerColor = "Black"<br>&nbsp; resp@gsMarkerIndex = 16 ; black dot<br>&nbsp; resp@gsMarkerSizeF = 0.01<br>&nbsp; gsn_polymarker(wks,plot,lon(0,loc1(1),loc1(0)),lat(0,loc1(1),loc1(0)),resp) ;<br>; Add the name of the station<br>&nbsp; &nbsp;txres= True<br>&nbsp; &nbsp;txres@txFontHeightF = 0.015 ; text font height<br>&nbsp; &nbsp;gsn_text(wks,plot,"Chanaral",lon(0,loc1(1),loc1(0)),lat(0,loc1(1),loc1(0)),txres)<br>&nbsp; &nbsp;draw(plot)<br>&nbsp; &nbsp;frame(wks)<br>end<br>&nbsp;&nbsp; &nbsp;</blockquote></div></div>
</div></blockquote><blockquote type="cite"><div><span>_______________________________________________</span><br><span>ncl-talk mailing list</span><br><span>List instructions, subscriber options, unsubscribe:</span><br><span><a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></span><br></div></blockquote></body></html>