Re: incorrect plotting of polymarkers

I think you want to do this:

nlat = dimsizes(lat)
lat_pts = new(nlat,typeof(lat2d))
lon_pts = new(nlat,typeof(lon2d))
do k=0,nlat-1
     n = nm(k,0)
     m = nm(k,1)
     lat_pts(k) = lat2d(n,m)
     lon_pts(k) = lon2d(n,m)
  end do

  dum = gsn_add_polymarker(wks,map,lon_pts,lat_pts, mres)

On Nov 26, 2012, at 10:57 AM, Kantave Greene wrote:

> Hello
>> I am having a problem with polymarkers being plotting incorrectly on a lambert conformal map although there indices are valid and correct. The attached image is the output from the script.
>> Based on the user specified lat/lon, getind_latlon2d returns the indices which polymarker then plots. Using the following statement:
>> print(n+" "+m+" "+lat(n,m)+" "+lon(n,m))
>> I confirmed that the indices are returning the nearest lat/lon to the indices.
>> Can you help me get the script to plot the indices correctly?
>>
>> Please find attached a copy of the code and output file I am using. Thanks for all your help.
>>
>
> --
> Kantave Greene
> Instructor
> Meteorology Department
> Jackson State University
> P.O. Box 17660
> Jackson, MS 39217
> (601) 979-1293
>
> <plot_polygon.ncl><polygon_orig.ps>_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Nov 26 14:25:47 2012