/[MITgcm]/MITgcm_contrib/gael/matlab_class/gcmfaces_IO/convert2gcmfaces.m
ViewVC logotype

Contents of /MITgcm_contrib/gael/matlab_class/gcmfaces_IO/convert2gcmfaces.m

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


Revision 1.6 - (show annotations) (download)
Fri Jul 1 15:18:18 2011 UTC (14 years ago) by gforget
Branch: MAIN
Changes since 1.5: +22 -22 lines
- add fifth dimension.

1 function [v1]=convert2gcmfaces(v0,varargin);
2 %object: converts model output to gcmfaces object
3 % or vice versa
4 %input: v0 model output array (resp. gcmfaces object)
5 %output: v1 gcmfaces object (model output array)
6 %
7 %note: global mygrid parameters (nFaces, fileFormat) are used
8
9 input_list_check('convert2gcmfaces',nargin);
10
11 aa=whos('v0'); doGcm2Faces=strcmp(aa.class,'double');
12
13 global mygrid;
14
15 if doGcm2Faces;
16
17 [n1,n2,n3,n4,n5]=size(v0);
18
19 if strcmp(mygrid.fileFormat,'straight');
20 v1={v0};
21 elseif strcmp(mygrid.fileFormat,'cube');
22 for ii=1:6; v1{ii}=v0(n2*(ii-1)+[1:n2],:,:,:,:); end;
23 elseif strcmp(mygrid.fileFormat,'compact');
24
25 nn=size(v0,1);
26 pp=size(v0,2)/nn;
27 mm=(pp+4-mygrid.nFaces)/4*nn;
28
29 v00=reshape(v0,[nn*nn*pp n3*n4*n5]);
30 i0=1; i1=nn*mm; v1{1}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]);
31 i0=i1+1; i1=i1+nn*mm; v1{2}=reshape(v00(i0:i1,:),[nn mm n3 n4 n5]);
32 i0=i1+1; i1=i1+nn*nn; v1{3}=reshape(v00(i0:i1,:),[nn nn n3 n4 n5]);
33 i0=i1+1; i1=i1+nn*mm; v1{4}=reshape(v00(i0:i1,:),[mm nn n3 n4 n5]);
34 i0=i1+1; i1=i1+nn*mm; v1{5}=reshape(v00(i0:i1,:),[mm nn n3 n4 n5]);
35 if mygrid.nFaces==6;
36 i0=i1+1; i1=i1+nn*nn; v1{6}=reshape(v00(i0:i1,:),[nn nn n3 n4 n5]);
37 end;
38
39 end;
40
41 if mygrid.nFaces==1; gridType='ll';
42 elseif mygrid.nFaces==5; gridType='llc';
43 elseif mygrid.nFaces==6; gridType='cube';
44 end;
45 v1=gcmfaces(v1,gridType);
46
47 else;
48
49 [n1,n2,n3,n4,n5]=size(v0{1});
50
51 if strcmp(mygrid.fileFormat,'straight');
52 v1=v0{1};
53 elseif strcmp(mygrid.fileFormat,'cube');
54 v1=zeros(n2*6,n2,n3,n4,n5);
55 for ii=1:6; v1([1:n2]+(ii-1)*n2,:,:,:,:)=v0{ii}; end;
56 elseif strcmp(mygrid.fileFormat,'compact');
57
58 v0_faces=v0; clear v0;
59 for iFace=1:mygrid.nFaces; eval(['v0{iFace}=get(v0_faces,''f' num2str(iFace) ''');']); end;
60
61 nn=size(v0{1},1); mm=size(v0{1},2);
62 pp=mm/nn*4+mygrid.nFaces-4;
63
64 n3=size(v0{1},3); n4=size(v0{1},4); n5=size(v0{1},5);
65 v1=NaN*zeros(nn,nn*pp,n3,n4,n5);
66
67 v11=NaN*zeros(nn*nn*pp,n3*n4*n5);
68
69 i0=1; i1=nn*mm; tmp1=reshape(v0{1},[nn*mm n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
70 i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{2},[nn*mm n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
71 i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{3},[nn*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
72 i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{4},[mm*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
73 i0=i1+1; i1=i1+nn*mm; tmp1=reshape(v0{5},[mm*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
74 if mygrid.nFaces==6;
75 i0=i1+1; i1=i1+nn*nn; tmp1=reshape(v0{6},[nn*nn n3*n4*n5]); v11(i0:i1,:)=tmp1(:,:);
76 end;
77
78 v1=reshape(v11,[nn nn*pp n3 n4 n5]);
79
80 end;
81
82 end;
83
84

  ViewVC Help
Powered by ViewVC 1.1.22