Can lammps calculate the shear modulus under the hydrostatic pressure?

Your fix npt command appears to have no options included which
allow the tilt params to vary. Hence they will remain 0,0.

Steve