Using the daily precipitation data at the 26 hydrological & meteorological stations, daily averaged air temperature data at the 11 stations, and daily pan evaporation data at the 14 stations in 2000, the distributed water-heat coupled model (DWHC) was calibrated. Using the daily averaged runoff data in 2000 at the Yingluoxia station where the runoff amount of Heihe mountainous river basin was measured to calibrate, 6 numerical tests that used different spatial interpolation methods to calculate the daily precipitation, daily averaged air temperature, and daily pan evaporation (E601) in each grid, were put up. Due to the spatial sparsity and asymmetry of the hydrological and meteorological stations, the results of the 6 numerical tests have little differences. The three-dimensional interpolation methods considering the altitudes are little different from the methods that take no account of the altitudes. The daily data measured at the 2 stations farside from the research basin were added for use in the model, and the model results were not better. At last the 2-dimensional interpolation method with a named as nearest, was used in the model. The Nash-Sutcliffe equation value is of 0.6101, and the balance error is of 0.0808% for the calibration processes. Using the daily data in 1999 to validate the model, the Nash-Sutcliffe equation value and balance error is about 0.6270 and 2.9824%, respectively. Because the soil and vegetation data are not detailed and very short, the model results are only qualitative to some extent. The water balance, and the spatial and temporal distributions of the hydrological factors, were both discussed. The model results showed that, in the high and cold regions, the runoff production processes mainly occurred at the land surfaces and in the low soil layers, the alpine meadow held up precipitation to a large extent, and there produced some ground ices due to the heat-water coupled processes in the high mountainous regions.