%omkar

nu_ft2ps = 1.764e-4;
rho_slugpft3 = 0.0020451;
c_ft = 10;
xtr_ft = .2*c_ft;

ft_2_in = 12;
Vinf_fps = 198;
q = 0.5 * rho_slugpft3 * (Vinf_fps^2);
Re_xtr = Vinf_fps * xtr_ft / nu_ft2ps;
delta_xtr_ft=( 5.2 * xtr_ft )/(Re_xtr)^0.5;
delta_l_ft = ((delta_xtr_ft*((Vinf_fps)^0.2))/(0.37*((nu_ft2ps)^0.2)))^(1/0.8);
lt_ft = c_ft - xtr_ft + delta_l_ft;
Re_lt = Vinf_fps * lt_ft / nu_ft2ps;
delta_ft = 0.37 * lt_ft /(Re_lt)^0.2


%         for k = 1:100
%             if (k <= 100*lam)
%                 Re_x = v_fps * x(k) / dynvisc;
%                 delta_x=( 5.2 * x(k) )/(Re_x)^0.5;
%                 wake(j,k)= delta_x;
%             else %if turbulent
%                 x_tr=x(100*lam);
%                 Re_xtr = v_fps * x_tr / dynvisc;
%                 delta_xtr=( 5.2 * x_tr )/(Re_xtr)^0.5;
%                 delta_l = ((delta_xtr*((v_fps)^0.2))/(0.37*((dynvisc)^0.2)))^(1/0.8);
%                 lt = x(k) - x_tr + delta_l;
%                 Re_lt = v_fps * lt / dynvisc;
%                 delta = 0.37 * lt /(Re_lt)^0.2;
%                 wake(j,k) = delta;
%             end
%
%         end