docs: def mcc_t(self,ifig=None,lims=[0,15,0,25],label=None,colour=None,
                mask=False,s2ms=False,dashes=None):
          """
          Plot mass of [oclark01@scandium 15M_led_f_print_nets]$ cp ../15M_led_f_ppcno/ective core as a function of time.
  
          Parameters
          ----------
          ifig : integer or string
              Figure label, if None the current figure is used
              The default value is None.
          lims : list [x_lower, x_upper, y_lower, y_upper]
          label : string
              Label for the model
              The default value is None
          colour : string
              The colour of the line
              The default value is None
          mask : boolean, optional
              Do you want to try to hide numerical spikes in the
              plot?
              The default is False
          s2ms : boolean, optional
              skip to main squence?
          dashes : list, optional
              Custom dashing style. If None, ignore.
              The default is None.
  
          """
  
          fsize=18
  
          params = {'axes.labelsize':  fsize,
          #    'font.family':       'serif',
          'font.family':        'Times New Roman',
          'figure.facecolor':  'white',
          'text.fontsize':     fsize,
          'legend.fontsize':   fsize,
          'xtick.labelsize':   fsize*0.8,
          'ytick.labelsize':   fsize*0.8,
          'text.usetex':       False}
  
          try:
              pl.rcParams.update(params)
          except:
              pass
  
          if ifig is not None:
              pl.figure(ifig)
  
          if s2ms:
              h1=self.get('center_h1')
              idx=np.where(h1[0]-h1>=3.e-3)[0][0]
              skip=idx
          else:
              skip=0
  
          age= self.get('star_age')
          x1 = old_div(age, 1.e6)
          x2 = old_div(age, 1.e6)
          y1 = self.get('mix_qtop_1')*self.get('star_mass')
          y2 = self.get('mix_qtop_2')*self.get('star_mass')
          mt1 = self.get('mix_type_1')
          mt2 = self.get('mix_type_2')
  
          x1 = x1[skip:]
          x2 = x2[skip:]
          y1 = y1[skip:]
          y2 = y2[skip:]
          mt1 = mt1[skip:]
          mt2 = mt2[skip:]
  
          # Mask spikes...
          if mask:
              x1 = np.ma.masked_where(mt1 != 1, x1)
              x2 = np.ma.masked_where(mt2 != 1, x2)
              y1 = np.ma.masked_where(mt1 != 1, y1)
              y2 = np.ma.masked_where(mt2 != 1, y2)
  
          if ifig is not None:
              pl.figure(ifig)
          if label is not None:
              if colour is not None:
                  line,=pl.plot(x1,y1,label=label,color=colour)
                  line,=pl.plot(x2,y2,color=colour)
              else:
                  line,=pl.plot(x1,y1,label=label)
                  line,=pl.plot(x2,y2)
          else:
              if colour is not None:
                  line,=pl.plot(x1,y1,color=colour)
                  line,=pl.plot(x2,y2,color=colour)
              else:
                  line,=pl.plot(x1,y1)
                  line,=pl.plot(x2,y2)
  
          if dashes is not None:
              line.set_dashes(dashes)
  
          if label is not None:
              pl.legend(loc='best').draw_frame(False)
  
          pl.xlim(lims[:2])
          pl.ylim(lims[2:])
          pl.ylabel('$M/M_\odot}$')
          pl.xlabel('$t/{\\rm Myr}$')