docs: def plot_teq(fignum, araiblock, s, pars):
      """
      plots directions  of pTRM steps and zero field steps
  
      Parameters
      __________
      fignum : figure number for matplotlib object
      araiblock : nested list of data from pmag.sortarai()
      s : specimen name
      pars : default is "",
          otherwise is dictionary with keys:
          'measurement_step_min' and 'measurement_step_max'
  
      Effects
      _______
      makes the equal area projection with color coded symbols
          red circles: ZI steps
          blue squares: IZ steps
          yellow : pTRM steps
  
  
      """
      first_Z, first_I = araiblock[0], araiblock[1]
  # make the stereonet
      plt.figure(num=fignum)
      plt.clf()
      ZIblock, IZblock, pTblock = [], [], []
      for zrec in first_Z:  # sort out the zerofield steps
          if zrec[4] == 1:  # this is a ZI step
              ZIblock.append([zrec[1], zrec[2]])
          else:
              IZblock.append([zrec[1], zrec[2]])
      plot_net(fignum)
      if pars != "":
          min, max = float(pars["measurement_step_min"]), float(
              pars["measurement_step_max"])
      else:
          min, max = first_I[0][0], first_I[-1][0]
      for irec in first_I:
          if irec[1] != 0 and irec[1] != 0 and irec[0] >= min and irec[0] <= max:
              pTblock.append([irec[1], irec[2]])
      if len(ZIblock) < 1 and len(IZblock) < 1 and len(pTblock) < 1:
          return
      if not isServer:
          plt.figtext(.02, .01, version_num)
  #
  #   put on the directions
  #
      sym = {'lower': ['o', 'r'], 'upper': ['o', 'm']}
      if len(ZIblock) > 0:
          plot_di_sym(fignum, ZIblock, sym)  # plot ZI directions
      sym = {'lower': ['s', 'b'], 'upper': ['s', 'c']}
      if len(IZblock) > 0:
          plot_di_sym(fignum, IZblock, sym)  # plot IZ directions
      sym = {'lower': ['^', 'g'], 'upper': ['^', 'y']}
      if len(pTblock) > 0:
          plot_di_sym(fignum, pTblock, sym)  # plot pTRM directions
      plt.axis("equal")
      plt.text(-1.1, 1.15, s)