matlab_class/gcmfaces_exch/exch_Z_llc.m

function [FLD]=exch_Z_llc(fld);
%[FLD]=exch_Z_llc(fld);
%adds vorticity points (to north and east of center points) for
%lat-lon cap 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)=NaN;
        FLD.f5(end,1:end-1,k3,k4)=NaN;
        
        %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));
        
        %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))/3;
        zzC(2)=(fld.f2(end,1,k3,k4)+fld.f4(end,1,k3,k4))/2;
        
        %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);
            warning('exch_Z_llc filled 2 missing corners assuming uncertain symmetries!');
        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);
        
    end;
end;

1;

1	function [FLD]=exch_Z_llc(fld);
2	%[FLD]=exch_Z_llc(fld);
3	%adds vorticity points (to north and east of center points) for
4	%lat-lon cap 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)=NaN;
35	FLD.f5(end,1:end-1,k3,k4)=NaN;
36
37	%overlap in j+1
38	FLD.f1(2:end,end,k3,k4)=fliplr(fld.f3(1,:,k3,k4));
39	FLD.f2(1:end-1,end,k3,k4)=fld.f3(:,1,k3,k4);
40	FLD.f3(2:end,end,k3,k4)=fliplr(fld.f5(1,:,k3,k4));
41	FLD.f4(1:end-1,end,k3,k4)=fld.f5(:,1,k3,k4);
42	FLD.f5(2:end,end,k3,k4)=fliplr(fld.f1(1,:,k3,k4));
43
44	%missing corners : use the average value (from the 3 or 2 neighbours)
45	zzC(1)=(fld.f1(1,end,k3,k4)+fld.f3(1,end,k3,k4)+fld.f5(1,end,k3,k4))/3;
46	zzC(2)=(fld.f2(end,1,k3,k4)+fld.f4(end,1,k3,k4))/2;
47
48	%missing corners : assume symmetry of the grid (that may be completly wrong!)
49	if exch_Z_assume_sym;
50	zzC(1)=fld.f4(1,1,k3,k4);
51	zzC(2)=fld.f1(1,1,k3,k4);
52	warning('exch_Z_llc filled 2 missing corners assuming uncertain symmetries!');
53	end
54
55	%- 1rst = N.W corner of face 1
56	FLD.f1(1,end,k3,k4)=zzC(1);
57	FLD.f3(1,end,k3,k4)=zzC(1);
58	FLD.f5(1,end,k3,k4)=zzC(1);
59	%- 2nd = S.E corner of face 2
60	FLD.f2(end,1,k3,k4)=zzC(2);
61	FLD.f4(end,1,k3,k4)=zzC(2);
62
63	end;
64	end;
65
66	1;
67