<div dir="ltr"><div><div><div>Hi Yun,<br></div><div>Please always include ncl-talk on your replies so that other ncl-talk folks can assist and also so the conversation gets saved in the ncl-talk archives to assist future NCL users.<br><br></div><div>As my initial suggestion didn't solve the issue I created a sample script that used your code snippet and tried to get it to do what you want. The script and resulting graphics are attached.<br></div><br></div><div>I saw two issues that needed addressing:<br></div><div>1) Unless you created your own custom colormap or mistyped, there isn't a colormap named WhAqGrYeOrReVi200. For my test case I switched to using BkBlAqGrYeOrReViWh200.<br></div><div>2) The GetFillColor documentation states that the function is deprecated, and to use either get_color_rgba or
get_color_index. GetFillColor assigns colors based on a scalar field, but uses the colors in the retrieved colormap sequentially. For example, if you set up your cnLevels array = (/5,10,15,20,25,30/), and are using the BkBlAqGrYeOrReViWh200 colormap, all points less than 5 will be colored index 2 here:<br><a href="http://www.ncl.ucar.edu/Document/Graphics/ColorTables/BkBlAqGrYeOrReViWh200.shtml">http://www.ncl.ucar.edu/Document/Graphics/ColorTables/BkBlAqGrYeOrReViWh200.shtml</a><br><br></div><div>All points greater than or equal to 5 and less than 10 will be colored index 3, and so on up to index 7. For this colormap, all those indices look black.<br><br></div><div>What I think you want to do is to use get_color_index instead of GetFillColor, as get_color_index assigns the colors in a fashion similar to the manner that NCL spans colormaps for color-filled graphics. See the documentation here:<br><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/get_color_index.shtml">http://www.ncl.ucar.edu/Document/Functions/Built-in/get_color_index.shtml</a> <br></div><br></div>Hope that helps!<br><div><div>Adam<br><br>begin<br><br>x = ispan(190,290,20)<br>y = ispan(25,67,6)<br>pmvalue = ispan(2,32,5)<br><br>wks = gsn_open_wks("png", "overlay")<br>gsn_define_colormap(wks,"BkBlAqGrYeOrReViWh200")<br>res = True<br>res@gsnDraw = False<br>res@gsnFrame = False<br>res@mpFillOn = False<br>res@mpOutlineOn = True<br>res@mpLimitMode = "LatLon"<br>res@mpMinLatF = 20.<br>res@mpMaxLatF = 75.<br>res@mpMinLonF = 180.<br>res@mpMaxLonF = 300.<br><br>res@mpCenterLonF = 240.<br>plot = gsn_csm_map(wks,res)<br>cmap = gsn_retrieve_colormap(wks)<br>cnLevels = (/5,10,15,20,25,30/)<br>res_poly = True<br>res_poly@gsMarkerIndex = 16<br>res_poly@gsMarkerSizeF = 10 <br>do i = 0, dimsizes(x)-1<br>res_poly@gsMarkerColor = get_color_index("BkBlAqGrYeOrReViWh200",cnLevels, pmvalue(i))<br>print("Value = "+pmvalue(i)+" "+res_poly@gsMarkerColor)<br>dum = gsn_add_polymarker(wks, plot, x(i), y(i), res_poly)<br>str1 = unique_string("dum")<br>plot@$str1$=dum<br>end do<br>draw (plot)<br>frame(wks)<br>end<br><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Mar 5, 2015 at 3:53 PM, Yun Yue <span dir="ltr"><<a href="mailto:yunyue@huskers.unl.edu" target="_blank">yunyue@huskers.unl.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif">
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Hi NCL users,<br>
</p>
<p style="margin-top:0px;margin-bottom:0px">I'm trying to plot observation data to model simulated data map. I think I should use gsn_add_polymarker to overlay them. But for the @gsMarkerColor, I want use the same color for plotting simulated data map,
so I tried to use code like the following. But I got black dots on the map. Can somebody help me to solve this? Thank you very much.<br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Yun<br>
</p>
<p style="margin-top:0px;margin-bottom:0px">wks = gsn_open_wks("x11", "overlay")<br>
</p>
<p style="margin-top:0px;margin-bottom:0px">gsn_define_colormap(wks,"WhAqGrYeOrReVi200")<br>
</p>
<p style="margin-top:0px;margin-bottom:0px">cmap = gsn_retrieve_colormap(wks)<br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;background-color:rgb(255,255,255)">cnLevels = (/5,10,15,20,25,30</span><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;background-color:rgb(255,255,255)">/</span><span style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:16px;background-color:rgb(255,255,255)">)</span><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">res_poly = True<br>
</p>
<p style="margin-top:0px;margin-bottom:0px">res_poly@gsMarkerIndex = 16<br>
</p>
<p style="margin-top:0px;margin-bottom:0px">res_poly@gsMarkerSizeF = 10 <br>
</p>
<p style="margin-top:0px;margin-bottom:0px">do i = 0, station_numb<br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(255,0,0)">res_poly@gsMarkerColor = GetFillColor(cnLevels, cmap, pmvalue(i))</span><br>
</p>
<div>
<div style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif">
<div name="divtagdefaultwrapper">
<div>dum = gsn_add_polymarker(wks, plot, x, y, res_poly)<br>
</div>
<div>str1 = unique_string("dum"<span style="font-size:12pt">)</span></div>
<div><span style="font-size:12pt">plot@$str1$=dum</span></div>
<div><span style="font-size:12pt">end do</span></div>
<div><span style="font-size:12pt">draw (plot)</span></div>
<div><span style="font-size:12pt"><br>
</span></div>
<div><span style="font-size:12pt"><br>
</span></div>
<div><br>
</div>
<div><span style="font-size:12pt"><br>
</span></div>
</div>
</div>
</div>
</div>
</div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist, </font></span><span><font color="#888888">Climate and Global Dynamics Division, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphilli/</a> </font></span><span><font color="#888888">303-497-1726 </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div>
</div>