I guess there is a bug in the diagnostics module in step3d_t.F. It forget to add horizontal advection term in diaTwrk when TS_MPDATA is defined.
After line 715, it may need to add:
# ifdef DIAGNOSTICS_TS
DiaTwrk(i,j,k,itrc,iThadv)=-cff1
# endif
Please take a check.
Search found 1 match
- Mon Mar 09, 2009 4:23 pm
- Forum: ROMS Problems
- Topic: passive tracer not conserved with TS_MPDATA in shallow water
- Replies: 5
- Views: 4721