def METHOD_NAME(self, ctx: np.recarray, imts, mean, sig, tau, phi):
    """
    See :meth:`superclass method
    <.base.GroundShakingIntensityModel.compute>`
    for spec of input and result values.
    """
    for m, imt in enumerate(imts):
        C = self.COEFFS[imt]
        mag = ctx.mag - 6
        d = np.sqrt(ctx.rjb ** 2 + C['c7'] ** 2)
        mean[m] += C['c1'] + C['c2'] * mag + C['c3'] * mag ** 2 + C['c6']
        idx = d <= 100.
        mean[m, idx] += C['c5'] * np.log10(d[idx])
        idx = d > 100.
        mean[m, idx] += (C['c5'] * np.log10(100.) -
                         np.log10(d[idx] / 100.) +
                         C['c4'] * (d[idx] - 100.))
        # convert from log10 to ln and from cm/s**2 to g
        mean[m] = np.log((10.0 ** (mean[m] - 2.0)) / g)
        sig[m] = C['sigma']