Documentation of remove_anncyc

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

Function Synopsis

[x, clim] = annave(sst, clim, plusorminus);

Help text

   [x, clim] = annave(y);

   x = matrix with annual average removed
   y = input matrix
   clim = climatology

Cross-Reference Information

Listing of function remove_anncyc

function [x, clim] = annave(sst, clim, plusorminus);
n = ndims(sst);

szsst = size(sst);
sst = reshape(sst, szsst(1), prod(szsst(2:n)));

[tmax,ngrid]=size(sst);

%  Get climatology
if nargin < 2; 
  plusorminus=-1;
  for m=1:12
    clim(m,:)=mean2(sst(m:12:(tmax),:));
  end
end;

%  Remove climatology from each time step
x=zeros(tmax,ngrid);
for m=1:tmax;
  l=rem(m-1,12)+1;
  x(m,:)=sst(m,:)+plusorminus*clim(l,:);
end

%  Put data back into output format
x = reshape(x, szsst);
clim = reshape(clim, [12 szsst(2:n)]);