<div dir="ltr">Hi Lyndz,<div>The wgt_runave documentation describes how the function works:</div><div><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/wgt_runave.shtml">https://www.ncl.ucar.edu/Document/Functions/Built-in/wgt_runave.shtml</a><br></div><div><br></div><div>As you are specifying an odd number of weights, the centered value is where the smoothed 3yr value will be placed. I find that doing a simple test in interactive NCL is a quick way to check how many NCL functions work:</div><div>ncl</div><div>a = (/2.3,1.2,5,-1.2/)</div><div>a@_FillValue = -9999.</div><div>b = wgt_runave(a,(/0.25,0.5,0.25/),0)</div><div>print(b)</div><div>;---- output-----</div><div><div>Variable: b</div><div>Type: float</div><div>Total Size: 16 bytes</div><div> 4 values</div><div>Number of Dimensions: 1</div><div>Dimensions and sizes:<span style="white-space:pre"> </span>[4]</div><div>Coordinates: </div><div>Number Of Attributes: 1</div><div> _FillValue :<span style="white-space:pre"> </span>-9999</div><div>(0)<span style="white-space:pre"> </span>-9999</div><div>(1)<span style="white-space:pre"> </span>2.425</div><div>(2)<span style="white-space:pre"> </span>2.5</div><div>(3)<span style="white-space:pre"> </span>-9999</div></div><div><br></div><div>The average of 2.3, 1.2 and 5 are placed in the (1) slot (=the middle). Values at the beginning and end are missing as there are no values at (-1) or (4). As the documentation states, you can set the option to 1 to get values at the end points:</div><div><div>c = wgt_runave(a,(/0.25,0.5,0.25/),1)</div><div>print(c)</div></div><div>;---output-----</div><div><div>Variable: c</div><div>Type: float</div><div>Total Size: 16 bytes</div><div> 4 values</div><div>Number of Dimensions: 1</div><div>Dimensions and sizes:<span style="white-space:pre"> </span>[4]</div><div>Coordinates: </div><div>Number Of Attributes: 1</div><div> _FillValue :<span style="white-space:pre"> </span>-9999</div><div>(0)<span style="white-space:pre"> </span>1.75</div><div>(1)<span style="white-space:pre"> </span>2.425</div><div>(2)<span style="white-space:pre"> </span>2.5</div><div>(3)<span style="white-space:pre"> </span>1.9</div></div><div><br></div><div>Hope that helps. If you have any further questions please reply to the ncl-talk email list.</div><div>Adam</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 5, 2017 at 1:22 AM, Lyndon Mark Olaguera <span dir="ltr"><<a href="mailto:olagueralyndonmark429@gmail.com" target="_blank">olagueralyndonmark429@gmail.com</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"><div><div class="m_1440626870289179716gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Dear ncl-experts,</div><div><br></div><div>I am just confused about the output of the 1-2-1 average in NCL.</div><div><br></div><div>I have a csv file containing 30 years of data (1979-2008). I attached the raw file (P28-P31.csv) and the smoothed file (P28-P31_filt.csv). After applying the 3 point 1-2-1 average, I got 29 points with -9999 (missing at both ends) after applying the script. </div><div><br></div><div>Question:</div><div><br></div><div>[1] I'm confused now which output point corresponds to a particular year (which point in the smoothed file corresponds to 1979, 1980,etc).</div><div>[2[ Is there a way to add a "filler" value to the points without data after smoothing automatically?</div><div><br></div><div><br></div><div><div>begin</div><div> diri = "./"</div><div> files = systemfunc("ls *.csv") </div><div> numfiles = dimsizes(files)</div><div> do ifil = 0,numfiles-1</div><div> lines = asciiread(files(ifil),-1,"<wbr>string")</div><div> data<span class="m_1440626870289179716gmail-Apple-tab-span" style="white-space:pre-wrap"> </span> = lines(1:)</div><div> rain<span class="m_1440626870289179716gmail-Apple-tab-span" style="white-space:pre-wrap"> </span> = tofloat(str_get_field(data,1,"<wbr>,"))</div><div>;;****************************<wbr>***********************</div><div>;;Apply 1-2-1 running average (three pentads)</div><div>;;****************************<wbr>************************</div><div> x = wgt_runave(rain, (/0.25, 0.50, 0.25/), 0)</div><div> x@_FillValue = -9999</div><div>;;****************************<wbr>************************</div><div>;;Save output to a csv file</div><div>;;****************************<wbr>************************</div><div> flname = str_get_cols(files(ifil),0,6)+<wbr>"_filt.csv"</div><div> asciiwrite(flname,x)</div><div> end do</div><div>end</div></div><div><br></div><div>Thank you in advance for the help.</div><div><br></div><div>Sincerely,</div><div><br></div><div>Lyndz</div><div><br></div></div></div></div></div></div></div></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><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist, </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphilli/</a> </font></span><span><font color="#888888">303-497-1726 </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</div>