<div dir="ltr"><div><div>FTR: <a href="https://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclStatements.shtml">https://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclStatements.shtml</a><br><br></div>See the 'do' section.<br><br></div>[snip]<br>The following form is also
necessary<b> to loop backwards</b>. With this form, if start is greater than end,
then the identifier is decremented by the stride value, which must always be
positive. If the start is less than the end, then the identifier is incremented
by the stride value.
<p>
</p><pre> do <i>loop_identifier</i> = <i>scalar_start_expr</i> , <i>scalar_end_expr</i> , <i> scalar_stride_expr</i>
<i> statement list </i>
end do
<br></pre><pre>[snip]<br><br></pre><pre>Given my fortran background, I found this a bit unusual but 'it-is-what-it-is'<br><br></pre><pre>D<br></pre></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Oct 6, 2017 at 11:24 AM, Dave Allured - NOAA Affiliate <span dir="ltr"><<a href="mailto:dave.allured@noaa.gov" target="_blank">dave.allured@noaa.gov</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">Oops. Please ignore my correction. Dennis's version was correct. The stride parameter is different between fortran and NCL.<span class="HOEnZb"><font color="#888888"><div><br></div><div>--Dave</div></font></span><div><div class="h5"><div class="gmail_extra"><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Oct 6, 2017 at 11:21 AM, Dave Allured - NOAA Affiliate <span dir="ltr"><<a href="mailto:dave.allured@noaa.gov" target="_blank">dave.allured@noaa.gov</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">Correction to Dennis's line 1:<br><div><br></div><div><font face="monospace, monospace"><span style="font-size:14px">Old: do m=4,0,1</span><br></font></div><div><div><font face="monospace, monospace"><span style="font-size:14px">New: do m=4,0,-1</span><br></font></div></div><div><span style="font-size:14px"><br></span></div><div><span style="font-size:14px">I can't confirm the rest of the code.</span></div><div><span style="font-size:14px"><br></span></div><div><span style="font-size:14px">--Dave</span></div><div><span style="font-size:14px"><br></span></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Oct 6, 2017 at 11:08 AM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.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">A person sent me the following offline:<br><br>do m=4,0,1<div> R(:,:,:,m) = dim_sum_n_Wrap(Q(:,:,:,m::,m::<wbr>,m::),(/3,4,5/))</div><div> if m.lt.4 then</div><div> R(:,:,:,m) = R(:,:,:,m) - dim_sum_n_Wrap(R(:,:,:,(m+1)::<wbr>),3)</div><div> end if</div><div>end do<br><br><br></div></div><div class="gmail_extra"><div class="gmail_quote">On Thu, Oct 5, 2017 at 5:37 PM, <span dir="ltr"><<a href="mailto:Arnold.Sullivan@csiro.au" target="_blank">Arnold.Sullivan@csiro.au</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div link="#0563C1" vlink="#954F72" lang="EN-AU">
<div class="m_3606792257068817146m_-2806852915113090976m_5214227225033462450m_-7047176370259980277m_4353728014666826418WordSection1">
<p class="MsoNormal">Dear all ncl users,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Could someone help me to convert Matlab code using sum to NCL<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><b><span style="color:#0070c0">R is [time] [lat] [lon] [m1]<u></u><u></u></span></b></p>
<p class="MsoNormal"><b><span style="color:#0070c0"><u></u> <u></u></span></b></p>
<p class="MsoNormal"><b><span style="color:#0070c0">Q is [time] [lat] [lon] [m1] [m2] [m3]<u></u><u></u></span></b></p>
<p class="MsoNormal"><b><span style="color:#0070c0"><u></u> <u></u></span></b></p>
<p class="MsoNormal"><b><span style="color:#0070c0">For m=5:-1:1<u></u><u></u></span></b></p>
<p class="MsoNormal"><b><span style="color:#0070c0"> R(:,:,:,m) = sum( sum( sum( Q(:,:,:,m:end,m:end,m:end), 4), 5), 6);<u></u><u></u></span></b></p>
<p class="MsoNormal"><b><span style="color:#0070c0"> if m < 5<u></u><u></u></span></b></p>
<p class="MsoNormal"><b><span style="color:#0070c0"> R(:,:,:,m) = R(:,:,:,m) - sum( R(:,:,:,(m+1):end), 4);<u></u><u></u></span></b></p>
<p class="MsoNormal" style="text-indent:9.6pt"><b><span style="color:#0070c0">end<u></u><u></u></span></b></p>
<p class="MsoNormal"><b><span style="color:#0070c0">end<u></u><u></u></span></b></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Is that the same way to use dim_sum_n_Wrap ?<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Regards,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Arnold</p></div></div></blockquote></div></div></blockquote></div></div></div></blockquote></div></div></div></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>