code: def get_magnitude_scaling_term(self, C, rup):
          
          if rup.mag <= self.CONSTANTS["m_c"]:
              return C["ccr"] * rup.mag
          else:
              return (C["ccr"] * self.CONSTANTS["m_c"]) +\
                  (C["dcr"] * (rup.mag - self.CONSTANTS["m_c"]))
docs: Returns the magnitude scaling term in equations 1 and 2