function epsilon=dissipation_calculation(u_star,z)

k=0.41;

epsilon = (u_star.^3)./(k*z);
