Supplementary material (part 1) for "Scale interactions of winds in a small Alpine valley: A modeling case study of Nafingalm"
收藏资源简介:
Part 1 of the supplementary material for the master's thesis "Scale interactions of winds in a small Alpine valley: A modeling case study of Nafingalm" (Freundorfer 2026, available at https://resolver.obvsg.at/urn:nbn:at:at-ubi:1-200524). Due to memory constraints, the supplementary material is split up into two parts: part 1 contains scripts, datasets, model setup files and figures (see description below), while part 2 includes the model output to create the figures of the thesis (available at https://zenodo.org/records/19371244). In order to reproduce the plots, the datasets and model output have to be unpacked and the location of the folder manually integrated into the scripts. Part of the scripts, mainly the cross sections, were designed to be run on an HPC, however it should be possible to execute the scripts normally as well. The figures that use observational data and tslist output were created locally in the jupyter notebook. Scripts This folder is divided into 3 types of scripts. All of them were used to process and plot model as well as observational data. Jupyter_Notebook Master_thesis_plotting.ipynb: Jupyter Notebook to plot figures using observational data as well as tslist output. Python_scripts This folder contains Python scripts that were run on the HPC to create the cross sections as well as additional miscellaneous scripts. A desciption can be found in the following. Plotting routines: plot_10m_ts_IBOX.py: script to plot the near-surface comparison at Kolsass (Fig. 3.8) plot_horizontal_cs_tiles.py: script to plot horizontal cross sections from wrfout files (Figs. 3.5, A.4) plot_vertical_cs_tiles.py: script to plot vertical cross sections from wrfout files (Figs. 3.6, 3.13, 3.14, 3.15, A.2, A.3) TKE_plotting.py: script to plot vertical cross sections from tslist output (Figs. 3.16d,e, 3.17d,e, 3.18d,e) windfield_plots.py: script to plot vertical cross sections from tslist output (Fig. 3.12d-l) plot_functions.py: script that contains the actual plotting functions for the scripts above. help_functions.py: script that contains various miscellaneous functions used in the plotting. Data processing: TKE_calculation_2.py: script to calculate the TKE from the tslist output. extract_QKE_from_wrfout.py: script to extract the SGS TKE from the MYNN BL scheme. extract_topo.py: script to extract the topography from the wrfout output of the simulations. smooth_topography.py: script to smooth the topography data of the models. A more detailed description can be found in the note "LES_ndown_notes.txt". check_wrf_slopes.py: script to check the slope angles in the topography data of the simulations. tslist_to_netcdf_optimized.py: script to convert the tslist output to .nc files. create_snow_in_glaciers.py: script to correct the snow field in the met_em files during the WPS preprocessing to avoid instabilities. From Paolo Deidda, more information can be found in his thesis, Chapter 2.2. MATLAB This folder contains MATLAB files to create Fig. 3.12a-c from the lidar data at Nafingalm. Main script: plot_sl88_rhi_user2_3panel.m Dependencies: cmap_bluered16.m, read_lidar_hpl_file.m Datasets This folder contains part of the dataset in order to plot the figures. The model output is located in another repository due to memory limitations. ERA5 This folder contains data from the ERA5 reanalysis. It consists of temperature, geopotential height and windfield data at the 925, 850, 800, 700 and 300 hPa pressure levels. It covers central Europe on a timeframe of 00 UTC 01 September - 12 UTC 02 September. iBox This folder contains AWS data from the iBox located at Kolsass. It features temperature and relative humitidy measurements 8.7 m AGL (taact3_avg and rhact3_avg) as well as measurements of wind speed and direction at 12 m AGL (wind_speed_4, avg_wdir4). Innsbruck This folder contains lidar data from the roof of the university of Innsbruck. It consists of a timeseries of vertical profiles of horizontal wind speeds for 01 and 02 September. Nafingalm This folder contains meansurements at Nafingalm in the Weer Valley from a test campaign previous to the TEAMx observational campaign, that were used in the creation of the figures for the thesis. The data consist of lidar measurements, near-surface temperature observations from an AWS and vertical temperature profiles from UAS measurements from 11 UTC 01 September to 15 UTC 02 September. The lidar data consists of a timeseries of vertical profiles of horizontal wind speeds for the whole measurement period, as well as selected .hpl files for the creation of a vertical cross section of radial wind velocities. The UAS data consists of two netcdf files, one for the vertical temperature structure and one for the drone trajectories. The AWS temperature time series was converted to a .mat file to be loaded into the Jupyter Notebook more easily. SRTM This folder contains SRTM-1 topography data used to compare to the model topographies. SGS TKE from MYNN scheme in mesoscale models This folder contains the subgrid-scale turbulent kinetic energy model data extracted from the wrfout files. WRF setup This folder contains files used for the WRF setup. It consists of: WPS namelist: Setup file for the WRF Preprocessing System (WPS). WRF namelists: Setup files for the WRF model itself. It consists of multiple files, as the thesis features 3 offline nested simulations. bash scripts: Scripts to run WPS and WRF (real.exe, wrf.exe and ndown.exe) on the HPC. tslist: Files to specify the locations of the tslist output. It features an extensive list of locations to verify the model (tslist_meso), a shorter version for relevant locations around the Alps (tslist_short, used in the mesoscale simulations) and an additional file that specifies a "curtain" of gridpoints from d03/D100 to create a vertical cross section (tslist_LES, used in the LES). METGRID.TBL: A modified table to use the ECMWF forecast data as input data for d01/D2500. GEOGRID.TBL: A modified table to use the SRTM-1 topography data as well as the CORINE land use dataset in the simulations. Note that the interpolation of the topography data onto the model grid, a different interpolation method was used for d03/D100 than for d01/D2500 and d02/D500. geo_em files: The model terrain for the three simulations, for d03/D100 an original and smoothed version. wrfinput d03: wrfinput files d03/D100, that were used directly for the simulation. After the usage of ndown.exe as well as an original file and a version with additonally smoothed terrain. The specifics can be found in this thesis, Chapter 2.2, with the python file "smooth_topography.py". WRF base code changes: This folder contains three files of the WRF source code that were changed to include the output of TKE, as well as a bug fix in the diffusion calculation in the file module_diffusion_em.F. The file Registry.EM_COMMON contains changes to the standard wrf history output (wrfout files) to include the SGS TKE from the MYNN scheme in d01/D2500 and d02/D500. The file wrf_timeseries.F contains a change to include the SGS TKE from the Deardorff TKE scheme in d03/D100 to be included in the tslist output. A more detailed description of the changes of the two latter files can be found in the note "TKE_output.txt". notes: Three text files to provide a description of the changes to the source code changes (TKE_output.txt), the nesting of the simulations (LES_ndown_notes.txt) as well as a file to change the order of the tslist output files for the along valley cross section of d03/D100 (V_TKE.zip in part 2 of the repository). Depending on the needs, this should be ideally specified correctly in the tslist before running the simulation. Figures This folder consists of the figures used in the thesis.



