Added calculated density field to ssu.R