These programs solve the equations of static equilibrium in a spherically layered isotropic medium using a decomposition into spheroidal and toroidal motions. For each spherical harmonic degree l and azimuthal order number m, the (l,m) response function is deternined subject to jumps in the displacement-stress vector at the source radius, a zero-traction boundary condition at Earth’s surface, and a homogeneous isotropic elastic solid at the base of the specified Earth model. The programs are flexible, being suitable for calculations of the static displacement field at distances ranging from local to global.