matlab_class/gcmfaces_calc/gcmfaces_lines_zonal.m

function []=gcmfaces_lines_zonal(varargin);
%object:    define the set of quasi longitudinal lines
%           along which transports will integrated;
%           LATS_MASKS that will be added to mygrid.
%optional input:
%           LATS_VAL is the latitudes vector ([-89:89]' by default)
%           MSKS_NAM is the list of fields to store in mygrid

global mygrid;

if nargin>0; LATS_VAL=varargin{1}; else; LATS_VAL=[-89:89]'; end;
if nargin>1; MSKS_NAM=varargin{2}; else; MSKS_NAM={'mskCint','mskCedge','mskWedge','mskSedge'}; end;

for iy=1:length(LATS_VAL);

    mskCint=1*(mygrid.YC>=LATS_VAL(iy));
    [mskCedge,mskWedge,mskSedge]=gcmfaces_edge_mask(mskCint);

    for im=1:length(MSKS_NAM);
      eval(['tmp1.' MSKS_NAM{im} '=' MSKS_NAM{im} ';']);
    end;
    tmp1.lat=LATS_VAL(iy);

    %store:
    if iy==1;
        LATS_MASKS=tmp1;
    else;
        LATS_MASKS(iy)=tmp1;
    end;
    
end;

mygrid.LATS_MASKS=LATS_MASKS;


1	gforget	1.1	function []=gcmfaces_lines_zonal(varargin);
2			%object: define the set of quasi longitudinal lines
3			% along which transports will integrated;
4			% LATS_MASKS that will be added to mygrid.
5			%optional input:
6			% LATS_VAL is the latitudes vector ([-89:89]' by default)
7	gforget	1.3	% MSKS_NAM is the list of fields to store in mygrid
8	gforget	1.1
9			global mygrid;
10
11			if nargin>0; LATS_VAL=varargin{1}; else; LATS_VAL=[-89:89]'; end;
12	gforget	1.3	if nargin>1; MSKS_NAM=varargin{2}; else; MSKS_NAM={'mskCint','mskCedge','mskWedge','mskSedge'}; end;
13	gforget	1.1
14			for iy=1:length(LATS_VAL);
15	gforget	1.2
16			mskCint=1*(mygrid.YC>=LATS_VAL(iy));
17			[mskCedge,mskWedge,mskSedge]=gcmfaces_edge_mask(mskCint);
18
19	gforget	1.3	for im=1:length(MSKS_NAM);
20			eval(['tmp1.' MSKS_NAM{im} '=' MSKS_NAM{im} ';']);
21			end;
22			tmp1.lat=LATS_VAL(iy);
23
24	gforget	1.1	%store:
25			if iy==1;
26	gforget	1.3	LATS_MASKS=tmp1;
27	gforget	1.1	else;
28	gforget	1.3	LATS_MASKS(iy)=tmp1;
29	gforget	1.1	end;
30
31			end;
32
33			mygrid.LATS_MASKS=LATS_MASKS;
34
35