function L = L_MO(ustar,T,w_T_)
% input: ustar [m/s], T [^oC], w_T_ [m/s K] (in this order)
% output L [m]
g = 9.81;
vonK = 0.4;

L = -ustar.^3 ./ (vonK * g ./ (T + 273.15) .* w_T_);