docs: def _get_a1(bbar, dbar, slip_moment, mmax):
          """
          Returns the A1 term (I.4 of Table 2 in Anderson & Luco)
          """
          return ((dbar - bbar) / dbar) * (slip_moment / _scale_moment(mmax))