2D single-phase flow case in MRST

[Nurzhaugan Omiraliyeva]

Hi, I am trying to reproduce a simple 2D single-phase flow case in MRST and compare it with COMSOL.

In COMSOL, the pressure distribution appears as a smooth diagonal gradient across the domain. However, in MRST I get a different result — the pressure distribution looks inconsistent, and values near the corners seem incorrect.

I am using hydrostatic pressure boundary conditions on left and right boundaries.  

Could someone help me understand what might be wrong in my setup?