: 3. 力学過程 : dcpam5 支配方程式系とその離散化 : 1. この文書について

2. 座標系・変換公式

2.1 はじめに

ここでは, 座標系および水平格子点, 鉛直レベルの取り方を記す. さらに, 力学過程の時間積分において使用する水平スペクトルを定義し, 格子点値とスペクトルの係数との変換則を記す.

2.2 座標系

座標系は, 水平方向には緯度 $\varphi$ , 経度 $\lambda$ を, 鉛直方向には $\sigma \equiv \frac{p}{p_s}$ をとる. ここでは気圧, は地表面気圧である.

座標の取り方に関する詳細は別紙『支配方程式系の導出に関する参考資料』の『座標系の取り方』を参照せよ.

2.3 水平格子点

水平方向の格子点の位置は, Gauss 緯度 (格子点数個 ^2.1), 等間隔の経度 (同個) である.

Gauss 緯度
Gauss 緯度を次のLegendre函数 $P_J(\sin \varphi)$ の零点 $\varphi_j (j=1,2,3,\cdots,J)$ として定義する. 順番としては, ${\displaystyle \frac{\pi}{2} > \varphi_1 > \varphi_2 > \cdots > \varphi_{J} > -\frac{\pi}{2} }$ とする ^2.2. なお以後, $\sin \varphi = \mu$ と書くことがある.
経度方向の格子点
経度方向の格子点の位置を

$\displaystyle \lambda_i = \frac{2 \pi (i-1)}{I} \ \ (i=1,2,\cdots,I)$ (2.2)

ととる.

2.4 鉛直レベル

Arakawa and Suarez (1983) のスキームを用いる. とり方は以下のとおりである ^2.3. 下の層から上へと層の番号をつける. 整数レベルと半整数レベルを定義する ^2.4. 半整数レベルでの $\sigma$ の値 $\sigma_{k-1/2} \ (k=1,2,\cdots,K)$ を定義する. ただし, レベル $\frac{1}{2}$ は下端( $\sigma=1$ ), レベル $K+\frac{1}{2}$ は上端( $\sigma=0$ )とする. 整数レベルの $\sigma$ の値 $\sigma_k \ (k=1,2,\ldots K)$ は次の式から求める.

$\displaystyle \sigma_k = \left\{ \frac{1}{1+\kappa} \left( \frac{ \sigma^{\kapp... ... +1}_{k+1/2} } { \sigma_{k-1/2} - \sigma_{k+1/2} } \right) \right\}^{1/\kappa}.$

(2.3)

ただし, $\kappa=\frac{R}{C_p}$ である. ここで,

は乾燥空気の気体定数,

は乾燥空気の定圧比熱である ^2.5. また, レベル加重 $\Delta \sigma$ は以下のように定義される.

$\begin{align*}\begin{split}\Delta \sigma_k &\equiv \sigma_{k-1/2} - \sigma_{k+1/... ...K+1/2} &\equiv \sigma_{K} - \sigma_{K+1/2} = \sigma_{K}. \end{split}\end{align*}$

(2.4)

$\begin{picture}(300,150)(50,10) \put(50,20){\line(1,0){220}} \put(50,40){\line... ...){\shortstack{$\sigma=1$}} \put(280,136){\shortstack{$\sigma=0$}} \end{picture}$

2.5 水平スペクトル

ここでは, 力学過程の時間積分での計算において用いるスペクトルを導入し, 格子点での値とスペクトルの係数とのやり取りの公式を示す.

2.5.1 水平スペクトルの基底の導入

格子点上の点で定義された物理量は, 格子点上でのみ値を持つ (以下このことを, 「離散化した」と呼ぶ) 球面調和函数の和の形で表現される. また, 各格子点における物理量の水平微分を評価するために, $(\lambda, \varphi)$ 面で定義された (以下, 「連続系の」と呼ぶ) 球面調和函数系で内挿して得られる関数を用いる. ここではその球面調和函数を導入する. なお, 簡単のために, 連続系の球面調和函数のみを陽に記す. 離散系の球面調和函数は連続系の球面調和函数に格子点の座標を代入したものから構成される.

$(\lambda, \varphi)$ 面において, 球面調和函数 $Y_n^m(\lambda,\varphi)$ は次のように定義される.

$\displaystyle Y_n^m(\lambda,\varphi) \equiv P_n^m(\sin \varphi) \exp(im \lambda),$

(2.5)

ただし,

は $\ 0 \le \vert m\vert \le n$ を満たす整数であり, $P_n^m(\sin \varphi)$ は 2で規格化されたLegendre函数・陪函数

	$\displaystyle P_n^m(\mu)\equiv \sqrt{\frac{(2n+1)(n-\vert m\vert)!}{(n+\vert m\... ...u^2)^{\frac{\vert m\vert}{2}} }{2^n n!} \DD[n+\vert m\vert]{}{\mu} (\mu^2-1)^n,$	(2.6)
	$\displaystyle \int_{-1}^1 P_n^m(\mu) P_{n'}^m(\mu) d \mu = 2 \delta_{nn'}$	(2.7)

である. なお,

を

とも書く. また $\sin \varphi = \mu$ であることを再掲しておく.

2.5.2 波数切断

波数切断は三角形切断 (T) または平行四辺形切断 (R) とする. , は三角形切断, 平行四辺形切断のときについてそれぞれ以下のとおりである. ただし, 切断波数を $N_{tr}$ とする.

三角形切断の場合
$M= N_{tr}, \ N =N_{tr}$ , $I \ge 3N_{tr}+1$ , かつ $J \ge \frac{3N_{tr}+1}{2}$ .
自由度は, $(N_{tr}+1)^2$ である.
平行四辺形切断の場合
$M= N_{tr}, \ N(m) =N_{tr}+\vert m\vert$ , $I \ge 3N_{tr}+1$ , かつ $J \ge 3N_{tr}+1$ .
自由度は, $(2N_{tr}+1) (N_{tr}+1)$ である.

よく用いられる値の例としては, T42 の場合 $I=128,\ J=64$ , R21 の場合 $I=64,\ J=64$ がある.

球面調和函数と波数切断に関する詳細は, 第A.1節および第A.8節を参照せよ.

2.5.3 離散化したスペクトルの基底の直交性

離散化したLegendre函数と三角関数は次の直交条件を満たす ^2.6.

$\displaystyle \sum_{j=1}^{J} P_n^m (\mu_j) P_{n'}^m (\mu_j) w_j$	$\displaystyle = \delta_{nn'},$	(2.8)
$\displaystyle \sum_{i=1}^{I} \exp(im \lambda_i) \exp(-im' \lambda_i)$	$\displaystyle = I \delta_{mm'}.$	(2.9)

ここで

は Gauss 荷重で, ${\displaystyle w_j \equiv \frac{(2J-1)(1-\sin^2 \varphi_j)} {\left\{J P_{J-1}(\sin \varphi_j)\right\}^2 } }$ である.

2.5.4 格子点値とスペクトルの係数との変換法

物理量の格子点 $(\lambda_i,\varphi_j)$ (ただし $i=1,2,\cdots,I. \quad j=1,2,\cdots,J$ ) での値 $A_{ij}=A(\lambda_i,\varphi_j)$ とスペクトル空間での (ただし $m=-M,\cdots,M. \quad n=\vert m\vert,\cdots,N(m)$ ) の係数 $\tilde{A}_n^m$ とは次の変換則に従う ^2.7.

$\displaystyle A_{ij}$	$\displaystyle \equiv \sum_{m=-M}^{M} \sum_{n=\vert m\vert}^{N} \tilde{A}_n^m Y_n^m (\lambda_i,\varphi_j),$	(2.10)
$\displaystyle \tilde{A}_n^m$	$\displaystyle = \frac{1}{I} \sum_{i=1}^{I} \sum_{j=1}^{J} A_{ij} Y_n^{m*} (\lambda_i, \varphi_j) w_j .$	(2.11)

が実数であることを用いると, ${\displaystyle \left\{\tilde{A}^m_n \exp(im\lambda) \right\}^* = \tilde{A}^{-m}_n \exp(-im\lambda) }$ なので,

については負でない整数の範囲で和をとることができる ^2.8. ここで, ``

'' は複素共役を表す. ただし,

の定義を以下のように修正していることに注意せよ.

$\displaystyle A_{ij}$	$\displaystyle = \sum_{m=0}^{M} \sum_{n=m}^{N} \Re \tilde{A}_n^m Y_n^m(\lambda_i, \varphi_j),$	(2.12)
$\displaystyle \tilde{A}_n^m$	$\displaystyle = \left\{ \begin{array}{ll} {\displaystyle \frac{1}{I} \sum_{i=1}... ...varphi_j) w_j , & \ \ \ 1 \le m \le M, \quad m \le n \le N. \end{array} \right.$	(2.13)

2.5.5 内挿公式

$(\lambda, \varphi)$ 空間で定義される物理量 $A(\lambda,\varphi)$ を格子点値 $A_{ij}$ をもとに内挿する場合には, 変換公式を用いて $A_{ij}$ から $\tilde{A}_n^m$ を求めた上で,

$\displaystyle A(\lambda,\varphi)$

$\displaystyle \equiv \sum_{m=-M}^{M} \sum_{n=\vert m\vert}^{N} \tilde{A}_n^m Y_n^m (\lambda, \varphi)$

(2.14)

として得る.

2.5.6 空間微分の評価

各格子点における空間微分値の評価は, 内挿公式を用いて得た連続関数の空間微分の格子点値で評価する.

$\lambda$ 微分

$\displaystyle \left( \DP{f}{\lambda} \right)_{ij}$ $\displaystyle \equiv \sum_{m=-M}^M \sum_{n=\vert m\vert}^N im \tilde{f}_n^m Y_n^m (\lambda_i,\varphi_j) ,$ (2.15)

$\displaystyle \widetilde{ \left( \DP{f}{\lambda} \right) }_n^m$ $\displaystyle = \frac{1}{I} \sum_{i=1}^I \sum_{j=1}^J im f_{ij} Y_n^{m*} (\lambda_i, \varphi_j) w_j .$ (2.16)
$\varphi$ 微分

$\displaystyle \left(\DP{f}{\varphi}\right)_{ij}$ $\displaystyle \equiv \sum_{m=-M}^M \sum_{n=\vert m\vert}^N \tilde{f}_n^m \left. \DD{P_n^m}{\varphi} \right\vert _j \exp(im \lambda_i) ,$ (2.17)

$\displaystyle \widetilde{\left(\DP{f}{\varphi}\right)}_n^m$ $\displaystyle = - \frac{1}{I} \sum_{i=1}^I \sum_{j=1}^J f_{ij} \left. \DD{P_n^m}{\varphi}\right\vert _j \exp(-im \lambda_i) w_j .$ (2.18)

2.6 参考文献

: Arakawa, A., Suarez, M. J., 1983: Vertical differencing of the primitive equations in sigma coordinates. Mon. Wea. Rev., 111, 34-35.
: 気象庁予報部, 1982 : スペクトル法による数値予報（その原理と実際）. 気象庁, 111pp.
: Haltiner, G.J., Williams, R.T., 1980: Numerical Prediction and Dynamic Meteorology (2nd ed.). John Wiley & Sons, 477pp.
: 森口, 宇田川, 一松編 ,1956 : 岩波数学公式I . 岩波書店, 318pp.
: 森口, 宇田川, 一松編 ,1960 : 岩波数学公式III . 岩波書店, 310pp.
: 一松信, 1982 : 数値解析. 朝倉書店, 163pp.
: 森正武, 1984 : 数値解析法. 朝倉書店, 202pp.
: 寺沢寛一, 1983 : 自然科学者のための数学概論（増訂版）. 岩波書店, 711pp.

... 個 ^2.1

以下,

は偶数とする. dcpam5では, (Gauss 緯度としてとる場合には)

は偶数でなければならない.

... とする ^2.2

次の Legendre 函数 $P_J (\mu)$ は

$\displaystyle \left[ \DD{}{\mu} \left\{ (1-\mu^2) \DD{}{\mu} \right\} + J(J+1) \right] P_J(\mu) = 0$

(2.1)

を満たす

次多項式であり, $P_J (\mu)$ の零点は全て $-1 < \mu < 1$ にある. なお, Gauss 緯度は近似的には ${\displaystyle \sin^{-1} \left( \cos \frac{j-1/2}{J}\pi \right) }$ で与えられる.

... とり方は以下のとおりである ^2.3

このスキームは次のような特徴をもつ. ( 2005/04/04 石渡: 始めの 4 つは不正確な表現かも )

全領域積分した質量を保存
全領域積分したエネルギーを保存
全領域積分の角運動量を保存
全質量積分した温位を保存
静水圧の式がlocalにきまる. (下層の高度が上層の温度に依存しない)
水平方向に一定の, ある特定の温度分布について, 静水圧の式が正確になり, 気圧傾度力が0になる.
等温位大気はいつまでも等温位に留まる

... 整数レベルと半整数レベルを定義する ^2.4

物理量により, 整数レベルで定義されるものと, 半整数レベルで定義されるものがある.

... は乾燥空気の定圧比熱である ^2.5

いずれも定数としている.

... 次の直交条件を満たす ^2.6

詳しくは第A.5節を参照せよ.

... とは次の変換則に従う ^2.7

正変換, 逆変換時の係数は consistent に与えてさえいれば問題がない.

... 和をとることができる ^2.8

さらに, 実際の計算手続きとしては, $P_n^m(\sin \varphi)$ が,

が偶数 (even) の時 $\varphi=0$ について対称,

が奇数 (odd) の時 $\varphi=0$ について反対称であることを考慮して演算回数を減らすことができる. すなわち, $A_{ij}$ の計算では北半球のみについて南北対昭�蓿繙就�粮㏍芍��轣蛹≒鳫�笏蜿遐�竚癈鷭∂焜聨纃瘟赧漓�籬�㏍聽轣蛹就箟号恒聿高⊂桿轣蛹Γ蔚飴頏阡繝�籟鹿畩� $A_{ij}^{even}$ と反対昭�蓿繙就�粮㏍芍��轣蛹≒鳫�笏蜿遐�竚癈鷭∂焜聨纃瘟赧漓�籬�㏍聽轣蛹就求唄堰郡碓⊂桿轣蛹Γ蔚飴頏阡繝�籟鹿畩� $A_{ij}^{odd}$ についてそれぞれ計算し, 南半球については $A{i,J-j}=A_{ij}^{even}-A_{ij}^{odd}$ とすればよい. また,

の計算においては, その対称性, 反対称性に基づいて $A_{i,j}+A_{i,J-j}$ または $A_{i,j}-A_{i,J-j}$ の一方を

について 1から

まで加えればよい.

: 3. 力学過程 : dcpam5 支配方程式系とその離散化 : 1. この文書について

Yasuhiro MORIKAWA 平成20年11月9日

$\displaystyle \left( \DP{f}{\lambda} \right)_{ij}$	$\displaystyle \equiv \sum_{m=-M}^M \sum_{n=\vert m\vert}^N im \tilde{f}_n^m Y_n^m (\lambda_i,\varphi_j) ,$	(2.15)
$\displaystyle \widetilde{ \left( \DP{f}{\lambda} \right) }_n^m$	$\displaystyle = \frac{1}{I} \sum_{i=1}^I \sum_{j=1}^J im f_{ij} Y_n^{m*} (\lambda_i, \varphi_j) w_j .$	(2.16)

$\displaystyle \left(\DP{f}{\varphi}\right)_{ij}$	$\displaystyle \equiv \sum_{m=-M}^M \sum_{n=\vert m\vert}^N \tilde{f}_n^m \left. \DD{P_n^m}{\varphi} \right\vert _j \exp(im \lambda_i) ,$	(2.17)
$\displaystyle \widetilde{\left(\DP{f}{\varphi}\right)}_n^m$	$\displaystyle = - \frac{1}{I} \sum_{i=1}^I \sum_{j=1}^J f_{ij} \left. \DD{P_n^m}{\varphi}\right\vert _j \exp(-im \lambda_i) w_j .$	(2.18)