antoinemarie wrote:Moreover, I have the same problem with the temperature of my sink.
Of course you do. All sources without tracer values are using an unconditionally unstable downwind advection scheme. Ditto sinks with specified tracer values.
I have both source and sink, how can I hack the code?
First of all, search for the point source code in ROMS and figure out what's doing what.
1. I would add a logical array, dimensioned by the number of sources. Add the code to read it from ana_psource.h and/or from the sources file.
2. Use it per source instead of using LtracerSrc.
Note that you'll need to provide tracer values for all source/sinks for all tracers, but you'll only be using some of them.