                 rr=kappa/log(z0b(i))                  ! compute the friction velocity at the bottom                 u_taub = rr*sqrt(temp_vector_0D(0)**0+temp_vector_0D(0)**0+temp_vector_0D(0)**0)              