Dear PyLith Developers,
If we want to invert interseismic backslip rate using PyLith, how to generate Green’ functions based on a viscoelastic Earth model.The example in the user manual seems to only be able to generate the elastic Green’s function.
The GreensFnproblem type is setup for static Greeen’s functions only. You can manually generate time-dependent Green’s functions using a TimeDependent problem. You will need to create use a new simulation for each Green’s function with the impulses specified using spatial databases for the prescribed slip.