_RL viscFac
viscFac=0.
IF (no_slip_bottom) viscFac=1.
& 2.*KappaRV(i,j,kDownC)*rdrckp1*viscFac
& 2.*KappaRU(i,j,kDownC)*rdrckp1*viscFac