Update diffusion_solvers.f90

added extra output line