[地球流体電脳倶楽部] [dcmodel | dcmodel-tools]

湿潤サーマル実験

目的

arare4 と arare5 で対流の立ち方が異なること (see 2012-12-27_sugiyama run1) をきっかけとして, arare4 と arare5 の違いを突き詰める. run21 までは, 乱流拡散係数を決めるのに用いる MixLen を dx, dy, dz の最小値にしていたが, 原論文では平均を使っているので, そのように変更する.

run23 までは, arare4 と arare5 で計算を揃えるために, 混合比の計算を arare4 に揃えていた. しかし, 問題があるので, arare5 のオリジナルに戻す. ファイルの接頭詞の sugiyama7 は計算しなおしたもの, sugiyama5 は run22 のもの.

• 音波減衰項の決め方を修正 (arare4 に合わせて, 水平方向・鉛直方向で同じ値を用いる)
• 乱流過程の MixLen の決め方が arare4 内で不整合を起こしていたので修正.
  • moist/moistbuoyancy.f90 と phys/turbulence.f90 の 2 カ所で別々の定式化によって MixLen を決めていた.
  • dx = dz ならば発現しないバグ
• 雨の落下を計算する部分で, 上部境界からのフラックスを明示的にゼロとした (arare4 のみ)
• 乱流拡散プロセスの MinLen の決め方を修正.

その他の設定など.

• 雲微物理有り
• 乱流過程有り
• 放射過程無し
• 地表面フラックス無し
• NDiffRatio = 50.0 => NDiffRatio = 1.0
• dy = 1000 m => dy = 400 m (dx と合わせる)
• 2 次の数値拡散を利用. AlphaNDiff2 = 1e-4
• 基本場は温位一定

arare5 のソース

• ソース: arare5 の cvs 版. 2012-12-25 に checkout.
  • 2 次の数値拡散を利用できるように修正.
  • 音波減衰項の係数の決め方を修正
• 設定ファイル

結果の比較

sugiyama7_Exner.nc_Exner_000000.png	sugiyama7_Exner.nc_Exner_000600.png	sugiyama7_Exner.nc_Exner_001200.png	sugiyama7_Exner.nc_Exner_001800.png	sugiyama7_Exner.nc_Exner_002400.png	sugiyama7_Exner.nc_Exner_003000.png	sugiyama7_Exner.nc_Exner_003600.png	sugiyama7_Exner.nc_Exner_004200.png
sugiyama5_Exner.nc_Exner_000000.png	sugiyama5_Exner.nc_Exner_000600.png	sugiyama5_Exner.nc_Exner_001200.png	sugiyama5_Exner.nc_Exner_001800.png	sugiyama5_Exner.nc_Exner_002400.png	sugiyama5_Exner.nc_Exner_003000.png	sugiyama5_Exner.nc_Exner_003600.png	sugiyama5_Exner.nc_Exner_004200.png
sugiyama7_Km.nc_Km_000000.png	sugiyama7_Km.nc_Km_000600.png	sugiyama7_Km.nc_Km_001200.png	sugiyama7_Km.nc_Km_001800.png	sugiyama7_Km.nc_Km_002400.png	sugiyama7_Km.nc_Km_003000.png	sugiyama7_Km.nc_Km_003600.png	sugiyama7_Km.nc_Km_004200.png
sugiyama5_Km.nc_Km_000000.png	sugiyama5_Km.nc_Km_000600.png	sugiyama5_Km.nc_Km_001200.png	sugiyama5_Km.nc_Km_001800.png	sugiyama5_Km.nc_Km_002400.png	sugiyama5_Km.nc_Km_003000.png	sugiyama5_Km.nc_Km_003600.png	sugiyama5_Km.nc_Km_004200.png
sugiyama7_PTemp.nc_PTemp_000000.png	sugiyama7_PTemp.nc_PTemp_000600.png	sugiyama7_PTemp.nc_PTemp_001200.png	sugiyama7_PTemp.nc_PTemp_001800.png	sugiyama7_PTemp.nc_PTemp_002400.png	sugiyama7_PTemp.nc_PTemp_003000.png	sugiyama7_PTemp.nc_PTemp_003600.png	sugiyama7_PTemp.nc_PTemp_004200.png
sugiyama5_PTemp.nc_PTemp_000000.png	sugiyama5_PTemp.nc_PTemp_000600.png	sugiyama5_PTemp.nc_PTemp_001200.png	sugiyama5_PTemp.nc_PTemp_001800.png	sugiyama5_PTemp.nc_PTemp_002400.png	sugiyama5_PTemp.nc_PTemp_003000.png	sugiyama5_PTemp.nc_PTemp_003600.png	sugiyama5_PTemp.nc_PTemp_004200.png
sugiyama7_VelX.nc_VelX_000000.png	sugiyama7_VelX.nc_VelX_000600.png	sugiyama7_VelX.nc_VelX_001200.png	sugiyama7_VelX.nc_VelX_001800.png	sugiyama7_VelX.nc_VelX_002400.png	sugiyama7_VelX.nc_VelX_003000.png	sugiyama7_VelX.nc_VelX_003600.png	sugiyama7_VelX.nc_VelX_004200.png
sugiyama5_VelX.nc_VelX_000000.png	sugiyama5_VelX.nc_VelX_000600.png	sugiyama5_VelX.nc_VelX_001200.png	sugiyama5_VelX.nc_VelX_001800.png	sugiyama5_VelX.nc_VelX_002400.png	sugiyama5_VelX.nc_VelX_003000.png	sugiyama5_VelX.nc_VelX_003600.png	sugiyama5_VelX.nc_VelX_004200.png
sugiyama7_VelZ.nc_VelZ_000000.png	sugiyama7_VelZ.nc_VelZ_000600.png	sugiyama7_VelZ.nc_VelZ_001200.png	sugiyama7_VelZ.nc_VelZ_001800.png	sugiyama7_VelZ.nc_VelZ_002400.png	sugiyama7_VelZ.nc_VelZ_003000.png	sugiyama7_VelZ.nc_VelZ_003600.png	sugiyama7_VelZ.nc_VelZ_004200.png
sugiyama5_VelZ.nc_VelZ_000000.png	sugiyama5_VelZ.nc_VelZ_000600.png	sugiyama5_VelZ.nc_VelZ_001200.png	sugiyama5_VelZ.nc_VelZ_001800.png	sugiyama5_VelZ.nc_VelZ_002400.png	sugiyama5_VelZ.nc_VelZ_003000.png	sugiyama5_VelZ.nc_VelZ_003600.png	sugiyama5_VelZ.nc_VelZ_004200.png
sugiyama7_H2O-g.nc_H2O-g_000000.png	sugiyama7_H2O-g.nc_H2O-g_000600.png	sugiyama7_H2O-g.nc_H2O-g_001200.png	sugiyama7_H2O-g.nc_H2O-g_001800.png	sugiyama7_H2O-g.nc_H2O-g_002400.png	sugiyama7_H2O-g.nc_H2O-g_003000.png	sugiyama7_H2O-g.nc_H2O-g_003600.png	sugiyama7_H2O-g.nc_H2O-g_004200.png
sugiyama5_H2O-g.nc_H2O-g_000000.png	sugiyama5_H2O-g.nc_H2O-g_000600.png	sugiyama5_H2O-g.nc_H2O-g_001200.png	sugiyama5_H2O-g.nc_H2O-g_001800.png	sugiyama5_H2O-g.nc_H2O-g_002400.png	sugiyama5_H2O-g.nc_H2O-g_003000.png	sugiyama5_H2O-g.nc_H2O-g_003600.png	sugiyama5_H2O-g.nc_H2O-g_004200.png
sugiyama7_H2O-l-Cloud.nc_H2O-l-Cloud_000000.png	sugiyama7_H2O-l-Cloud.nc_H2O-l-Cloud_000600.png	sugiyama7_H2O-l-Cloud.nc_H2O-l-Cloud_001200.png	sugiyama7_H2O-l-Cloud.nc_H2O-l-Cloud_001800.png	sugiyama7_H2O-l-Cloud.nc_H2O-l-Cloud_002400.png	sugiyama7_H2O-l-Cloud.nc_H2O-l-Cloud_003000.png	sugiyama7_H2O-l-Cloud.nc_H2O-l-Cloud_003600.png	sugiyama7_H2O-l-Cloud.nc_H2O-l-Cloud_004200.png
sugiyama5_H2O-l-Cloud.nc_H2O-l-Cloud_000000.png	sugiyama5_H2O-l-Cloud.nc_H2O-l-Cloud_000600.png	sugiyama5_H2O-l-Cloud.nc_H2O-l-Cloud_001200.png	sugiyama5_H2O-l-Cloud.nc_H2O-l-Cloud_001800.png	sugiyama5_H2O-l-Cloud.nc_H2O-l-Cloud_002400.png	sugiyama5_H2O-l-Cloud.nc_H2O-l-Cloud_003000.png	sugiyama5_H2O-l-Cloud.nc_H2O-l-Cloud_003600.png	sugiyama5_H2O-l-Cloud.nc_H2O-l-Cloud_004200.png
sugiyama7_H2O-l-Rain.nc_H2O-l-Rain_000000.png	sugiyama7_H2O-l-Rain.nc_H2O-l-Rain_000600.png	sugiyama7_H2O-l-Rain.nc_H2O-l-Rain_001200.png	sugiyama7_H2O-l-Rain.nc_H2O-l-Rain_001800.png	sugiyama7_H2O-l-Rain.nc_H2O-l-Rain_002400.png	sugiyama7_H2O-l-Rain.nc_H2O-l-Rain_003000.png	sugiyama7_H2O-l-Rain.nc_H2O-l-Rain_003600.png	sugiyama7_H2O-l-Rain.nc_H2O-l-Rain_004200.png
sugiyama5_H2O-l-Rain.nc_H2O-l-Rain_000000.png	sugiyama5_H2O-l-Rain.nc_H2O-l-Rain_000600.png	sugiyama5_H2O-l-Rain.nc_H2O-l-Rain_001200.png	sugiyama5_H2O-l-Rain.nc_H2O-l-Rain_001800.png	sugiyama5_H2O-l-Rain.nc_H2O-l-Rain_002400.png	sugiyama5_H2O-l-Rain.nc_H2O-l-Rain_003000.png	sugiyama5_H2O-l-Rain.nc_H2O-l-Rain_003600.png	sugiyama5_H2O-l-Rain.nc_H2O-l-Rain_004200.png