code: def _get_sof_terms(self, C, rake):
          
          if rake >= 45.0 and rake <= 135.0:
              # Reverse faulting
              return C["b0R"], C["b1R"]
          elif rake <= -45. and rake >= -135.0:
              # Normal faulting
              return C["b0N"], C["b1N"]
          else:
              # Strike slip
              return C["b0SS"], C["b1SS"]
docs: Returns the style-of-faulting scaling parameters