docs: def create_singleplots(plotman, cov, mag, pha, pha_fpi, alpha, options):
      '''Plot the data of the tomodir in individual plots.
      '''
      magunit = 'log_rho'
      if not pha == []:
          [real, imag] = calc_complex(mag, pha)
          if not pha_fpi == []:
              [real_fpi, imag_fpi] = calc_complex(mag, pha_fpi)
              if options.cmaglin:
                  mag = np.power(10, mag)
                  magunit = 'rho'
              data = np.column_stack((mag, cov, pha, real, imag,
                                      pha_fpi, real_fpi, imag_fpi))
              titles = ['Magnitude', 'Coverage',
                        'Phase', 'Real Part', 'Imaginary Part',
                        'FPI Phase', 'FPI Real Part', 'FPI Imaginary Part']
              unites = [
                  magunit, 'cov',
                  'phi', 'log_real', 'log_imag',
                  'phi', 'log_real', 'log_imag'
              ]
              vmins = [options.mag_vmin, options.cov_vmin,
                       options.pha_vmin, options.real_vmin, options.imag_vmin,
                       options.pha_vmin, options.real_vmin, options.imag_vmin]
              vmaxs = [options.mag_vmax, options.cov_vmax,
                       options.pha_vmax, options.real_vmax, options.imag_vmax,
                       options.pha_vmax, options.real_vmax, options.imag_vmax]
              cmaps = ['jet', 'GnBu',
                       'jet_r', 'jet_r', 'plasma_r',
                       'plasma', 'jet_r', 'plasma_r']
              saves = ['rho', 'cov',
                       'phi', 'real', 'imag',
                       'fpi_phi', 'fpi_real', 'fpi_imag']
          else:
              if options.cmaglin:
                  mag = np.power(10, mag)
                  magunit = 'rho'
              data = np.column_stack((mag, cov, pha, real, imag))
              titles = ['Magnitude', 'Coverage',
                        'Phase', 'Real Part', 'Imaginary Part']
              unites = [magunit, 'cov',
                        'phi', 'log_real', 'log_imag']
              vmins = [options.mag_vmin, options.cov_vmin,
                       options.pha_vmin, options.real_vmin, options.imag_vmin]
              vmaxs = [options.mag_vmax, options.cov_vmax,
                       options.pha_vmax, options.real_vmax, options.imag_vmax]
              cmaps = ['jet', 'GnBu',
                       'jet_r', 'jet_r', 'plasma_r']
              saves = ['rho', 'cov',
                       'phi', 'real', 'imag']
      else:
          data = np.column_stack((mag, cov))
          titles = ['Magnitude', 'Coverage']
          unites = [magunit, 'cov']
          vmins = [options.mag_vmin, options.cov_vmin]
          vmaxs = [options.mag_vmax, options.cov_vmax]
          cmaps = ['jet', 'GnBu']
          saves = ['rho', 'cov']
      try:
          mod_rho = np.genfromtxt('rho/rho.dat', skip_header=1, usecols=([0]))
          mod_pha = np.genfromtxt('rho/rho.dat', skip_header=1, usecols=([1]))
          data = np.column_stack((data, mod_rho, mod_pha))
          titles.append('Model')
          titles.append('Model')
          unites.append('rho')
          unites.append('phi')
          vmins.append(options.mag_vmin)
          vmins.append(options.pha_vmin)
          vmaxs.append(options.mag_vmax)
          vmaxs.append(options.pha_vmax)
          cmaps.append('jet')
          cmaps.append('plasma')
          saves.append('rhomod')
          saves.append('phamod')
      except:
          pass
      for datum, title, unit, vmin, vmax, cm, save in zip(
              np.transpose(data), titles, unites, vmins, vmaxs, cmaps, saves):
          sizex, sizez = getfigsize(plotman)
          f, ax = plt.subplots(1, figsize=(sizex, sizez))
          cid = plotman.parman.add_data(datum)
          # handle options
          cblabel = units.get_label(unit)
          if options.title is not None:
              title = options.title
          zlabel = 'z [' + options.unit + ']'
          xlabel = 'x [' + options.unit + ']'
          xmin, xmax, zmin, zmax, vmin, vmax = check_minmax(
                  plotman,
                  cid,
                  options.xmin, options.xmax,
                  options.zmin, options.zmax,
                  vmin, vmax
                  )
          # plot
          cmap = mpl_cm.get_cmap(cm)
          fig, ax, cnorm, cmap, cb, scalarMap = plotman.plot_elements_to_ax(
                  cid=cid,
                  cid_alpha=alpha,
                  ax=ax,
                  xmin=xmin,
                  xmax=xmax,
                  zmin=zmin,
                  zmax=zmax,
                  cblabel=cblabel,
                  title=title,
                  zlabel=zlabel,
                  xlabel=xlabel,
                  plot_colorbar=True,
                  cmap_name=cm,
                  over=cmap(1.0),
              under=cmap(0.0),
                  no_elecs=options.no_elecs,
                  cbmin=vmin,
                  cbmax=vmax,
                  )
          f.tight_layout()
          f.savefig(save + '.png', dpi=300)