code: def plot_di_mean(dec, inc, a95, color='k', marker='o', markersize=20, label='', legend='no'):
      
      DI_dimap = pmag.dimap(dec, inc)
      if inc < 0:
          plt.scatter(DI_dimap[0], DI_dimap[1],
                      edgecolors=color, facecolors='white',
                      marker=marker, s=markersize, label=label)
      if inc >= 0:
          plt.scatter(DI_dimap[0], DI_dimap[1],
                      edgecolors=color, facecolors=color,
                      marker=marker, s=markersize, label=label)
      Xcirc, Ycirc = [], []
      Da95, Ia95 = pmag.circ(dec, inc, a95)
      if legend == 'yes':
          plt.legend(loc=2)
      for k in range(len(Da95)):
          XY = pmag.dimap(Da95[k], Ia95[k])
          Xcirc.append(XY[0])
          Ycirc.append(XY[1])
      plt.plot(Xcirc, Ycirc, c=color)
      plt.tight_layout()
docs: Plot a mean direction (declination, inclination) with alpha_95 ellipse on
      an equal area plot.
  
      Before this function is called, a plot needs to be initialized with code
      that looks something like:
      >fignum = 1
      >plt.figure(num=fignum,figsize=(10,10),dpi=160)
      >ipmag.plot_net(fignum)
  
      Required Parameters
      -----------
      dec : declination of mean being plotted
      inc : inclination of mean being plotted
      a95 : a95 confidence ellipse of mean being plotted
  
      Optional Parameters (defaults are used if not specified)
      -----------
      color : the default color is black. Other colors can be chosen (e.g. 'r').
      marker : the default is a circle. Other symbols can be chosen (e.g. 's').
      markersize : the default is 20. Other sizes can be chosen.
      label : the default is no label. Labels can be assigned.
      legend : the default is no legend ('no'). Putting 'yes' will plot a legend.