The capability to evaluate city-scale outdoor thermal comfort is crucial for understanding the public's experience of heat stress, especially during heat waves. Nevertheless, there remains a methodological gap in accurately mapping thermal comfort with high spatial and temporal resolutions in complex urban environment and in quantifying the contributions of various urban climate drivers, such as sky view factors and tree canopy fraction. To bridge this gap, we propose an innovative WRF-UCM-SOLWEIG framework that couples a mesoscale model with an urban microclimate model. Specifically, it integrates the Weather Research and Forecasting model coupled with the urban canopy model (WRF-UCM) and the Solar and Longwave Environmental Irradiance Geometry model (SOLWEIG). This framework is then applied to a densely populated tropical city in China, to quantify the impact of urban climate drivers on the Universal Thermal Climate Index (UTCI). The results reveal a significant diurnal variation in the impact of urban morphology. Notably, the most significant drivers affecting daytime UTCI are the impervious surface fraction and the tree canopy fraction, with maximum Pearson's correlation coefficients of 0.80 and -0.70, respectively. These findings suggest that urban heat mitigation strategies should prioritize on reducing impervious surfaces and enhancing shade management, alongside expanding urban forestry measures.