---Added density output option 'bes'