% Create Masked Regions
% Hannah Ross
% 10/23/2019

% Create masked regions for calculating wake extent from TKE.

function [ymask1,ymask2,zmask1,zmask2] = maskRegions(posGrid)

% UW plane 1 (X/D = 0.75)
ymask1.MongoUW{1,1} = posGrid.MongoUW{1,1}(:,:,2) < -31 | posGrid.MongoUW{1,1}(:,:,2) > -23;
ymask2.MongoUW{1,1} = posGrid.MongoUW{1,1}(:,:,2) < 27 | posGrid.MongoUW{1,1}(:,:,2) > 28;

zmask1.MongoUW{1,1} = posGrid.MongoUW{1,1}(:,:,3) < -15 | posGrid.MongoUW{1,1}(:,:,3) > -10;
zmask2.MongoUW{1,1} = posGrid.MongoUW{1,1}(:,:,3) < 11 | posGrid.MongoUW{1,1}(:,:,3) > 15.1;

% UW plane 2 (X/D = 1.25)
ymask1.MongoUW{2,1} = posGrid.MongoUW{2,1}(:,:,2) < -31 | posGrid.MongoUW{2,1}(:,:,2) > -22;
ymask2.MongoUW{2,1} = posGrid.MongoUW{2,1}(:,:,2) < 25 | posGrid.MongoUW{2,1}(:,:,2) > 31;

zmask1.MongoUW{2,1} = posGrid.MongoUW{2,1}(:,:,3) < -18 | posGrid.MongoUW{2,1}(:,:,3) > -8;
zmask2.MongoUW{2,1} = posGrid.MongoUW{2,1}(:,:,3) < 11 | posGrid.MongoUW{2,1}(:,:,3) > 16;

% UW plane 3 (X/D = 1.75)
ymask1.MongoUW{3,1} = posGrid.MongoUW{3,1}(:,:,2) < -30 | posGrid.MongoUW{3,1}(:,:,2) > -25.5;
ymask2.MongoUW{3,1} = posGrid.MongoUW{3,1}(:,:,2) < 25.5 | posGrid.MongoUW{3,1}(:,:,2) > 31;

zmask1.MongoUW{3,1} = posGrid.MongoUW{3,1}(:,:,3) < -16 | posGrid.MongoUW{3,1}(:,:,3) > -4;
zmask2.MongoUW{3,1} = posGrid.MongoUW{3,1}(:,:,3) < 8 | posGrid.MongoUW{3,1}(:,:,3) > 16;

% UW plane 4 (X/D = 2.25)
ymask1.MongoUW{4,1} = posGrid.MongoUW{4,1}(:,:,2) < -31 | posGrid.MongoUW{4,1}(:,:,2) > -25.5;
ymask2.MongoUW{4,1} = posGrid.MongoUW{4,1}(:,:,2) < 26 | posGrid.MongoUW{4,1}(:,:,2) > 27;

zmask1.MongoUW{4,1} = posGrid.MongoUW{4,1}(:,:,3) < -15 | posGrid.MongoUW{4,1}(:,:,3) > -2;
zmask2.MongoUW{4,1} = posGrid.MongoUW{4,1}(:,:,3) < 4 | posGrid.MongoUW{4,1}(:,:,3) > 13.4;

% BMSC plane 1 (X/D = 0.75)
ymask1.MongoBMSC{1,1} = posGrid.MongoBMSC{1,1}(:,:,2) < -33 | posGrid.MongoBMSC{1,1}(:,:,2) > -25;
ymask2.MongoBMSC{1,1} = posGrid.MongoBMSC{1,1}(:,:,2) < 28 | posGrid.MongoBMSC{1,1}(:,:,2) > 33;

zmask1.MongoBMSC{1,1} = posGrid.MongoBMSC{1,1}(:,:,3) < -16.4 | posGrid.MongoBMSC{1,1}(:,:,3) > -10;
zmask2.MongoBMSC{1,1} = posGrid.MongoBMSC{1,1}(:,:,3) < 10 | posGrid.MongoBMSC{1,1}(:,:,3) > 19;

% BMSC plane 2 (X/D = 1.25)
ymask1.MongoBMSC{2,1} = posGrid.MongoBMSC{2,1}(:,:,2) < -36 | posGrid.MongoBMSC{2,1}(:,:,2) > -29;
ymask2.MongoBMSC{2,1} = posGrid.MongoBMSC{2,1}(:,:,2) < 28 | posGrid.MongoBMSC{2,1}(:,:,2) > 35;

zmask1.MongoBMSC{2,1} = posGrid.MongoBMSC{2,1}(:,:,3) < -18 | posGrid.MongoBMSC{2,1}(:,:,3) > -8;
zmask2.MongoBMSC{2,1} = posGrid.MongoBMSC{2,1}(:,:,3) < 8 | posGrid.MongoBMSC{2,1}(:,:,3) > 20;

% BMSC plane 3 (X/D = 1.75)
ymask1.MongoBMSC{3,1} = posGrid.MongoBMSC{3,1}(:,:,2) < -31.4 | posGrid.MongoBMSC{3,1}(:,:,2) > -25;
ymask2.MongoBMSC{3,1} = posGrid.MongoBMSC{3,1}(:,:,2) < 26 | posGrid.MongoBMSC{3,1}(:,:,2) > 32;

zmask1.MongoBMSC{3,1} = posGrid.MongoBMSC{3,1}(:,:,3) < -14.5 | posGrid.MongoBMSC{3,1}(:,:,3) > -6.9;
zmask2.MongoBMSC{3,1} = posGrid.MongoBMSC{3,1}(:,:,3) < 7.8 | posGrid.MongoBMSC{3,1}(:,:,3) > 17;

% BMSC plane 4 (X/D = 2.25)
ymask1.MongoBMSC{4,1} = posGrid.MongoBMSC{4,1}(:,:,2) < -36 | posGrid.MongoBMSC{4,1}(:,:,2) > -26;
ymask2.MongoBMSC{4,1} = posGrid.MongoBMSC{4,1}(:,:,2) < 27 | posGrid.MongoBMSC{4,1}(:,:,2) > 36;

zmask1.MongoBMSC{4,1} = posGrid.MongoBMSC{4,1}(:,:,3) < -18 | posGrid.MongoBMSC{4,1}(:,:,3) > -3;
zmask2.MongoBMSC{4,1} = posGrid.MongoBMSC{4,1}(:,:,3) < 5.5 | posGrid.MongoBMSC{4,1}(:,:,3) > 17;

for i = 1:size(ymask1.MongoUW,1)
    ymask1.MongoUW{i,1} = ~ymask1.MongoUW{i,1};
    zmask1.MongoUW{i,1} = ~zmask1.MongoUW{i,1};
    ymask1.MongoBMSC{i,1} = ~ymask1.MongoBMSC{i,1};
    zmask1.MongoBMSC{i,1} = ~zmask1.MongoBMSC{i,1};
    
    ymask2.MongoUW{i,1} = ~ymask2.MongoUW{i,1};
    zmask2.MongoUW{i,1} = ~zmask2.MongoUW{i,1};
    ymask2.MongoBMSC{i,1} = ~ymask2.MongoBMSC{i,1};
    zmask2.MongoBMSC{i,1} = ~zmask2.MongoBMSC{i,1};
end

end