open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

I need to access the TA and TRN (tracer concentrations) at the time when the source file is read and interpolated on line 434 of telemac3d.f. The reason being is that I'm modifying the code to allow for re-circulation between an intake and outfall whereby I want to modify the tracer concentration (TA_SRC) read from the sources file with the modelled value at the intake location.

The issue is that these variables are split into various subdomains when running in parallel and this section of the code does not appear to be looping over the subdomains. I understand how to convert a global to local subdomain number and check whether the current subdomain has the global node of interest. The issue is that this section of the code is not looping over the subdomains.

Does anybody know how I might be able to get around this? Perhaps its better to modify the TA_SRC value within another area of the source?

Thanks,
Toby

So I managed to resolve the issue and was initially mistaken that the telemac3d.f source was not being looped over each subdomain. Once I took a look into the log files for the other processor cores I could indeed see that the local node was correctly detected and the tracer value read.

I think the thread can be closed for now, but perhaps if the community is interested in incorporating functionality to link intake/outfall sources together for re-circulation studies I'd be happy to share my modifications. I'm sure these could be implemented in a more elegant way.