Re: daily mean

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Wed, 1 Mar 2006 06:35:24 -0700 (MST)

> I have 6 hourly data for a month. I want to get daily mean of each day at each grid point. What is the function I can use.
>

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/4

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

         do nt=0,ntim-1,4
            xday(nt/4,:,:,:) = dim_avg(x(lev|:,lat|:,lon|:,time|nt:nt+3))
         end do
         xday_at_long_name = "Daily Mean: "+x_at_long_name

Good luck
D
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Mar 01 2006 - 06:35:24 MST

This archive was generated by hypermail 2.2.0 : Wed Mar 01 2006 - 10:20:53 MST