matlab_class/gcmfaces_exch/exch_Z_cube.m

function [FLD]=exch_Z_cub(fld);
%[FLD]=exch_Z_cub(fld, [flag]);
%adds vorticity points (to north and east of center points) for
%cubed-sphere grid

gcmfaces_global;

%missing corners : assume symmetry of the grid (that may be completly wrong!)
global exch_Z_assume_sym;
if isempty(exch_Z_assume_sym);
    exch_Z_assume_sym=0;
end;

%determine vertical and/or time dimensions
n3=max(size(fld.f1,3),1); n4=max(size(fld.f1,4),1);

%initialize FLD
FLD=gcmfaces;
for iF=1:mygrid.nFaces
    [n1,n2]=size(mygrid.XC{iF});
    FLD{iF}=NaN*zeros([n1+1 n2+1 n3 n4]);
end;

for k3=1:n3; for k4=1:n4;
        
        for iF=1:mygrid.nFaces
            FLD{iF}(1:end-1,1:end-1,k3,k4)=fld{iF}(:,:,k3,k4);
        end
        
        %overlap in i+1
        FLD.f1(end,1:end-1,k3,k4)=fld.f2(1,:,k3,k4);
        FLD.f2(end,2:end,k3,k4)=flipud(fld.f4(:,1,k3,k4));
        FLD.f3(end,1:end-1,k3,k4)=fld.f4(1,:,k3,k4);
        FLD.f4(end,2:end,k3,k4)=flipud(fld.f6(:,1,k3,k4));
        FLD.f5(end,1:end-1,k3,k4)=fld.f6(1,:,k3,k4);
        FLD.f6(end,2:end,k3,k4)=flipud(fld.f2(:,1,k3,k4));
        
        %overlap in j+1
        FLD.f1(2:end,end,k3,k4)=fliplr(fld.f3(1,:,k3,k4));
        FLD.f2(1:end-1,end,k3,k4)=fld.f3(:,1,k3,k4);
        FLD.f3(2:end,end,k3,k4)=fliplr(fld.f5(1,:,k3,k4));
        FLD.f4(1:end-1,end,k3,k4)=fld.f5(:,1,k3,k4);
        FLD.f5(2:end,end,k3,k4)=fliplr(fld.f1(1,:,k3,k4));
        FLD.f6(1:end-1,end,k3,k4)=fld.f1(:,1,k3,k4);
        
        %missing corners : use the average value (from the 3 or 2 neighbours)
        zzC(1)=fld.f1(1,end,k3,k4)+fld.f3(1,end,k3,k4)+fld.f5(1,end,k3,k4);
        zzC(2)=fld.f2(end,1,k3,k4)+fld.f4(end,1,k3,k4)+fld.f6(end,1,k3,k4);
        zzC=zzC/3;
        
        %missing corners : assume symmetry of the grid (that may be completly wrong!)
        if exch_Z_assume_sym;
            zzC(1)=fld.f4(1,1,k3,k4);
            zzC(2)=fld.f1(1,1,k3,k4);
        end
        
        %- 1rst = N.W corner of face 1
        FLD.f1(1,end,k3,k4)=zzC(1);
        FLD.f3(1,end,k3,k4)=zzC(1);
        FLD.f5(1,end,k3,k4)=zzC(1);
        %- 2nd  = S.E corner of face 2
        FLD.f2(end,1,k3,k4)=zzC(2);
        FLD.f4(end,1,k3,k4)=zzC(2);
        FLD.f6(end,1,k3,k4)=zzC(2);
        
    end; 
end;

1	function [FLD]=exch_Z_cub(fld);
2	%[FLD]=exch_Z_cub(fld, [flag]);
3	%adds vorticity points (to north and east of center points) for
4	%cubed-sphere grid
5
6	gcmfaces_global;
7
8	%missing corners : assume symmetry of the grid (that may be completly wrong!)
9	global exch_Z_assume_sym;
10	if isempty(exch_Z_assume_sym);
11	exch_Z_assume_sym=0;
12	end;
13
14	%determine vertical and/or time dimensions
15	n3=max(size(fld.f1,3),1); n4=max(size(fld.f1,4),1);
16
17	%initialize FLD
18	FLD=gcmfaces;
19	for iF=1:mygrid.nFaces
20	[n1,n2]=size(mygrid.XC{iF});
21	FLD{iF}=NaN*zeros([n1+1 n2+1 n3 n4]);
22	end;
23
24	for k3=1:n3; for k4=1:n4;
25
26	for iF=1:mygrid.nFaces
27	FLD{iF}(1:end-1,1:end-1,k3,k4)=fld{iF}(:,:,k3,k4);
28	end
29
30	%overlap in i+1
31	FLD.f1(end,1:end-1,k3,k4)=fld.f2(1,:,k3,k4);
32	FLD.f2(end,2:end,k3,k4)=flipud(fld.f4(:,1,k3,k4));
33	FLD.f3(end,1:end-1,k3,k4)=fld.f4(1,:,k3,k4);
34	FLD.f4(end,2:end,k3,k4)=flipud(fld.f6(:,1,k3,k4));
35	FLD.f5(end,1:end-1,k3,k4)=fld.f6(1,:,k3,k4);
36	FLD.f6(end,2:end,k3,k4)=flipud(fld.f2(:,1,k3,k4));
37
38	%overlap in j+1
39	FLD.f1(2:end,end,k3,k4)=fliplr(fld.f3(1,:,k3,k4));
40	FLD.f2(1:end-1,end,k3,k4)=fld.f3(:,1,k3,k4);
41	FLD.f3(2:end,end,k3,k4)=fliplr(fld.f5(1,:,k3,k4));
42	FLD.f4(1:end-1,end,k3,k4)=fld.f5(:,1,k3,k4);
43	FLD.f5(2:end,end,k3,k4)=fliplr(fld.f1(1,:,k3,k4));
44	FLD.f6(1:end-1,end,k3,k4)=fld.f1(:,1,k3,k4);
45
46	%missing corners : use the average value (from the 3 or 2 neighbours)
47	zzC(1)=fld.f1(1,end,k3,k4)+fld.f3(1,end,k3,k4)+fld.f5(1,end,k3,k4);
48	zzC(2)=fld.f2(end,1,k3,k4)+fld.f4(end,1,k3,k4)+fld.f6(end,1,k3,k4);
49	zzC=zzC/3;
50
51	%missing corners : assume symmetry of the grid (that may be completly wrong!)
52	if exch_Z_assume_sym;
53	zzC(1)=fld.f4(1,1,k3,k4);
54	zzC(2)=fld.f1(1,1,k3,k4);
55	end
56
57	%- 1rst = N.W corner of face 1
58	FLD.f1(1,end,k3,k4)=zzC(1);
59	FLD.f3(1,end,k3,k4)=zzC(1);
60	FLD.f5(1,end,k3,k4)=zzC(1);
61	%- 2nd = S.E corner of face 2
62	FLD.f2(end,1,k3,k4)=zzC(2);
63	FLD.f4(end,1,k3,k4)=zzC(2);
64	FLD.f6(end,1,k3,k4)=zzC(2);
65
66	end;
67	end;
68