docs: def view_hmap(token, dstore):
      """
      Display the highest 20 points of the mean hazard map. Called as
      $ oq show hmap:0.1  # 10% PoE
      """
      try:
          poe = valid.probability(token.split(':')[1])
      except IndexError:
          poe = 0.1
      mean = dict(extract(dstore, 'hcurves?kind=mean'))['mean']
      oq = dstore['oqparam']
      hmap = calc.make_hmap_array(mean, oq.imtls, [poe], len(mean))
      dt = numpy.dtype([('sid', U32)] + [(imt, F32) for imt in oq.imtls])
      array = numpy.zeros(len(hmap), dt)
      for i, vals in enumerate(hmap):
          array[i] = (i, ) + tuple(vals)
      array.sort(order=list(oq.imtls)[0])
      return rst_table(array[:20])