Matlab files used to produce the figures in the publication: On the Derivation of Closed-Form Expressions for Displacements, Strains and Stresses Inside Poroelastic Reservoirs

MATLAB code to recreate figures from the paper:
Cornelissen, P., Meulenbroek, B.J., Jansen, J.D., (2023). On the Derivation of Closed-Form Expressions for Displacements, Strains and Stresses Inside Poroelastic Reservoirs. Submitted to JGR Solid Earth.
The MATLAB scripts are used to compute the analytical solution for the poroelastic strain and stress field for rectangular and triangular reservoirs. The analytical solutions have been derived using inclusion theory (Equations 40-42 and 45-47 of the associated manuscript). The purpose of these scripts is to verify the analytical solutions by computing the strain discontinuity across the inclusion boundary and computing the volumetric strain in and outside the inclusion, and to show an example of an application of inclusion theory to reservoir depletion. These scripts produce Figures 5, 6, and 9 of the associated manuscript.