matlab_class/gcmfaces_convert/convert2arctic.m

function [fld0]=convert2arctic(fld00,varargin);

if fld00.nFaces==1; fld00=convert2cube(fld00); end;

n1=size(fld00{1},1);
n2=size(fld00{1},2);
n3=size(fld00{1},3);
n4=size(fld00{1},4);

f1=flipdim(permute(fld00{1},[2 1 3 4]),2);
f2=fld00{2};                       
f3=fld00{3};
f4=fld00{4};
f5=flipdim(permute(fld00{5},[2 1 3 4]),1);
if isempty(fld00{3}); f3=NaN*zeros(n1,n1,n3,n4); end;

fld0=NaN*zeros(n2*2+n1,n2*2+n1,n3,n4);
fld0(1:n2,n2+1:n2+n1,:,:)=f1;
fld0(n2+1:n2+n1,n2+1:n2+n1,:,:)=f3;
fld0(n2+n1+1:end,n2+1:n2+n1,:,:)=f4;
fld0(n2+1:n2+n1,1:n2,:,:)=f2;
fld0(n2+1:n2+n1,n2+n1+1:end,:,:)=f5;

nn=round(n2/2); 
nn2=n2-nn+1; nn1=n2*2+n1-nn2+1;
fld0=fld0(nn2:nn1,nn2:nn1);

if nargin==2; doFill=varargin{1}; else; doFill=1; end;
if doFill;

if n3>1; fprintf('not implemented yet\n'); return; end;
%fill the corners of fld0:
warning('off','MATLAB:interp1:NaNinY');
for pp=1:4;
%do the interpolation in polar coordinates...
  ii=[0:nn]'*ones(1,nn+1); jj=ones(nn+1,1)*[0:nn];
  aa=angle(ii+i*jj); bb=min(abs(ii+i*jj),nn);
%
  v0=fld0(end-nn:end,end-nn:end);
  k1=find(aa==0); k2=find(aa>0&aa<=pi/4); v0(k2)=interp1(bb(k1),v0(k1),bb(k2));
  k1=find(aa==pi/2); k2=find(aa>pi/4&aa<=pi/2); v0(k2)=interp1(bb(k1),v0(k1),bb(k2));
  fld0(end-nn+1:end,end-nn+1:end)=v0(2:end,2:end);
  fld0=sym_g(fld0,5,0);
%
end;
warning('on','MATLAB:interp1:NaNinY');

end;%if doFill;

1	function [fld0]=convert2arctic(fld00,varargin);
2
3	if fld00.nFaces==1; fld00=convert2cube(fld00); end;
4
5	n1=size(fld00{1},1);
6	n2=size(fld00{1},2);
7	n3=size(fld00{1},3);
8	n4=size(fld00{1},4);
9
10	f1=flipdim(permute(fld00{1},[2 1 3 4]),2);
11	f2=fld00{2};
12	f3=fld00{3};
13	f4=fld00{4};
14	f5=flipdim(permute(fld00{5},[2 1 3 4]),1);
15	if isempty(fld00{3}); f3=NaN*zeros(n1,n1,n3,n4); end;
16
17	fld0=NaNzeros(n22+n1,n2*2+n1,n3,n4);
18	fld0(1:n2,n2+1:n2+n1,:,:)=f1;
19	fld0(n2+1:n2+n1,n2+1:n2+n1,:,:)=f3;
20	fld0(n2+n1+1:end,n2+1:n2+n1,:,:)=f4;
21	fld0(n2+1:n2+n1,1:n2,:,:)=f2;
22	fld0(n2+1:n2+n1,n2+n1+1:end,:,:)=f5;
23
24	nn=round(n2/2);
25	nn2=n2-nn+1; nn1=n2*2+n1-nn2+1;
26	fld0=fld0(nn2:nn1,nn2:nn1);
27
28	if nargin==2; doFill=varargin{1}; else; doFill=1; end;
29	if doFill;
30
31	if n3>1; fprintf('not implemented yet\n'); return; end;
32	%fill the corners of fld0:
33	warning('off','MATLAB:interp1:NaNinY');
34	for pp=1:4;
35	%do the interpolation in polar coordinates...
36	ii=[0:nn]'ones(1,nn+1); jj=ones(nn+1,1)[0:nn];
37	aa=angle(ii+ijj); bb=min(abs(ii+ijj),nn);
38	%
39	v0=fld0(end-nn:end,end-nn:end);
40	k1=find(aa==0); k2=find(aa>0&aa<=pi/4); v0(k2)=interp1(bb(k1),v0(k1),bb(k2));
41	k1=find(aa==pi/2); k2=find(aa>pi/4&aa<=pi/2); v0(k2)=interp1(bb(k1),v0(k1),bb(k2));
42	fld0(end-nn+1:end,end-nn+1:end)=v0(2:end,2:end);
43	fld0=sym_g(fld0,5,0);
44	%
45	end;
46	warning('on','MATLAB:interp1:NaNinY');
47
48	end;%if doFill;
49