pro lit_des,fich,logp,logt

ct1=['M/Mtot','R','logP','logT','logRho','logL','Krad',$
	'S']
car='rien'
close,1
openr,1,fich
for i=1,6 do readf,1,car
readf,1,n
print,'Nombre de couches: ',n
a=fltarr(8,n)
readf,1,a
close,1

logp=fltarr(n)
logp(*)=a(2,*)
logt=fltarr(n)
logt(*)=a(3,*)

end

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
pro struct

fic1='results/hd209458/1400k/alex/str-std-1,35rj.dat'
fic2='results/hd209458/1400k/alex-1000/str-1000-1,35rj.dat'

lit_des,fic1,logp1,logt1
lit_des,fic2,logp2,logt2

!x.style=1
!y.style=1
!x.range=[3,5]
!y.range=[13,5]
!p.thick=4
!p.charthick=3
!p.charsize=1.4
!x.title='LogT [K]'
!y.title='LogP [dyn cm!E-2!N]'
plot,logt1,logp1,/nodata

oplot,logt1,logp1,color=255
oplot,logt2,logp2,color=124

end
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
pro ps
set_plot,'ps'
device,/color,file='struct.ps'
loadct,13
struct
device,/close
set_plot,'x'
end