docs: def get_kinds(self, kind, R):
          """
          Yield 'rlz-000', 'rlz-001', ...', 'mean', 'quantile-0.1', ...
          """
          stats = self.hazard_stats()
          if kind == 'stats':
              yield from stats
              return
          elif kind == 'rlzs':
              for r in range(R):
                  yield 'rlz-%d' % r
              return
          elif kind:
              yield kind
              return
          # default: yield stats (and realizations if required)
          if R > 1 and self.individual_curves or not stats:
              for r in range(R):
                  yield 'rlz-%03d' % r
          yield from stats