Request to accommodate for an initial porosity profile instead of a flat (i.e. uniform) initial porosity distribution. The profile is described [here](https://github.com/MindTheGap-ERC/LHeureuxEqs/blob/main/Eqs.html). This will also be implemented in the [Matlab code](https://github.com/MindTheGap-ERC/LMA-Matlab).