Update diffusionsolver.f90