「投稿論文 revise 用の図集」

● 成長率の解像度依存性1(-1≦logE≦1 の範囲の全部の最大不安定モードを描いた図)

  
 
 左図: y方向64grid, 波数0.01刻みで計算した場合の結果       右図: y方向128grid, 波数0.01刻みで計算した場合の結果

● 成長率の解像度依存性2(-0.5≦logE≦0 の範囲の全部の最大不安定モードを描いた図)

  
 
 左図: y方向64grid, 波数0.01刻みで計算した場合の結果       右図: y方向128grid, 波数0.01刻みで計算した場合の結果


 ※計算する波数の刻み幅を細かくすると(下右図), 階段状の不連続は小さくなっているように見える)

 
  
 
 左図: y方向128grid, 波数0.01刻みで計算した場合の結果       右図: y方向128grid, 波数0.002刻みで計算した場合の結果

● 成長率の解像度依存性3(-2.5≦logE≦7.5 の範囲の全部の最大不安定モードを描いた図)

 ※全部の最大不安定モードを描くと, logE=-0.5 付近の飛びが良くわかるようになった

  
 
 左図: y方向64grid, 波数0.01刻みで計算した場合の結果       右図: y方向128grid, 波数0.01刻みで計算した場合の結果

● 最大不安定モードの構造の logE 依存性


● 最大不安定モードが存在する分散曲線上の不安定モードの構造の波数依存性

 (下図の青枠で囲まれた部分の不安定モードの構造を低波数側から並べた図)



● 連続モードの構造一覧(k を固定して見た場合)

 (下図の青線上(k=0.10)の連続モードの構造図)


● 連続モードの構造一覧(n を固定して見た場合)

 (下図の青線上(n=64,81,87,100)の連続モードの構造図)


● 一様Γ平面近似解で得られた分散曲線(ロスビーモードの場所を探してみた)

 ※一様Γ平面近似で表した Rossby モードの分散関係が 正しければ, 下図の2本の黄色の曲線の間に64本(n=1〜64)の Rossby モードの 分散曲線が収まっていると予想される. logE が大きくなると, Rossby モード の分散曲線は c=2.5 付近に集まる様子が見られる. Rossby モードの位相速度 は, 力学的赤道での基本場の流速値 2.5 を差し引くと, (無風の赤道波の場合 と同様に)負の値しか持たないように見える.

  
 
左図: (logE=-0.90 の場合の一様Γ平面近似解の分散曲線)    右図: (logE=-0.40 の場合一様Γ平面近似解の分散曲線)

 ※ 分散曲線の色の内訳は次の通り. ピンク(n=1 の E-IG), オレンジ(n=1 の W-IG) , 黄色(上: n=64 / 下: n=1 の Rossby モード).
その他の色は, n=0 の曲線(投稿論文準拠).
  
 
左図: (logE=-0.20 の場合の一様Γ平面近似解の分散曲線)    右図: (logE=-0.10 の場合の一様Γ平面近似解の分散曲線)

 ※ 分散曲線の色の内訳は次の通り. ピンク(n=1 の E-IG), オレンジ(n=1 の W-IG) , 黄色(上: n=64 / 下: n=1 の Rossby モード).
その他の色は, n=0 の曲線(投稿論文準拠).
  
 
左図: (logE=+1.00 の場合の一様Γ平面近似解の分散曲線)    右図: (logE=+1.50 の場合の一様Γ平面近似解の分散曲線)

 ※ 分散曲線の色の内訳は次の通り. ピンク(n=1 の E-IG), オレンジ(n=1 の W-IG) , 黄色(上: n=64 / 下: n=1 の Rossby モード).
その他の色は, n=0 の曲線(投稿論文準拠).
  
 
左図: (logE=+1.80 の場合の一様Γ平面近似解の分散曲線)    右図: (logE=+2.00 の場合の一様Γ平面近似解の分散曲線)

 ※ 分散曲線の色の内訳は次の通り. ピンク&紺(n=1 の E-IG), オレンジ(n=1 の W-IG) , 黄色(上: n=64 / 下: n=1 の Rossby モード). 淡い紫(n=1 の Rossby と W-IG の共鳴による不安定モード). その他の色は, n=0 の曲線(投稿論文準拠).

● 一様Γ平面近似解の不安定モードの範囲

  
 
左図: 影のある部分 K+WMRG の共鳴 / 影の無い部分 K+C の共鳴   右図: 影のある部分 R+WIG の共鳴 / 影の無い部分不安定モード無し


● 固有値計算で得られた中立波の構造1(MRG)

  
 
左図: 東進混合ロスビー重力波(logE=1.00, k=0.05, n=63)    右図: 西進混合ロスビー重力波(logE=1.00, k=0.05, n=129)
(それぞれ下図の★印のモードの構造を描いてある)

● Figure 1

  
 
左図: 基本場(-2 ≦ y ≦ 3)          右図: 基本場(-3 ≦ y ≦ 3)

● Figure 2

 
 
  左図: 基本場(-2 ≦ y ≦ 3 / 南北 grid 数 64)     右図: 基本場(-2 ≦ y ≦ 3 / 南北 grid 数 128)

 
基本場(-3 ≦ y ≦ 3 / 南北 grid 数 78)




● Figure 7 (アニメーション版1− grid 数の比較)

● Figure 7 (アニメーション版2− 南北計算領域の比較)

● Figure 8(a) (Kelvin mode − k を変えた場合の構造比較)
 
  左図: k=0.01                  右図: k=0.10
 
  左図: k=0.20                  右図: k=0.30
 
  左図: k=0.40                  右図: k=0.50
 
左図: k=0.60



● Figure 8(b) (continuous mode − k を変えた場合の構造比較)
 
  左図: k=0.01                  右図: k=0.10
 
  左図: k=0.20                  右図: k=0.30
 
  左図: k=0.40                  右図: k=0.50
 
左図: k=0.60



● Figure 14 (grid 数の比較 − k=0.01 の場合)
  基本場(-2 ≦ y ≦ 3 / 南北 grid 数 64)
  基本場(-2 ≦ y ≦ 3 / 南北 grid 数 128)



● Figure 14 (南北計算領域の比較 − k=0.01 の場合)
  基本場(-2 ≦ y ≦ 3 / 南北 grid 数 64)
  基本場(-3 ≦ y ≦ 3 / 南北 grid 数 78)


※ 北側の計算領域を広げて北側の境界ケルビン波が不安定になるかどうかのチェック(下図)

  基本場(-2 ≦ y ≦ 5 / 南北 grid 数 89)


● Figure 14 (南北計算領域の比較 − k > 0.01 の場合)

※ 北側の計算領域を広げて北側の境界ケルビン波が不安定になるかどう かのチェック
(-2≦y≦3 の領域では全波数で不安定になる場合があったので, 他の波数 でも見てみることにした)

  基本場(-2 ≦ y ≦ 5 / 南北 grid 数 89)
  基本場(-2 ≦ y ≦ 5 / 南北 grid 数 89)
  基本場(-2 ≦ y ≦ 5 / 南北 grid 数 89)
  基本場(-2 ≦ y ≦ 5 / 南北 grid 数 89)