docs: def plot_np(fignum, indata, s, units):
      """
      makes plot of de(re)magnetization data for Thellier-Thellier type experiment
  
      Parameters
      __________
      fignum : matplotlib figure number
      indata :  araiblock from, e.g., pmag.sortarai()
      s : specimen name
      units : [K, J, ""] (kelvin, joules, unknown)
  
      Effect
      _______
      Makes a plot
      """
      global globals
      first_Z, first_I, ptrm_check, ptrm_tail = indata[0], indata[1], indata[2], indata[3]
      plt.figure(num=fignum)
      plt.clf()
      if not isServer:
          plt.figtext(.02, .01, version_num)
      X, Y, recnum = [], [], 0
  #
      for rec in first_Z:
          if units == "K":
              if rec[0] != 0:
                  X.append(rec[0] - 273.)
              else:
                  X.append(rec[0])
          if (units == "J") or (not units) or (units == "T"):
              X.append(rec[0])
          Y.append(old_div(rec[3], first_Z[0][3]))
          delta = .02 * Y[0]
          if recnum % 2 == 0:
              plt.text(X[-1] - delta, Y[-1] + delta,
                       (' ' + str(recnum)), fontsize=9)
          recnum += 1
      plt.plot(X, Y)
      plt.scatter(X, Y, marker='o', color='b')
      X, Y = [], []
      for rec in first_I:
          if units == "K":
              if rec[0] != 0:
                  X.append(rec[0] - 273)
              else:
                  X.append(rec[0])
          if (units == "J") or (not units) or (units == "T"):
              X.append(rec[0])
          Y.append(old_div(rec[3], first_Z[0][3]))
      if globals != 0:
          globals.DIlist = X
          globals.DIlisty = Y
      plt.plot(X, Y)
      plt.scatter(X, Y, marker='s', color='r')
      plt.ylabel("Circles: NRM; Squares: pTRM")
      if units == "K":
          plt.xlabel("Temperature (C)")
      elif units == "J":
          plt.xlabel("Microwave Energy (J)")
      else:
          plt.xlabel("")
      title = s + ": NRM = " + '%9.2e' % (first_Z[0][3])
      plt.title(title)
      plt.axhline(y=0, xmin=0, xmax=1, color='k')
      plt.axvline(x=0, ymin=0, ymax=1, color='k')