@common_multi
default,"MgI_marcs_sun_ETL"
mulrd
cntrb
kr=18
freq=500

plot,taulg,cntrbf(*,freq,kr),xtitle='log !4s!3!d5000',ytitle='dF/d(log !4s!3!d5000!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="cntrb.eps",/color,/encap
plot,taulg,cntrbf(*,freq,kr),xtitle='log !4s!3!d5000',ytitle='dF/d(log !4s!3!d5000!n)'
device,/close

set_plot,'x'