docs: def plot_slnp(fignum, SiteRec, datablock, key):
      """
      plots lines and planes on a great  circle with alpha 95 and mean
      deprecated (used in pmagplotlib)
      """
  # make the stereonet
      plt.figure(num=fignum)
      plot_net(fignum)
      s = SiteRec['er_site_name']
  #
  #   plot on the data
  #
      coord = SiteRec['site_tilt_correction']
      title = ''
      if coord == '-1':
          title = s + ": specimen coordinates"
      if coord == '0':
          title = s + ": geographic coordinates"
      if coord == '100':
          title = s + ": tilt corrected coordinates"
      DIblock, GCblock = [], []
      for plotrec in datablock:
          if plotrec[key + '_direction_type'] == 'p':  # direction is pole to plane
              GCblock.append(
                  (float(plotrec[key + "_dec"]), float(plotrec[key + "_inc"])))
          else:  # assume direction is a directed line
              DIblock.append(
                  (float(plotrec[key + "_dec"]), float(plotrec[key + "_inc"])))
      if len(DIblock) > 0:
          plot_di(fignum, DIblock)  # plot directed lines
      if len(GCblock) > 0:
          for pole in GCblock:
              plot_circ(fignum, pole, 90., 'g')  # plot directed lines
  #
  # put on the mean direction
  #
      x, y = [], []
      XY = pmag.dimap(float(SiteRec["site_dec"]), float(SiteRec["site_inc"]))
      x.append(XY[0])
      y.append(XY[1])
      plt.scatter(x, y, marker='d', s=80, c='g')
      plt.title(title)
  #
  # get the alpha95
  #
      Xcirc, Ycirc = [], []
      Da95, Ia95 = pmag.circ(float(SiteRec["site_dec"]), float(
          SiteRec["site_inc"]), float(SiteRec["site_alpha95"]))
      for k in range(len(Da95)):
          XY = pmag.dimap(Da95[k], Ia95[k])
          Xcirc.append(XY[0])
          Ycirc.append(XY[1])
      plt.plot(Xcirc, Ycirc, 'g')