Separating horizontal XI-advection term for tracer diagnostics in the source code (in step3d_t.F)

[wavewatch]

Hi all,
There is a horizontal XI-advection variable in the DIAGNOSTIC output file, i.e temp_xadv (d(u*T)/dx = u dT/dx+T du/dx).Since the two properties (u dT/dx and T du/dx) of the tracer concentration are both valuable, the advective form of tracer concentration equation is changed into a conservation form.
But now I want to see what u dT/dx and T du/dx do, how should it be implemented in code step3d_t.F ? (The convection scheme is TS_U3HADVECTION)
Any suggestion is appreciated!! thanks

[fionanxhjh]

Hi, nice to see your question. I have thought of the output of temperature advection, i think the temp_xadv is the Tdu/dx, not the d(Tu)/dx.