<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;"><div>刘振</div>
<div>中山大学<br>环境科学与工程学院<br>大气科学系<br>Phone: +86-15013246049</div>
<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">liuzhen9@mail2.sysu.edu.cn</font> </div></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;">------------------ 原始邮件 ------------------</div><div style="font-size: 12px;background:#efefef;padding:8px;"><div><b>发件人:</b> "Mary Haley";<haley@ucar.edu>;</div><div><b>发送时间:</b> 2017年11月14日(星期二) 凌晨5:44</div><div><b>收件人:</b> "刘振"<286909655@qq.com>;<wbr></div><div><b>抄送:</b> "ncl-talk"<ncl-talk@ucar.edu>; <wbr></div><div><b>主题:</b> Re: [ncl-talk] 回复：Re: didn't show all markers</div></div><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">http://www.ncl.ucar.edu/ftp_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 <wbr>cleaning tips (I have changed <wbr>the script related with tip 1 <wbr>in the last email based on <wbr>your tips before), which are <wbr>very useful. Sometimes I <wbr>indeed find the script is <wbr>quite slow. Yes, I do still <wbr>have problems with the dots. <wbr>The following figure shows <wbr>parts of the dots are not <wbr>shown on the plot (print("lat <wbr>/ plev = "+lat(indices(ip,1))+<wbr>"/"+plev(indices(ip,0)))). <wbr>The dots are not out of the <wbr>range of the X and Y axes. I <wbr>still cannot figure out what <wbr>is going on and need your <wbr>further 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.<wbr>00000000000001/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 class=""><br> 中山大学<br>环境科学与工程学院<br>大气科学系<br>Phone: <a value="+8615013246049" target="_blank">+86-15013246049</a><br> Liu Zhen<br> Department of Atmospheric <wbr>Science<br>School of Environmental <wbr>Science 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.<wbr>edu</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 <wbr>show all markers<br><br><br>Zhen,<br><br>Which markers in particular <wbr>are not showing up? It's <wbr>too hard for us to look at <wbr>your output and try to figure <wbr>out which dots are on the <wbr>plot and which are not.<br><br>Also, it's very unlikely <wbr>that NCL would just drop some <wbr>of those dots.  There <wbr>has to be something else <wbr>going on, like the dots are <wbr>out-of-range of the limits of <wbr>your X and/or Y axes.<span class=""><br><br>[1] Clean coding tip: you do <wbr>not need to loop across each <wbr>index of the array to add <wbr>each dot, if you are drawing <wbr>them all with the same color <wbr>and size. Instead of this <wbr>code:<br><br>dot=new((/dimsizes(forc),<wbr>dimsizes(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) = <wbr>gsn_add_polymarker(wks,plot(<wbr>iex),lat(indices(ip,1)),plev(<wbr>indices(ip,0)),resmarker)<br>  end do<span class=""><br><br>you can have this code:<br><br>dot=new((/dimsizes(forc),<wbr>dimsizes(lat)/),graphic)<br>do iex=0,dimsizes(forc)-1<br></span> . . .<br>  dot(iex) = gsn_add_<wbr>polymarker(wks,plot(iex),lat(<wbr>indices(:,1)),plev(indices(:,<wbr>0)),resmarker)<span class=""><br><br><br>This will run faster, because <wbr>each call to gsn_add_<wbr>polymarker causes the <wbr>creation of an internal <wbr>object, and if you do this <wbr>for each single dot, you are <wbr>going to end up using more <wbr>memory and taking more time.<br><br>[2] Clean coding tip: if you <wbr>find yourself calling "<wbr>dimsizes" on the same <wbr>variable multiple times, I <wbr>suggest saving it to a <wbr>variable so it can save you <wbr>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((/<wbr>dimsizes(forc),dimsizes(ens),<wbr>23,90/),"float") ; 10 years&#<wbr>xA0;       &#<wbr>xA0;       &#<wbr>xA0;       &#<wbr>xA0;    <span class=""><br>rc_m = new((/dimsizes(forc),<wbr>23,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 class=""><br>nyear = dimsizes(year)<br>nvars = dimsizes(vars)<br>var_run = new((/nforc,nens,<wbr>nyear),23,90/),"float")<br></span>var_djf_clim =  new((/<wbr>nforc,nens),23,90/),"float") ;<wbr> 10 years      <wbr>        &#<wbr>xA0;       &#<wbr>xA0;      <span class=""><br>rc_m = new((/nforc,23,90/),"<wbr>float")<br>do iv=0,nvars-1<br></span>I don't know what "23" <wbr>and "90" are supposed to <wbr>represent, but it would be <wbr>good to also assign these to <wbr>variables as well, so if you <wbr>need to change the 23 or 90 <wbr>values later, you can just do <wbr>it in one place.<br>[3] Clean coding tip: you don&<wbr>#39;t need to use "new" to <wbr>preallocate an array that is <wbr>getting calculated by a <wbr>function.  You have this <wbr>code:<span class=""><br>rc_m = new((/dimsizes(forc),<wbr>23,90/),"float")<br></span>rc_m = regCoef_n(year,<wbr>var_ens,0,1)<br>Remove the "new" line as it <wbr>is not needed.  The <wbr>regCoef_n function will <wbr>allocate the memory for you, <wbr>as does every single NCL <wbr>function.<span class=""><br>If you continue to have <wbr>problems with the dots, <wbr>please be more specific about <wbr>which dots are not showing up.<br>--Mary<br><br><br>On Thu, Nov 9, 2017 at 2:56 <wbr>PM, 刘振 <<a href="mailto:286909655@qq.com" target="_blank">286909655@qq.com</a>> <wbr>wrote:<br><br><br>Dear all,<br><br></span>  I am using ncl 6.4 to <wbr>plot and want to print dot <wbr>when value pass student t <wbr>test using function gsn_add_<wbr>polymaker (script and figure <wbr>attached, i.e. 850hPa). But <wbr>the dots are only printed at <wbr>some levels (not all levels) <wbr>in vertical direction. <wbr>Actually, I got 23 pressure <wbr>levels for my data and indeed <wbr>some dot are not printed on <wbr>the plot when I print the <wbr>indices passing significance <wbr>level (see log files attached)<wbr>. Could you please help me to <wbr>solve this problem. Thanks.<span class=""><br><br><br>zhen<br><br>------------------<br>刘振<br> 中山大学<br>环境科学与工程学院<br>大气科学系<br>Phone: <a value="+8615013246049" target="_blank">+86-15013246049</a><br></span> Liu  Zhen<span class=""><br> Department of Atmospheric <wbr>Science<br>School of Environmental <wbr>Science 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="HOEnZb"><div class="h5"><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, <wbr>subscriber options, <wbr>unsubscribe:<br> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/<wbr>mailman/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">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/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div></div>