docs: def _compute_magnitude_term(self, C, dc1, mag):
          """
          Computes the magnitude scaling term given by equation (2)
          """
          base = C['theta1'] + (self.CONSTS['theta4'] * dc1)
          dmag = self.CONSTS["C1"] + dc1
          if mag > dmag:
              f_mag = (self.CONSTS['theta5'] * (mag - dmag)) +\
                  C['theta13'] * ((10. - mag) ** 2.)
  
          else:
              f_mag = (self.CONSTS['theta4'] * (mag - dmag)) +\
                  C['theta13'] * ((10. - mag) ** 2.)
  
          return base + f_mag