Diff of /MITgcm_contrib/enderton/Diagnostics/DiagAverage.m

Parent Directory | Revision Log | Revision Graph | Patch

--- MITgcm_contrib/enderton/Diagnostics/DiagAverage.m	2005/01/31 15:43:26	1.1
+++ MITgcm_contrib/enderton/Diagnostics/DiagAverage.m	2005/06/28 21:33:51	1.2
@@ -32,7 +32,6 @@
 % Calcuate months of year.
 months = mod(absmonths-1,12)+1;
 
-
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %                       Take monthly average of data                      %
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -58,13 +57,15 @@
             if isempty(find(months == imon))
                 dataMonAvg(:,:,imon) = NaN * zeros(size(data(:,:,1)));
             else
-                dataMonAvg(:,:,imon) = meanovernan(data(:,:,months == imon),3);
+                inds=find(months==imon);
+                dataMonAvg(:,:,imon) = meanovernan(data(:,:,absmonths(inds)),3);
             end
         elseif ismember(fln,fields3D) || isequal(Dim,3)
             if isempty(find(months == imon))
                 dataMonAvg(:,:,:,imon) = NaN * zeros(size(data(:,:,:,1)));
             else
-                dataMonAvg(:,:,:,imon) = meanovernan(data(:,:,:,months == imon),4);
+                inds=find(months==imon);
+                dataMonAvg(:,:,:,imon) = meanovernan(data(:,:,:,absmonths(inds)),4);
             end
         else
             error('Field not accounted for in ''fields2D'' of ''fields3D''');
@@ -149,4 +150,4 @@
 %                           Prepare for output                            %
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
-data = temp;
\ No newline at end of file
+data = temp;