<div>Thanks for all of the responses. I have tried adding in the code you mentioned (see below). I am not sure if I am putting it in the correct place? and I am now receiving another error:</div><div>&quot;UserWarning: Warning: converting a masked element to nan.&quot;</div>
<div>Not sure if that is bringing me any closer? Any feedback will be greatly appreciated. </div><div><br></div><div>from netCDF4 import Dataset</div><div>import matplotlib.pyplot as plt</div><div>import numpy as N</div><div>
from mpl_toolkits.basemap import Basemap</div><div>import os</div><div><br></div><div>MainFolder=r&quot;E:/DSE_BushfireClimatologyProject/griddeddatasamples/GriddedData/T_SFC/&quot;</div><div>all_TSFC=[] </div><div>for (path, dirs, files) in os.walk(MainFolder):</div>
<div>    for dir in dirs:</div><div>        print dir</div><div>    path=path+&#39;/&#39;</div><div>    for ncfile in files:</div><div>        if ncfile[-3:]==&#39;.nc&#39;:</div><div>            ncfile=os.path.join(path,ncfile)</div>
<div>            ncfile=Dataset(ncfile, &#39;r+&#39;, &#39;NETCDF4&#39;)</div><div>            TSFC=ncfile.variables[&#39;T_SFC&#39;][4::24,:,:]</div><div>            LAT=ncfile.variables[&#39;latitude&#39;][:]</div><div>
            LON=ncfile.variables[&#39;longitude&#39;][:]</div><div>            #TIME=ncfile.variables[&#39;time&#39;][:]</div><div>            fillvalue=ncfile.variables[&#39;T_SFC&#39;]._FillValue</div><div>            ncfile.close()</div>
<div><br></div><div>            array=N.true_divide(TSFC[0],len(TSFC))</div><div><br></div><div>            for i in xrange(1, len(TSFC)-1,1):</div><div>                            array=N.add(array, N.true_divide(array[i],len(TSFC)))</div>
<div><br></div><div>#plot output summary stats</div><div>map = Basemap(projection=&#39;merc&#39;,llcrnrlat=-40,urcrnrlat=-33,</div><div>              llcrnrlon=139.0,urcrnrlon=151.0,lat_ts=0,resolution=&#39;i&#39;)</div><div>
<br></div><div>x,y=map(*N.meshgrid(LON,LAT))</div><div>plt.title(&#39;TSFC Mean at 3pm&#39;)</div><div>ticks=[-5,0,5,10,15,20,25,30,35,40,45,50]</div><div>CS = map.contourf(x,y,array, cmap=plt.cm.jet)</div><div>l,b,w,h =0.1,0.1,0.8,0.8</div>
<div>cax = plt.axes([l+w+0.025, b, 0.025, h])</div><div>plt.colorbar(CS,cax=cax, drawedges=True)</div><div><br></div><div>plt.savefig((os.path.join(MainFolder, &#39;Mean.png&#39;)))</div><div>plt.show()</div><div>plt.close()</div>
<div><br></div><div>  </div><div><br></div><div><br></div><br><div class="gmail_quote">On Sat, Aug 27, 2011 at 10:54 AM, srean <span dir="ltr">&lt;<a href="mailto:srean.list@gmail.com">srean.list@gmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><br><div class="gmail_quote"><div class="im">On Fri, Aug 26, 2011 at 2:33 PM, Phil Morefield <span dir="ltr">&lt;<a href="mailto:philmorefield@yahoo.com" target="_blank">philmorefield@yahoo.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204, 204, 204);padding-left:1ex">
<div><div style="color:rgb(0, 0, 0);background-color:rgb(255, 255, 255);font-family:tahoma,new york,times,serif;font-size:12pt"><div><div><span><br></span></div></div>
<div><span>The formula you have written looks like you&#39;re collapsing everything into a single value. I think he&#39;s trying to average a bunch of 2D arrays into a single 2D array.</span></div></div></div></blockquote>

</div><div><br>You are correct, the form that I posted can be read as if it is  for updating single mean vector \mu, but you can use the same for an nd-array trivially. Just have \mu and t as nd-arrays. m can be one too.  Numpy broadcasting will take care of the rest.<br>

<br>One advantage is that it requires only a constant amount of memory for the computation, you can even read the data in from an infinite pipe or generator that yields a single vector or a matrix at a time (or bundles them up m at a time). It will always be uptodate with the current estimate of the means. In fact will work for any moment too.<br>
<font color="#888888">
<br>--srean<br></font></div></div><br>
<br>_______________________________________________<br>
SciPy-User mailing list<br>
<a href="mailto:SciPy-User@scipy.org">SciPy-User@scipy.org</a><br>
<a href="http://mail.scipy.org/mailman/listinfo/scipy-user" target="_blank">http://mail.scipy.org/mailman/listinfo/scipy-user</a><br>
<br></blockquote></div><br>