1 |
function [v1]=convert2gcmfaces(v0,varargin); |
function [v1]=convert2gcmfaces(v0,varargin); |
2 |
%object: converts model output to gcmfaces object |
%object: converts model output to gcmfaces object |
3 |
% or vice versa |
% or vice versa |
4 |
%input: v0 model output array (resp. gcmfaces object) |
%input: v0 model output array (resp. gcmfaces object) |
5 |
%output: v1 gcmfaces object (model output array) |
%output: v1 gcmfaces object (model output array) |
13 |
global mygrid; |
global mygrid; |
14 |
|
|
15 |
if doGcm2Faces&mygrid.gcm2facesFast; |
if doGcm2Faces&mygrid.gcm2facesFast; |
16 |
|
|
17 |
[n1,n2,n3,n4,n5]=size(v0); |
[n1,n2,n3,n4,n5]=size(v0); |
18 |
|
|
19 |
v0=reshape(v0,[n1*n2 n3*n4*n5]); |
v0=reshape(v0,[n1*n2 n3*n4*n5]); |
20 |
for iFace=1:mygrid.nFaces; |
for iFace=1:mygrid.nFaces; |
21 |
v1{iFace}=reshape(mygrid.gcm2faces{iFace}*v0,[size(mygrid.XC{iFace}) n3 n4 n5]); |
v1{iFace}=reshape(mygrid.gcm2faces{iFace}*v0,[size(mygrid.XC{iFace}) n3 n4 n5]); |
22 |
end; |
end; |
23 |
|
|
24 |
if mygrid.nFaces==1; gridType='ll'; |
v1=gcmfaces(v1); |
25 |
elseif mygrid.nFaces==5; gridType='llc'; |
|
|
elseif mygrid.nFaces==6; gridType='cube'; |
|
|
end; |
|
|
v1=gcmfaces(v1,gridType); |
|
|
|
|
26 |
elseif ~doGcm2Faces&mygrid.gcm2facesFast; |
elseif ~doGcm2Faces&mygrid.gcm2facesFast; |
27 |
|
|
28 |
[n1,n2,n3,n4,n5]=size(v0{1}); |
[n1,n2,n3,n4,n5]=size(v0{1}); |
29 |
|
|
30 |
v1=zeros(mygrid.faces2gcmSize(1)*mygrid.faces2gcmSize(2),n3*n4*n5); |
v1=zeros(mygrid.faces2gcmSize(1)*mygrid.faces2gcmSize(2),n3*n4*n5); |
31 |
for iFace=1:mygrid.nFaces; |
for iFace=1:mygrid.nFaces; |
32 |
nn=size(mygrid.XC{iFace}); |
nn=size(mygrid.XC{iFace}); |
33 |
v1(mygrid.faces2gcm{iFace},:)=reshape(v0{iFace},[nn(1)*nn(2) n3*n4*n5]); |
v1(mygrid.faces2gcm{iFace},:)=reshape(v0{iFace},[nn(1)*nn(2) n3*n4*n5]); |
34 |
end; |
end; |
35 |
v1=reshape(v1,[mygrid.faces2gcmSize n3 n4 n5]); |
v1=reshape(v1,[mygrid.faces2gcmSize n3 n4 n5]); |
36 |
|
|
37 |
elseif doGcm2Faces; |
elseif doGcm2Faces; |
38 |
|
|
39 |
[n1,n2,n3,n4,n5]=size(v0); |
[n1,n2,n3,n4,n5]=size(v0); |
40 |
|
|
41 |
if strcmp(mygrid.fileFormat,'straight'); |
if strcmp(mygrid.fileFormat,'straight'); |
42 |
v1={v0}; |
v1={v0}; |
43 |
elseif strcmp(mygrid.fileFormat,'cube'); |
elseif strcmp(mygrid.fileFormat,'cube'); |
44 |
for ii=1:6; v1{ii}=v0(n2*(ii-1)+[1:n2],:,:,:,:); end; |
for ii=1:6; v1{ii}=v0(n2*(ii-1)+[1:n2],:,:,:,:); end; |
45 |
elseif strcmp(mygrid.fileFormat,'compact'); |
elseif strcmp(mygrid.fileFormat,'compact')|strcmp(mygrid.fileFormat,'nctiles'); |
46 |
|
v00=reshape(v0,[n1*n2 n3*n4*n5]); |
47 |
nn=size(v0,1); |
i0=0; i1=0; |
48 |
pp=size(v0,2)/nn; |
for iFace=1:mygrid.nFaces; |
49 |
mm=(pp+4-mygrid.nFaces)/4*nn; |
i0=i1+1; |
50 |
|
nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2); |
51 |
v00=reshape(v0,[nn*nn*pp n3*n4*n5]); |
i1=i1+nn*mm; |
52 |
i0=1; i1=nn*mm; v1{1}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]); |
v1{iFace}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]); |
53 |
i0=i1+1; i1=i1+nn*mm; v1{2}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]); |
end; |
|
i0=i1+1; i1=i1+nn*nn; v1{3}=reshape(v00(i0:i1,:),[nn nn n3 n4 n5]); |
|
|
i0=i1+1; i1=i1+nn*mm; v1{4}=reshape(v00(i0:i1,:),[mm nn n3 n4 n5]); |
|
|
i0=i1+1; i1=i1+nn*mm; v1{5}=reshape(v00(i0:i1,:),[mm nn n3 n4 n5]); |
|
|
if mygrid.nFaces==6; |
|
|
i0=i1+1; i1=i1+nn*nn; v1{6}=reshape(v00(i0:i1,:),[nn nn n3 n4 n5]); |
|
54 |
end; |
end; |
55 |
|
|
56 |
end; |
v1=gcmfaces(v1); |
57 |
|
|
58 |
if mygrid.nFaces==1; gridType='ll'; |
if ~isempty(mygrid.facesExpand); v1=convert2widefaces(v1); end; |
59 |
elseif mygrid.nFaces==5; gridType='llc'; |
|
|
elseif mygrid.nFaces==6; gridType='cube'; |
|
|
end; |
|
|
v1=gcmfaces(v1,gridType); |
|
|
|
|
60 |
else; |
else; |
61 |
|
|
62 |
[n1,n2,n3,n4,n5]=size(v0{1}); |
[n1,n2,n3,n4,n5]=size(v0{1}); |
63 |
|
|
64 |
if strcmp(mygrid.fileFormat,'straight'); |
if strcmp(mygrid.fileFormat,'straight'); |
65 |
v1=v0{1}; |
v1=v0{1}; |
66 |
elseif strcmp(mygrid.fileFormat,'cube'); |
elseif strcmp(mygrid.fileFormat,'cube'); |
67 |
v1=zeros(n2*6,n2,n3,n4,n5); |
v1=zeros(n2*6,n2,n3,n4,n5); |
68 |
for ii=1:6; v1([1:n2]+(ii-1)*n2,:,:,:,:)=v0{ii}; end; |
for ii=1:6; v1([1:n2]+(ii-1)*n2,:,:,:,:)=v0{ii}; end; |
69 |
elseif strcmp(mygrid.fileFormat,'compact'); |
elseif strcmp(mygrid.fileFormat,'compact')|strcmp(mygrid.fileFormat,'nctiles'); |
70 |
|
|
71 |
v0_faces=v0; clear v0; |
% is there a reason for this? |
72 |
for iFace=1:mygrid.nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end; |
% v0_faces=v0; clear v0; |
73 |
|
% for iFace=1:mygrid.nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end; |
74 |
nn=size(v0{1},1); mm=size(v0{1},2); |
|
75 |
pp=mm/nn*4+mygrid.nFaces-4; |
if ~isempty(mygrid.facesExpand); v0=convert2widefaces(v0); end; |
76 |
|
|
77 |
n3=size(v0{1},3); n4=size(v0{1},4); n5=size(v0{1},5); |
aa=0;bb=0; |
78 |
v1=NaN*zeros(nn,nn*pp,n3,n4,n5); |
for iFace=1:mygrid.nFaces; |
79 |
|
aa=aa+prod(mygrid.facesSize(iFace,:)); |
80 |
v11=NaN*zeros(nn*nn*pp,n3*n4*n5); |
bb=bb+prod(size(v0{iFace})); |
81 |
|
end; |
82 |
i0=1; i1=nn*mm; tmp1=reshape(v0{1},[nn*mm n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:); |
bb=bb/aa; |
83 |
i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{2},[nn*mm n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:); |
v11=NaN*zeros(aa,bb); |
84 |
i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{3},[nn*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:); |
|
85 |
i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{4},[mm*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:); |
i0=0; i1=0; |
86 |
i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{5},[mm*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:); |
for iFace=1:mygrid.nFaces; |
87 |
if mygrid.nFaces==6; |
i0=i1+1; |
88 |
i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{6},[nn*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:); |
nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2); |
89 |
end; |
i1=i1+nn*mm; |
90 |
|
v11(i0:i1,:)=reshape(v0{iFace},[nn*mm bb]); |
91 |
v1=reshape(v11,[nn nn*pp n3 n4 n5]); |
end; |
92 |
|
nn=mygrid.ioSize(1); mm=mygrid.ioSize(2); |
93 |
end; |
v1=reshape(v11,[nn mm n3 n4 n5]); |
94 |
|
|
95 |
|
end; |
96 |
|
|
97 |
end; |
end; |
98 |
|
|
99 |
|
|