pro dhfe
;This IDL program plots dh vs dfe with error bars (from dhp to dhm) and
;(from dfep to dfem) and a weighted least squares linear fit to the data
; The plot is in file idl.ps
set_plot,'ps'
device,/color
name = ['36 Oph', 'Eps Eri', '31 Com',  'HZ 43', 'Procyon', 'Beta Cas',$ 
'G191-B2B', 'Capella', 'Beta Gem', 'Alpha Tri', 'WD0621-376', 'WD2211-495',$
'WD1634-573', 'IX Vel', 'Gam2 Vel', 'REJ1738+665', 'Gam Cas', 'Iota Ori',$
'Del Ori', 'Eps Ori', 'PG0038+199', 'JL 9', 'HD 195965', 'LSS 1274',$
'HD 191877', 'HD 90087']

sym = [4, 4, 4, 6, 4, 4,$ 
6, 4, 4, 4, 6, 6,$
6, 6, 5, 6, 2, 2,$
5, 2, 6, 6, 6, 6,$
6, 6]

dh = [15.0, 14.0, 20.0, 16.6, 16.0, 17.0,$ 
16.6, 16.0, 14.7, 13.2, 14.1, 15.1,$
16.0, 12.3, 24.6, 20.0, 10.5, 13.8,$
6.5, 5.9, 24.3, 9.8, 8.4, 7.5,$
7.7, 9.7]

dhp = [17.5, 16.0, 22.0, 18.0, 18.0, 18.5,$ 
17.05, 17.4, 16.7, 15.2, 19.7, 21.1,$
18.5, 21.4, 27.8, 22.6, 7.5, 19.4,$
7.8, 9.1, 31.5, 11.8, 10.1, 9.3,$
10.3, 11.7]

dhm = [12.5, 12.0, 18.0, 15.2, 14.0, 15.5,$ 
16.3, 14.1, 12.7, 10.2, 8.5, 9.1,$
13.5, 6.6, 21.9, 17.4, 13.5, 2.5,$
5.5, 2.7, 18.9, 7.8, 7.1, 5.7,$
6.4, 7.7]

dfe = [-0.71, -0.83, -0.84, -1.03, -1.04, -1.17,$ 
-0.64, -0.94, -1.09, -1.03, -1.33, -1.02,$
-1.28, -1.19, -1.14, -1.34, -1.45, -1.47,$
-1.65, -1.71, -1.57, -1.60, -1.65, -1.68,$
-1.61, -1.51]

dfep = [-0.46, -0.67, -0.74, -0.93, -0.99, -1.07,$ 
-0.55, -0.89, -0.95, -0.83, -1.02, -0.87,$
-1.15, -1.01, -1.01, -1.24, -1.25, -1.27,$
-1.50, -1.61, -1.49, -1.43, -1.61, -1.53,$
-1.56, -1.43]

dfem = [-0.96, -0.99, -0.94, -1.13, -1.09, -1.27,$ 
-0.73, -0.98, -1.23, -1.23, -1.64, -1.27,$
-1.41, -1.37, -1.27, -1.44, -1.65, -1.67,$
-1.80, -1.81, -1.65, -1.77, -1.69, -1.83,$
-1.67, -1.59]

For k=0,25 do begin
print,name(k),dh(k),dhp(k), dhm(k), dfe(k),dfep(k),dfem(k)
ENDFOR

!xtitle = '!6Log D(Fe)'
!ytitle = '!6D/H (ppm) Inside and Beyond the Local Bubble'
!ymin=0.
!ymax=30.
!xmin=0.
!xmax=-2.0

A=linfit(dfe, dh, measure_errors=dhp-dh)
print,A
x=fltarr(2)
y=x
x(0)=0.
y(0)=A(0)
x(1)=-2.
y(1)=A(0) - A(1)*2.
plot,x,y

For k=0,25 do begin
oplot,[dfe(k)],[dh(k)],psym=sym(k)
plots,[dfe(k),dfe(k)],[dhm(k),dhp(k)]
plots,[dfem(k),dfep(k)],[dh(k),dh(k)]
ENDFOR

device,/close
end