
% time history of v_3 velocity at spanwise cells 2, 6, ...
w_k2 = load("w_k2_time.dat");
w_k6 = load("w_k6_time.dat");
w_k10 =load("w_k10_time.dat");
w_k14 =load("w_k14_time.dat");
w_k18 =load("w_k18_time.dat");
w_k22 =load("w_k22_time.dat");


nt=length(w_k2)
nk=6

% put all files in an array
w=zeros(nt,nk);
w(:,1) =w_k2;
w(:,2) =w_k6;
w(:,3) =w_k10;
w(:,4) =w_k14;
w(:,5) =w_k18;
w(:,6) =w_k22;

% compute the two=point correlation at node 1

