docs: def tcrhoc(self,ifig=None,lims=[3.,10.,8.,10.],label=None,colour=None,
                 dashes=None):
          """
          Central temperature again central density plot
  
          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
          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 label is not None:
              if colour is not None:
                  line,=pl.plot(self.get('log_center_Rho'),self.get('log_center_T'),label=label,
                          color=colour)
              else:
                  line,=pl.plot(self.get('log_center_Rho'),self.get('log_center_T'),label=label)
  
          else:
              if colour is not None:
                  line,=pl.plot(self.get('log_center_Rho'),self.get('log_center_T'),
                          color=colour)
              else:
                  line,=pl.plot(self.get('log_center_Rho'),self.get('log_center_T'))
          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.xlabel('log $\\rho_{\\rm c}$')
          pl.ylabel('log $T_{\\rm c}$')