IDL>restore,'siess_tracks25m.sav'IDL>plot,alog10(teff1d0),alog10(lumi1d0),psym=-3,xrange=[4.5,3.4],yrange=[-3,4],xtitle='Teff',ytitle='L',$> ystyle=1,xstyle=1IDL>oplot,alog10(teff0d10),alog10(lumi0d10),psym=-3IDL>oplot,alog10(teff0d13),alog10(lumi0d13),psym=-3IDL>oplot,alog10(teff0d16),alog10(lumi0d16),psym=-3IDL>oplot,alog10(teff0d20),alog10(lumi0d20),psym=-3IDL>oplot,alog10(teff0d25),alog10(lumi0d25),psym=-3IDL>oplot,alog10(teff0d30),alog10(lumi0d30),psym=-3IDL>oplot,alog10(teff0d40),alog10(lumi0d40),psym=-3IDL>oplot,alog10(teff0d50),alog10(lumi0d50),psym=-3IDL>oplot,alog10(teff0d60),alog10(lumi0d60),psym=-3IDL>oplot,alog10(teff0d70),alog10(lumi0d70),psym=-3IDL>oplot,alog10(teff0d80),alog10(lumi0d80),psym=-3IDL>oplot,alog10(teff0d90),alog10(lumi0d90),psym=-3IDL>oplot,alog10(teff1d1),alog10(lumi1d1),psym=-3IDL>oplot,alog10(teff1d2),alog10(lumi1d2),psym=-3IDL>oplot,alog10(teff1d3),alog10(lumi1d3),psym=-3IDL>oplot,alog10(teff1d4),alog10(lumi1d4),psym=-3IDL>oplot,alog10(teff1d5),alog10(lumi1d5),psym=-3IDL>oplot,alog10(teff1d6),alog10(lumi1d6),psym=-3IDL>oplot,alog10(teff1d7),alog10(lumi1d7),psym=-3IDL>oplot,alog10(teff1d8),alog10(lumi1d8),psym=-3IDL>oplot,alog10(teff1d9),alog10(lumi1d9),psym=-3IDL>oplot,alog10(teff2d0),alog10(lumi2d0),psym=-3IDL>oplot,alog10(teff2d2),alog10(lumi2d2),psym=-3IDL>oplot,alog10(teff2d5),alog10(lumi2d5),psym=-3IDL>oplot,alog10(teff2d7),alog10(lumi2d7),psym=-3IDL>oplot,alog10(teff3d0),alog10(lumi3d0),psym=-3IDL>oplot,alog10(teff3d5),alog10(lumi3d5),psym=-3IDL>oplot,alog10(teff4d0),alog10(lumi4d0),psym=-3IDL>oplot,alog10(teff5d0),alog10(lumi5d0),psym=-3IDL>oplot,alog10(teff6d0),alog10(lumi6d0),psym=-3IDL>oplot,alog10(teff7d0),alog10(lumi7d0),psym=-3IDL>fori=0,31dooplot,alog10(teff[i,*]),alog10(lbol[i,*]),line=2IDL>