using MOP for Pzz but IK for Pxx and Pyy
implemented Schofield and Henderson scheme for local pressure