# Documentation of cosweight

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

## Function Synopsis

x=cosweight(sst, nya);

## Help text

x=cosweight(sst, nya);

This function weights data in matrix 'sst' by the cosine of
latitudes in matrix 'nya'.

## Cross-Reference Information

This function is called by

## Listing of function cosweight

function x=cosweight(sst, nya);

sst = squeeze(sst);

ndim = sum(size(size(sst)))-1;
if ndim == 3;
[ntim, nlat, nlon] = size(sst);

% Assume the second dimension is y

sst = reshape(shiftdim(sst, 1), nlat, nlon*ntim);
[m,n] = size(nya);
if m == 1; nya = nya'; end;
nya = nya .* (3.1415927)/180;
sst = sst .* sqrt(abs(cos(nya * ones(1, nlon*ntim))));
sst = shiftdim(reshape(sst, nlat, nlon, ntim), 2);

else

%  [tmax,ngrid]=size(sst);
%  ny=length(nya);
%  nx=ngrid/ny;
%  for i=1:ny
%    sst(:,(nx*(i-1)+1):(nx*i))=sst(:,(nx*(i-1)+1):(nx*i))*cos(nya(i)*pi/180);
%  end
%  error('2D matrices not done correctly, use 3D');

[nlat, nlon] = size(sst);

%  Assume the first dimension is y

[m,n] = size(nya);
if m == 1; nya = nya'; end;
nya = nya.*(3.1415927)/180;
sst = sst .* sqrt(abs(cos(nya * ones(1, nlon))));

end

x=sst;