Generate stress in 2D subduction

[baagaard]

For one-point quadrature the point is in the center of the triangle (note that the basis functions are all equal at the quadrature point).

Yes, you can add functions to the script as necessary.