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.