message:l,m:{ [ 地球流体電脳倶楽部 / DCPAM / 計算結果 / Held and Suarez (1994) テスト / SIGEN.htm ] } message:l,m,rd:{ = Held and Suarez (1994) 実験 # * 2009/09/21 (高橋 芳幸) 2007-12-10_morikawa 以下のファイルをコピーして作成 == 概要 Held and Suarez (1994) の力学コアベンチマークテストを行った. 環境, およびパラメータは以下の通り. } message:l,m,rd:{ == 実験環境 } message:l,m:{
システム (uname -a) Linux karei 2.6.18-6-686 #1 SMP
コンパイラ frt, Fujitsu Fortran Compiler Driver Version 5.0 (Nov 12 2003 20:08:33)
依存するライブラリとバージョン
  • netcdf ver 3.6.1 compiled with frt and gcc
  • ispack ver 0.71 compiled with Fujitsu compilers
  • gtool5 ver 20090324 compiled with frt
  • spml ver 0.5.0 compiled with frt
} message:l,m,rd:{ == 実験設定 } message:l,m:{ =begin * 解像度: T63L20 * 鉛直層分布: sigma = P / Ps を等間隔で分割 * 水平拡散: * ラプラシアンの 4 乗 * 最大全波数, 63, において時定数 0.1 日 * 時間ステップ: 15 分 =end } message:l,m:{ # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
プログラムのソースコード # dcpam4-20071012 に若干手を加えたもの # [TAR.GZ]
# ビルド時の設定ファイル #
NAMELIST ファイル # info00/dcpam_hs94_T42L20_1200day.nml #
初期値 # 渦度, 発散, 比湿はゼロ. # 地表面気圧は一定 (1000 hPa). # 温度は基本場が 250 K で微笑擾乱が与えられている. # 擾乱の分布については # 初期値および実験結果 の項目を参照のこと. #
 
計算領域 # 水平(λ,φ) : 全球 (λ=[0, 360], φ=[-90, 90])
# 鉛直(σ=p/p0) : 大気下端〜上端 (σ = [0, 1]) #
時間間隔30 分
積分時間1200 日
境界条件 # 水平(λ,φ) : 周期境界条件
# 鉛直(σ) : dσ/dt = 0 (σ = 0, 1) #
数値解法 #
    #
  • 水平(λ,φ) : スペクトル法 (球面調和関数展開・三角形切断)
    #
  • 鉛直(σ) : 差分法 # (Arakawa and Suarez, 1983)
    #
  • 時間(t) : リープフロッグスキーム. #
      #
    • ただし重力波の項にはセミインプリシット法を用いる. # 時間フィルター (Asselin, 1972). # を用いる. フィルター係数は 0.05 #
    #
#
 
# ## # # # # # # # # # # # # ## ##
空間解像度 # 東西格子点数 : 192
# 南北格子点数 : 96
# 鉛直格子点数 : 20 #
最大波数63
鉛直層分布sigma = P / Ps を等間隔で分割
 
# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
惑星半径6.371e+6 [m]
回転角速度7.292e-5 [1/s]
重力加速度 9.8 [m/s^2]
大気定圧比熱1004.6 [J K^-1 kg^-1]
大気気体定数287.04 [J K^-1 kg^-1]
最大波数に対する e-folding time8640.0 [s]
超粘性の次数8
} #message:l,m,rd:{ #== 実行結果に関する付加情報 #} #message:l,m:{ # # # # # # # # # # # # # # # #
標準出力 # info00/dcpam_hs94_T42L20_1200day_stdout.txt #
標準エラー出力 # info00/dcpam_hs94_T42L20_1200day_stderr.txt #
プロファイラによる速度計測 # info00/ftrace.txt #
#} message:: message:l,m,rd:{ =begin == 参考文献 #* Arakawa, A., Suarez, M. J., 1983: # Vertical differencing of the primitive equations in sigma coordinates. # ((|Mon. Wea. Rev.|)), ((*111*)), 34--35. #* Asselin, R. A., 1972: # Frequency filter for time integrations. # ((|Mon. Wea. Rev.|)), ((*100*)), 487--490. * Held, I. M., and Suarez, M. J., 1994: A proposal for the intercomparison of the dynamical cores of atmospheric general circuation models. ((|Bull. Am. Meteor. Soc.|)), ((* 75 *)), 1825--1830. ((<(abstract)|URL:http://ams.allenpress.com/perlserv/?request=get-abstract&doi=10.1175%2F1520-0477%281994%29075%3C1825%3AAPFTIO%3E2.0.CO%3B2>)) ((<(pdf)|URL:http://ams.allenpress.com/archive/1520-0477/75/10/pdf/i1520-0477-75-10-1825.pdf>)) =end } #message:m,rd:{ #== 初期値および実験結果 #} ## ## Following parts are Thumbnail comments ## #label:m:{ #
#初期温度場 #
#} #:: #:: #:: #init_Temp_0.97Sig.png:: σ=0.97 #init_Temp_0.52Sig.png:: σ=0.52 #label:m:{ #最上端である σ=0.02 では 250K で一定. #} #:: label:m:{
実験結果
} :: :: :: ta.gif:m,<:{ 200 日目から 1200 日目までの時間・東西平均温度 (bar(T)) } ua.gif:m,<:{ 200 日目から 1200 日目までの時間・東西平均東西風 (bar(U)) } ta2.gif:m,<:{ 200 日目から 1200 日目までの時間・東西平均温度分散 (bar(T'^2)) } uaspe.gif:m,<:{ 200 日目から 1200 日目までの時間・東西・鉛直平均東西風東西スペクトル } vt.gif:m,<:{ 200 日目から 1200 日目までの時間・東西平均南北風-温度共分散 (bar(v'T')) }