<div dir="ltr"><div><div><div><div><div><div><div>re: runave<br></div>The answer is no. It uses centered values. For a 3 pt average"<br><br></div> x3 = (x(n-1) + x(n) + x(n+1))/3<br><br></div>To me, you want<br><br></div> x3 = (x(n-2) + x(n-1) +x(n))/3<br></div>or<br> x3 = (x(n-3) + x(n-2) +x(n-1))/3<br><br></div>which one?<br><br>===<br>I'm assuming 'x' is type float or double.<br></div>This is for a 1D time series<br><div><br>undef("jenny")<br>function jenny(x[*]:numeric, nsum[1]:integer)<br>local nx, xavg, nStrt, nLast<br>begin<br> nx = dimsizes(x)<br> xavg = x <br> xavg = x@_FillValue<br><br> nStrt = 0<br> nLast = nsum-1<br> do n=nLast,nx-1<br> xavg(n) = avg(x(nStrt:nLast))<br> print("n="+n+" nStrt="+nStrt+" nLast="+nLast+" xavg(n)="+xavg(n))<br> nStrt = nStrt+1<br> nLast = nLast+1<br> end do<br><br> return(xavg)<br>end<br><br> N = 100<br> x = conform_dims( N, 1.0, -1) ; all 1s<br> x@_FillValue = -999.0<br><br> k = 10<br><br> z = jenny(x,k)<br> print(z) all 1s<br><br><div><div><div><div><div><br> <br><br><br></div></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Feb 2, 2017 at 4:11 PM, Andrea Jenney <span dir="ltr"><<a href="mailto:andrea@atmos.colostate.edu" target="_blank">andrea@atmos.colostate.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Hi NCL Community,<div><br></div><div>I am hoping you can provide me with some clarification on the documentation for the ‘runave’ function. I am trying to compute a running mean for a time series but using only points from previous (and not future) time steps to compute the mean for each point. Is this something I can do with the ‘runave’ function? If so, which option do I select?</div><div><br></div><div>Example, for clarity:</div><div>- x is my array of values</div><div>- running average window of size 10</div><div>- I would like an array, y, of size(x), at which every point, the value at y(i) is computed using the mean of x(i-9:i), but do not want to utilize cyclic conditions (runave opt -1)</div><div><br></div><div>Thank you</div><div><div>
<div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div><span style="background-color:rgb(255,255,255)">---</span></div><div><span style="background-color:rgb(255,255,255)">Andrea Jenney</span><br><span style="background-color:rgb(255,255,255)">Graduate Research Assistant | Randall Group</span></div><div><span style="background-color:rgb(255,255,255)">Colorado State University | Department of Atmospheric Science </span></div><div><br></div></div><br class="m_-6220714160037038896Apple-interchange-newline"><br class="m_-6220714160037038896Apple-interchange-newline">
</div>
<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>