[ncl-talk] Running average

Dennis Shea shea at ucar.edu
Mon Dec 3 15:19:09 MST 2018


Hi Soma,

What you want are *block* 5-day and 7-day averages. EG: *"73* times for a
*365* day series."
It is your responsibility to be clear on your objective.
Otherwise, people spend their valuable time trying to figure out what you
want.

x(time, lat, lon)     .... (ntim,nlat,mlon)
     0      1    2     <=== dimension numbers [ Please read NCL
documentation]

Two approaches:

[1]

   endopt = ... ; -1, 0, 1   PLEASE read the documentation
   ndim_time = 0

   navg = 5
   xRunAvg = *runave_n_Wrap
<http://www.ncl.ucar.edu/Document/Functions/Contributed/runave_n_Wrap.shtml>(*x,
navg, endopt, ndim_time)  ; classic running average (ntim,nlat,mlon)
   xBlock_5 = xRunAvg(::navg,:,:)
   printVarSummary(xBlock_5)

   navg = 7
   xRunAvg = *runave_n_Wrap
<http://www.ncl.ucar.edu/Document/Functions/Contributed/runave_n_Wrap.shtml>(*x,
navg, endopt, ndim_time)  ; classic running average (ntim,nlat,mlon)
   xBlock_7 = xRunAvg(::navg,:,:)
   printVarSummary(xBlock_7)

[2]
*http://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_segment_values.shtml*
<http://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_segment_values.shtml>

This requires that 'time' have CF conforming units. I suggest you use [1].

On Mon, Dec 3, 2018 at 2:27 PM Soma Roy <somaroy892 at gmail.com> wrote:

> Hello,
>
> I used the example 2 from the below link;
> https://www.ncl.ucar.edu/Document/Functions/Built-in/runave_n.shtml
>
> But there is no change in output time dimension.
> I used 365 days data over nlat x nlon as input.
>
> I set nave=5 with expectation to get nlat x nlon for 73 times in output.
>
> But there is no change in output dimensions. Then how to extract 5 days
> running avg. over 365 days dataset?
>
> Kindly inform me.
>
> Thanking you,
> Soma
>
> On Mon, Dec 3, 2018, 21:10 Soma Roy <somaroy892 at gmail.com wrote:
>
>> Thanks for the hints
>>
>> Soma
>>
>> On Mon, Dec 3, 2018, 19:21 Appo derbetini <appopson4 at gmail.com wrote:
>>
>>> Hi,
>>> I think that you easily do it with
>>> *https://www.ncl.ucar.edu/Document/Functions/Built-in/runave_n.shtml
>>> <https://www.ncl.ucar.edu/Document/Functions/Built-in/runave_n.shtml>*
>>>
>>> Regards
>>>
>>> Le lun. 3 déc. 2018 à 10:54, Soma Roy <somaroy892 at gmail.com> a écrit :
>>>
>>>> Thank you.
>>>> Efficient way to do.
>>>>
>>>> But I would like to know also the steps for this by using ncl.
>>>>
>>>> Thanks,
>>>> Soma
>>>>
>>>> On Mon, Dec 3, 2018, 15:14 Oluwaseun Ilori <wilseun at gmail.com wrote:
>>>>
>>>>> Cdo will do it easily and better
>>>>>
>>>>> cdo -timselmean,5 infile outfile
>>>>>
>>>>>
>>>>> On Mon, Dec 3, 2018, 10:36 AM Soma Roy <somaroy892 at gmail.com wrote:
>>>>>
>>>>>> Hello,
>>>>>>
>>>>>> I am using daily gridded rainfall datasets.
>>>>>> I want to compute 5 days and 7 days running avg. of rainfall from the
>>>>>> datasets over each grid points.
>>>>>>
>>>>>> Structure of the input file is like below;
>>>>>>
>>>>>> Rain (730, 129, 135)  ; (time (365*2), lat, lon)
>>>>>>
>>>>>> Please kindly inform me which function should I use in this case.
>>>>>>
>>>>>> Thanking you,
>>>>>> Soma
>>>>>>
>>>>>>
>>>>>> _______________________________________________
>>>>>> ncl-talk mailing list
>>>>>> ncl-talk at ucar.edu
>>>>>> List instructions, subscriber options, unsubscribe:
>>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>>
>>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20181203/c781f5e4/attachment.html>


More information about the ncl-talk mailing list