| 1 |
gforget |
1.1 |
function [X,Y,FLD]=convert2pcol(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 |
|
|
input_list_check('convert2pcol',nargin); |
| 13 |
|
|
|
| 14 |
gforget |
1.2 |
gcmfaces_global; |
| 15 |
|
|
|
| 16 |
gforget |
1.1 |
c=varargin{1}; |
| 17 |
|
|
|
| 18 |
gforget |
1.2 |
if isfield(mygrid,'xtrct'); |
| 19 |
|
|
[X,Y,FLD]=convert2pcol_xtrct(varargin{:}); |
| 20 |
|
|
elseif strcmp(c.gridType,'llc'); |
| 21 |
gforget |
1.1 |
[X,Y,FLD]=convert2pcol_llc(varargin{:}); |
| 22 |
|
|
elseif strcmp(c.gridType,'cube'); |
| 23 |
|
|
[X,Y,FLD]=convert2pcol_cube(varargin{:}); |
| 24 |
|
|
elseif strcmp(c.gridType,'llpc'); |
| 25 |
|
|
[X,Y,FLD]=convert2pcol_llpc(varargin{:}); |
| 26 |
|
|
elseif strcmp(c.gridType,'ll'); |
| 27 |
|
|
[X,Y,FLD]=convert2pcol_ll(varargin{:}); |
| 28 |
|
|
else; |
| 29 |
|
|
error(['convert2pcol not implemented for ' c.gridType '!?']); |
| 30 |
|
|
end; |
| 31 |
|
|
|