This a simple stokes solver to calculate normal stress at the surface and topography from the flow in the mantle. It is planned to be coupled with LitMod2D_2.0. It reads post_processing.out from LitMod2D_2.0 and takes temperature and density distribution to solve for the flow.
Note: still under development stage.