| 1 |
gforget |
1.1 |
function [v1]=convert2gcmfaces(v0,varargin); |
| 2 |
|
|
|
| 3 |
|
|
aa=whos('v0'); doGcm2Faces=strcmp(aa.class,'double'); |
| 4 |
|
|
|
| 5 |
|
|
if doGcm2Faces; |
| 6 |
|
|
|
| 7 |
|
|
if nargin==1; nFaces=5; else; nFaces=varargin{1}; end; |
| 8 |
|
|
|
| 9 |
|
|
if nFaces>1; |
| 10 |
|
|
|
| 11 |
|
|
nn=size(v0,1); |
| 12 |
|
|
pp=size(v0,2)/nn; |
| 13 |
|
|
mm=(pp+4-nFaces)/4*nn; |
| 14 |
|
|
|
| 15 |
|
|
n3=size(v0,3); n4=size(v0,4); |
| 16 |
|
|
for i3=1:n3; for i4=1:n4; |
| 17 |
|
|
v00=v0(:,:,i3,i4); |
| 18 |
|
|
|
| 19 |
|
|
v00=v00(:); |
| 20 |
|
|
i0=1; i1=nn*mm; v1{1}(:,:,i3,i4)=reshape(v00(i0:i1),nn,mm); |
| 21 |
|
|
i0=i1+1; i1=i1+nn*mm; v1{2}(:,:,i3,i4)=reshape(v00(i0:i1),nn,mm); |
| 22 |
|
|
i0=i1+1; i1=i1+nn*nn; v1{3}(:,:,i3,i4)=reshape(v00(i0:i1),nn,nn); |
| 23 |
|
|
i0=i1+1; i1=i1+nn*mm; v1{4}(:,:,i3,i4)=reshape(v00(i0:i1),mm,nn); |
| 24 |
|
|
i0=i1+1; i1=i1+nn*mm; v1{5}(:,:,i3,i4)=reshape(v00(i0:i1),mm,nn); |
| 25 |
|
|
if nFaces==6; |
| 26 |
|
|
i0=i1+1; i1=i1+nn*nn; |
| 27 |
|
|
end; |
| 28 |
|
|
|
| 29 |
|
|
end; end; |
| 30 |
|
|
else; |
| 31 |
|
|
v1={v0}; |
| 32 |
|
|
end; |
| 33 |
|
|
|
| 34 |
|
|
if nFaces==1; gridType='ll'; elseif nFaces==5; gridType='llc'; elseif nFaces==6; gridType='cube'; end; |
| 35 |
|
|
v1=gcmfaces(v1,gridType); |
| 36 |
|
|
|
| 37 |
|
|
else; |
| 38 |
|
|
|
| 39 |
|
|
nFaces=get(v0,'nFaces'); |
| 40 |
|
|
|
| 41 |
|
|
if nFaces>1; |
| 42 |
|
|
|
| 43 |
|
|
v0_faces=v0; clear v0; |
| 44 |
|
|
for iFace=1:nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end; |
| 45 |
|
|
|
| 46 |
|
|
nn=size(v0{1},1); mm=size(v0{1},2); |
| 47 |
|
|
pp=mm/nn*4+nFaces-4; |
| 48 |
|
|
|
| 49 |
|
|
n3=size(v0{1},3); n4=size(v0{1},4); |
| 50 |
|
|
v1=NaN*zeros(nn,nn*pp,n3,n4); |
| 51 |
|
|
for i3=1:n3; for i4=1:n4; |
| 52 |
|
|
|
| 53 |
|
|
v11=NaN*zeros(nn*nn*pp,1); |
| 54 |
|
|
|
| 55 |
|
|
i0=1; i1=nn*mm; tmp1=v0{1}(:,:,i3,i4); v11(i0:i1)=tmp1(:); |
| 56 |
|
|
i0=i1+1; i1=i1+nn*mm; tmp1=v0{2}(:,:,i3,i4); v11(i0:i1)=tmp1(:); |
| 57 |
|
|
i0=i1+1; i1=i1+nn*nn; tmp1=v0{3}(:,:,i3,i4); v11(i0:i1)=tmp1(:); |
| 58 |
|
|
i0=i1+1; i1=i1+nn*mm; tmp1=v0{4}(:,:,i3,i4); v11(i0:i1)=tmp1(:); |
| 59 |
|
|
i0=i1+1; i1=i1+nn*mm; tmp1=v0{5}(:,:,i3,i4); v11(i0:i1)=tmp1(:); |
| 60 |
|
|
if nFaces==6; |
| 61 |
|
|
i0=i1+1; i1=i1+nn*nn; |
| 62 |
|
|
end; |
| 63 |
|
|
|
| 64 |
|
|
v1(:,:,i3,i4)=reshape(v11,nn,nn*pp); |
| 65 |
|
|
end; end; |
| 66 |
|
|
|
| 67 |
|
|
else; |
| 68 |
|
|
v1=v0{1}; |
| 69 |
|
|
end; |
| 70 |
|
|
|
| 71 |
|
|
end; |
| 72 |
|
|
|
| 73 |
|
|
|
Parent Directory
| 
Revision Log
| 
Revision Graph