code: def _get_a2_value(bbar, dbar, slip, beta, mmax):
          
          return ((dbar - bbar) / bbar) * (slip / beta) *\
              np.exp(-(dbar / 2.) * mmax)
docs: Returns the A2 value defined in II.8 (Table 3)