The mesoscale numerical weather prediction system (Advanced Research WRF /ARW/ Version 3.3.1) was initialized with US National Center for Environmental Prediction Final Analyses (FNL) of 1° × 1° spatial and 6 h temporal resolution. Two-way nesting was applied on three model domains of grid size 25, 5, and 1 km and number of grid points 26 × 21, 36 × 36, 111 × 111, respectively. A Lambert projection with truelat1 = 30°, truelat2 = 60°, and coordinates of the center of the outermost domain 42.084°N, 27.951°E was used. 24 categories for land surface according to USGS (US Geological Survey) were taken.

The choice of the physical processes’ parameterizations is listed in Table. The upper boundary of the model domains was set at 50 hPa and 43 vertical levels were defined, 30 of which below 2000 m.

For each case, the first 12 h were used for model spin-up. The model data were extracted with 10 min temporal resolution to correspond to the time resolution of the sodar data. The ability of WRF model to represent the observed evolution of the vertical structure of the wind field at a coastal site for one of the most challenging situations of high wind speed and sea breeze development (17-18 September 2016) was tested.