| 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) |
| 6 |
% |
% |
| 7 |
%note: global mygrid parameters (nFaces, fileFormat) are used |
%note: global mygrid parameters (nFaces, fileFormat) are used |
| 8 |
|
|
| 9 |
input_list_check('convert2gcmfaces',nargin); |
input_list_check('convert2gcmfaces',nargin); |
| 10 |
|
|
| 11 |
aa=whos('v0'); doGcm2Faces=strcmp(aa.class,'double'); |
aa=whos('v0'); |
| 12 |
|
if strcmp(aa.class,'single'); v0=double(v0); end; |
| 13 |
global mygrid; |
aa=whos('v0'); |
| 14 |
|
doGcm2Faces=strcmp(aa.class,'double'); |
| 15 |
if doGcm2Faces&mygrid.gcm2facesFast; |
|
| 16 |
|
global mygrid; |
| 17 |
[n1,n2,n3,n4,n5]=size(v0); |
|
| 18 |
|
if doGcm2Faces&mygrid.gcm2facesFast; |
| 19 |
v0=reshape(v0,[n1*n2 n3*n4*n5]); |
|
| 20 |
for iFace=1:mygrid.nFaces; |
[n1,n2,n3,n4,n5]=size(v0); |
| 21 |
v1{iFace}=reshape(mygrid.gcm2faces{iFace}*v0,[size(mygrid.XC{iFace}) n3 n4 n5]); |
|
| 22 |
end; |
v0=reshape(v0,[n1*n2 n3*n4*n5]); |
| 23 |
|
for iFace=1:mygrid.nFaces; |
| 24 |
v1=gcmfaces(v1); |
v1{iFace}=reshape(mygrid.gcm2faces{iFace}*v0,[size(mygrid.XC{iFace}) n3 n4 n5]); |
| 25 |
|
end; |
| 26 |
elseif ~doGcm2Faces&mygrid.gcm2facesFast; |
|
| 27 |
|
v1=gcmfaces(v1); |
| 28 |
[n1,n2,n3,n4,n5]=size(v0{1}); |
|
| 29 |
|
elseif ~doGcm2Faces&mygrid.gcm2facesFast; |
| 30 |
v1=zeros(mygrid.faces2gcmSize(1)*mygrid.faces2gcmSize(2),n3*n4*n5); |
|
| 31 |
for iFace=1:mygrid.nFaces; |
[n1,n2,n3,n4,n5]=size(v0{1}); |
| 32 |
nn=size(mygrid.XC{iFace}); |
|
| 33 |
v1(mygrid.faces2gcm{iFace},:)=reshape(v0{iFace},[nn(1)*nn(2) n3*n4*n5]); |
v1=zeros(mygrid.faces2gcmSize(1)*mygrid.faces2gcmSize(2),n3*n4*n5); |
| 34 |
end; |
for iFace=1:mygrid.nFaces; |
| 35 |
v1=reshape(v1,[mygrid.faces2gcmSize n3 n4 n5]); |
nn=size(mygrid.XC{iFace}); |
| 36 |
|
v1(mygrid.faces2gcm{iFace},:)=reshape(v0{iFace},[nn(1)*nn(2) n3*n4*n5]); |
| 37 |
elseif doGcm2Faces; |
end; |
| 38 |
|
v1=reshape(v1,[mygrid.faces2gcmSize n3 n4 n5]); |
| 39 |
[n1,n2,n3,n4,n5]=size(v0); |
|
| 40 |
|
elseif doGcm2Faces; |
| 41 |
if strcmp(mygrid.fileFormat,'straight'); |
|
| 42 |
v1={v0}; |
[n1,n2,n3,n4,n5]=size(v0); |
| 43 |
elseif strcmp(mygrid.fileFormat,'cube'); |
|
| 44 |
for ii=1:6; v1{ii}=v0(n2*(ii-1)+[1:n2],:,:,:,:); end; |
if strcmp(mygrid.fileFormat,'straight'); |
| 45 |
elseif strcmp(mygrid.fileFormat,'compact')|strcmp(mygrid.fileFormat,'nctiles'); |
v1={v0}; |
| 46 |
v00=reshape(v0,[n1*n2 n3*n4*n5]); |
elseif strcmp(mygrid.fileFormat,'cube'); |
| 47 |
i0=0; i1=0; |
for ii=1:6; v1{ii}=v0(n2*(ii-1)+[1:n2],:,:,:,:); end; |
| 48 |
for iFace=1:mygrid.nFaces; |
elseif strcmp(mygrid.fileFormat,'compact')|strcmp(mygrid.fileFormat,'nctiles'); |
| 49 |
i0=i1+1; |
v00=reshape(v0,[n1*n2 n3*n4*n5]); |
| 50 |
nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2); |
i0=0; i1=0; |
| 51 |
i1=i1+nn*mm; |
for iFace=1:mygrid.nFaces; |
| 52 |
v1{iFace}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]); |
i0=i1+1; |
| 53 |
end; |
nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2); |
| 54 |
end; |
i1=i1+nn*mm; |
| 55 |
|
v1{iFace}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]); |
| 56 |
v1=gcmfaces(v1); |
end; |
| 57 |
|
end; |
| 58 |
if ~isempty(mygrid.facesExpand); v1=convert2widefaces(v1); end; |
|
| 59 |
|
v1=gcmfaces(v1); |
| 60 |
else; |
|
| 61 |
|
if ~isempty(mygrid.facesExpand); v1=convert2widefaces(v1); end; |
| 62 |
[n1,n2,n3,n4,n5]=size(v0{1}); |
|
| 63 |
|
else; |
| 64 |
if strcmp(mygrid.fileFormat,'straight'); |
|
| 65 |
v1=v0{1}; |
[n1,n2,n3,n4,n5]=size(v0{1}); |
| 66 |
elseif strcmp(mygrid.fileFormat,'cube'); |
|
| 67 |
v1=zeros(n2*6,n2,n3,n4,n5); |
if strcmp(mygrid.fileFormat,'straight'); |
| 68 |
for ii=1:6; v1([1:n2]+(ii-1)*n2,:,:,:,:)=v0{ii}; end; |
v1=v0{1}; |
| 69 |
elseif strcmp(mygrid.fileFormat,'compact')|strcmp(mygrid.fileFormat,'nctiles'); |
elseif strcmp(mygrid.fileFormat,'cube'); |
| 70 |
|
v1=zeros(n2*6,n2,n3,n4,n5); |
| 71 |
% is there a reason for this? |
for ii=1:6; v1([1:n2]+(ii-1)*n2,:,:,:,:)=v0{ii}; end; |
| 72 |
% v0_faces=v0; clear v0; |
elseif strcmp(mygrid.fileFormat,'compact')|strcmp(mygrid.fileFormat,'nctiles'); |
| 73 |
% for iFace=1:mygrid.nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end; |
|
| 74 |
|
% is there a reason for this? |
| 75 |
if ~isempty(mygrid.facesExpand); v0=convert2widefaces(v0); end; |
% v0_faces=v0; clear v0; |
| 76 |
|
% for iFace=1:mygrid.nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end; |
| 77 |
aa=0;bb=0; |
|
| 78 |
for iFace=1:mygrid.nFaces; |
if ~isempty(mygrid.facesExpand); v0=convert2widefaces(v0); end; |
| 79 |
aa=aa+prod(mygrid.facesSize(iFace,:)); |
|
| 80 |
bb=bb+prod(size(v0{iFace})); |
aa=0;bb=0; |
| 81 |
end; |
for iFace=1:mygrid.nFaces; |
| 82 |
bb=bb/aa; |
aa=aa+prod(mygrid.facesSize(iFace,:)); |
| 83 |
v11=NaN*zeros(aa,bb); |
bb=bb+prod(size(v0{iFace})); |
| 84 |
|
end; |
| 85 |
i0=0; i1=0; |
bb=bb/aa; |
| 86 |
for iFace=1:mygrid.nFaces; |
v11=NaN*zeros(aa,bb); |
| 87 |
i0=i1+1; |
|
| 88 |
nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2); |
i0=0; i1=0; |
| 89 |
i1=i1+nn*mm; |
for iFace=1:mygrid.nFaces; |
| 90 |
v11(i0:i1,:)=reshape(v0{iFace},[nn*mm bb]); |
i0=i1+1; |
| 91 |
end; |
nn=mygrid.facesSize(iFace,1); mm=mygrid.facesSize(iFace,2); |
| 92 |
nn=mygrid.ioSize(1); mm=mygrid.ioSize(2); |
i1=i1+nn*mm; |
| 93 |
v1=reshape(v11,[nn mm n3 n4 n5]); |
v11(i0:i1,:)=reshape(v0{iFace},[nn*mm bb]); |
| 94 |
|
end; |
| 95 |
end; |
nn=mygrid.ioSize(1); mm=mygrid.ioSize(2); |
| 96 |
|
v1=reshape(v11,[nn mm n3 n4 n5]); |
| 97 |
end; |
|
| 98 |
|
end; |
| 99 |
|
|
| 100 |
|
end; |
| 101 |
|
|
| 102 |
|
|
Parent Directory
| 
Revision Log
| 
Revision Graph
| 
Patch