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

湿潤サーマル実験

目的

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

run26 に比べて, 雲微物理過程の計算の順序を変更. 雲 => 雨, 雨 => 蒸気, 蒸気 <=> 雲. 湿度はほぼ 100%.

gpmaxmin sugiyama10_H2O-g_Hum.nc@H2O-g_Hum
sugiyama10_H2O-g_Hum.nc@H2O-g_Hum : max=100.009948730469, at, 
    ...(total 1 points)...

sugiyama10_H2O-g_Hum.nc@H2O-g_Hum : min=28.4048194885254, at,  
    ...(total 1 points)...

• 音波減衰項の決め方を修正 (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 次の数値拡散を利用できるように修正.
  • 音波減衰項の係数の決め方を修正
• 設定ファイル

結果の比較

sugiyama10_Exner.nc_Exner_000000.png	sugiyama10_Exner.nc_Exner_000600.png	sugiyama10_Exner.nc_Exner_001200.png	sugiyama10_Exner.nc_Exner_001800.png	sugiyama10_Exner.nc_Exner_002400.png	sugiyama10_Exner.nc_Exner_003000.png	sugiyama10_Exner.nc_Exner_003600.png	sugiyama10_Exner.nc_Exner_004200.png
sugiyama10_H2O-g.nc_H2O-g_000000.png	sugiyama10_H2O-g.nc_H2O-g_000600.png	sugiyama10_H2O-g.nc_H2O-g_001200.png	sugiyama10_H2O-g.nc_H2O-g_001800.png	sugiyama10_H2O-g.nc_H2O-g_002400.png	sugiyama10_H2O-g.nc_H2O-g_003000.png	sugiyama10_H2O-g.nc_H2O-g_003600.png	sugiyama10_H2O-g.nc_H2O-g_004200.png
sugiyama10_H2O-g_Hum.nc_H2O-g_Hum_000000.png	sugiyama10_H2O-g_Hum.nc_H2O-g_Hum_000600.png	sugiyama10_H2O-g_Hum.nc_H2O-g_Hum_001200.png	sugiyama10_H2O-g_Hum.nc_H2O-g_Hum_001800.png	sugiyama10_H2O-g_Hum.nc_H2O-g_Hum_002400.png	sugiyama10_H2O-g_Hum.nc_H2O-g_Hum_003000.png	sugiyama10_H2O-g_Hum.nc_H2O-g_Hum_003600.png	sugiyama10_H2O-g_Hum.nc_H2O-g_Hum_004200.png
sugiyama10_H2O-l-Cloud.nc_H2O-l-Cloud_000000.png	sugiyama10_H2O-l-Cloud.nc_H2O-l-Cloud_000600.png	sugiyama10_H2O-l-Cloud.nc_H2O-l-Cloud_001200.png	sugiyama10_H2O-l-Cloud.nc_H2O-l-Cloud_001800.png	sugiyama10_H2O-l-Cloud.nc_H2O-l-Cloud_002400.png	sugiyama10_H2O-l-Cloud.nc_H2O-l-Cloud_003000.png	sugiyama10_H2O-l-Cloud.nc_H2O-l-Cloud_003600.png	sugiyama10_H2O-l-Cloud.nc_H2O-l-Cloud_004200.png
sugiyama10_H2O-l-Rain.nc_H2O-l-Rain_000000.png	sugiyama10_H2O-l-Rain.nc_H2O-l-Rain_000600.png	sugiyama10_H2O-l-Rain.nc_H2O-l-Rain_001200.png	sugiyama10_H2O-l-Rain.nc_H2O-l-Rain_001800.png	sugiyama10_H2O-l-Rain.nc_H2O-l-Rain_002400.png	sugiyama10_H2O-l-Rain.nc_H2O-l-Rain_003000.png	sugiyama10_H2O-l-Rain.nc_H2O-l-Rain_003600.png	sugiyama10_H2O-l-Rain.nc_H2O-l-Rain_004200.png
sugiyama10_Km.nc_Km_000000.png	sugiyama10_Km.nc_Km_000600.png	sugiyama10_Km.nc_Km_001200.png	sugiyama10_Km.nc_Km_001800.png	sugiyama10_Km.nc_Km_002400.png	sugiyama10_Km.nc_Km_003000.png	sugiyama10_Km.nc_Km_003600.png	sugiyama10_Km.nc_Km_004200.png
sugiyama10_PTemp.nc_PTemp_000000.png	sugiyama10_PTemp.nc_PTemp_000600.png	sugiyama10_PTemp.nc_PTemp_001200.png	sugiyama10_PTemp.nc_PTemp_001800.png	sugiyama10_PTemp.nc_PTemp_002400.png	sugiyama10_PTemp.nc_PTemp_003000.png	sugiyama10_PTemp.nc_PTemp_003600.png	sugiyama10_PTemp.nc_PTemp_004200.png
sugiyama10_VelX.nc_VelX_000000.png	sugiyama10_VelX.nc_VelX_000600.png	sugiyama10_VelX.nc_VelX_001200.png	sugiyama10_VelX.nc_VelX_001800.png	sugiyama10_VelX.nc_VelX_002400.png	sugiyama10_VelX.nc_VelX_003000.png	sugiyama10_VelX.nc_VelX_003600.png	sugiyama10_VelX.nc_VelX_004200.png
sugiyama10_VelZ.nc_VelZ_000000.png	sugiyama10_VelZ.nc_VelZ_000600.png	sugiyama10_VelZ.nc_VelZ_001200.png	sugiyama10_VelZ.nc_VelZ_001800.png	sugiyama10_VelZ.nc_VelZ_002400.png	sugiyama10_VelZ.nc_VelZ_003000.png	sugiyama10_VelZ.nc_VelZ_003600.png	sugiyama10_VelZ.nc_VelZ_004200.png