A modeling of the Black Sea circulation for 2016 was carried out with different sets of atmospheric forcing data to determine the optimal atmospheric forcing for retrospective analysis of hydrophysical fields. An eddy-resolving z-model of Marine Hydrophysical Institute with a resolution of 1.6 km was used for the calculations. Differences in the circulation structure for the two experiments were revealed. It was shown that in the SKIRON experiment compared to ERA5, the cyclonic circulation of the Black Sea was weakened, the isopycnic surfaces were aligned, and the cold intermediate layer was not determined by the 8°C isotherm due to the underestimation of the solar radiation flux and weakening of the wind influence. A comparison of the model thermohaline characteristics calculated using ERA5 and SKIRON atmospheric forcings and measurement data of temperature and salinity obtained by ARGO profiling floats and onboard equipment in 87, 89, 91 cruises of R/V Professor Vodyanitsky was carried out. According to the validation results, it was obtained that in the upper 300-meter layer, for all measurement stations the mean RMSE of temperature and salinity in the ERA5 experiment were 28 and 17% lower, respectively, than the RMSE calculated from the SKIRON data.