MITgcm_contrib/gag_subduction_scripts/cshift.m

function [out] = cshift(in,DIM,shift)
%function [out] = cshift(in,DIM,shift)
%
% Replicate the CSHIFT function in F90 (?). 
%
% G. Gebbie, MIT-WHOI, Dec 2003.

 totaldims = ndims(in);
 index = 1: totaldims;
 index(index==DIM) = [];
 index = [DIM index];
 sizin = size(in);
 in = permute(in,index);
 in = reshape(in,sizin(DIM),prod(sizin)./sizin(DIM));

 if shift>=0
   shift = shift - size(in,1);
 end  

 out = [in(shift+1+size(in,1):size(in,1),:);in(1:size(in,1)+shift,:)];
 out = reshape(out,sizin(index));
 out = ipermute(out,index);
1	gebbie	1.1	function [out] = cshift(in,DIM,shift)
2			%function [out] = cshift(in,DIM,shift)
3			%
4			% Replicate the CSHIFT function in F90 (?).
5			%
6			% G. Gebbie, MIT-WHOI, Dec 2003.
7
8			totaldims = ndims(in);
9			index = 1: totaldims;
10			index(index==DIM) = [];
11			index = [DIM index];
12			sizin = size(in);
13			in = permute(in,index);
14			in = reshape(in,sizin(DIM),prod(sizin)./sizin(DIM));
15
16			if shift>=0
17			shift = shift - size(in,1);
18			end
19
20			out = [in(shift+1+size(in,1):size(in,1),:);in(1:size(in,1)+shift,:)];
21			out = reshape(out,sizin(index));
22			out = ipermute(out,index);