AS I remember the theory, division with H results from depth-averaging.
Correct.
Then the coefficient -0.0001 in the given example is multiplied by Hprop because subroutine CVDFTR divides source term with Hprop.
Correct.
Variable TIMP is not listed in the table of variable names description but is the function argument.
Yes. Unfortunately, the tables of variable names are out of date in many areas of the TELEMAC source code. Without wanting to sound like an advertisement, we're also working on an automated documentation system that ought to avoid this problem in the future.
It is also for the NTRAC that I presume is total number of tracers. Right?
.
Correct.
With regard to the keyword, I agree with Christophe (telemac) that adding a new keyword is a difficult thing to do without a fairly intimate knowledge of the code. To expand on Christophe's point about external files, are you aware of the
FORTRAN FILE keyword? Adding something like
FORTRAN FILE : './princi.f'
to your steering file means that TELEMAC looks in
princi.f and uses any subroutines contained therein to overwrite the standard TELEMAC subroutines. Many of the validation cases use such a file to modify
CONDIN to specify particular initial conditions.
In your case, putting an appropriately modified version of
DIFSOU in
princi.f (or whatever else you call your Fortran file) would mean that TELEMAC needs only to recompile this subroutine, and repeat the linking procedure, each time you change the decay law. The computational cost of this is quite small, and all of the recompilation steps are hidden from view in the telemac2d command.
I hope that's useful, and that I didn't just tell you things you already knew!
