code: def _get_magnitude_scaling_term(self, C, mag):
          
          dmag = mag - self.CONSTS["Mh"]
          if mag < self.CONSTS["Mh"]:
              return C["e1"] + (C["b1"] * dmag) + (C["b2"] * (dmag ** 2.0))
          else:
              return C["e1"] + (C["b3"] * dmag)
docs: Returns the magnitude scaling term of the GMPE described in
          equation 3