INTEGER I,J,kUp,kUpC,kTop
kUpC = k
kUpC = kUp
rdrckp1=recip_drC(kUpC)
& - 2.*KappaRU(i,j,kUpC)*rdrckp1*viscFac
& - 2.*KappaRV(i,j,kUpC)*rdrckp1*viscFac