<div dir="ltr"><div><div>
Hi Hesham,<br></div>For simple non-weighted running averages you have
two choices when it comes to the class of functions to use: the runave*
suite of functions and the dim_avg suite of functions. Missing values
are allowed and ignored in the dim_avg* functions, while missing data
results in missing averages in the runave suite of functions. I would
suggest you use dim_avg_n and run the data through a do loop. <br><br><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/dim_avg.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Built-in/dim_avg.shtml</a><br><br></div><div>For example, you can do something like this:<br></div><div>int = 7 ; 7 day average<br></div><div>data_7day = new(dimsizes(data),typeof(data)<br></div><div>do gg = 3,dimsizes(data)-4<br></div><div> data_7day(gg) = (/ dim_avg(data(gg-3:gg+3)) /)<br></div><div>end do<br><br></div><div>data_7day
would hold 7 day running averages at each timestep, and would compute
it if there is 1 day out of 7 that has data. The above can be
generalized without the do loop as well using syntax from the script I
sent you:<br><br> stime = 19991001 ; plot start time in YYYYMMDD<br> etime = 19991031 ; plot end time in YYYYMMDD<br> si = ind(time.eq.stime)<br> ei = ind(time.eq.etime)<br><br></div><div>data_avg = dim_avg_Wrap(data(si:ei))<br><br></div><div>data_avg would hold the October 1999 average.<br></div><div><br></div><div>Finally, if you would like to interpolate to fill in the missing data, I would recommend using the linmsg function:<br><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/linmsg.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Built-in/linmsg.shtml</a><br><br></div><div>Hope all that helps. If not, please respond to the ncl-talk email list and not to
me personally.<br></div>Adam</div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Apr 6, 2015 at 1:00 PM, Hesham Afify <span dir="ltr"><<a href="mailto:algocomp@gmail.com" target="_blank">algocomp@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 class="gmail_default" style="font-size:large">Hi Adam</div><div class="gmail_default" style="font-size:large">Thanks for your kind attention and reply, I used the script and it's really helpful and hits the point.</div><div class="gmail_default" style="font-size:large">The aspect facing me now is that the missing data not temporally evenly distributed, so I'm wondering if NCL has some function that could make weekly or monthly average with minimal error, on other words, Is there a way to relax missing values in data to some ideal function (or climatological values)</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">thanks in advance</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">Hesham</div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Mar 27, 2015 at 3:28 PM, Adam Phillips <span dir="ltr"><<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@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"><div>Hi Hesham,<br></div><div>Due to the fact that you are reading in an
ascii file, are working on daily data, and wish to plot customized time
ranges, this is a case where conceptually what you ask for sounds
simple but in reality it is not.<br><br></div><div>I have come up with an example code that (somewhat) does what you want. You will have to modify the attached script to suit your own needs. I cannot spend any more time on this. You will also have to take a close look at the coding to verify that the data is being read in and used correctly. <br></div><div><br></div><div>I
used a number of functions to parse and prepare the ascii data for
plotting: str_squeeze, str_split, str_split_csv, tofloat, toint and
sprinti. I also used cd_inv_calendar to prepare a time array that could
be used with the time_axis_labels procedure so that you would get nicer
tick mark labels. Note that you will likely need to play with
time_axis_labels somewhat to get the tick mark labels where you want
them:<br><br><a href="http://www.ncl.ucar.edu/Document/Functions/User_contributed/time_axis_labels.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/User_contributed/time_axis_labels.shtml</a><br><br></div><div>Finally, note that you can set the range of times to be plotted by setting stime/etime. You can plot any range of times from a few days to a few months.<br></div><div><br></div><div>I
have attached the code and an example graphic here. If you have any
further questions please respond to the ncl-talk email list and not to
me personally. That way others can see the dialog and can assist.<br></div><div>Good luck,<br></div>Adam </div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Tue, Mar 24, 2015 at 8:45 AM, Hesham Afify <span dir="ltr"><<a href="mailto:algocomp@gmail.com" target="_blank">algocomp@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr"><div style="font-size:large">Hi All</div><div style="font-size:large">I have some observational daily data for a couple of years, this data have many changes in regularity by the mean it toke on daily, one per three days, and one per six days patterns, have some missed periods.</div><div style="font-size:large">the point here that I want to plot this data and my inquiry about that is how to plot this data in daily, weekly, and monthly basis.</div><div style="font-size:large">a sample of this data is attached</div><div style="font-size:large"><br></div><div style="font-size:large">best regards</div><span><font color="#888888"><div style="font-size:large">Hesham </div><div><br></div>
</font></span></div>
<br></div></div>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><span><font color="#888888"><br><br clear="all"><span class="HOEnZb"><font color="#888888"><br>-- <br><div><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"><a href="tel:303-497-1726" value="+13034971726" target="_blank">303-497-1726</a> </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>
</font></span></font></span></div><span class="HOEnZb"><font color="#888888">
</font></span></blockquote></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div><div dir="ltr"><p> </p><p style="background-image:initial;background-repeat:initial"><br></p></div></div>
</font></span></div></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="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>