D. 人工的な散逸項

計算の安定性のために, 以下に示す人工的な散逸項を加える.

音波モードを安定して時間積分するために, 音波減衰項を運動方程式 (A.4), (A.5) に加える ^[16].

但し,

ここで Δ τ は短い時間ステップである. 減衰係数は &alpha = 5 × 10^-9 とする.

移流項を中心差分で離散化したことによって生じる数値ノイズを抑制するため, すべての予報方程式に数値粘性項を追加する.

ここで φ は各予報変数であり,

ここで Δ t は長い時間ステップである. 減衰係数は α_h = α_v = 1 × 10^-4 とする.

上部境界付近での波の鉛直伝播と波の反射を抑制するため, 運動方程式 (A.4), (A.5) にレイリー摩擦を加える.

ここで γ は以下のように与える.

z_b = 2.0 × 10⁶ m, z_t = 3.0 × 10⁶ m, α_γ = 1/500 s^-1 とする.