Rocks are one of the most abundant and economical building materials available to humans, which are used in many construction projects. The flow in the coarse porous media becomes turbulent due to their large particles and pores size. As a results of this fact Darcy law is not valid in such media, and alternatively should utilize a set of relations called non-Darcy relations for the investigation of the hydraulic flow in porous medium. One of the hydraulic studies in the coarse porous media is the determining the profile of the surface water inside this media, which is important in designing rockfill structures. In general, non-Darcy relations are presented in forms of either power or quadratic between Reynolds number and friction coefficient or velocity and hydraulic gradient. These relations have coefficients that are depend on physical properties of rock and fluid. Hence, any relation is achieved based on a certain range of particle diameter. In the present study, four coarse porous media with average diameters of 1.16, 1.83, 3.85 and 6.4 cm and dimensions of 50×40×80 cm in a laboratory flume have been tested and the power relations coefficients were obtained using laboratory data and optimization method for the mentioned media. By modifying the gradually varied flow equation, which is used to determine the flow profile in the open channels, an equation was derived to determine the water surface profile inside the coarse porous media. Hydraulic characteristics of the flow in porous media were analyzed using the calibrated relations obtained from the experimental media. The results of this study showed the inverse relationship between the Reynolds number and the coefficient of friction and the direct and nonlinear relation between the velocity and hydraulic gradient flow in these media. Also, the study of the changes of the hydraulic characteristics with the diameter of the gravel particles shows that the velocity, the Reynolds and the Froude number are directly related to the diameter of the rocks, but the hydraulic gradient decreases with increase of the materials diameter. The study of laboratory profiles showed that the finer the media and the lower the flow rate, the fluctuations of the surface water profiles in the coarse porous media became lower. The comparison results between the laboratory profiles with the computational profiles derived from the numerical solution of the modified gradually varied flow equation showed a very low difference between the experimental and the computational values, in such a way that the NRMSE error rate for all flows was less than 10%. Keywords : Coarse porous media, Non-Darcy relations, Water surface profile, Gradually varied flow