Correct a small density calculation bug

Previously, disk.f90 evaluate the cell density at the center.  Now evaluate
based on an average over the whole area, to be consistent with the treatment in
grid.f90.