  The computation of the Love numbers for a spherically symmetric
self-gravitating viscoelastic Earth is a classical problem in global
geodynamics. Here we revisit the problem of the numerical evaluation of loading
and tidal Love numbers in the static limit for an incompressible planetary
body, adopting a Laplace inversion scheme based upon the Post-Widder formula as
an alternative to the {traditional viscoelastic normal modes method. We also
consider, whithin the same framework, complex-valued, frequency-dependent Love
numbers that describe the response to a periodic forcing, which are paramount
in the study of the tidal deformation of planets. Furthermore, we numerically
obtain the time-derivatives of Love numbers, suitable for modeling geodetic
signals in response to surface loads variations. A number of examples are
shown, in which time and frequency-dependent Love numbers are evaluated for the
Earth and planets adopting realistic rheological profiles. The numerical
solution scheme is implemented in ALMA${}^3$ (the plAnetary Love nuMbers
cAlculator, version 3), an upgraded open-source Fortran 90 program that
computes the Love numbers for radially layered planetary bodies with a wide
range of rheologies, including transient laws like Andrade or Burgers.
