docs: def get_fault_term(self, rake):
          """
          Returns coefficient for faulting style (pg 156)
          """
          rake = rake + 360 if rake < 0 else rake
  
          if (rake >= 45) & (rake <= 135):
              f = 1.
          elif (rake >= 225) & (rake <= 315):
              f = 0.5
          else:
              f = 0.
          return f