<div dir="ltr"><div><br></div><div>Hi all, </div><div><br></div><div>I&#39;m trying to add the location of a city and it&#39;s name in the same plot of a contour-terrain heights drawn over a map.</div><div><br></div><div>I&#39;m editing the  <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&#39;t shows me the polymarker and the name of the city. </div><div><br></div><div>If someone see my mistake, I would really aprecciate his help.<br><br>Best regards,<br>Sebastián. </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 &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br>load &quot;$NCARG_ROOT<i>/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl&quot;</i><br>begin<br> a=addfile(&quot;/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>&quot;,&quot;r&quot;)<br><br>  type = &quot;pdf&quot;<br>  wks = gsn_open_wks(type,&quot;Norte&quot;)         ; Create a plot workstation<br><br>  opts = True                                  ; Set some Basic Plot options<br>  opts@MainTitle = &quot;GEOGRID FIELDS&quot;<br>  opts@InitTime = False                        ; Do not plot time or footers<br>  opts@Footer = False<br>  ter = wrf_user_getvar(a,&quot;HGT_M&quot;,0)           ; Read the variable to memory<br>  lat  = wrf_user_getvar(a,&quot;lat&quot;,-1)<br>  lon  = wrf_user_getvar(a,&quot;lon&quot;,-1)<br>  opt = True<br>  loc1 = wrf_user_ll_to_ij(a,-70.63929,-33.35079,opt)  ; Chanaral<br>  loc1 = loc1-1;<br>  res = opts            <br> res@cnFillOn = True                          ; Create a color fill plot<br>  res@ContourParameters = (/ 0., 5000., 50. /) ; Set the levels<br>  <br>  contour = wrf_contour(a,wks,ter,res)<br>  pltres = True<br>  pltres@PanelPlot = True                             ; Set plot options<br>  mpres = True                                 ; Set map options<br>  mpres@mpGeophysicalLineColor      = &quot;Black&quot;  ; Overwrite basic map settings<br>  mpres@mpGridLineColor             = &quot;Black&quot;<br>  mpres@mpLimbLineColor             = &quot;Black&quot;<br>  mpres@mpNationalLineColor         = &quot;Black&quot;<br>  mpres@mpPerimLineColor            = &quot;Black&quot;<br>  mpres@mpUSStateLineColor          = &quot;Black&quot;<br>  mpres@gsnDraw = False<br>  mpres@gsnFrame =False<br>  plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres) ; Plot field over map background<br> ; Poly res desired<br>  resp = True<br>  resp@gsMarkerColor = &quot;Black&quot;<br>  resp@gsMarkerIndex = 16 ; black dot<br>  resp@gsMarkerSizeF = 0.01<br>  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>   txres= True<br>   txres@txFontHeightF = 0.015 ; text font height<br>   gsn_text(wks,plot,&quot;Chanaral&quot;,lon(0,loc1(1),loc1(0)),lat(0,loc1(1),loc1(0)),txres)<br>   draw(plot)<br>   frame(wks)<br>end<br>    </blockquote></div></div>