code: def get_magnitude_scaling_term(self, C, rup):
          
          m_c = self.CONSTANTS["m_c"]
          if rup.mag <= m_c:
              return C["cSL"] * rup.mag +\
                  C["cSL2"] * ((rup.mag - self.CONSTANTS["m_sc"]) ** 2.)
          else:
              return C["cSL"] * m_c +\
                 C["cSL2"] * ((m_c - self.CONSTANTS["m_sc"]) ** 2.) +\
                 C["dSL"] * (rup.mag - m_c)
docs: Returns the magnitude scaling defined in equation 1