Implementation of a Full P1 Method in the Diffusion Code DONJON/NDF