![Page 1: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/1.jpg)
1
Reformulation of the LM fast-waves equation part
including a radiative upper boundary condition
Almut Gassmann and Hans-Joachim Herzog
(Meteorological Institute of Bonn University, DWD Potsdam)
![Page 2: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/2.jpg)
2
Time integration is divided into 2 parts
1. Fast waves (gravity and sound waves)
2. Slow tendencies (including advection)
Short review…
n+1n*
F(n*)
nF(n)
Fast waves and slow tendencies
improper mode separation improper combination in
case of the Runge-Kutta-variants
Further numerical shortcomings
divergence damping vertical implicit weights:
symmetry: buoyancy term <-> other implicit terms
lower boundary condition
![Page 3: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/3.jpg)
3
2-time-level scheme (KW-RK2-short)
Comments on the Murthy-Nanundiah-test (Baldauf 2004)
The test relies on the equation
whose stationary solution is known as g . Baldauf claimed that the KW-RK2-short-scheme was not suitable since the splitting into fast forcing and slow relaxation does not yield the correct stationary solution. BUT: The splitting into slow forcing and fast relaxation does.What is our approch alike? -> slow forcing and fast relaxation Forcing: physical and advective processes, nonlinear onesRelaxation: wave processes, linear ones
gdtd
![Page 4: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/4.jpg)
4
2-dimensional linear analysis of the fast wave part
Vertical advection of background pressure and temperature
These terms are essential for wave propagation and energy
consistency
•Which is the state to linearize around?
•LM basic state (current)
or•State at timestep „n“, slow mode backgroundBrunt-Vaisala-Frequency for
the isothermal atmosphere
scale height
variables are scaled to get rid of the density
![Page 5: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/5.jpg)
5
Time scheme for fast waveshorizontal explicit – vertical implicit
divergence damping
symmetric implicitness(treatment as in other
implicit terms)vertical temperature advection
Remark: Acoustic and gravity waves are not neatly separable!
![Page 6: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/6.jpg)
6
Divergence dampingRelative phase change
Phase speeds of gravity waves are distorted.
With divergence damping Without divergence damping
![Page 7: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/7.jpg)
7
Symmetric implicitnessAmplification factor
unsymmetric
symmetric
![Page 8: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/8.jpg)
8
Vertical advection of temperatureRelative phase change
Phase speeds are incorrect. The impact in forecasts can hardly be estimated.
Without T-advection (nonisothermal atmosphere) With T-advection
![Page 9: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/9.jpg)
9
Conclusions from linear analysis No divergence damping! Symmetric implicit formulation! Vertical temperature advection
belongs to fast waves as well as vertical pressure advection! Further conclusion: state to linearize
around is state at time step „n“ and not the LM base state!
![Page 10: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/10.jpg)
10
Appropriate splitting
with
slow tendencies fast wavesin vertical advection for
perturbation pressure or temperature
![Page 11: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/11.jpg)
11
•prescribe Neumann boundary conditions
with access to
which is also used to derive surface pressure
„fast“ LBC•prescribe w(ke1) via
•prescribe metrical term
in momentum equation via
(Almut Gassmann, COSMONewsletter 4, 2004, 155-158)
Lower boundary conditionfast waves
„slow and fast“ LBC
In that way we avoid any computational boundary condition.
slow tendencies
„slow“ LBC
•prescribe
via
out of fast waves
![Page 12: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/12.jpg)
12
• Gassmann, Meteorol Atmos Phys (2004):“An improved two-time-level split-explicit integration scheme for non-hydrostatic compressible models“ • Crank-Nicolson-method is used for vertical advection.• Runge-Kutta-method RK3/2 is used for horizontal advection only and should not be mixed up with the fast waves part. • Splitting errors of mixed methods (Wicker-Skamarock-type) are larger.
Splitting slow modes and fast waves
Gain of efficiency: • No mixing of slow-tendency computation with fast waves• No mixing of vertical advection with Runge-Kutta-steps
![Page 13: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/13.jpg)
K42
K5.33
K0
Background profile
![Page 14: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/14.jpg)
14
Mountain wave with RUBC
w-field
isothermalbackgroundand base state
![Page 15: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/15.jpg)
15
Strong sensitvity of surface pressure at the lee side of the Alps, if different formulations of metric terms in the wind divergence are used
Conservation form (not used in the default LM version), Direct control over in- and outflow across the edges
Alternative representation (used in the default LM version)
Divergence and metric terms
GvGa
uGaG
D )cos()cos(
1)cos(
11
w
Gv
Ga
JuGa
Jva
ua
D 1)cos(
)cos()cos(
1)cos(
1
p
u
u
G
![Page 16: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/16.jpg)
16
Southerly flow over the Alps12UTC, 3. April 2005, Analysis
![Page 17: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/17.jpg)
17
Pressure problem at the lee side of the Alps – Reference LM
![Page 18: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/18.jpg)
18
Pressure problem at the lee side of the Alps – 2tls ALM
![Page 19: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/19.jpg)
19
Cross section: pressure problem
![Page 20: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/20.jpg)
20
Significance?
![Page 21: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/21.jpg)
21
Potential temperature
![Page 22: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/22.jpg)
22
Potential temperature
![Page 23: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/23.jpg)
23
Northerly wind over the Alps
![Page 24: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/24.jpg)
24
Northerly wind over the Alps
![Page 25: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/25.jpg)
25
Moisture profiles in Lindenberg with different LM-Versions (Thanks to Gerd Vogel)
7-day mean with oper. LM version and new version, dx=2.8 km
![Page 26: 1 Reformulation of the LM fast- waves equation part including a radiative upper boundary condition Almut Gassmann and Hans-Joachim Herzog (Meteorological](https://reader033.vdocument.in/reader033/viewer/2022052706/5a4d1b2b7f8b9ab0599991db/html5/thumbnails/26.jpg)
28
Conclusions and plans Conclusions
The presented split-explicit algorithm is fully consistent and proven by linear analysis.
It needs no artificial assumptions and thus overcomes intuitive ad hoc methods.
Divergence formulation in terrain following is a very crucial point.
Plans Higher order advection and completetion for more
prognostic variables Further realistic testing Comparison with Lindenberg data