code: def _get_fault_type_dummy_variables(self, rup):
          
          U, SS, NS, RS = 0, 0, 0, 0
          if np.abs(rup.rake) <= 30.0 or (180.0 - np.abs(rup.rake)) <= 30.0:
              # strike-slip
              SS = 1
          elif rup.rake > 30.0 and rup.rake < 150.0:
              # reverse
              RS = 1
          else:
              # normal
              NS = 1
          return U, SS, NS, RS
docs: 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,
          because rake is always given.