<div dir="ltr"><div>Hi,</div><div><br></div><div>I have used gsn-coordinates to successfully draw a lon/lat grid over an island. </div><div><br></div><div>However, when I am trying to mark a station location using gsn_polymarker, the marker is placed elsewhere. I get the same result using gsn_add_polymarker. There is no error message, script runs well - it&#39;s just that the resultant image is not what I was trying to get. </div><div>Could you please let me know why this is happening, and how the problem can be fixed?</div><div><br></div><div>My script is copied below, and the resultant image (plus data file) are attached. </div><div><br></div><div>(NCL 6.1.0 on Ubuntu 12.0.4)</div><div><br></div><div>;--------------------------------------------------</div><div>; mapgrid_CCAM.ncl</div><div>;--------------------------------------------------</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;  </div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;  </div><div><br></div><div>begin</div><div><br></div><div>  inFile=&quot;pr_Fiji_1980.nc&quot; ; Read netCDF file</div><div>  </div><div>  a = addfile(inFile,&quot;r&quot;)</div><div>  prc  = a-&gt;rnd(0,:,:) </div><div>   </div><div>  lat = a-&gt;lat</div><div>  lon = a-&gt;lon</div><div>   </div><div>  lon1 = 177.26</div><div>  lat1 = -17.45</div><div>  </div><div>;---Start the graphics.</div><div>  wks = gsn_open_wks(&quot;eps&quot;,&quot;Fijigrid&quot;)           ; open an eps file</div><div><br></div><div>  res               = True                        ; plot mods desired</div><div><br></div><div>  res@gsnDraw       = False                       ; don&#39;t draw yet</div><div>  res@gsnFrame      = False                       ; don&#39;t advance frame</div><div>  res@tiMainString  = &quot;CCAM 8 km grid&quot;            ; Title</div><div>  </div><div>  res@mpMinLonF            =  177.1               ; set maximum and minimum lon lat values     </div><div>  res@mpMaxLonF            =  178.9</div><div>  res@mpMinLatF            =  -18.4</div><div>  res@mpMaxLatF            =  -17.10</div><div> </div><div>  mres = True                                      ; set marker resources</div><div>  mres@gsMarkerIndex = 16</div><div>  mres@gsMarkerColor = &quot;red&quot;</div><div><br></div><div>  res@pmTickMarkDisplayMode = &quot;Always&quot;            ;add tickmark labels</div><div><br></div><div>;---Create map (don&#39;t draw it yet)</div><div>  map = gsn_csm_map_ce(wks,res)</div><div><br></div><div>;---Add marker</div><div>  gsn_polymarker(wks,map,lon1,lat1,mres)</div><div><br></div><div>;---Call function to attach lat/lon lines.</div><div>  </div><div>  pres                   = True</div><div>  pres@gsnCoordsAsLines  = True        ; Default is points</div><div>  gsn_coordinates(wks,map,prc,pres)</div><div>  </div><div>end</div><div><br></div><div><br></div><div>Thanks,</div><div>Vandhna K.</div><div><br></div></div>