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.
%           LONS_MASKS that will be added to mygrid.
%optional input:
%           LATS_VAL is the latitudes vector ([-89:89]' by default)
%           LONS_VAL is the latitudes vector ([-179.5:179.5]' by default)

gcmfaces_global;

if nargin>0; LATS_VAL=varargin{1}; else; LATS_VAL=[-89:89]'; end;
if nargin>1; LONS_VAL=varargin{2}; else; LONS_VAL=[-179.5:179.5]'; end;

if iscell(LONS_VAL);
    warning('gcmfaces_lines_zonal now takes LONS_VAL as 2nd argument');
    LONS_VAL=[-179.5:179.5]';
end;

MSKS_NAM={'mskCint','mskCedge','mskWedge','mskSedge'};    

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;

d2r=pi/180;
for iy=1:length(LONS_VAL);

    mskCint=1*(sin(d2r*(mygrid.XC-LONS_VAL(iy)))>=0);
    [mskCedge,mskWedge,mskSedge]=gcmfaces_edge_mask(mskCint);
    tmp1=1*(cos(d2r*(mygrid.XC-LONS_VAL(iy)))>0);
    mskCedge=mskCedge.*tmp1;
    tmp1=1*(cos(d2r*(mygrid.XG-LONS_VAL(iy)))>0);
    mskWedge=mskWedge.*tmp1;
    tmp1=1*(cos(d2r*(mygrid.XC-LONS_VAL(iy)))>0);
    mskSedge=mskSedge.*tmp1;

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

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

mygrid.LATS_MASKS=LATS_MASKS;
mygrid.LONS_MASKS=LONS_MASKS;
mygrid.LATS=[mygrid.LATS_MASKS.lat]';
mygrid.LONS=[mygrid.LONS_MASKS.lon]';

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	% LONS_MASKS that will be added to mygrid.
6	%optional input:
7	% LATS_VAL is the latitudes vector ([-89:89]' by default)
8	% LONS_VAL is the latitudes vector ([-179.5:179.5]' by default)
9
10	gcmfaces_global;
11
12	if nargin>0; LATS_VAL=varargin{1}; else; LATS_VAL=[-89:89]'; end;
13	if nargin>1; LONS_VAL=varargin{2}; else; LONS_VAL=[-179.5:179.5]'; end;
14
15	if iscell(LONS_VAL);
16	warning('gcmfaces_lines_zonal now takes LONS_VAL as 2nd argument');
17	LONS_VAL=[-179.5:179.5]';
18	end;
19
20	MSKS_NAM={'mskCint','mskCedge','mskWedge','mskSedge'};
21
22	for iy=1:length(LATS_VAL);
23
24	mskCint=1*(mygrid.YC>=LATS_VAL(iy));
25	[mskCedge,mskWedge,mskSedge]=gcmfaces_edge_mask(mskCint);
26
27	for im=1:length(MSKS_NAM);
28	eval(['tmp1.' MSKS_NAM{im} '=' MSKS_NAM{im} ';']);
29	end;
30	tmp1.lat=LATS_VAL(iy);
31
32	%store:
33	if iy==1;
34	LATS_MASKS=tmp1;
35	else;
36	LATS_MASKS(iy)=tmp1;
37	end;
38
39	end;
40
41	d2r=pi/180;
42	for iy=1:length(LONS_VAL);
43
44	mskCint=1(sin(d2r(mygrid.XC-LONS_VAL(iy)))>=0);
45	[mskCedge,mskWedge,mskSedge]=gcmfaces_edge_mask(mskCint);
46	tmp1=1(cos(d2r(mygrid.XC-LONS_VAL(iy)))>0);
47	mskCedge=mskCedge.*tmp1;
48	tmp1=1(cos(d2r(mygrid.XG-LONS_VAL(iy)))>0);
49	mskWedge=mskWedge.*tmp1;
50	tmp1=1(cos(d2r(mygrid.XC-LONS_VAL(iy)))>0);
51	mskSedge=mskSedge.*tmp1;
52
53	clear tmp1;
54	for im=1:length(MSKS_NAM);
55	eval(['tmp1.' MSKS_NAM{im} '=' MSKS_NAM{im} ';']);
56	end;
57	tmp1.lon=LONS_VAL(iy);
58
59	%store:
60	if iy==1;
61	LONS_MASKS=tmp1;
62	else;
63	LONS_MASKS(iy)=tmp1;
64	end;
65
66	end;
67
68	mygrid.LATS_MASKS=LATS_MASKS;
69	mygrid.LONS_MASKS=LONS_MASKS;
70	mygrid.LATS=[mygrid.LATS_MASKS.lat]';
71	mygrid.LONS=[mygrid.LONS_MASKS.lon]';
72