Limiter in gamma scheme
 

Hi everyone,

I wrote a little post-processing tool to evaluate in which areas of my domain the upwind scheme is used for interpolation when using the gamma scheme.
The limiter / blending factor in the gamma scheme is calculated by

For a phict larger than k_, the limiter is 1 and the scheme is 100% central differencing.
However, I don't understand what happens if phict is larger than +1. The scheme should then be 100% upwind, as stated in the gamma-paper by Prof. Jasak. However, the limiter still is 1 instead of 0. What am I missing here?

Many thanks in advance,
Oliver