Hi all!
The cause of the problem is that there is a big difference in density law when I modify it (see the attachment). In my opinion, the difference in density is too large. I am very grateful to Chi-Tuan for giving me an idea to modify drsurr.f, but I still have too many problems that I do not know how to solve.
According to THE subroutine drsurr,DELTAR=(Rho-Rho0)/RHO0, where RHO is WATER DENSITY and RHO0 is AVERAGE WATER DENSITY IN THE DOMAIN, I know that RHO is a pointer variable after understanding. It doesn't seem to be directly defined and output in a file, so my first question is is RHO an array variable? Secondly, the project I want to use is DENLAW=3. When I check its calculation formula, I find that the reference tracer parameter is Y=TA%ADR(IND_T)%P. TA%ADR(1)%P%R(J) is more like an element of TA%ADR(IND_S)%P than TA%ADR(1)%P%R(J) in the subroutine USER_CONDI3D_TRAC, so is the RHO calculated as a node or as a PLANE or in other ways? What should I do if I want to calculate the density of nodes?
The last thing I want to conclude is, I now have a temperature and salinity and sediment concentration that is changing in time and space, how do I get a density field that is changing in time? Can this be done with telemac software?
Look forward to hearing from you all!
Xinjie Zhang