The dataset contains model output files and examples of scripts for depicting figures related to the article "Influence of thermospheric impacts of solar activity on the general circulation and stationary planetary waves in the middle atmosphere" by A.V. Koval, N. M. Gavrilov, A. I. Pogoreltsev, N. O. Shevchuk. Contents: Output data from model simulations averaged over 16 pairs of model runs under high and low solar activity: gh_m1_hsa5.dx, gh_m1_lsa5.dx – amplitudes of the geopotential height variations in g.p.m. at high and low solar activity caused by SPW modes with zonal wavenumber 1 averaged over 80 time subintervals selected from pairs of the MUAM runs (structure described in w1_tp1_hsa.ctl, w1_tp1_lsa.ctl). gh_m2_hsa5.dx, gh_m2_lsa5.dx, gh_m3_hsa5.dx, gh_m3_lsa5.dx, gh_m4_hsa5.dx, gh_m4_lsa5.dx – the same but for zonal wavenumbers 2,3,4. gh_m1_disp_dif5.dx – dispersion of differences in corresponding wave amplitudes between high and low solar activitiy (w1_tp1_disp.ctl). gh_m2_disp_dif5.dx, gh_m3_disp_dif5.dx, gh_m4_disp_dif5.dx – the same but for zonal wavenumbers 2,3,4. index_refr1_1_5_hsa.dx, index_refr1_2_5_hsa.dx, index_refr1_3_5_hsa.dx, index_refr1_4_5_hsa.dx – the zonal-mean quasi-geostrophic complex refractivity index squared (RI2) for PW modes with zonal wavenumbers 1-4 calculated on the basis of 80 time subintervals selected from the MUAM runs (ind_refr1_1.ctl). High SA. index_refr1_1_5_lsa.dx, index_refr1_2_5_lsa.dx, index_refr1_3_5_lsa.dx, index_refr1_4_5_lsa.dx – the same but for the low solar activity uq_vwres1_pv_Jan_1_5_hsa.dx, uq_vwres1_pv_Jan_2_5_hsa.dx, uq_vwres1_pv_Jan_3_5_hsa.dx, uq_vwres1_pv_Jan_4_5_hsa.dx – the Eliassen-Palm flux for PW modes with zonal wavenumbers 1-4 calculated on the basis of 80 time subintervals selected from the MUAM runs (uqep1_vwres_1.ctl). uq_vwres1_pv_Jan_1_5_lsa.dx, uq_vwres1_pv_Jan_2_5_lsa.dx, uq_vwres1_pv_Jan_3_5_lsa.dx, uq_vwres1_pv_Jan_4_5_lsa.dx – the same but for the low solar activity zw_a0_grads5_hsa.dx, tp_a0_grads5_hsa.dx – zonal-mean zonal wind in m/s and temperature in K for December – February averaged over 16-member ensemble of model runs for high solar activity (m0_zw5_hsa.ctl, m0_tp5_hsa.ctl). zw_a0_grads5_lsa.dx, tp_a0_grads5_lsa.dx – the same but for the low solar activity disp_U_dif5_.dx, disp_T_dif5_.dx – dispersion of differences in zonal wind and temperature between high and low solar activitiy for December – February averaged over 16-member ensemble of model runs (tp1_disp.ctl). uq_vwres_pv_Jan_1_ave_hsa.dx, uq_vwres_pv_Jan_1_ave_lsa.dx - The Eliassen-Palm flux and its divergence calculated for December – February averaged over 16-member ensemble of model runs for high solar activity (uqep1_vwres_ave.ctl). Additional data including outputs from separate model runs are available from the authors upon request.