Fixed mistake in limiting pressure. We were using density_average instead of cell_average. Needs to be tested.