docs: def get_magnitude_scaling(self, C, mag):
          """
          Returns the magnitude scaling term
          """
          d_m = mag - self.CONSTANTS["Mh"]
          if mag < self.CONSTANTS["Mh"]:
              return C["e1"] + C["b1"] * d_m + C["b2"] * (d_m ** 2.0)
          else:
              return C["e1"] + C["b3"] * d_m