<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p style="margin-top: 0px; margin-bottom: 0px;">Thank you Adam,<br>
</p>
<p style="margin-top: 0px; margin-bottom: 0px;">I found another way to get the plot. But I appreciate your ideas here. Thanks very much.<br>
</p>
<p style="margin-top: 0px; margin-bottom: 0px;">Yun<br>
</p>
<div style="color: rgb(33, 33, 33);">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Adam Phillips <asphilli@ucar.edu><br>
<b>Sent:</b> Friday, March 6, 2015 12:09<br>
<b>To:</b> Yun Yue<br>
<b>Cc:</b> ncl-talk@ucar.edu<br>
<b>Subject:</b> Re: [ncl-talk] add color dots to map</font>
<div> </div>
</div>
<div>
<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>
</div>
</div>
</div>
</body>
</html>