## 5 Conclusions  
In this study, numerical simulation was used to identify the patterns of surface wind pressure and wind load shape factors for five sizes of solar greenhouses at 0\({}^{\circ}\) and 180\({}^{\circ}\) wind angles. Various zoning methods for the front roof of solar greenhouses were discussed, and an improved zoning method was proposed, accounting for the effect of greenhouse size on the wind load shape factor of the front roof. The main conclusions are as follows:  
(1) At a coming wind angle of 0\({}^{\circ}\), with spans growing from 8 m to 16 m, the value of the positive wind pressure coefficient of the north wall increases with the increase in size. The shape factors of wind load within the B3-B4 region of the north wall also display an increase as the size length grows. Simultaneously, the region of negative wind pressure coefficients on the front roof, with an absolute value less than 0.4, diminishes as size length increases. Conversely, the region, with an absolute between 0.4 and 0.6, expands with increasing size. Separately, the region, with an absolute value exceeding 0.6, is almost unaffected by changes in size length. (2) At a coming wind angle of 180\({}^{\circ}\), the size exhibits minimal influence on the distribution of wind pressure coefficients along the north wall and the rear roof. Particularly, the width of the positive wind pressure coefficient region on the front roof expands proportionally with the size length. In contrast, the shape factor within the central region (H3) of the rear roof remains relatively constant, irrespective of size length. However, the regions at both ends (H1, H2, H4, and H5) experience a reduction in shape factor as the size increases, with a decrease of approximately 22%, when comparing the 16 m span to the 8 m span. (3) Regarding the shape factors of wind load, at a coming wind angle of 0\({}^{\circ}\), the absolute value of the localized shape factors along the 1.6 m width of the front roof exhibits a substantial increase, ranging from 28% to 41%, compared to the shape factors of the upper portion of the front roof determined using Method 1. (4) The proposed zoning approach for the front roof involves the following criteria: at a coming wind angle of 0\({}^{\circ}\), Method 1 is applied, dividing the entire arc length of the front roof into two equal segments; at a coming wind angle of 180\({}^{\circ}\), Method 3 is employed, where the boundary for the front roof's wind pressure coefficient was set at 0. In summary, variations in the dimensions of the greenhouse have a discernible effect on its surface wind pressure and shape factor. As size increases, the wind load per unit area shows a consistent upward trend. Wind-induced suction is particularly concentrated along the roof ridge, with greater intensity observed on the leeward side in comparison to the windward side. Furthermore, this suction effect becomes more pronounced closer to the roof ridge. Special attention should be paid to the wind resistance design of greenhouses in areas with significant wind suction and wind pressure. This paper focuses on two wind directions, highlighting the necessity to investigate the effects of additional wind directions and dimensional variations on surface wind pressures and the shape factor of wind load in groups of solar greenhouses. Conceptualization, Z.L. and J.X. ; software, Z.G. ; investigation, Y.L. ; writing--original draft preparation, Z.G. ; writing--review and editing, Z.L. and J.X. ; visualization, Z.G. and R.W. ; supervision, S.Z. All authors have read and agreed to the published version of the manuscript. This research was funded by the National Natural Science Foundation of China (U20A2020), the Beijing Innovation Consortium of Agriculture Research System (BAIC01-2023-20), the Cangnan