add funtion to compute superficial density