the important parameter to change is the relaxation factor for p_rgh - if you put it to 1.0 your simulation will run. We’ve seen many users running into troubles with this parameter and we’re planning on some changes to make this less error prone.
In addition, I’d recommend some more changes for speed + efficiency - change the linear system solvers to the following:
Last but not least, you’re running a steady state analysis with 4m cells but you want to write out 5 steps, which will give you quite some result data, where most of it is not really interesting, is it? So I’d recommend only writing out the last time step but running more iterations. I ran your 200 iterations, but the solution is not yet converged - so 1000 seems a good next try.