Documentation of area_mean

Global Index (all files) (short | long) | Local Index (files in subdir) (short | long)

Function Synopsis

[a2,lat2,lon2] = area_mean(a, inc, lat, lon);

Help text

  [a2, lat2, lon2] = area_mean(in1, inc, lat, lon);

  Thin will condense the matrix in1 by averaging
  over every inc element in each direction.  This
  program assumes the input data is three-
  dimensional, with time, lat, lon, as the dimensions.

Cross-Reference Information

Listing of function area_mean

function [a2,lat2,lon2] = area_mean(a, inc, lat, lon);

[ntim, nlat, nlon] = size(a);
nlat2 = floor(nlat/inc(1));
nlon2 = floor(nlon/inc(2));
a2 = repmat(NaN, [ntim, nlat2, nlon2]);

%  Average the data
for i = 1:nlat2;
  ind1 = inc(1)*(i-1)+[1:inc(1)];
  for j = 1:nlon2;
    ind2 = inc(2)*(j-1)+[1:inc(2)];
    a2(:,i,j) = squeeze(mean2(mean2(shiftdim(a(:,ind1,ind2),1))))';
  end
end

%  Do the same for lat and lon
if nargin == 4;
  lat2 = repmat(NaN, [nlat2, 1]);
  for i = 1:nlat2;
    ind1 = inc(1)*(i-1)+[1:inc(1)];
    lat2(i) = mean(lat(ind1));
  end

  lon2 = repmat(NaN, [nlon2, 1]);
  for i = 1:nlon2;
    ind1 = inc(2)*(i-1)+[1:inc(2)];
    lon2(i) = mean(lon(ind1));
  end
end