<div dir="ltr">Oops.  Please ignore my correction.  Dennis's version was correct.  The stride parameter is different between fortran and NCL.<div><br></div><div>--Dave</div><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_-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>