Diurnal variation of the flow over a forest canopy on a mountain slope is simulated numerically. In the daytime, the earth surface is heated by the solar radiation and the flow goes up the mountain due to the buoyancy force, and during the night, the air is drained downward along the slope owing to the cooling of the surface by radiation. In this flow process the forest canopy that consists of leaf region and the trunk region plays a dominant role as a momentum sink to the flow, thus the modeling of the leaf area region and trunk region is critical to the successful flow simulation. In the present study, a field measurement in an experimental forest in the State of Oregon in the United States is numerically analyzed. The resistance to the flow in the leaf region is directly related to the leaf area density (LAD), and the trunk is modeled as a cylinder.