_RL rdrckp1,maskUp
rdrckp1=recip_drC(kUpC)
& - 2.*KappaRU(i,j,kUpC)*rdrckp1*viscFac
& - 2.*KappaRV(i,j,kUpC)*rdrckp1*viscFac