Today we need to run an initialization model and save sig0.
We must be able to assign stress to each material following a law.
The gravity part (sigzz) must be derived from the density
Horizontal stresses must be derived by a multiplier or a difference from that, or from a fracture gradient