This page describes a numerical simulation of gravity waves over the Rocky
Mountains.

Computational Domain

The layout of the of the simulation is shown in the figure below.
The origin is located at Mount Blue Sky 39.6^{o} N 105.7^{o} W.
Note that the computational domain is rotated 7^{o} clockwise with
respect to lines of constant latitude.

The mesh is clustered in both horizontal directions in order to achieve
250 x 250 meter spacing over the mountain range in the region shown by
the black rectangle. Weak stretching of ~1.5% is used to the edges. The domain
extends to an altitude of 140 km and uses uniform vertical spacing of 250 m.
A total of 768 x 460 x 560 mesh points are used. Characteristic boundary
conditions are used on all sides except the surface.

Wind and Thermodynamic Profiles

The mean winds and temperature profile are taken from radiosonde data on
February 19th, 2016, from launch site in Grand Junction.
These profiles extend to an altitude of 31 km.
MERRA wind profiles are used above the radiosonde measuerments to an altitude of 80 km.
A third order interpolating polynomial is then used to smoothly extend the
winds above this altitude using the condition that U=V=0 at the upper
boundary.
The temperature is extended using the NRLMSIS Atmosphere Model (Composition).
Plots of various profiles are shown below.

Wind Condition

The mean winds near the surface are increased in time.
Forcing terms gradually introduce winds near the surface with the
objective of achieving the wind profile within a two hour period.
A hyperbolic tangent function is used in order to produce gentle
acceleration of the wind near the beginning and end of the forcing period.
The maximum forcing rate is equivalent to that of a linear ramp with a
duration of thirty minutes.

Cropped Movies

Animation of w' in the xz plane at the position y = -20

Animation of w' in the xy plane at altitude of 21 km

Animation of w' in the xy plane at altitude of 110 km

Long Results

Animation of w' in the xz plane at the position y = -20

Animation of w' in the xy plane at altitude of 21 km

Animation of w' in the xy plane at altitude of 110 km