docs: def _compute_style_of_faulting_term(self, C, rup):
          """
          Returns the style of faulting factor, depending on the mechanism (rake)
          and top of rupture depth (equations (4) and (5), pages 145 - 146)
          """
          frv, fnm = self._get_fault_type_dummy_variables(rup.rake)
  
          if frv > 0.:
              # Top of rupture depth term only applies to reverse faults
              if rup.ztor < 1.:
                  ffltz = rup.ztor
              else:
                  ffltz = 1.
          else:
              ffltz = 0.
          return (C['c7'] * frv * ffltz) + (C['c8'] * fnm)