Model a density bump following Bai 2023 prescription

We want to create a density bump in the disk with no planet and no midification of the viscosity