code: def _compute_faulting_mechanism(self, C, rake, dip):
          
          # flag for reverse faulting
          frv = float((dip > 45) and (22.5 <= rake <= 157.5))
          # flag for thrust faulting
          fth = float((dip <= 45) and (22.5 <= rake <= 157.5))
  
          return C['c10'] * frv + C['c11'] * fth
docs: Compute faulting mechanism term (see eq. 5, page 319).
  
          Reverse faulting is defined as occurring on steep faults (dip > 45)
          and rake in (22.5, 157.5).
  
          Thrust faulting is defined as occurring on shallow dipping faults
          (dip <=45) and rake in (22.5, 157.5)