import scipy.io as sio
import numpy as np
import sys
import matplotlib.pyplot as plt
plt.rcParams.update({'font.size': 22})


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


nt=len(w_k2)
nk=6

# put all files in an array
w=np.zeros((nt,nk))
w[:,0] =w_k2
w[:,1] =w_k6
w[:,2] =w_k10
w[:,3] =w_k14
w[:,4] =w_k18
w[:,5] =w_k22

# compute the two=point correlation at node [0]