pro lit_evol,fich,nbx,age,ray,lum,teff,rhp common lsj,lsj car='rien' unit=1 close,unit openr,unit,fich readf,unit,car readf,unit,car readf,unit,nbx print,'dimension du tableau:',nbx age=dblarr(nbx) ray=dblarr(nbx) lum=dblarr(nbx) teff=dblarr(nbx) rhp=dblarr(nbx) a=dblarr(6) for i=0,nbx-1 do begin readf,unit,a age(i)=a(0) ray(i)=a(1) lum(i)=alog10(a(2)*lsj) teff(i)=a(3) rhp(i)=a(5) endfor close,unit end ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; pro lit_bur,fich,nbx,age,ray,lum,teff,rhp common lsj,lsj car='rien' unit=1 close,unit openr,unit,fich for i=1,45 do readf,unit,car pn=700 nbx=0 age=dblarr(pn) ray=dblarr(pn) lum=dblarr(pn) teff=dblarr(pn) rhp=dblarr(pn) a=dblarr(11) repeat begin readf,unit,a age(nbx)=a(1)*1d3 ;Myrs ray(nbx)=a(4)/6.9894 ;Rj lum(nbx)=alog10(a(3)) ;Lsol teff(nbx)=a(2) rhp(nbx)=a(5) ;Tc nbx=nbx+1 endrep until (teff(nbx-1) le 960.) print,'dimension du tableau:',pn age=age(0:nbx-1) ray=ray(0:nbx-1) lum=lum(0:nbx-1) teff=teff(0:nbx-1) rhp=rhp(0:nbx-1) close,unit end ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; pro diaghr ;colors: Black,white,red,red,green,blue,gb,rg,rb,orange (original) ; obs lowt medt ht - - obs ; 0 1 2 3 4 5 6 7 8 9 ;red =[0,1,1,1,0,0,0,1,1,1] ;green=[0,1,0,0,1,0,1,1,0,0.87] ;blue =[0,1,0,0,0,1,1,0,1,0] ;tvlct,255*red,255*green,255*blue loadct,13 common lsj,lsj a=findgen(4)*(!pi *2/4.) usersym,cos(a),sin(a),/fill lsj=3.35d24/3.826d33 xrg=[2000,10] yrg=[-10,-2] ww=(xrg(0)-xrg(1))*3/100. ;lit_evol,'evol.des',nbx1,age1,ray1,lum1,teff1,rhp plot,[0,0],[0,0],xtitle='!3 Effective Temperature (K)',$ ytitle='Log Luminosity (L!D!9n!3!N)',$ ; title='Newly discovered planets',$ xrange=xrg,yrange=yrg,thick=3,charthick=3,linestyle=0,$ charsize=1.8,xstyle=1,ystyle=1,xminor=5,yminor=4,/nodata plot,[0,0],[0,0],xtitle='!3 Effective Temperature (K)',$ ytitle='Log Luminosity (L!D!9n!3!N)',$ ; title='Newly discovered planets',$ xrange=xrg,yrange=yrg,thick=3,charthick=3,linestyle=0,$ charsize=1.8,xstyle=1,ystyle=1,xminor=5,yminor=4,/nodata,/noerase ;findage,nbx1,age1,teff1,lum1,ww lit_evol,'Exo/ev51peg-0.46-1190.des',nbx,age,ray,lum,teff,rhp n=0 while ((n lt nbx-1)and(age(n) lt 1d4)) do n=n+1 oplot,teff(0:n-1),lum(0:n-1),linestyle=0,thick=3,col=255 oplot,[teff(n-1)],[lum(n-1)],thick=2,psym=8,symsize=1.5,col=255 lit_evol,'Exo/evtauboo-3.87-1300.des',nbx,age,ray,lum,teff,rhp n=0 while ((n lt nbx-1)and(age(n) lt 6d3)) do n=n+1 oplot,teff(0:n-1),lum(0:n-1),linestyle=0,thick=3,col=63 oplot,[teff(n-1)],[lum(n-1)],thick=2,psym=8,symsize=1.5,col=63 lit_evol,'Exo/ev55cnc-0.8-637.des',nbx,age,ray,lum,teff,rhp n=0 while ((n lt nbx-1)and(age(n) lt 1d4)) do n=n+1 oplot,teff(0:n-1),lum(0:n-1),linestyle=0,thick=3,col=111 oplot,[teff(n-1)],[lum(n-1)],thick=2,psym=8,symsize=1.5,col=111 lit_evol,'Exo/evlnl-0.9-42.des',nbx,age,ray,lum,teff,rhp n=0 while ((n lt nbx-1)and(age(n) lt 1d4)) do n=n+1 oplot,teff(0:n-1),lum(0:n-1),linestyle=0,thick=2;,col=6 oplot,[teff(n-1)],[lum(n-1)],thick=2,psym=8,symsize=1.5;,col=6 lit_evol,'Exo/ev47uma-2.39-173.des',nbx,age,ray,lum,teff,rhp n=0 while ((n lt nbx-1)and(age(n) lt 8d3)) do n=n+1 oplot,teff(0:n-1),lum(0:n-1),linestyle=0,thick=3,col=159 oplot,[teff(n-1)],[lum(n-1)],thick=2,psym=8,symsize=1.5,col=159 lit_evol,'Exo/ev70vir-6.6-335.des',nbx,age,ray,lum,teff,rhp n=0 while ((n lt nbx-1)and(age(n) lt 1d4)) do n=n+1 oplot,teff(0:n-1),lum(0:n-1),linestyle=0,thick=3,col=255 oplot,[teff(n-1)],[lum(n-1)],thick=2,psym=8,symsize=1.5,col=255 lit_evol,'Exo/evhd-10-413.des',nbx,age,ray,lum,teff,rhp n=0 while ((n lt nbx-1)and(age(n) lt 6d3)) do n=n+1 oplot,teff(0:n-1),lum(0:n-1),linestyle=0,thick=3,col=111 oplot,[teff(n-1)],[lum(n-1)],thick=2,psym=8,symsize=1.5,col=111 lit_evol,'Exo/evupsand-0.6-1140.des',nbx,age,ray,lum,teff,rhp n=0 while ((n lt nbx-1)and(age(n) lt 6d3)) do n=n+1 oplot,teff(0:n-1),lum(0:n-1),linestyle=0,thick=3,col=207 oplot,[teff(n-1)],[lum(n-1)],thick=2,psym=8,symsize=1.5,col=207 lit_bur,'/export/nanosun/burrows/mbd1.0/mbd1.05',nbx,age,ray,lum,teff,rhp oplot,teff,lum,linestyle=0,thick=3,col=207 oplot,[teff(nbx-1)],[lum(nbx-1)],thick=2,psym=8,symsize=1.5,col=207 ;carre,'Exo/ev51peg.txt',3 ;carre,'Exo/evtauboo.txt',4 ;carre,'Exo/ev55cnc.txt',5 ;carre,'Exo/evlnl21185.txt',6 ;carre,'Exo/ev47uma.txt',7 ;carre,'Exo/ev70vir.txt',8 ;carre,'Exo/evhd114762.txt',0 noms,1180,-4.9,'!351 Peg b',1.0,1.6,255 noms,1370,-4.5,'!4s!3 Boo b',1.0,1.6,63 noms,1100,-4.6,'!4t!3 And b',0.0,1.6,207 noms,1000,-5.4,'Gl 229 B',1.0,1.6,207 noms,700,-5.9,'55 Cnc b',1.0,1.6,111 noms,220,-9.3,'Lal 21185 b',1.0,1.6,0 noms,230,-8.0,'47 UMa b',1.0,1.6,159 noms,400,-7.0,'70 Vir b',1.0,1.6,255 noms,500,-6.5,'HD 114762 b',1.0,1.6,111 a=findgen(16)*(!pi *2/16.) usersym,cos(a),sin(a),/fill oplot,[124.4],[-9.062],psym=8,symsize=1.5,thick=3 oplot,[95.0],[-9.651],psym=8,symsize=1.5,thick=3 xyouts,200,-9.062-.1,'!20s!3',align=0.0,charsize=1.6,charthick=2 xyouts,180,-9.651-.1,'!20T!3',align=0.0,charsize=1.6,charthick=2 ;findage,nbx,age,teff,lum,ww ;stop ;read,car ;plot,mass,rayon,xstyle=1,ystyle=1,xrange=[0,3.5],yrange=[0.6,1.8] ;oplot,mass,rayon,psym=1 ;plot,age2,ray2,xrange=[7,10],yrange=[0,16],ystyle=1,charsize=1.6,thick=2,$ ; xtitle='!3 log(t) (yrs)',ytitle='R/R!DJup!N',charthick=2 ;oplot,age3,ray3,linestyle=2,thick=2 ;oplot,age4,ray4 end ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; pro noms,xx,yy,nom,al,chrsz,cl if (cl ne 0) then xyouts,xx-2.5,yy-0.02,nom,align=al,charsize=chrsz,charthick=2,col=cl xyouts,xx,yy,nom,align=al,charsize=chrsz,charthick=2 end ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; pro carre,fic,cl common lsj,lsj car='rien' close,1 openr,1,fic readf,1,car a=fltarr(8) l=a t=a for i=0,3 do begin readf,1,car readf,1,a l(2*i+0)=a(2) t(2*i+0)=a(3) l(2*i+1)=a(6) t(2*i+1)=a(7) endfor close,1 l=alog10(l*lsj) for i=0,5 do begin for j=i+1,6 do begin for k=j+1,7 do begin oplot,[t(i),t(j),t(k)],[l(i),l(j),l(k)],thick=2 polyfill,[t(i),t(j),t(k)],[l(i),l(j),l(k)],col=cl;,col=0.50*!d.n_colors endfor endfor endfor for i=0,5 do begin for j=i+1,6 do begin for k=j+1,7 do begin polyfill,[t(i),t(j),t(k)],[l(i),l(j),l(k)],col=cl;,col=0.50*!d.n_colors endfor endfor endfor end ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; pro verif lit_evol,'evol.des',nbx,age,ray,lum,teff,rhp age=age-min(age)+1d-6 age=alog10(age) plot,age,ray oplot,age,1*(teff/teff(0)-1)+1,linestyle=2 end ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; pro findage,nbx,age,x,y,ww,coul lumm=7d-4 ;7d-4 ;7d-3 y=10^y if (age(nbx-1) lt 0) then begin bid=age(nbx-1) age=age-bid ia=1 while (y(ia) lt lumm) do ia=ia+1 bid=age(ia-1)+(lumm-y(ia-1))*(age(ia)-age(ia-1))/(y(ia)-y(ia-1)) age(ia)=bid y(ia)=lumm age=age-bid for i=ia+1,nbx-1 do begin age(i)=0 x(i)=x(ia) y(i)=y(ia) endfor ageb=age xb=x yb=y for i=0,nbx-1 do begin age(i)=ageb(nbx-1-i) x(i)=xb(nbx-1-i) y(i)=yb(nbx-1-i) ; print,i,age(i),x(i),y(i) endfor endif for i=0,nbx-1 do age(i)=6+alog10(max([1.d-33,age(i)])) y=alog10(y) logt=[-100,7,7.5,8,8.5,9,9.5,10,11] ;c=['0','7','7.5','8','8.5','9','9.5','10'] c=[' ','7',' ','8',' ','9',' ','10'] nba=8 tx=dblarr(nba) ty=dblarr(nba) tx(0)=x(0) ty(0)=y(0) ia=1 while (logt(ia) lt age(1)) do ia=ia+1 for i=1,nbx-2 do begin if ((age(i) gt logt(ia)) and (age(i+1) gt logt(ia))) then ia=ia+1 if ((age(i) le logt(ia)) and (age(i+1) gt logt(ia))) then begin tx(ia)=x(i);+(logt(ia)-age(i))*(x(i+1)-x(i))/(age(i+1)-age(i)) ty(ia)=y(i);+(logt(ia)-age(i))*(y(i+1)-y(i))/(age(i+1)-age(i)) print,ia,age(i),tx(ia),ty(ia) ia=ia+1 endif endfor if (age(nbx-1) eq logt(ia)) then begin tx(ia)=x(nbx-1) ty(ia)=y(nbx-1) ia=ia+1 endif ia=nba oplot,tx(1:nba-1),ty(1:nba-1),psym=8,symsize=0.8,thick=2;,col=coul for i=0,nba-1 do begin w=ww*(((i mod 2)-1)*1. + .5) xyouts,tx(i)+w,ty(i)+alog10(.95),c(i),alignment=(i mod 2),charsize=.9;,col=coul endfor end ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; pro agebis,nbx,age,x,y,ww,car,bid lumm=8d-3 ;7d-4 ;7d-3 y=10^y if (age(nbx-1) lt 0) then begin bid=age(nbx-1) age=age-bid ia=1 while (y(ia) lt lumm) do ia=ia+1 bid=age(ia-1)+(lumm-y(ia-1))*(age(ia)-age(ia-1))/(y(ia)-y(ia-1)) age(ia)=bid y(ia)=lumm age=age-bid for i=ia+1,nbx-1 do begin age(i)=0 x(i)=x(ia) y(i)=y(ia) endfor ageb=age xb=x yb=y for i=0,nbx-1 do begin age(i)=ageb(nbx-1-i) x(i)=xb(nbx-1-i) y(i)=yb(nbx-1-i) ; print,i,age(i),x(i),y(i) endfor endif for i=0,nbx-1 do age(i)=6+alog10(max([1.d-33,age(i)])) y=alog10(y) logt=[-100,7,7.5,8,8.5,9,9.5,10,11] ;c=['0','7','7.5','8','8.5','9','9.5','10'] c=[' ','7',' ','8',' ','9',' ','10'] nba=8 tx=dblarr(nba) ty=dblarr(nba) tx(0)=x(0) ty(0)=y(0) ia=1 while (logt(ia) lt age(1)) do ia=ia+1 for i=1,nbx-2 do begin if ((age(i) gt logt(ia)) and (age(i+1) gt logt(ia))) then ia=ia+1 if ((age(i) le logt(ia)) and (age(i+1) gt logt(ia))) then begin tx(ia)=x(i);+(logt(ia)-age(i))*(x(i+1)-x(i))/(age(i+1)-age(i)) ty(ia)=y(i);+(logt(ia)-age(i))*(y(i+1)-y(i))/(age(i+1)-age(i)) print,ia,age(i),tx(ia),ty(ia) ia=ia+1 endif endfor if (age(nbx-1) eq logt(ia)) then begin tx(ia)=x(nbx-1) ty(ia)=y(nbx-1) ia=ia+1 endif ia=nba oplot,tx(1:nba-1),ty(1:nba-1),psym=8,symsize=0.8,thick=2;,col=coul for i=0,nba-1 do begin w=ww*(((i mod 2)-1)*1. + .5) xyouts,tx(i)+w,ty(i)+alog10(.95),c(i),alignment=(i mod 2),charsize=.9;,col=coul endfor end ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; pro evorad age=[2.7,4.76,10,30,70,100,300,1000,3000,10000] ;prad=[7.395,7.187,6.823,6.894,7.364,7.446,7.939,8.299,8.748,9.315] prad=[7.395,7.187,7.0,7.0,7.364,7.446,7.939,8.299,8.748,9.315] pc=[9.763,9.858,10.22,13.06,13.17,13.21,13.29,13.35,13.40,13.45] age=alog10(age)+6 prad=prad-6 pc=pc-6 plot,age,prad,yrange=[7.5,0],xstyle=1,ystyle=1,charsize=1.3,thick=2,$ xtitle='!17log t (Yrs)',ytitle='log P (bar)',charthick=2 ; title='Evolution of a Jupiter at 0.05 A.U. from a G2V star' pxage=[age(0),age,age(9)] pypr=[1,prad,1] pypc=[7.5,pc,7.5] polyfill,[age(0),age(9),age(9),age(0),age(0)],[7.5,7.5,0,0,7.5],col=255 polyfill,pxage,pypr,col=0.75*!d.n_colors ;207 polyfill,pxage,pypc,col=0.50*!d.n_colors ;0 plot,age,prad,yrange=[7.5,0],xstyle=1,ystyle=1,charsize=1.3,thick=2,$ xtitle='log t (Yrs)',ytitle='log P (bar)',charthick=2,$ ; title='Evolution of a Jupiter at 0.05 A.U. from a G2V star',$ /noerase,/nodata oplot,age,prad,thick=2,col=255 oplot,age,pc,thick=2 ;oplot,[6,10],[1,1],thick=2 xyouts,8.2,0.8,'"Atmosphere"',charsize=1.6,charthick=3,alignment=0.5 xyouts,9.0,1.5,'Radiative Region',charsize=1.6,charthick=3,alignment=0.5 xyouts,8.5,4.3,'Convective Region',charsize=1.6,charthick=3,alignment=0.5 xyouts,7.27,5.53,'Central!C Pressure',charsize=1.6,charthick=3,alignment=0.5,$ orientation=-65;,col=255 end ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; pro hpt common lsj,lsj lsj=3.35d24/3.826d33 lit_evol,'../EGP/ev050_100.des',nbx,age,ray,lum,teff,rhp set_viewport,.1,.95,.6,.95 plot_io,ray,rhp,ytitle='R/H!DP!N(10bar)',$ title='1M!DJ!N at 5.2, 0.1, 0.05 and 0.02 A.U.',yrange=[1,3000],$ ystyle=1,charsize=1.6,xrange=[0,15] lit_evol,'../EGP/ev100_100.des',nbx,age,ray,lum,teff,rhp oplot,ray,rhp,linestyle=2 lit_evol,'../EGP/ev520_100r.des',nbx,age,ray,lum,teff,rhp oplot,ray,rhp,linestyle=1,thick=3 lit_evol,'../EGP/ev020_100.des',nbx,age,ray,lum,teff,rhp oplot,ray,rhp,linestyle=3 set_viewport,.1,.95,.15,.50 lit_evol,'../EGP/ev050_100.des',nbx,age,ray,lum,teff,rhp times,tau,nbx,age,ray plot,ray,tau,xtitle='!17 R/R!DJup!N',ytitle='Log(Tau!Dcontraction!N) (Gyrs)',$ yrange=[-5,2],ystyle=1,xrange=[0,15],charsize=1.6,/noerase lit_evol,'../EGP/ev100_100.des',nbx,age,ray,lum,teff,rhp times,tau,nbx,age,ray oplot,ray,tau,linestyle=2 lit_evol,'../EGP/ev520_100r.des',nbx,age,ray,lum,teff,rhp oplot,ray,tau,linestyle=1,thick=3 times,tau,nbx,age,ray lit_evol,'../EGP/ev020_100.des',nbx,age,ray,lum,teff,rhp times,tau,nbx,age,ray oplot,ray,tau,linestyle=3 end ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; pro times,tau,nbx,age,ray if (age(nbx-1) lt 0) then begin age1=age ray1=ray for i=0,nbx-1 do begin age(i)=age1(nbx-1-i) ray(i)=ray1(nbx-1-i) endfor endif tau=dblarr(nbx) for i=nbx-1,1,-1 do begin j=i-1 while (j gt 0) and (ray(j) lt 2*ray(i)) do j=j-1 tau(i)=-33 if (j gt 0) then tau(i)=alog10((age(i)-age(j))/1000.) ; print,i,j,tau(i),ray(i),ray(j) endfor tau(0)=-33 end ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; pro diag loadct,13 peg=255 jup=63 tid=192 hay=223 common lsj,lsj a=findgen(16)*(!pi *2/16.) usersym,cos(a),sin(a),/fill lsj=3.35d24/3.826d33 lxx=1/2.18636d-9 alx=alog10(lxx) xrg=[2200,100] yrg=[-9,-2]+alx ww=(xrg(0)-xrg(1))*3/100. lit_evol,'../EGP/ev520_100m.des',nbx1,age1,ray1,lum1,teff1,rhp lum1=lum1+alx plot,teff1,lum1,xtitle='!17 Surface Temperature (Kelvins)',ytitle='Luminosity',$ xrange=xrg,yrange=yrg,thick=2,yminor=1,xminor=1,$ ytickname=['1',' ','10!E2!N',' ','10!E4!N',' ','10!E6'],$ charsize=1.6,xstyle=1,ystyle=1,linestyle=2 ;--- Limite de Roche ---; mpla=1 ;rox=[xrg(0),1.33818D+03,1.22634D+03,9.6d2,1d2] ;roy=[5.98d5,5.98d5,7.23735D+05,1.5d6,2d7] rox=[xrg(0),xrg(1)] roy=[5.98d5,5.98d5] roy=alog10(roy*lsj*lxx*mpla^(2/3.)) rochex=[rox(0),rox,rox(1)] rochey=[yrg(1),roy,yrg(1)] polyfill,rochex,rochey,col=tid oplot,rox,roy,linestyle=1,thick=2,col=tid xyouts,2000,-2.6+alx,'!17Planet is torn apart!Cby star`s gravity',charsize=1.6,charthick=1 ;--- Limite de Hayashi ---; hayx=[xrg(1),teff1(0),teff1,teff1(nbx1-1),xrg(1)] hayy=[yrg(1),lum1(0),lum1,lum1(nbx1-1),yrg(0)] ;for i=0,nbx1+1 do begin ; if (hayy(i) gt yrg(1)) then begin ; hayy(i)=yrg(1) ; hayx(i)=xrg(1) ; endif ;endfor polyfill,hayx,hayy,col=hay,clip=[1500,1e-5,1000,1e-2] xyouts,700,-4+alx,'Forbidden!CRegion',charsize=1.6,charthick=1 ;--------------------------------------------------------------------; plot,teff1,lum1,xtitle='!17 Surface Temperature (Kelvins)',ytitle='Luminosity',$ xrange=xrg,yrange=yrg,thick=2,yminor=1,xminor=1,$ ytickname=['1',' ','10!E2!N',' ','10!E4!N',' ','10!E6'],$ charsize=1.6,xstyle=1,ystyle=1,linestyle=2,/noerase te=dblarr(5) lu=dblarr(5) lit_evol,'../EGP/ev520_100r.des',nbx,age,ray,lum,teff,rhp lum=lum+alx oplot,teff,lum,linestyle=0,thick=10,col=jup agebis,nbx,age,teff,lum,ww,$ [' ','10 million yrs ','100 million yrs ','1 billion yrs ',' '],jup oplot,[124.4],[0],psym=4,symsize=1.5,thick=4,col=jup xyouts,[124.4],[0],'!17Jupiter ',alignment=1.0,$ charsize=1.6,charthick=3 oplot,[430,140],[-0.1,-0.1],thick=4,col=jup lit_evol,'../EGP/ev050_100.des',nbx,age,ray,lum,teff,rhp lum=lum+alx age1=age teff1=teff lum1=lum agebis,nbx,age,teff,lum,ww,$ [' ','10 million yrs ','100 million yrs ','1 billion yrs ','10 billion yrs '],peg ; [' ','10!E7!N yrs ','10!E8!N yrs ','10!E9!N yrs ','10!E10!N yrs '],peg lum=lum-.05 oplot,teff,lum,linestyle=0,thick=10,col=peg agebis,nbx,age1,teff1,lum1,ww,$ [' ','10 million yrs ','100 million yrs ','1 billion yrs ','10 billion yrs '],peg te(1)=teff(nbx-1) oplot,[te(1)],[lum(nbx-1)+0.03],psym=5,symsize=1.5,thick=4,col=peg xyouts,te(1),lum(nbx-1)-0.04,' Case a',alignment=0.0,charsize=1.2 oplot,[1230,1050],[lum(nbx-1)-0.07,lum(nbx-1)-0.07],col=peg xyouts,te(1),-6.5+alx,'!1751PegB',charsize=1.6,charthick=3,alignment=0.5 oplot,[1400,1110],[-6.6,-6.6]+alx,thick=4,col=peg rtell=0.35 lu=alog10(4*!pi*(rtell*7d9)^2*5.67051d-5*te^4/3.826d33*lxx) oplot,te,lu+0.03,psym=5,symsize=1.5,thick=4,col=peg xyouts,te(1),lu(1)-0.03,' Case b',alignment=0.0,charsize=1.2 oplot,[1230,1050],[lu(1)-0.06,lu(1)-0.06],col=peg x=xrg(0)+findgen(100)/99.*(xrg(1)-xrg(0)) cte=9.126d-18*lxx y=alog10(cte/16*x^4) ;L pour 1/4 Rj cc=['1/4 R!DJ!N','1/2 R!DJ!N','1 R!DJ!N','2 R!DJ!N',' '] for i=1,5 do begin oplot,x,y,linestyle=1 xyouts,2100,y(0),cc(i-1),charthick=1,charsize=1.2,orientation=-12 y=y+alog10(4) endfor end ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; pro hr ;loadct,13 peg=255 jup=63 tid=192 hay=223 common lsj,lsj a=findgen(16)*(!pi *2/16.) usersym,cos(a),sin(a),/fill lsj=3.35d24/3.826d33 lxx=1/2.18636d-9 alx=alog10(lxx) xrg=[2200,100] yrg=[-9,-2]+alx ww=(xrg(0)-xrg(1))*3/100. lit_evol,'evol.des',nbx1,age1,ray1,lum1,teff1,rhp lum1=lum1+alx plot,teff1,lum1,xtitle='!17 Surface Temperature (Kelvins)',ytitle='Luminosity',$ xrange=xrg,yrange=yrg,thick=2,yminor=1,xminor=1,$ ytickname=['1',' ','10!E2!N',' ','10!E4!N',' ','10!E6'],$ charsize=1.6,xstyle=1,ystyle=1,linestyle=1,/nodata ;--------------------------------------------------------------------; plot,teff1,lum1,xtitle='!17 Surface Temperature (Kelvins)',ytitle='Luminosity',$ xrange=xrg,yrange=yrg,thick=2,yminor=1,xminor=1,$ ytickname=['1',' ','10!E2!N',' ','10!E4!N',' ','10!E6'],$ charsize=1.6,xstyle=1,ystyle=1,linestyle=1,/noerase,/nodata te=dblarr(5) lu=dblarr(5) ;lit_evol,'evol.des',nbx,age,ray,lum,teff,rhp ;lum=lum+alx ;oplot,teff,lum,linestyle=0,thick=2;,col=jup ;findage,nbx,age,teff,lum,ww,0 cmas='2.0' ctem=['0','100','200','300','400','500','600','700','800','900','1000',$ '1100','1200','1300','1400','1500'] for i=0,15 do begin fich='MV/ev'+cmas+'-'+ctem(i)+'.des' lit_evol,fich,nbx,age,ray,lum,teff,rhp lum=lum+alx oplot,teff,lum,linestyle=0,thick=2;,col=jup ;findage,nbx,age,teff,lum,ww,0 endfor oplot,[124.4],[0],psym=4,symsize=1.5,thick=2;,col=jup x=xrg(0)+findgen(100)/99.*(xrg(1)-xrg(0)) cte=9.126d-18*lxx y=alog10(cte/16*x^4) ;L pour 1/4 Rj cc=['1/4 R!DJ!N','1/2 R!DJ!N','1 R!DJ!N','2 R!DJ!N','4 R!DJ!N',' ',' ',' '] for i=1,8 do begin oplot,x,y,linestyle=1 xyouts,xrg(0)-100,y(0),cc(i-1),charthick=1,charsize=1.2,orientation=-12 y=y+alog10(4) endfor xyouts,0.90,0.85,cmas+' M!DJ!N',charsize=1.8,charthick=2,/normal,align=1.0 end ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; pro psdiag set_plot,'ps' device,/landscape,/color,/bold loadct,13 ;evorad diaghr device,/close set_plot,'x' loadct,0 end