Using numerical models for the purpose of channel-routing calculation has been well accepted in engineering practice. However, most traditional models fail to predict the transcritical flows because of numerical instability. This paper presents two high-r