code: def _get_style_of_faulting_term(self, C, rup):
          
          SS, NS, RS = 0.0, 0.0, 0.0
          if np.abs(rup.rake) <= 30.0 or (180.0 - np.abs(rup.rake)) <= 30.0:
              # strike-slip
              SS = 1.0
          elif rup.rake > 30.0 and rup.rake < 150.0:
              # reverse
              RS = 1.0
          else:
              # normal
              NS = 1.0
          return (C["sofN"] * NS) + (C["sofR"] * RS) + (C["sofS"] * SS)
docs: Returns the style-of-faulting term.
          Fault type (Strike-slip, Normal, Thrust/reverse) is
          derived from rake angle.
          Rakes angles within 30 of horizontal are strike-slip,
          angles from 30 to 150 are reverse, and angles from
          -30 to -150 are normal.
          Note that the 'Unspecified' case is not considered in this class
          as rake is required as an input variable