matlab_class/gcmfaces_convert/convert2pcol_llpc.m

function [X,Y,FLD]=convert2pcol_llpc(x,y,fld,varargin);
%object:    gcmfaces to 'pcolor format' conversion
%inputs:    x is longitude (e.g. mygrid.XC)
%           y is latitude (e.g. mygrid.YC)
%           fld is the 2D field of interest (e.g. mygrid.hFacC(:,:,1))
%outputs:   X,Y,FLD are array versions of x,y,fld
%
%note:      this function is designed so that one may readily
%           plot the output in geographic coordinates
%           using e.g. 'figure; pcolor(X,Y,FLD);'

k3=1; k4=1;

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

%FIX DATE CHANGE LINE IN LATLON PART:
%--------------------------------------
X=circshift(X,[55 0]);
Y=circshift(Y,[55 0]);
FLD=circshift(FLD,[55 0]);

%ADD POINTS TO FIX DATE CHANGE LINE ELSEWHERE:
%---------------------------------------------
X=[X(353:end,:);X;X(1:32,:)];
Y=[Y(353:end,:);Y;Y(1:32,:)];
FLD=[FLD(353:end,:);FLD;FLD(1:32,:)];

tmp1=X(1:64,:); tmp1(tmp1>0)=tmp1(tmp1>0)-360; X(1:64,:)=tmp1;
tmp1=X(353:end,:); tmp1(tmp1<0)=tmp1(tmp1<0)+360; X(353:end,:)=tmp1;


1	function [X,Y,FLD]=convert2pcol_llpc(x,y,fld,varargin);
2	%object: gcmfaces to 'pcolor format' conversion
3	%inputs: x is longitude (e.g. mygrid.XC)
4	% y is latitude (e.g. mygrid.YC)
5	% fld is the 2D field of interest (e.g. mygrid.hFacC(:,:,1))
6	%outputs: X,Y,FLD are array versions of x,y,fld
7	%
8	%note: this function is designed so that one may readily
9	% plot the output in geographic coordinates
10	% using e.g. 'figure; pcolor(X,Y,FLD);'
11
12	k3=1; k4=1;
13
14	%ASSEMBLE "LATLON" FACES:
15	%----------------------
16	X=[x{1};x{2};sym_g(x{3},7,0); sym_g(x{4},7,0)];
17	Y=[y{1};y{2};sym_g(y{3},7,0); sym_g(y{4},7,0)];
18	FLD=[fld{1}(:,:,k3,k4);fld{2}(:,:,k3,k4);sym_g(fld{3}(:,:,k3,k4),7,0);sym_g(fld{4}(:,:,k3,k4),7,0)];
19
20	%FIX DATE CHANGE LINE IN LATLON PART:
21	%--------------------------------------
22	X=circshift(X,[55 0]);
23	Y=circshift(Y,[55 0]);
24	FLD=circshift(FLD,[55 0]);
25
26	%ADD POINTS TO FIX DATE CHANGE LINE ELSEWHERE:
27	%---------------------------------------------
28	X=[X(353:end,:);X;X(1:32,:)];
29	Y=[Y(353:end,:);Y;Y(1:32,:)];
30	FLD=[FLD(353:end,:);FLD;FLD(1:32,:)];
31
32	tmp1=X(1:64,:); tmp1(tmp1>0)=tmp1(tmp1>0)-360; X(1:64,:)=tmp1;
33	tmp1=X(353:end,:); tmp1(tmp1<0)=tmp1(tmp1<0)+360; X(353:end,:)=tmp1;
34
35