& ( (wBarYp*vZp+wBarYm*vZm)
& +(ABS(wBarYp)*vZp-ABS(wBarYm)*vZm)
vShearTerm(I,J)=-0.5*(wBarYp*vZp+wBarYm*vZm)
_RL mask_Kp1,mask_Km1,wBarYm,wBarYp
wBarYm=0.5*(