<div dir="ltr"><div class="gmail_default" style="font-size:small">Zhen,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The reason why this is so hard to debug is that I can't run your program in order to print out values before they are passed to the plotting functions.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The marker values you printed all have the level equal to 85000, but the range of your data is 100 to 1000. If you really are trying to plot the marker at value 85,000, then this is definitely out-of-range. This value needs to be converted to hPa (multiply by 0.01) before you try to add these markers.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The thing I don't understand is why you have other markers on the plot that are showing up. Are these markers in hPa units?</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 13, 2017 at 2:47 PM, 刘振 <span dir="ltr"><<a href="mailto:286909655@qq.com" target="_blank">286909655@qq.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>Hi Mary,</div><div><br></div><div> Thanks for your reply. But I am sorry I don't have the authority to copy the data. Maybe I have to leave this behind. Thanks for your really helpful tips and advices. I really appreciate it. </div><div><br></div><div>Best regards,</div><div>zhen</div><div><br></div><div><div style="color:#909090;font-family:Arial Narrow;font-size:12px">------------------</div><div style="font-size:14px;font-family:Verdana;color:#000"><span class=""><div>刘振</div>
<div>中山大学<br>环境科学与工程学院<br>大气科学系<br>Phone: <a href="tel:+86%20150%201324%206049" value="+8615013246049" target="_blank">+86-15013246049</a></div>
</span><span class=""><div>Liu Zhen</div>
<div>Department of Atmospheric Science<br>School of Environmental Science and Engineering</div>
<div>Sun Yat-sen University<br>Email address: <font color="#7f7f7f"><a href="mailto:liuzhen9@mail2.sysu.edu.cn" target="_blank">liuzhen9@mail2.sysu.edu.cn</a></font> </div></span></div></div><div> </div><div><div><br></div><div><br></div><div style="font-size:12px;font-family:Arial Narrow;padding:2px 0 2px 0">------------------ 原始邮件 ------<wbr>------------</div><div style="font-size:12px;background:#efefef;padding:8px"><div><b>发件人:</b> "Mary Haley";<<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>>;</div><div><b>发送时间:</b> 2017年11月14日(星期二) 凌晨5:44</div><div><b>收件人:</b> "刘振"<<a href="mailto:286909655@qq.com" target="_blank">286909655@qq.com</a>>;</div><div><b>抄送:</b> "ncl-talk"<<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.<wbr>edu</a>>; </div><div><b>主题:</b> Re: [ncl-talk] 回复：Re: didn't show all markers</div></div><div><div class="h5"><div><br></div><div dir="ltr"><div class="gmail_default" style="font-size:small">Zhen,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The only way I can look at this is if you provide your full script and data, so I can run it here.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If the data file is not large, then can use our ftp:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/ftp_files.shtml" target="_blank">http://www.ncl.ucar.edu/ftp_<wbr>files.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">If it is large (say > 500 MB), then please provide the file somewhere that I can download it.</div><div class="gmail_default"><br></div><div class="gmail_default">Thank you,</div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div><div class="gmail_default"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Nov 10, 2017 at 12:23 PM, 刘振 <span dir="ltr"><<a href="mailto:286909655@qq.com" target="_blank">286909655@qq.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Mary,<br><br> Thank you so much for your c<wbr>leaning tips (I have changed t<wbr>he script related with tip 1 i<wbr>n the last email based on your<wbr> tips before), which are very <wbr>useful. Sometimes I indeed <wbr>find the script is quite slow.<wbr> Yes, I do still have <wbr>problems with the dots. The <wbr>following figure shows parts <wbr>of the dots are not shown on <wbr>the plot (print("lat / plev = <wbr>"+lat(indices(ip,1))+"/"+plev(<wbr>indices(ip,0)))). The dots <wbr>are not out of the range of <wbr>the X and Y axes. I still <wbr>cannot figure out what is <wbr>going on and need your further<wbr> help. Thanks.<br>(0) lat / plev = -25/85000<br>(0) lat / plev = -23/85000<br>(0) lat / plev = -21/85000<br>(0) lat / plev = -19/85000<br>(0) lat / plev = -17/85000<br>(0) lat / plev = -15/85000<br>(0) lat / plev = -11/85000<br>(0) lat / plev = -3/85000<br>(0) lat / plev = -1/85000<br>(0) lat / plev = 3/85000<br>(0) lat / plev = 5/85000<br>(0) lat / plev = 7/85000<br>(0) lat / plev = 13/85000<br>(0) lat / plev = 15/85000<br>(0) lat / plev = 27/85000<br>(0) lat / plev = 45/85000<br>(0) lat / plev = 47/85000<br>(0) lat / plev = 59.000000<wbr>00000001/85000<br>(0) lat / plev = 61/85000<br>(0) lat / plev = 69/85000<br>(0) lat / plev = 73/85000<br>(0) lat / plev = 75/85000<br>(0) lat / plev = 77/85000<br><br>Best regards,<br>zhen<br><br><br>------------------刘振<span><br> 中山大学<br>环境科学与工程学院<br>大气科学系<br>Phone: <a value="+8615013246049">+86-15013246049</a><br> Liu Zhen<br> Department of Atmospheric Sci<wbr>ence<br>School of Environmental Scienc<wbr>e and Engineering<br> Sun Yat-sen University<br>Email address: <a href="mailto:liuzhen9@mail2.sysu.edu.cn" target="_blank">liuzhen9@mail2.<wbr>sysu.edu.cn</a><br><br></span>---原始邮件---<br>发件人:"Mary Haley"<<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.ed<wbr>u</a>><br>发送时间:2017年11月11日(星期六) 凌晨1:33<br>收件人:"刘振"<<a href="mailto:286909655@qq.com" target="_blank">286909655@qq.com</a>>;<br>主题: [ncl-talk] didn't show<wbr> all markers<br><br><br>Zhen,<br><br>Which markers in particular ar<wbr>e not showing up? It's too<wbr> hard for us to look at your <wbr>output and try to figure out <wbr>which dots are on the plot <wbr>and which are not.<br><br>Also, it's very unlikely t<wbr>hat NCL would just drop some o<wbr>f those dots.  There has <wbr>to be something else going on,<wbr> like the dots are out-of-<wbr>range of the limits of your X <wbr>and/or Y axes.<span><br><br>[1] Clean coding tip: you do n<wbr>ot need to loop across each in<wbr>dex of the array to add each <wbr>dot, if you are drawing them <wbr>all with the same color and <wbr>size. Instead of this code:<br><br>dot=new((/dimsizes(forc),dimsi<wbr>zes(lat)*dimsizes(plev)/),<wbr>graphic)<br>do iex=0,dimsizes(forc)-1<br></span> . . .<br><br>  do ip=0,npts-1<br>    dot(iex,ip) = gs<wbr>n_add_polymarker(wks,plot(iex)<wbr>,lat(indices(ip,1)),plev(indic<wbr>es(ip,0)),resmarker)<br>  end do<span><br><br>you can have this code:<br><br>dot=new((/dimsizes(forc),dimsi<wbr>zes(lat)/),graphic)<br>do iex=0,dimsizes(forc)-1<br></span> . . .<br>  dot(iex) = gsn_add_poly<wbr>marker(wks,plot(iex),lat(indic<wbr>es(:,1)),plev(indices(:,0)),<wbr>resmarker)<span><br><br><br>This will run faster, because <wbr>each call to gsn_add_polymarke<wbr>r causes the creation of an <wbr>internal object, and if you <wbr>do this for each single dot, <wbr>you are going to end up using <wbr>more memory and taking more <wbr>time.<br><br>[2] Clean coding tip: if you f<wbr>ind yourself calling "dimsizes<wbr>" on the same variable <wbr>multiple times, I suggest <wbr>saving it to a variable so it <wbr>can save you some typing:<br><br>For example:<br><br>var_run = new((/dimsizes(forc)<wbr>,dimsizes(ens),dimsizes(year),<wbr>23,90/),"float")<br></span>var_djf_clim =  new((/dim<wbr>sizes(forc),dimsizes(ens),23,<wbr>90/),"float") ; 10 years <wbr>         &<wbr>#xA0;       &#<wbr>xA0;       &#<wbr>xA0;  <span><br>rc_m = new((/dimsizes(forc),23<wbr>,90/),"float")<br>do iv=0,dimsizes(vars)-1<br><br>could be changed to:<br><br>nforc = dimsizes(forc)<br></span>nens  = dimsizes(ens)<span><br>nyear = dimsizes(year)<br>nvars = dimsizes(vars)<br>var_run = new((/nforc,nens,nye<wbr>ar),23,90/),"float")<br></span>var_djf_clim =  new((/nfo<wbr>rc,nens),23,90/),"float") ; <wbr>10 years      &<wbr>#xA0;       &#x<wbr>A0;        <wbr>;      <span><br>rc_m = new((/nforc,23,90/),"fl<wbr>oat")<br>do iv=0,nvars-1<br></span>I don't know what "23" and<wbr> "90" are supposed to represen<wbr>t, but it would be good to <wbr>also assign these to variables<wbr> as well, so if you need to <wbr>change the 23 or 90 values <wbr>later, you can just do it in <wbr>one place.<br>[3] Clean coding tip: you don&<wbr>#39;t need to use "new" to pre<wbr>allocate an array that is gett<wbr>ing calculated by a function.&<wbr>#xA0; You have this code:<span><br>rc_m = new((/dimsizes(forc),23<wbr>,90/),"float")<br></span>rc_m = regCoef_n(year,var<wbr>_ens,0,1)<br>Remove the "new" line as it is<wbr> not needed.  The regCoef<wbr>_n function will allocate the <wbr>memory for you, as does every <wbr>single NCL function.<span><br>If you continue to have proble<wbr>ms with the dots, please be <wbr>more specific about which <wbr>dots are not showing up.<br>--Mary<br><br><br>On Thu, Nov 9, 2017 at 2:56 PM<wbr>, 刘振 <<a href="mailto:286909655@qq.com" target="_blank">286909655@qq.com</a>> wrote:<br><br><br>Dear all,<br><br></span>  I am using ncl 6.4 to p<wbr>lot and want to print dot when<wbr> value pass student t test <wbr>using function gsn_add_polymak<wbr>er (script and figure attached<wbr>, i.e. 850hPa). But the dots <wbr>are only printed at some <wbr>levels (not all levels) in <wbr>vertical direction. Actually, <wbr>I got 23 pressure levels for <wbr>my data and indeed some dot <wbr>are not printed on the plot <wbr>when I print the indices <wbr>passing significance level (<wbr>see log files attached). <wbr>Could you please help me to so<wbr>lve this problem. Thanks.<span><br><br><br>zhen<br><br>------------------<br>刘振<br> 中山大学<br>环境科学与工程学院<br>大气科学系<br>Phone: <a value="+8615013246049">+86-15013246049</a><br></span> Liu  Zhen<span><br> Department of Atmospheric Sci<wbr>ence<br>School of Environmental Scienc<wbr>e and Engineering<br> Sun Yat-sen University<br></span>Email address: <a href="mailto:liuzhen9@mail2.sysu.edu.cn" target="_blank">liuzhen9@mail2.<wbr>sysu.edu.cn</a> <br><br><br> <div class="m_-5422580639854949200HOEnZb"><div class="m_-5422580639854949200h5"><br>______________________________<wbr>_________________<br> ncl-talk mailing list<br> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br> List instructions, subscriber<wbr> options, unsubscribe:<br> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailm<wbr>an/listinfo/ncl-talk</a><br> <br><br><br></div></div><br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div></div></div></div></blockquote></div><br></div>