Hello Gao,
Can you upload:
- the 3 steering files (T3D, GAIA, WAQTEL),
- your modified user_t3d_trsce.f subroutine,
- SOURCES FILE (with extension .txt e.g. so that it can be read).
Could you also try to modify t3d_trsce (you can do it in your FORTRAN FILE) by moving the
CALL USER_T3D_TRSCE(T3D_TRSCE, TIME, I, ITRAC)
outside the condition
IF(.NOT.OKTRSCE(I,ITRAC).OR.
& T3D_FILES(T3DVEF)%NAME(1:1).EQ.' ') THEN
e.g. moving it 3 lines after after the last ENDIF but before RETURN.
With the current implementation of T3D_TRSCE, USER_T3D_TRSCE is only called if you do not use a SOURCE FILE or if there is no information for the specific tracer in the SOURCE FILE.
If you want to apply your modifications implemented in USER_T3D_TRSCE for tracers at sources read in a SOURCE FILE, it should be modified that way (I think it is a bug).
Chi-Tuan