The evaluation of thermal bioclimate can be conducted employing either observational or modeling techniques. The advantage of the numerical modeling approach lies in that it can be applied in areas where there is lack of observational data, providing a detailed insight on the prevailing thermal bioclimatic conditions. However, this approach should be exploited carefully since model simulations can be frequently biased. The aim of this paper is to examine the suitability of a mesoscale atmospheric model in terms of evaluating thermal bioclimate. For this, the numerical weather prediction Weather Research and Forecasting (WRF) model and the radiation RayMan model are employed for simulating thermal bioclimatic conditions in Greece during a 1-year time period. The physiologically equivalent temperature (PET) is selected as an index for evaluating thermal bioclimate, while synoptic weather station data are exploited for verifying model performance. The results of the present study shed light on the strengths and weaknesses of the numerical modeling approach. Overall, it is shown that model simulations can provide a useful alternative tool for studying thermal bioclimate. Specifically for Greece, the WRF/RayMan modeling system was found to perform adequately well in reproducing the spatial and temporal variations of PET.