function [datmat] = array2matrix(dat)

[nlin,ncol,nlayer] = size(dat);


datNaN = NaN(nlin,nlayer*ncol);
datNaN(:,1:3:end) = dat(:,:,1);
datNaN(:,2:3:end) = dat(:,:,2);
datNaN(:,3:3:end) = dat(:,:,3);

% datmat = [[1:nlin]',[(0:nlin-1)/freq]',datNaN];
datmat = datNaN;
end