In cases that use surface momentum flux parametrization = C|U|u, distribute momentum flux using SMG scheme from libmpdata++ (set cdrag).
In cases that use other parametreizations (e.g. DYCOMS), keep the explicit momentum fluxes (with exponential vertical redistribution) or, better, find a way to use SGS to do this.