code: def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
          
          # This is just used for testing purposes
          if len(stddev_types) == 0:
              stddev_types = [StdDev.TOTAL]
          mean, stds = self._get_mean_and_stddevs(sites, rup, dists, imt,
                                                  stddev_types)
          stddevs = [np.ones(len(dists.repi))*get_sigma(imt)]
          delta = self._get_delta(stds, dists)
          mean = mean + stds + delta
          mean = np.squeeze(mean)
          return mean, stddevs
docs: See documentation for method `GroundShakingIntensityModel` in
          :class:~`openquake.hazardlib.gsim.base.GSIM`