Heat waves are increasing in frequency and also exhibit high spatial variability in its distribution over India. There are limited studies focused on the weather related human thermal comfort over India due to non-availability of high resolution (HR) climate data. Here we develop dynamically downscaled HR (4x4 km) daily climate information for the months of April to June during 2001-2016 using a regional climate model called Weather Research and Forecasting (WRF) Model, which are validated with station observations. The thermal comfort and its spatio-temporal variability over India are quantified in terms of indices like Excessive Heat Factor (EHF), Heat Index (HI), Humidex, Apparent Temperature (AT) and Wet Bulb Globe Temperature (WBGT). The daily surface air temperature and thermal comfort indices of HR WRF model simulations are in good agreement with station observations. The results show that there is an increasing trend in annual heat waves coverage (22240km2/year), annual frequency (0.07 days/year) and average intensity (0.04 °C/year) during 2001-2016. The distributions of indices have spatial and temporal variability. The days with severe discomfort are significantly increasing (99% significance level) over north India and it is quantified with increase of extreme category of indices at the rate of EHF (15.9%), HI (14.9%), Humidex (15.9%), AT (13.4%) and WBGT (13.8%). During heat waves, prolonged exposure or physical activity under sun will led to adverse health impacts and it is mostly observed over northwest and south eastern states. These findings stress the need for developing suitable mitigation strategies for a sustainable ecosystem