code: def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
          
          assert all(stddev_type in self.DEFINED_FOR_STANDARD_DEVIATION_TYPES
                     for stddev_type in stddev_types)
  
          C = self.COEFFS[imt]
          mean = self._get_mean(
              C, rup.mag, rup.rake, rup.dip, dists.rrup, dists.rjb
          )
          stddevs = self._get_stddevs(C, rup.mag, stddev_types, dists.rrup.size)
  
          return mean, stddevs
docs: See :meth:`superclass method
          <.base.GroundShakingIntensityModel.get_mean_and_stddevs>`
          for spec of input and result values.