docs: def _get_magnitude_term(self, C, mag):
          """
          Returns the magnitude scaling term defined in equation 2
          """
          f_mag = C["c0"] + C["c1"] * mag
          if (mag > 4.5) and (mag <= 5.5):
              return f_mag + (C["c2"] * (mag - 4.5))
          elif (mag > 5.5) and (mag <= 6.5):
              return f_mag + (C["c2"] * (mag - 4.5)) + (C["c3"] * (mag - 5.5))
          elif mag > 6.5:
              return f_mag + (C["c2"] * (mag - 4.5)) + (C["c3"] * (mag - 5.5)) +\
                  (C["c4"] * (mag - 6.5))
          else:
              return f_mag