Fixed a few compaction problems: Deallocating certain osolve fields after compaction; Storing densityp=density in initialize_compaction; Several small bugs related to calculating the compaction velocity