& ( (wBarXp*uZp+wBarXm*uZm)
& +(ABS(wBarXp)*uZp-ABS(wBarXm)*uZm)
uShearTerm(I,J)=-0.5*(wBarXp*uZp+wBarXm*uZm)
_RL mask_Kp1,mask_Km1,wBarXm,wBarXp
wBarXp=0.5*(