docs: def plot_eqarea_mean(self, meanpars, fig):
          """
          Given a dictionary of parameters from pmag.dofisher, pmag.dolnp, or
          pmag.dobingham (meanpars) plots parameters to fig
          """
          mpars_to_plot = []
          if meanpars == {}:
              return
          if meanpars['calculation_type'] == 'Fisher by polarity':
              for mode in list(meanpars.keys()):
                  if type(meanpars[mode]) == dict and meanpars[mode] != {}:
                      mpars_to_plot.append(meanpars[mode])
          else:
              mpars_to_plot.append(meanpars)
          ymin, ymax = fig.get_ylim()
          xmin, xmax = fig.get_xlim()
          if 'color' in meanpars:
              color = meanpars['color']
          else:
              color = 'black'
          size, alpha = 30, 1.
          # put on the mean direction
          for mpars in mpars_to_plot:
              XYM = pmag.dimap(float(mpars["dec"]), float(mpars["inc"]))
              if float(mpars["inc"]) > 0:
                  FC = color
                  EC = 'black'
              else:
                  FC = 'white'
                  EC = color
              self.displayed_means.append(fig.scatter([XYM[0]], [
                                          XYM[1]], marker='o', edgecolor=EC, facecolor=FC, s=size, lw=1, clip_on=False, alpha=alpha))
  
              if "alpha95" in list(mpars.keys()):
                  # get the alpha95
                  Xcirc, Ycirc = [], []
                  Da95, Ia95 = pmag.circ(float(mpars["dec"]), float(
                      mpars["inc"]), float(mpars["alpha95"]))
                  for k in range(len(Da95)):
                      XY = pmag.dimap(Da95[k], Ia95[k])
                      Xcirc.append(XY[0])
                      Ycirc.append(XY[1])
                  self.displayed_means.append(
                      fig.plot(Xcirc, Ycirc, color, alpha=alpha))
  
              if self.ie_open:
                  self.displayed_means.append(self.ie.scatter([XYM[0]], [
                                              XYM[1]], marker='o', edgecolor=EC, facecolor=FC, s=size, lw=1, clip_on=False, alpha=alpha))
                  if "alpha95" in list(mpars.keys()):
                      self.displayed_means.append(
                          self.ie.plot(Xcirc, Ycirc, color, alpha=alpha))
                  self.ie.eqarea.set_xlim(xmin, xmax)
                  self.ie.eqarea.set_ylim(ymin, ymax)
  
          fig.set_xlim(xmin, xmax)
          fig.set_ylim(ymin, ymax)