

dx=0.007769
dz= 0.00234375

% load grid
viscos=1.0684e-06
nj=162
load exam_data.txt

data=exam_data;

yc=data(:,1);      % coordinate of the control volume faces
u=data(:,2);       % velocity, v_1
k=data(:,3);       % modeled turbulent kinetic energy
omega=data(:,4);   % omega
uv=data(:,5);      % resolved shear stress


for j=2:nj-1
   y(j)=0.5*(yc(j)+yc(j-1));  % cell center
end

y(nj)=yc(nj-1);




