MPDATA tiling bug

[jcwarner]

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 afew 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 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 to = 0.0, then the method is purely upwind.

-john

[jcwarner]

mpdata_adiff corrected ranges