[ncl-talk] Time series script

[Alan Brammer]

Copying a similar response from a previous question on NCL talk.

google: hourly to daily ncl
will give you many queries and responses along similar line.

 For a more specific answer, a more specific question is going to be
needed.

let: x(time,lev,lat,lon)

         dimx = dimsizes(x)
         ntim = dimx(0) ; number of time steps
         klev = dimx(1)
         nlat = dimx(2)
         mlon = dimx(3)

         nday = ntim/24

         xDay = x(::24,:,:,:) ; 'trick' to allocate memory + retain meta
data

         do nt=0,ntim-1,24
            xday(nt/24,:,:,:) = dim_avg_n( x(nt:nt+23,:,:,:),0 )  ;

         end do
         xday_at_long_name = "Daily Mean: "+x_at_long_name

Alan.


On Thu, Sep 24, 2015 at 11:21 AM, Jesús Garcia Rosales <jesus21gr at gmail.com>
wrote:

> Hi users
>
> The thing that I want to do is the next:
>
> I have many wrf files for february ( example: wrfout_2000_feb,
> wrfout_2001_feb,.... wrfout_2012_feb), and in each file has hourly data
> (example: to wrfout_2000_feb it has day 1: 00z, 01z, 02z,....day 28: 00z,
> 01z.....). And I have daily data from my stations, so I want to compare
> both each station with each wrf file in a time-serie during a february
> month, but the data from station is daily and the data from wrf is hourly.
> Can someone help me please?.
>
>
>  Thanks for your time.
> Jesús
>
> _______________________________________________
> 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/20150924/a2efc1cf/attachment.html