The climatological dataset was produced using the Weather and Research Forecasting (WRF) model, version 4.2.2, configured with two nested domains at 10 km (D1) and 2 km (D2) horizontal grid spacing. It covers most of the South Island of New Zealand and is centered over Brewster Glacier in the Southern Alps. The model was forced every three hours by ERA5 reanalysis data at its outer lateral boundaries. The dataset spans the period of 1 January 2005 to 31 December 2020, providing daily output in the outer domain (D1) and 3-hourly output in the innermost domain (D2).

                The data provided here are a selection of daily averages from the inner WRF domain (D2; 2-km grid spacing). They are distributed among three different file types containing 4-dimensional, 3-dimensional and time-invariant output variables, respectively. For the 4-dimensional fields, perturbation and base-state atmospheric pressure (WRF variables P and PB) and geopotential (PH and PHB) were combined to produce full model fields (PRES and GEOPT). Perturbation potential temperature (T) was converted to total potential temperature (THETA). Wind vectors (U,V, and W) were converted to mass points and rotated to earth coordinates.

                -------
                Acknowledgements: The modeling and related research was supported by the German Research Foundation (DFG) grant no. 453305163. The authors gratefully acknowledge the scientific support and HPC resources provided by the Erlangen National High Performance Computing Center (NHR@FAU) of the Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) under the NHR project b128dc / ATMOS ("Numerical atmospheric modeling for the attribution of climate change and for model improvement"). NHR funding is provided by federal and Bavarian state authorities. NHR@FAU hardware is partially funded by the German Research Foundation (DFG) – 440719683.