/[MITgcm]/MITgcm/utils/matlab/merccube.m
ViewVC logotype

Diff of /MITgcm/utils/matlab/merccube.m

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

revision 1.4 by adcroft, Thu Aug 30 14:00:21 2001 UTC revision 1.5 by adcroft, Fri Jun 4 17:03:50 2004 UTC
# Line 33  else Line 33  else
33  end  end
34  Q=C;  Q=C;
35    
36  if ndims(X)==2 & size(X,1)==6*size(X,2)  if ndims(Q)==2 & size(Q,1)==6*size(Q,2)
37   [nx ny nt]=size(X);   [nx ny nt]=size(X);
38   X=permute( reshape(X,[nx/6 6 ny]),[1 3 2]);   X=permute( reshape(X,[nx/6 6 ny]),[1 3 2]);
39   Y=permute( reshape(Y,[nx/6 6 ny]),[1 3 2]);   Y=permute( reshape(Y,[nx/6 6 ny]),[1 3 2]);
40   Q=permute( reshape(Q,[nx/6 6 ny]),[1 3 2]);   Q=permute( reshape(Q,[nx/6 6 ny]),[1 3 2]);
41  elseif ndims(X)==3 & size(X,2)==6  elseif ndims(Q)==3 & size(Q,2)==6
42   X=permute( X,[1 3 2]);   X=permute( X,[1 3 2]);
43   Y=permute( Y,[1 3 2]);   Y=permute( Y,[1 3 2]);
44   Q=permute( Q,[1 3 2]);   Q=permute( Q,[1 3 2]);
45  elseif ndims(X)==3 & size(X,3)==6  elseif ndims(Q)==3 & size(Q,3)==6
46   [nx ny nt]=size(X);   [nx ny nt]=size(X);
47  else  else
48   size(XX)   size(XX)
# Line 67  if size(X,1)==size(Q,1) Line 67  if size(X,1)==size(Q,1)
67  end  end
68  [nx ny nt]=size(X);  [nx ny nt]=size(X);
69    
70  Q(end+1,:,:)=0;  Q(end+1,:,:)=NaN;
71  Q(:,end+1,:)=0;  Q(:,end+1,:)=NaN;
72  Q(end,:,[1 3 5])=Q(1,:,[2 4 6]);  Q(end,:,[1 3 5])=Q(1,:,[2 4 6]);
73  Q(:,end,[2 4 6])=Q(:,1,[3 5 1]);  Q(:,end,[2 4 6])=Q(:,1,[3 5 1]);
74  Q(:,end,[1 3 5])=squeeze(Q(1,end:-1:1,[3 5 1]));  Q(:,end,[1 3 5])=squeeze(Q(1,end:-1:1,[3 5 1]));
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22