Skip to content

wind

This app creates a meteorological data file with synthetic wind fields suitable for testing advection using the method of Williamson et al. (1992, Sect. 3). The synthetic wind field describes a global rotation around the Earth with a given rotation axis. The output is written in form of MPTRAC's netCDF meteo file format.

The generated file now also contains a configurable set of auxiliary meteorological fields (PS, BLH, ISHF, IEWS, INSS, Q, etc.) written through MPTRAC's standard write_met_nc() path, so the output can be used directly in idealized diffusion and PBL test cases without an additional post-processing step. WIND_BLH is converted internally to boundary-layer top pressure, while WIND_Q and WIND_O3 are still specified as mass mixing ratios and converted internally to the representation used by MPTRAC's meteo writer.

  • Williamson, D. L., Drake, J. B., Hack, J. J., Jakob, R., & Swarztrauber, P. N. (1992). A standard test set for numerical approximations to the shallow water equations in spherical geometry. Journal of computational physics, 102(1), 211-224.
# Calling wind
$ ./wind <ctl> <metbase>

The required parameters are: * ctl: the control parameter file * metbase: the basename of the output file

Optional control parameters: * WIND_T0: time step of meteo data (default: 0 s) * WIND_NX: number of longitudes (default: 360, minimum: 2) * WIND_NY: number of latitudes (default: 181, minimum: 2) * WIND_NZ: number of pressure levels (default: 61, minimum: 2) * WIND_Z0: log-pressure height of lowermost pressure level (default: 0 km) * WIND_Z1: log-pressure height of uppermost pressure level (default: 60 km) * WIND_U0: horizontal wind speed at lowermost level (default: 38.587660177302 m/s) * WIND_U1: horizontal wind speed at uppermost level (default: 38.587660177302 m/s) * WIND_W0: vertical velocity (default: 0 m/s) * WIND_ALPHA: rotation angle (default: 0 deg) * WIND_LAT_REVERSE: latitude direction flag (default: 0, set to 1 for reversed latitude direction) * WIND_TEMP0: temperature at the lowermost level (default: 280 K) * WIND_TEMP1: temperature at the uppermost level (default: 280 K) * WIND_PS: surface pressure (default: 100000 Pa) * WIND_ZS: surface height (default: 0 km) * WIND_T2M: 2 m temperature (default: 280 K) * WIND_IEWS: eastward turbulent surface stress (default: 0 N m^-2) * WIND_INSS: northward turbulent surface stress (default: 0 N m^-2) * WIND_ISHF: surface sensible heat flux in ECMWF convention (positive downward; default: 0 W m^-2) * WIND_LSM: land-sea mask (default: 1) * WIND_SST: sea surface temperature (default: 280 K) * WIND_BLH: boundary layer height (default: 1000 m) * WIND_Q: specific humidity (default: 0 kg kg^-1) * WIND_O3: ozone mass mixing ratio (default: 0 kg kg^-1)