Hello,
As far as I can see in the subroutine 'SOUKEP.f' (TELEMAC 3D v7p2r1), the buoyancy term G in the k-epsilon turbulence model only takes into account the derivative of traceur N°1 with respect to Z, and not the derivative of the density with respect to Z. Is it right? Do I miss something?
In my case, I have two traceurs : salinity and temperature, and the density evolves with both. Does the buoyancy term evolve only with the salinity?
Thanks in advance for your help with this matter.
Best regards,
Sophie