dfer/matlab_stuff/CreateDDrmask.m

function [msk_SB,msk_LB] = CreateDDrmask(grdO);

% function [msk_SB,msk_LB] = CreateDDrmask(grdO);
% Create mask for Small and Large Basin of Double Drake 

nc=24;
nr=length(grdO.drF);

maskW=ceil(grdO.HFacW(1:6*nc,1:nc,:));
maskS=ceil(grdO.HFacS(1:6*nc,1:nc,:));
maskC=ceil(grdO.HFacC(1:6*nc,1:nc,:));

%%%% Southern most latitudes of walls
%min(grdO.YC( find(grdO.HFacC(1:6*nc,1:nc,1) == 0) ))     %% C-points
%min(grdO.YC( find(grdO.HFacW(1:6*nc,1:nc,1) == 0) ))     %% U-points
%min(grdO.YG( find(grdO.HFacS(1:6*nc,1:nc,1) == 0) ))     %% V-points

%%%% Small-Basin Masks

mskSB_G=ones(6*nc,nc); mskSB_G(grdO.XG(1:6*nc,1:nc) < 90 )=0;
mskSB_C=ones(6*nc,nc); mskSB_C(grdO.XC(1:6*nc,1:nc) < 90 )=0;

maskW_SB = maskW .* repmat(mskSB_G,[1 1 nr]);
maskS_SB = maskS .* repmat(mskSB_C,[1 1 nr]);
maskC_SB = maskC .* repmat(mskSB_C,[1 1 nr]);

msk_SB.maskW=maskW_SB;
msk_SB.maskS=maskS_SB;
msk_SB.maskC=maskC_SB;

%%%% Large-Basin Masks

mskLB_G=ones(6*nc,nc); mskLB_G(grdO.XG(1:6*nc,1:nc) > 90 )=0;
mskLB_C=ones(6*nc,nc); mskLB_C(grdO.XC(1:6*nc,1:nc) > 90 )=0;

maskW_LB = maskW .* repmat(mskLB_G,[1 1 nr]);
maskS_LB = maskS .* repmat(mskLB_C,[1 1 nr]);
maskC_LB = maskC .* repmat(mskLB_C,[1 1 nr]);

msk_LB.maskW=maskW_LB;
msk_LB.maskS=maskS_LB;
msk_LB.maskC=maskC_LB;

return
1	function [msk_SB,msk_LB] = CreateDDrmask(grdO);
2
3	% function [msk_SB,msk_LB] = CreateDDrmask(grdO);
4	% Create mask for Small and Large Basin of Double Drake
5
6	nc=24;
7	nr=length(grdO.drF);
8
9	maskW=ceil(grdO.HFacW(1:6*nc,1:nc,:));
10	maskS=ceil(grdO.HFacS(1:6*nc,1:nc,:));
11	maskC=ceil(grdO.HFacC(1:6*nc,1:nc,:));
12
13	%%%% Southern most latitudes of walls
14	%min(grdO.YC( find(grdO.HFacC(1:6*nc,1:nc,1) == 0) )) %% C-points
15	%min(grdO.YC( find(grdO.HFacW(1:6*nc,1:nc,1) == 0) )) %% U-points
16	%min(grdO.YG( find(grdO.HFacS(1:6*nc,1:nc,1) == 0) )) %% V-points
17
18	%%%% Small-Basin Masks
19
20	mskSB_G=ones(6nc,nc); mskSB_G(grdO.XG(1:6nc,1:nc) < 90 )=0;
21	mskSB_C=ones(6nc,nc); mskSB_C(grdO.XC(1:6nc,1:nc) < 90 )=0;
22
23	maskW_SB = maskW .* repmat(mskSB_G,[1 1 nr]);
24	maskS_SB = maskS .* repmat(mskSB_C,[1 1 nr]);
25	maskC_SB = maskC .* repmat(mskSB_C,[1 1 nr]);
26
27	msk_SB.maskW=maskW_SB;
28	msk_SB.maskS=maskS_SB;
29	msk_SB.maskC=maskC_SB;
30
31	%%%% Large-Basin Masks
32
33	mskLB_G=ones(6nc,nc); mskLB_G(grdO.XG(1:6nc,1:nc) > 90 )=0;
34	mskLB_C=ones(6nc,nc); mskLB_C(grdO.XC(1:6nc,1:nc) > 90 )=0;
35
36	maskW_LB = maskW .* repmat(mskLB_G,[1 1 nr]);
37	maskS_LB = maskS .* repmat(mskLB_C,[1 1 nr]);
38	maskC_LB = maskC .* repmat(mskLB_C,[1 1 nr]);
39
40	msk_LB.maskW=maskW_LB;
41	msk_LB.maskS=maskS_LB;
42	msk_LB.maskC=maskC_LB;
43
44	return