[MITgcm-support] Question about internal tide energy flux balance

[Shiqiang Hu (ENV - Postgraduate Researcher)]

Hi MITgcm commuity,

Hope you all well.

I built an idealized 3D internal tide model based on internal_wave experiment<https://github.com/MITgcm/MITgcm/blob/master/verification/internal_wave> to calculate the energy flux F=<u'p'>(Nash et al., 2005<https://doi.org/10.1175/JTECH1784.1>). Periodic 1-mode internal tides propagate to the right and interact with a simple topography. The model works but I got stuck at the energy flux balance part. Sometimes when I chose a small area within the topography to calculate the net energy flux across its boundary, the result is positive, which means there are rather dissipation but genaration of enerny. So I wanna ask if someone has experience with this.
Basically, I use RHOAnoma+rhoConst to get original density rho, cause the paper gets density anomaly from rho(z,t)-rho_mean(z), where rho_mean(z) is the time mean vertical density profile (I guess I could also just use the Rhoref(z) as the reference profile?). The question is the C-grids, should I calculate these values on T-points or U-V points so the energy could be more balenced? Or maybe some other factors I didn't realize?
By the way, the non-hydrostatic pressure is quite small, so I ignored it.

Thanks for any suggestions.
Kind regards,

Sean
University of East Anglia
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20250718/6e4a46fc/attachment.html>