/[MITgcm]/MITgcm_contrib/enderton/Diagnostics/DiagAverage.m
ViewVC logotype

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

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by enderton, Mon Jan 31 15:43:26 2005 UTC revision 1.2 by molod, Tue Jun 28 21:33:51 2005 UTC
# Line 32  DiagGenParam; Line 32  DiagGenParam;
32  % Calcuate months of year.  % Calcuate months of year.
33  months = mod(absmonths-1,12)+1;  months = mod(absmonths-1,12)+1;
34    
   
35  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
36  %                       Take monthly average of data                      %  %                       Take monthly average of data                      %
37  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# Line 58  else Line 57  else
57              if isempty(find(months == imon))              if isempty(find(months == imon))
58                  dataMonAvg(:,:,imon) = NaN * zeros(size(data(:,:,1)));                  dataMonAvg(:,:,imon) = NaN * zeros(size(data(:,:,1)));
59              else              else
60                  dataMonAvg(:,:,imon) = meanovernan(data(:,:,months == imon),3);                  inds=find(months==imon);
61                    dataMonAvg(:,:,imon) = meanovernan(data(:,:,absmonths(inds)),3);
62              end              end
63          elseif ismember(fln,fields3D) || isequal(Dim,3)          elseif ismember(fln,fields3D) || isequal(Dim,3)
64              if isempty(find(months == imon))              if isempty(find(months == imon))
65                  dataMonAvg(:,:,:,imon) = NaN * zeros(size(data(:,:,:,1)));                  dataMonAvg(:,:,:,imon) = NaN * zeros(size(data(:,:,:,1)));
66              else              else
67                  dataMonAvg(:,:,:,imon) = meanovernan(data(:,:,:,months == imon),4);                  inds=find(months==imon);
68                    dataMonAvg(:,:,:,imon) = meanovernan(data(:,:,:,absmonths(inds)),4);
69              end              end
70          else          else
71              error('Field not accounted for in ''fields2D'' of ''fields3D''');              error('Field not accounted for in ''fields2D'' of ''fields3D''');
# Line 149  end Line 150  end
150  %                           Prepare for output                            %  %                           Prepare for output                            %
151  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
152    
 data = temp;  
153    data = temp;
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22