Re: how to calculate average of seasonal means

Oh, you could use the addfiles function with the join option (see
example 3) and dim_avg_n.

http://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml

Something like:

   fils = systemfunc ("ls year*.nc") ; file paths
   f = addfiles (fils, "r")

   ListSetType (f, "cat") ; concatenate (=default)
   w = f[:]->w

   w_seas01 = dim_avg_n(w(:,0,:,:),0)
   w_seas02 = dim_avg_n(w(:,1,:,:),0)
   w_seas03 = dim_avg_n(w(:,2,:,:),0)
   w_seas04 = dim_avg_n(w(:,3,:,:),0)

On Tue, Oct 19, 2010 at 6:54 PM, Vishali P. <vishali.prat@gmail.com> wrote:
> Hi Ryan
>
> Thanks. I already used that command. I'll explain you in steps what I did
> and what I want
>
> 1. Firstly, my data are separated in files from 1998 to 2004, 2 files per
> year(one for component u and another for v)
>
> 2. I calculated (with ncl) wind speed with formula
>
>         w = sqrt(u^2+v^2)
>
> 3. Created outputs of w for each year (year file is aprox 950MB)
>
> 4. Used the command 'cdo yseasmean' from CDO for each year (7 files)
>
> 5. NOW, I want to 'pick' all first season of each year to calculate average,
> then the 2nd season , then 3rd and then 4th.
>
> 6. And finaly...plot 4 graphics of sazonal mean of wind speed(1998-2004)
>
> If something is not clear...please tell.
>
> Vp
>

-- 
------------------------------------------
Ryan Pavlick, PhD Student
Biospheric Theory and Modellling Group
Max Planck Institute for Biogeochemistry
phone: +49.36.41.57.6262
web: http://gaia.mpg.de
------------------------------------------
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk