As noticed by Nicole, the recent change of intent(inout) to intent(in) was
not correct, even if our compilers did not noticed it on the build farm.
I changed the code to keep the intent(in), as this routine is doing a dump of the density, and therefore should not be allowed to change it.