matlab_class/gcmfaces_convert/convert2pcol_llc.m

function [X,Y,FLD]=convert2pcol_llc(x,y,fld,varargin);

if nargin>3; k3=varargin{1}; else; k3=1; end;
if nargin>4; k4=varargin{2}; else; k4=1; end;

%ASSEMBLE "LATLON" FACES:
%----------------------
X=[x{1};x{2};sym_g(x{4},7,0); sym_g(x{5},7,0)];
Y=[y{1};y{2};sym_g(y{4},7,0); sym_g(y{5},7,0)];
FLD=[fld{1}(:,:,k3,k4);fld{2}(:,:,k3,k4);sym_g(fld{4}(:,:,k3,k4),7,0);sym_g(fld{5}(:,:,k3,k4),7,0)];

%ADD POLAR CAP:
%--------------
pp=x{3}; 
M1=pp; M1(M1>M1(1,1))=NaN; M1(M1<M1(1,end))=NaN; M1(~isnan(M1))=1;
%M1=NaN*pp;
%M2=pp; M2(M2<M2(1,1))=NaN; M2(M2>M2(end,1))=NaN; M2(~isnan(M2))=1;
M2=pp; M2(pp<0)=NaN; M2(pp>=0)=1;
M3=pp; M3(M3<M3(end,1)&M3>M3(end,end))=NaN; M3(~isnan(M3))=1;
%M3=NaN*pp;
%M4=pp; M4(M4<M4(end,end))=NaN; M4(M4>M4(1,end))=NaN; M4(~isnan(M4))=1;
M4=pp; M4(pp>=0)=NaN; M4(pp<0)=1;

pp=x{3}; Xp=[sym_g(pp.*M1,5,0);pp.*M2;sym_g(pp.*M3,7,0);sym_g(pp.*M4,6,0)]; X=[X Xp];
pp=y{3}; Yp=[sym_g(pp.*M1,5,0);pp.*M2;sym_g(pp.*M3,7,0);sym_g(pp.*M4,6,0)]; Y=[Y Yp];
pp=fld{3}(:,:,k3,k4); FLDp=[sym_g(pp.*M1,5,0);pp.*M2;sym_g(pp.*M3,7,0);sym_g(pp.*M4,6,0)]; FLD=[FLD FLDp];

%FIX DATE CHANGE LINE IN LATLON PART:
%--------------------------------------
s1=size(X); s1=round(s1(2)/2); s2=find(diff(X(:,s1))<-180); 
X=circshift(X,[-s2 0]);
Y=circshift(Y,[-s2 0]);
FLD=circshift(FLD,[-s2 0]);

%ADD POINTS TO FIX DATE CHANGE LINE ELSEWHERE:
%---------------------------------------------
s1=size(X); s1=round(s1(1)/12); 
X=[X(end-s1+1:end,:);X;X(1:s1,:)];
Y=[Y(end-s1+1:end,:);Y;Y(1:s1,:)];
FLD=[FLD(end-s1+1:end,:);FLD;FLD(1:s1,:)];

tmp1=X(1:2*s1,:); tmp1(tmp1>0)=tmp1(tmp1>0)-360; X(1:2*s1,:)=tmp1;
tmp1=X(end-2*s1+1:end,:); tmp1(tmp1<0)=tmp1(tmp1<0)+360; X(end-2*s1+1:end,:)=tmp1;


1	gforget	1.1	function [X,Y,FLD]=convert2pcol_llc(x,y,fld,varargin);
2
3			if nargin>3; k3=varargin{1}; else; k3=1; end;
4			if nargin>4; k4=varargin{2}; else; k4=1; end;
5
6			%ASSEMBLE "LATLON" FACES:
7			%----------------------
8			X=[x{1};x{2};sym_g(x{4},7,0); sym_g(x{5},7,0)];
9			Y=[y{1};y{2};sym_g(y{4},7,0); sym_g(y{5},7,0)];
10			FLD=[fld{1}(:,:,k3,k4);fld{2}(:,:,k3,k4);sym_g(fld{4}(:,:,k3,k4),7,0);sym_g(fld{5}(:,:,k3,k4),7,0)];
11
12			%ADD POLAR CAP:
13			%--------------
14			pp=x{3};
15			M1=pp; M1(M1>M1(1,1))=NaN; M1(M1<M1(1,end))=NaN; M1(~isnan(M1))=1;
16			%M1=NaN*pp;
17			%M2=pp; M2(M2<M2(1,1))=NaN; M2(M2>M2(end,1))=NaN; M2(~isnan(M2))=1;
18			M2=pp; M2(pp<0)=NaN; M2(pp>=0)=1;
19			M3=pp; M3(M3<M3(end,1)&M3>M3(end,end))=NaN; M3(~isnan(M3))=1;
20			%M3=NaN*pp;
21			%M4=pp; M4(M4<M4(end,end))=NaN; M4(M4>M4(1,end))=NaN; M4(~isnan(M4))=1;
22			M4=pp; M4(pp>=0)=NaN; M4(pp<0)=1;
23
24			pp=x{3}; Xp=[sym_g(pp.M1,5,0);pp.M2;sym_g(pp.M3,7,0);sym_g(pp.M4,6,0)]; X=[X Xp];
25			pp=y{3}; Yp=[sym_g(pp.M1,5,0);pp.M2;sym_g(pp.M3,7,0);sym_g(pp.M4,6,0)]; Y=[Y Yp];
26			pp=fld{3}(:,:,k3,k4); FLDp=[sym_g(pp.M1,5,0);pp.M2;sym_g(pp.M3,7,0);sym_g(pp.M4,6,0)]; FLD=[FLD FLDp];
27
28			%FIX DATE CHANGE LINE IN LATLON PART:
29			%--------------------------------------
30			s1=size(X); s1=round(s1(2)/2); s2=find(diff(X(:,s1))<-180);
31			X=circshift(X,[-s2 0]);
32			Y=circshift(Y,[-s2 0]);
33			FLD=circshift(FLD,[-s2 0]);
34
35			%ADD POINTS TO FIX DATE CHANGE LINE ELSEWHERE:
36			%---------------------------------------------
37			s1=size(X); s1=round(s1(1)/12);
38			X=[X(end-s1+1:end,:);X;X(1:s1,:)];
39			Y=[Y(end-s1+1:end,:);Y;Y(1:s1,:)];
40			FLD=[FLD(end-s1+1:end,:);FLD;FLD(1:s1,:)];
41
42			tmp1=X(1:2s1,:); tmp1(tmp1>0)=tmp1(tmp1>0)-360; X(1:2s1,:)=tmp1;
43			tmp1=X(end-2s1+1:end,:); tmp1(tmp1<0)=tmp1(tmp1<0)+360; X(end-2s1+1:end,:)=tmp1;
44
45