Re: how to calculate average of seasonal means

Oops, that should have been:

ListSetType (f, "join")

On Tue, Oct 19, 2010 at 7:10 PM, Ryan Pavlick <ryan.pavlick@gmail.com> wrote:
> 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
> ------------------------------------------
>

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