docs: def _compute_finite_fault_correction(self, mag):
          """
          Compute finite fault correction term as geometric mean of correction
          terms obtained from Mw values calculated with Johnston 1996 and
          Atkinson and Boore 1987 conversion equations.
  
          Implement equations as in lines 1653 - 1658 in hazgridXnga2.f
          """
          mw_j96 = mblg_to_mw_johnston_96(mag)
          mw_ab87 = mblg_to_mw_atkinson_boore_87(mag)
  
          t1 = np.exp(-1.25 + 0.227 * mw_j96)
          t2 = np.exp(-1.25 + 0.227 * mw_ab87)
  
          return np.sqrt(t1 * t2)