Update diffusion_solvers.f90

Added two comments, added time as an argument to solvers.crank and solvers.super (for use in boundary conditions)