Rather than solving a linear system with the full diffusion matrix (implicit treatment of horizontal and vertical diffusion), only a tridiagonal solver in the SET_DIF subroutine is used leading to a substantial speedup, in particular when advecting a large number of tracers. No horizontal diffusion is considered.
Keywords: SCHEME FOR DIFFUSION OF TRACERS (in TELEMAC-3D) and SCHEME FOR DIFFUSION OF SUSPENDED SEDIMENTS IN 3D (in GAIA) with 1 value per tracer/sediment (new since release 9.1)