docs: def plot_selected_meas(self):
          """
          Goes through all measurements selected in logger and draws darker
          marker over all specimen plots to display which measurements have
          been selected
          """
          self.close_warning = True
          # set hex colors for cover and size of selected meas marker
          blue_cover = "#9999FF"
          red_cover = "#FF9999"
          eqarea_outline = "#FF0000"
          MS_selected = 40
  
          # remove old selected points
          for a in self.selected_meas_artists:
              if a in self.zijplot.collections:
                  self.zijplot.collections.remove(a)
              if a in self.specimen_eqarea.collections:
                  self.specimen_eqarea.collections.remove(a)
              if a in self.mplot.collections:
                  self.mplot.collections.remove(a)
              if a in self.mplot_af.collections:
                  self.mplot_af.collections.remove(a)
  
          # do zijderveld plot
          self.selected_meas_artists = []
          x, y, z = self.CART_rot[self.selected_meas,
                                  0], self.CART_rot[self.selected_meas, 1], self.CART_rot[self.selected_meas, 2]
          self.selected_meas_artists.append(self.zijplot.scatter(
              x, -1*y, c=red_cover, marker='o', s=MS_selected, zorder=2))
          self.selected_meas_artists.append(self.zijplot.scatter(
              x, -1*z, c=blue_cover, marker='s', s=MS_selected, zorder=2))
  
          # do down data for eqarea
          x_eq = array([row[0] for i, row in enumerate(self.zij_norm)
                        if i in self.selected_meas and row[2] > 0])
          y_eq = array([row[1] for i, row in enumerate(self.zij_norm)
                        if i in self.selected_meas and row[2] > 0])
          z_eq = abs(array([row[2] for i, row in enumerate(
              self.zij_norm) if i in self.selected_meas and row[2] > 0]))
          if len(x_eq) > 0:
              # from Collinson 1983
              R = array(sqrt(1-z_eq)/sqrt(x_eq**2+y_eq**2))
              eqarea_data_x = y_eq*R
              eqarea_data_y = x_eq*R
              self.selected_meas_artists.append(self.specimen_eqarea.scatter([eqarea_data_x], [
                                                eqarea_data_y], marker='o', edgecolor=eqarea_outline, facecolor="#808080", s=15*self.GUI_RESOLUTION, lw=1, clip_on=False))
  
          # do up data for eqarea
          x_eq = array([row[0] for i, row in enumerate(self.zij_norm)
                        if i in self.selected_meas and row[2] < 0])
          y_eq = array([row[1] for i, row in enumerate(self.zij_norm)
                        if i in self.selected_meas and row[2] < 0])
          z_eq = abs(array([row[2] for i, row in enumerate(
              self.zij_norm) if i in self.selected_meas and row[2] < 0]))
          if len(x_eq) > 0:
              # from Collinson 1983
              R = array(sqrt(1-z_eq)/sqrt(x_eq**2+y_eq**2))
              eqarea_data_x = y_eq*R
              eqarea_data_y = x_eq*R
              self.selected_meas_artists.append(self.specimen_eqarea.scatter([eqarea_data_x], [
                                                eqarea_data_y], marker='o', edgecolor=eqarea_outline, facecolor="#FFFFFF", s=15*self.GUI_RESOLUTION, lw=1, clip_on=False))
  
          # do M/M0 plot
          steps = self.Data[self.s]['zijdblock_steps']
          flags = self.Data[self.s]['measurement_flag']
          selected_af_meas = [
              i for i in self.selected_meas if "T" in steps[i] or steps[i] == "0" and flags[i] != "b"]
          selected_T_meas = [i for i in self.selected_meas if "C" in steps[i]
                             or steps[i] == "0" and flags[i] != "b"]
          data = array(self.Data[self.s]['zijdblock'])
          af_x = array(list(map(float, data[selected_af_meas, 0])))
          af_y = array(
              list(map(float, data[selected_af_meas, 3])))/float(data[0, 3])
          T_x = array(list(map(float, data[selected_T_meas, 0])))
          T_y = array(
              list(map(float, data[selected_T_meas, 3])))/float(data[0, 3])
  
          xmin, xmax = self.mplot.get_xlim()
          ymin, ymax = self.mplot.get_ylim()
          if T_x.astype(float).any() or T_y.astype(float).any():
              self.selected_meas_artists.append(self.mplot.scatter(
                  T_x, T_y, facecolor=red_cover, edgecolor="#000000", marker='o', s=30, lw=1, clip_on=False, zorder=3))
          self.mplot.set_xlim([xmin, xmax])
          self.mplot.set_ylim([ymin, ymax])
  
          xmin, xmax = self.mplot_af.get_xlim()
          ymin, ymax = self.mplot_af.get_ylim()
          if af_x.astype(float).any() or af_y.astype(float).any():
              self.selected_meas_artists.append(self.mplot_af.scatter(
                  af_x, af_y, facecolor=blue_cover, edgecolor="#000000", marker='o', s=30, lw=1, clip_on=False, zorder=3))
          self.mplot_af.set_xlim([xmin, xmax])
          self.mplot_af.set_ylim([ymin, ymax])
  
          self.canvas1.draw()
          self.canvas2.draw()
          self.canvas3.draw()