This paper deals with applying a numerical finite element model of water dynamics to sub-ice flow in a retention reservoir. On the basis of earlier experiences, the AdH model is chosen for this task. The previous usage of this model is outlined along with the physical basics used by this algorithm. The application procedure is then described and finally a few simulation results obtained for different conditions are shown. The results comparison and analysis indicate that the chosen model is capable of performing the required tasks.