_RL smooth3D_KzMax
smooth3D_KzMax=drC(k)
smooth3D_KzMax=smooth3D_KzMax*smooth3D_KzMax/smooth3DtotTime/2
if (smooth3D_kappaR(i,j,k,bi,bj).GT.smooth3D_KzMax) then
smooth3D_kappaR(i,j,k,bi,bj)=smooth3D_KzMax