@common_multi
default,"MgI_marcs_sun_NETL"
mulrd
kr=18

plot,taulg,alog10(n(1,*)/nstar(1,*)),xtitle='log !4s!3!d5000',ytitle='b!di!n'

set_plot,'ps'
!P.thick=3
!P.charsize=2
!P.charthick=3
angstrom = "!3"+STRING(197B)+"!X"
delta    = "!4D!3";"!4"+String(144B)+"!X"
tek_color

device,filename="bi.eps",/color,/encap
plot,taulg,alog10(n(1,*)/nstar(1,*)),xtitle='log !4s!3!d5000',ytitle='log b!di!n'
device,/close

set_plot,'x'