A MATLAB function to compute the elastic response to a disk load
The diskload MATLAB function computes the response to a uniform surface pressure load imposed in a disc of a given angular radius and height. The elastic response is found at one or more points located on the surface of the Earth at specified angular distance(s) from the center of the disc load.
The elastic response is computed using user-supplied elastic loading Love numbers (h,k,l) generated using a specific elastic structure model for the Earth.
This code is described in the following publication:
Bevis, M., Melini, D., Spada, G., 2016. On computing the geoelastic response to a disk load, Geophys. J. Int., 205 (3), 1,804-1,812, doi:10.1093/gji/ggw115.