_RL rdrckp1
rdrckp1=recip_drC(kUpC)
& * kappaRU(i,j,kUpC)*rdrckp1*viscFac
& * kappaRV(i,j,kUpC)*rdrckp1*viscFac