[Tutor] numpy.mean across multiple netcdf files

questions anon questions.anon at gmail.com
Thu Aug 11 03:57:02 CEST 2011

I have many ncfiles each containing one month of hourly temperature data.
I have worked out how to loop through a number of ncfiles and calculate the
mean for each file at a particular time and even plot this and output as a
What I am unsure of is how to calculate the mean at a particular time for
all files.
Do I somehow output the mean for each folder and then calculate the mean of
all the means? Or is there some way to loop through and calculate it
Below is the code I have working that calculates the mean for each file. Any
help will be greatly appreciated!!

from netCDF4 import Dataset
import numpy as N

print MainFolder
for (path, dirs, files) in os.walk(MainFolder):
    for dir in dirs:
        print dir

    for ncfile in files:
        if ncfile[-3:]=='.nc':
            print "dealing with ncfiles:", ncfile
            ncfile=Dataset(ncfile, 'r+', 'NETCDF4')
            TSFC=ncfile.variables['T_SFC'][4::24,:,:] #this chooses a
particular time for each day for the whole month

            #calculate summary stats
            TSFCmean=N.mean(TSFC, axis=0)
            print "The TSFCmean array is:", TSFCmean

# there is a lot more code that follows this to plot the mean for each
ncfile but I have cut that out as I think it is irrelevant to my question
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/tutor/attachments/20110811/c3d61683/attachment.html>

More information about the Tutor mailing list