ToDos:
OGSModel.from_feflow()
run_model
sim_result = model.run_model(logfile=f"{out_dir}/out.txt", args=f"-o {out_dir} -m .")
ms = ot.MeshSeries(sim_result)
ogs/ogs!5217 (diffs)