So, I assume you have the following
X = f1->X
a1 = dim_avg_Wrap(X(south_north|:,west_east|:, Time|:) )
;detete(X) ; only if each T has different Time dimension sizes
X = f2->X
b1 = dim_avg_Wrap(X(south_north|:,west_east|:, Time|:) )
etc,
Hence a1, b1, ... will have the time average for a specific period
They will also be 2 dimensional sime the time timension is degenerate
eg
printVarSummary( a1 )
To get the overall time average
xAVG = a1 + b1 + .... ; array syntax
------------
An alternative is to use addfiles
http://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml
diri = "..." ; input directory
fils = systemfunc ("ls "+diri+"wrf*.nc") ; file paths
f = addfiles (fils, "r")
X = f[:]->X ; read X from all files
printVarSummary( X )
xAVG = dim_avg_Wrap( X )
printVarSummary( xAVG )
---------------
good luck
On Wed, 1 Aug 2007, Sean Heuser wrote:
> Hello NCL Users,
> I'm am working with a 2-month WRF simulation and I'm trying to get the
> average temperature for the entire time frame. Since computer time was
> scarce, I have 9 input files(call them a,b,c,d,e,f,g,h,i) that span the
> time frame. I called the temperature variable for each(so I have 9 temp
> variables called a1,b1,c1,etc.). I tried using the dim_Avg function and
> was successful with one variable. For ex:
>
> Avg = dim_Avg(a1(south_north|:, west_east|:, Time|:)
>
> works. However, when I try it with more than one, it reads an error saying
> that it only expects 1 subscript. Does anyone know how I can possibly do
> an average with all of my temperature variables? Any help would be
> appreciated. Thank you.
>
>
>
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Aug 01 2007 - 12:03:50 MDT
This archive was generated by hypermail 2.2.0 : Mon Aug 06 2007 - 08:08:45 MDT