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


dx=0.007769
dz= 0.00234375

# load grid
viscos=1.0684e-06
nj=162
data= np.loadtxt("exam_data.txt")

yc=data[:,0]      # coordinate of the control volume faces
u=data[:,1]       # velocity, v_1
k=data[:,2]       # modeled turbulent kinetic energy
omega=data[:,3]   # omega
uv=data[:,4]      # resolved shear stress

y=np.zeros(nj)
for j in range (0,nj-1):
   y[j]=0.5*(yc[j]+yc[j-1])   # cell center

y[-1]=yc[-1-1]