Hello Ziyou,
Looking at your steering file, I have a few remarks:
- please use GAIA rather than sediment features from TELEMAC-3D, it is now the official module to deal with sediment transport both in 2D and 3D,
- I would not use an advection scheme for the advection of sediment which is not mass-conservative as the characteristics you use
(SCHEME FOR ADVECTION OF TRACERS = 1).
Try the LIPS:
SCHEME FOR ADVECTION OF TRACERS = 5 + SCHEME OPTION FOR ADVECTION OF TRACERS = 4
or NERD:
SCHEME FOR ADVECTION OF TRACERS = 13 + SCHEME OPTION FOR ADVECTION OF TRACERS = 1
if you have dry areas in your domain (otherwise you can try PSI or N schemes).
Anyway, in your implementation of TRISOU, you can see for other parts of this subroutine that there is an operation to activate the treatment of the computation of source terms (see type TYPR = 'Q' to be activated, not '0' which means nothing). Your adding would be OK if before CV1 and CV2 are already computed in this subroutine, otherwise it should be first initialised at 0.
Hope this helps,
Chi-Tuan