docs: def _get_a1_value(bbar, dbar, slip, beta, mmax):
          """
          Returns the A1 value defined in I.9 (Table 2)
          """
          return ((dbar - bbar) / dbar) * (slip / beta) *\
              np.exp(-(dbar / 2.) * mmax)