code: def get_magnitude_scaling_term(self, C, rup):
          
          if rup.ztor > 25.0:
              # Deep interface events
              c_int = C["cint"]
          else:
              c_int = C["cintS"]
  
          if rup.mag <= self.CONSTANTS["m_c"]:
              return c_int * rup.mag
          else:
              return (c_int * self.CONSTANTS["m_c"]) +\
                  (C["dint"] * (rup.mag - self.CONSTANTS["m_c"]))
docs: Returns magnitude scaling term, which is dependent on top of rupture
          depth - as described in equations 1 and 2