: 10. 飽和比湿 : dcpam5 支配方程式系とその離散化 : 8. 蒸発・凝結による地表面気圧変化

9. 鉛直フィルター

9.1 背景と目的

太陽定数を増加させた計算を行ったところ 2-grid noise が生じた. このノイズの振幅は非常に大きくなり計算が破綻する.

このノイズを消すためにフィルターをかける. 物理的実体は無い. 強いて言えば, 重力波でつぶれるかもね, 程度である.

9.2 基本的な手続き

調節のための基本温度場を決める.
調節を行う部分を決める. 具体的には, 全層か一部か. 一部としたらどの領域か
温度構造を基本場を調節する. その後, 誤差の補正を行う.

9.3 調節のための基本温度場

ギザギザ成分を取り除いた分布を基本場としたい. そこで「基本場」は次で与えられるとする.

$\displaystyle T_{Bk} = \frac{T_{k+1/2} + T_{k-1/2}}{2}$

(9.1)

ここで, $1 \leq k \leq KMAX$

9.4 調節する部分の決定

全層の場合
何も考えずに全部やる
ギザギザ部分のみ調節する場合
以下の式を満たす level はギザギザな点と判定する.

$\displaystyle (T_{k} - T_{k-1}) \cdot (T_{k+1} - T_{k}) < 0$ (9.2)

実際には, こんな凝ったやってないと思う.

9.5 調節および誤差の補正

節のタイトルが正しくないような気がする. 「誤差の補正」ではないなあ.

全層で調節を行い, 全層に誤差をばらまく場合.

これについては未実装である.
これはもっとも安直な調節方法である. 調節前の温度の値を $T_{k}$ , 調節後を $\hat{T}_{k}$ とする.
何も考えずに調節するなら

$\displaystyle \hat{T}_{k} = T_{k} + S_{grst} (T_{Bk} - T_{k})$ (9.3)

温度変化量は

$\displaystyle \Delta T_{k} \equiv \hat{T}_{k} - T_{k} = S_{grst} (T_{Bk} - T_{k})$ (9.4)

この場合, 全体で

$\displaystyle \sum^{KMAX}_{k=1} c_v \Delta T_{k} \Delta p_{k} = \sum^{KMAX}_{k=1} c_v S_{grst} (T_{Bk} - T_{k}) (p_s \Delta \sigma_{k})$ (9.5)

の内部エネルギーをコラムに与えてしまうことになる. これが誤差になる.
調節の前後で全体の内部エネルギー量に変化がないようにするためには

$\displaystyle \sum^{KMAX}_{k=1} \Delta T_{k} \Delta p_{k} = 0$ (9.6)

としてやらなければいけない. ただし, ここで比熱が一定の場合を仮定した.
これを解決するもっとも安直な方法は上の誤差を全層にばらまくこと. この場合, 温度の調節量は以下の式で与えられる.

$\displaystyle \hat{T}_{k} = T_{k} + S_{grst} (T_{Bk} - T_{k}) - \frac{ \sum^{KM... ...k=1} S_{grst} (T_{Bk} - T_{k}) \Delta p_{k} }{ \sum^{KMAX}_{k=1} \Delta p_{k} }$ (9.7)

こうすると

$\displaystyle \sum^{KMAX}_{k=1} \hat{T}_{k} \Delta p_{k} = \sum^{KMAX}_{k=1} \l... ...- T_{k}) \Delta p_{k} }{ \sum^{KMAX}_{k=1} \Delta p_{k} } \Delta p_{k} \right\}$ (9.8)

となる.
全層で調節し, 誤差を局所的に解消していく場合
AGCM5 の場合だと, p2grstA.F
3 点トリオで考え, そこで生じた誤差をその 3 層にばらまくことにする. 下から順に, 以下の式で温度の調節量を計算していく.

$\displaystyle \hat{T}_{k-1}$ $\displaystyle =$ $\displaystyle T_{k-1} + S_{grst} (T_{Bk-1} - T_{k-1}) - \frac{ \sum^{k+1}_{k=k-1} S_{grst} (T_{Bk} - T_{k}) \Delta p_{k} }{ \sum^{k+1}_{k=k-1} \Delta p_{k} },$ (9.9)

$\displaystyle \hat{T}_{k}$ $\displaystyle =$ $\displaystyle T_{k} + S_{grst} (T_{Bk} - T_{k}) - \frac{ \sum^{k+1}_{k=k-1} S_{grst} (T_{Bk} - T_{k}) \Delta p_{k} }{ \sum^{k+1}_{k=k-1} \Delta p_{k} },$ (9.10)

$\displaystyle \hat{T}_{k+1}$ $\displaystyle =$ $\displaystyle T_{k+1} + S_{grst} (T_{Bk+1} - T_{k+1}) - \frac{ \sum^{k+1}_{k=k-1} S_{grst} (T_{Bk} - T_{k}) \Delta p_{k} }{ \sum^{k+1}_{k=k-1} \Delta p_{k} }$ (9.11)

この $\hat{T}$ を用いて, 1 つ上に上がり,

$\displaystyle \hat{\hat{T}}_{k}$ $\displaystyle =$ $\displaystyle \hat{T}_{k} + S_{grst} (T_{Bk} - \hat{T}_{k}) - \frac{ \sum^{k+2}... ..._{grst} (T_{Bk} - \hat{T}_{k}) \Delta p_{k} }{ \sum^{k+2}_{k=k} \Delta p_{k} },$ (9.12)

$\displaystyle \hat{\hat{T}}_{k+1}$ $\displaystyle =$ $\displaystyle \hat{T}_{k+1} + S_{grst} (T_{Bk+1} - \hat{T}_{k+1}) - \frac{ \sum... ..._{grst} (T_{Bk} - \hat{T}_{k}) \Delta p_{k} }{ \sum^{k+2}_{k=k} \Delta p_{k} },$ (9.13)

$\displaystyle \hat{T}_{k+2}$ $\displaystyle =$ $\displaystyle T_{k+1} + S_{grst} (T_{Bk+2} - T_{k+2}) - \frac{ \sum^{k+2}_{k=k} S_{grst} (T_{Bk} - \hat{T}_{k}) \Delta p_{k} }{ \sum^{k+2}_{k=k} \Delta p_{k} }$ (9.14)

ただし, 右辺 3 項目の分子の和においてがの場合は, $\hat{T}_{k+2}$ ではなく $T_{k+2}$ である(まだ調節を行っていないので).
以上をlevel を 1 つづつ上がりながら順番に行う.
ギザギザ部分だけを調節し誤差を 3 点ごとにばらまく
下から順に

$\displaystyle (T_{k} - T_{k-1}) \cdot (T_{k+1} - T_{k}) < 0$ (9.15)

となる部分の 3 点トリオ出 (2) と同じことを行う.
ギザギザ部分だけを調節し, 誤差をギザギザ部分全体に均等にばらまく場合
連続してギザギザ部分になっているところを判定し (1) と同じことをする. 下記のように判定する.
1. 下から順に登っていって初めて
  
  $\displaystyle (T_{k} - T_{k-1}) \cdot (T_{k+1} - T_{k}) < 0$ (9.16)
  
  となったところで, level がギザギザ部分の底と判定する.
2. ギザギザ部分で登っていって
  
  $\displaystyle (T_{k} - T_{k-1}) \cdot (T_{k+1} - T_{k}) \geq 0$ (9.17)
  
  となったら, level がギザギザ部分の天井.
このように判定されたギザギザ部分全体に誤差をばらまく.

: 10. 飽和比湿 : dcpam5 支配方程式系とその離散化 : 8. 蒸発・凝結による地表面気圧変化

Yasuhiro MORIKAWA 平成20年8月9日

$\displaystyle \hat{T}_{k-1}$	$\displaystyle =$	$\displaystyle T_{k-1} + S_{grst} (T_{Bk-1} - T_{k-1}) - \frac{ \sum^{k+1}_{k=k-1} S_{grst} (T_{Bk} - T_{k}) \Delta p_{k} }{ \sum^{k+1}_{k=k-1} \Delta p_{k} },$	(9.9)
$\displaystyle \hat{T}_{k}$	$\displaystyle =$	$\displaystyle T_{k} + S_{grst} (T_{Bk} - T_{k}) - \frac{ \sum^{k+1}_{k=k-1} S_{grst} (T_{Bk} - T_{k}) \Delta p_{k} }{ \sum^{k+1}_{k=k-1} \Delta p_{k} },$	(9.10)
$\displaystyle \hat{T}_{k+1}$	$\displaystyle =$	$\displaystyle T_{k+1} + S_{grst} (T_{Bk+1} - T_{k+1}) - \frac{ \sum^{k+1}_{k=k-1} S_{grst} (T_{Bk} - T_{k}) \Delta p_{k} }{ \sum^{k+1}_{k=k-1} \Delta p_{k} }$	(9.11)

$\displaystyle \hat{\hat{T}}_{k}$	$\displaystyle =$	$\displaystyle \hat{T}_{k} + S_{grst} (T_{Bk} - \hat{T}_{k}) - \frac{ \sum^{k+2}... ..._{grst} (T_{Bk} - \hat{T}_{k}) \Delta p_{k} }{ \sum^{k+2}_{k=k} \Delta p_{k} },$	(9.12)
$\displaystyle \hat{\hat{T}}_{k+1}$	$\displaystyle =$	$\displaystyle \hat{T}_{k+1} + S_{grst} (T_{Bk+1} - \hat{T}_{k+1}) - \frac{ \sum... ..._{grst} (T_{Bk} - \hat{T}_{k}) \Delta p_{k} }{ \sum^{k+2}_{k=k} \Delta p_{k} },$	(9.13)
$\displaystyle \hat{T}_{k+2}$	$\displaystyle =$	$\displaystyle T_{k+1} + S_{grst} (T_{Bk+2} - T_{k+2}) - \frac{ \sum^{k+2}_{k=k} S_{grst} (T_{Bk} - \hat{T}_{k}) \Delta p_{k} }{ \sum^{k+2}_{k=k} \Delta p_{k} }$	(9.14)