;$Id: prho.pro,v 1.2 2006/08/03 04:28:07 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=18,yoffset=3
  !p.charthick=4 & !p.thick=4 & !x.thick=4 & !y.thick=4 
end
;
;  mv idl.ps ~/tex/gal/galti/fig/prho.ps
;
file='rho=0.7.sav'
restore,file
;
!p.charsize=1.8
!p.multi=[0,1,2]
!x.range=[0.,.16]
!y.margin=[3.3,1.7]
;
!x.title='!6'
!y.title='!8x!6 [kpc]'
!p.title='!6ln!8T!6'
!y.range=[-1,1]*.1
lev=grange(3.5,10.,12)
contour,shift(transpose(map),0,450),tvid,zvid,/fil,lev=lev
;
!x.title='!8t!6 [Gyr]'
!y.title='!8p!6 [10!u-14!n dyn]'
!y.range=0.
!p.title='!6'
nevery=10
;
;plot_io,tt,ppm,yr=[10,300]
plot_io,every(tt,nevery),every(ppm,nevery),yr=[10,300]
;
;  calculate pressure from mean values
;
Rgas=8.314d7
mu=.62
unit_rho=1d-24
unit_pp=1d-14
unit_TT=1.
fact=Rgas*unit_rho*unit_TT/(unit_pp*mu)
;oplot,tt,fact*rhom*TTm,li=0
oplot,every(tt,nevery),every(fact*rhom*TTm,nevery),li=0
;
siz=1.8
xyouts,.05,100,'!6(!13R/!7l!6)<!7q!6><!8T!6>',siz=siz
xyouts,.08,35,'!6(!13R/!7l!6)<!7q!8T!6>=<!8p!6>',siz=siz
;
END