docs: def create_anisophaplot(plotman, x, y, z, alpha, options):
      '''Plot the data of the tomodir in one overview plot.
      '''
      sizex, sizez = getfigsize(plotman)
      # create figure
      f, ax = plt.subplots(2, 3, figsize=(3 * sizex, 2 * sizez))
      if options.title is not None:
          plt.suptitle(options.title, fontsize=18)
          plt.subplots_adjust(wspace=1, top=0.8)
      # plot phase
      cidx = plotman.parman.add_data(x)
      cidy = plotman.parman.add_data(y)
      cidz = plotman.parman.add_data(z)
      cidxy = plotman.parman.add_data(np.subtract(x, y))
      cidyz = plotman.parman.add_data(np.subtract(y, z))
      cidzx = plotman.parman.add_data(np.subtract(z, x))
      plot_pha(cidx, ax[0, 0], plotman, 'x', alpha,
               options.pha_vmin, options.pha_vmax,
               options.xmin, options.xmax, options.zmin, options.zmax,
               options.unit, options.pha_cbtiks, options.no_elecs,
               )
      plot_pha(cidy, ax[0, 1], plotman, 'y', alpha,
               options.pha_vmin, options.pha_vmax,
               options.xmin, options.xmax, options.zmin, options.zmax,
               options.unit, options.pha_cbtiks, options.no_elecs,
               )
      plot_pha(cidz, ax[0, 2], plotman, 'z', alpha,
               options.pha_vmin, options.pha_vmax,
               options.xmin, options.xmax, options.zmin, options.zmax,
               options.unit, options.pha_cbtiks, options.no_elecs,
               )
      plot_ratio(cidxy, ax[1, 0], plotman, 'x-y', alpha,
                 options.rat_vmin, options.rat_vmax,
                 options.xmin, options.xmax, options.zmin, options.zmax,
                 options.unit, options.mag_cbtiks, options.no_elecs,
                 )
      plot_ratio(cidyz, ax[1, 1], plotman, 'y-z', alpha,
                 options.rat_vmin, options.rat_vmax,
                 options.xmin, options.xmax, options.zmin, options.zmax,
                 options.unit, options.mag_cbtiks, options.no_elecs,
                 )
      plot_ratio(cidzx, ax[1, 2], plotman, 'z-x', alpha,
                 options.rat_vmin, options.rat_vmax,
                 options.xmin, options.xmax, options.zmin, options.zmax,
                 options.unit, options.mag_cbtiks, options.no_elecs,
                 )
      f.tight_layout()
      f.savefig('pha_aniso.png', dpi=300)
      return f, ax