function [data1, data2, data3, NewVar1, NewVar2, NewVar3] = extr_by_varA2(data, AllVar, column)
% divides the data in a 2d matrix (such as an RMS array) into 3 separate 
% datastacks, each corresponding to three variable values. 
% Also divides the variable matrix into 3. 
data1 = data;
data2 = data;
data3 = data;
NewVar1 = AllVar;
NewVar2 = AllVar;
NewVar3 = AllVar;
for i = length(AllVar(:,column)):-1:1
    if (NewVar1(i,column) == 2 | NewVar1(i,column) == 3)
        NewVar1(i,:) = [];
        data1(:,i) = [];
    end
end
for i = length(AllVar(:,column)):-1:1
    if (NewVar2(i,column) == 1 | NewVar2(i,column) == 3)
        NewVar2(i,:) = [];
        data2(:,i) = [];
    end
end
for i = length(AllVar(:,column)):-1:1
    if (NewVar3(i,column) == 1 | NewVar3(i,column) == 2)
        NewVar3(i,:) = [];
        data3(:,i) = [];
    end
end
end