SLOPE_LIMITER = max(0.D0,max(min(1.D0,2.D0*Cr),
_RL FUNCTION SLOPE_LIMITER (Cr)
phi = SLOPE_LIMITER (Cr)
_RL SLOPE_LIMITER