The thermal environment in a small city rapidly deteriorates due to the urbanization and overpopulation. It is important to understand and predict the thermal environment in a city area. The thermal environment is highly affected by the solar radiation and temperature distributions changing over time periodically. To predict the thermal environment precisely, the solar radiation calculation including radiation strength, incidence angle, and thermal radiation between building surface and ground should be considered. In this study, the computational domain includes various artificial structures such as building, ground, asphalt, brick and grass. To consider the solar radiation, the unsteady state numerical calculation is performed from sun rise to mid-day (2:00pm). The numerical methods consist of solar load and one dimensional heat conduction through the boundaries to reduce the computational load and improve the flexibility of the calculation.