_RL ws(imt), wm(imt)
O wm, ws, myThid )
wm(i) = sign(eins,wm(i))*max(phepsi,abs(wm(i)))
gat1m(i) = visch / hbl(i) / wm(i)
dat1m(i) = -viscp / wm(i) + f1 * visch
blmc(i,ki,1) = hbl(i) * wm(i) * sig * (1. + sig * Gm)
dkm1(i,1) = hbl(i) * wm(i) * sig * (1. + sig * Gm)
_RL wm(imt), ws(imt)
O wm, ws,
wm(i) = (1.-ufrac) * wbm + ufrac * wam
wm(i) = vonk * ustar(i) * u3 / tempVar
ws(i) = wm(i)