id	summary	reporter	owner	description	type	status	priority	milestone	component	version	resolution	keywords	cc
649	MPDATA tiling bug	jcwarner	arango	"'''TS_MPDATA''' developed a tiling bug in the '''mpdata_adiff.F''' routine when all the ranges were recoded to be more general.  There were only a few places, but it was critical. I found the offending lines and provided a routine (attached) with correct ranges.

Also, i added a new flag '''TS_MPDATA_LIMIT'''. When this option is activated, the upwind corrector fluxes are further limited to be a factor of '''fac = 0.25''' times the computed corrector flux. We have found that in some limited applications with rapidly varying topography and strong tracer gradients, that the upwind flux was creating some oscillations. This limits that behavior.  Of course, if you set the '''fac = 1.0''', then the method is purely upwind.


"	bug	closed	major	Release ROMS/TOMS 3.7	Nonlinear	3.7	Fixed