def METHOD_NAME(model_horizontal_two_layers):
    import subsurface
    m = model_horizontal_two_layers
    vals = gp.compute_model(m, set_solutions=True)
    unstruct = m.solutions.set_meshes(m.surfaces)
    ts = subsurface.TriSurf(unstruct)
    s  = subsurface.visualization.to_pyvista_mesh(ts)
    subsurface.visualization.pv_plot([s], image_2d=True)
    print(unstruct)