code: def plot_arai_zij(ZED, araiblock, zijdblock, s, units):
      
      angle = zijdblock[0][1]
      norm = 1
      if units == "U":
          norm = 0
      plot_arai(ZED['arai'], araiblock, s, units)
      plot_teq(ZED['eqarea'], araiblock, s, "")
      plot_zij(ZED['zijd'], zijdblock, angle, s, norm)
      plot_np(ZED['deremag'], araiblock, s, units)
      return ZED
docs: calls the four plotting programs for Thellier-Thellier experiments
  
      Parameters
      __________
      ZED : dictionary with plotting figure keys:
          deremag : figure for de (re) magnezation plots
          arai : figure for the Arai diagram
          eqarea : equal area projection of data, color coded by
              red circles: ZI steps
              blue squares: IZ steps
              yellow triangles : pTRM steps
          zijd : Zijderveld diagram color coded by ZI, IZ steps
          deremag : demagnetization and remagnetization versus temperature
      araiblock : nested list of required data from Arai plots
      zijdblock : nested list of required data for Zijderveld plots
      s : specimen name
      units : units for the arai and zijderveld plots
  
      Effects
      ________
      Makes four plots from the data by calling
      plot_arai : Arai plots
      plot_teq : equal area projection for Thellier data
      plotZ : Zijderveld diagram
      plot_np : de (re) magnetization diagram