message:l,m:{
<small>
  [ <a href="/">地球流体電脳倶楽部</a> /
    <a href="/library/dcpam">DCPAM</a> /
    <a href="/library/dcpam/sample">計算結果</a> /
    <a href="/library/dcpam/sample/held-suarez-1994.htm">Held and Suarez (1994) テスト</a> /
    <a href="SIGEN.htm">SIGEN.htm</a> ] 
</small>
}
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:{
<table border="1" cellspacing="1" cellpadding="5" align="center">
  <tr>
    <td width="20%">システム (uname -a)</td>
    <td>
      Linux karei 2.6.18-6-686 #1 SMP
    </td>
  </tr>
  
  <tr>
    <td width="20%">コンパイラ</td>
    <td>
      frt, Fujitsu Fortran Compiler Driver Version 5.0 (Nov 12 2003 20:08:33)
    </td>
  </tr>

  <tr>
    <td width="20%">依存するライブラリとバージョン</td>
    <td>
      <ul>
	<li> netcdf ver 3.6.1 compiled with frt and gcc
	<li> ispack ver 0.71 compiled with Fujitsu compilers
	<li> gtool5 ver 20090324 compiled with frt
	<li> spml   ver 0.5.0 compiled with frt
      </ul>
    </td>
  </tr>

</table>
}
message:l,m,rd:{
== 実験設定
}
message:l,m:{
=begin
* 解像度: T63L20
  * 鉛直層分布: sigma = P / Ps を等間隔で分割
* 水平拡散: 
  * ラプラシアンの 4 乗
  * 最大全波数, 63, において時定数 0.1 日
* 時間ステップ: 15 分
=end
}
message:l,m:{
#<table border="1" cellspacing="1" cellpadding="5" align="center">
#  <tr>
#    <td width="20%">プログラムのソースコード</td>
#    <td>
#      <a href="dcpam4-20071012-custom">dcpam4-20071012 に若干手を加えたもの</a>
#      [<a href="dcpam4-20071012-custom.tar.gz">TAR.GZ</a>] <br>
#      <a href="dcpam4-20071012-custom/Config.mk">ビルド時の設定ファイル</a>
#    </td>
#  </tr>
#
#  <tr>
#    <td width="20%">NAMELIST ファイル</td>
#    <td>
#      <a href="info00/dcpam_hs94_T42L20_1200day.nml">info00/dcpam_hs94_T42L20_1200day.nml</a>
#    </td>
#  </tr>
#
#  <tr>
#    <td width="20%">初期値</td>
#    <td>
#      渦度, 発散, 比湿はゼロ.
#      地表面気圧は一定 (1000 hPa).
#      温度は基本場が 250 K で微笑擾乱が与えられている.
#      擾乱の分布については
#      <a href="#label-6">初期値および実験結果</a> の項目を参照のこと. 
#    </td>
#  </tr>
#
#  <tr><td colspan="2">&nbsp;<td></tr>
#
#  <tr>
#    <td>計算領域</td>
#    <td>
#      水平(λ,φ)      : 全球 (λ=[0, 360], φ=[-90, 90])<br>
#      鉛直(σ=p/p<sub>0</sub>) : 大気下端〜上端 (σ = [0, 1])
#    </td>
#  </tr>
#  <tr>
#    <td>時間間隔</td>
#    <td>30 分</td>
#  </tr>
#  <tr>
#    <td>積分時間</td>
#    <td>1200 日</td>
#  </tr>
#  <tr>
#    <td>境界条件</td>
#    <td>
#      水平(λ,φ) : 周期境界条件<br>
#      鉛直(σ)    : dσ/dt = 0 (σ = 0, 1)
#    </td>
#  </tr>
#  <tr>
#    <td>数値解法</td>
#    <td>
#      <ul>
#        <li>水平(λ,φ) : スペクトル法 (球面調和関数展開・三角形切断)<br>
#        <li>鉛直(σ)    : 差分法
#                          (<a href="#References">Arakawa and Suarez, 1983</a>)<br>
#        <li>時間(t)     : リープフロッグスキーム.
#	  <ul>
#		     <li> ただし重力波の項にはセミインプリシット法を用いる. 
#                          時間フィルター (<a href="#References">Asselin, 1972</a>).
#		          を用いる. フィルター係数は 0.05
#	  </ul>
#      </ul>
#    </td>
#  </tr>
#  <tr><td colspan="2">&nbsp;<td></tr>
#</table>
#
##<table border="1" cellspacing="1" cellpadding="5" align="center">
#  <tr>
#    <td>空間解像度</td>
#    <td>
#      東西格子点数 : 192<br>
#      南北格子点数 : 96<br>
#      鉛直格子点数 : 20
#    </td>
#  </tr>
#  <tr>
#    <td>最大波数</td>
#    <td>63</td>
#  </tr>
#  <tr>
#    <td>鉛直層分布</td>
#    <td>sigma = P / Ps を等間隔で分割</td>
#  </tr>
##  <tr><td colspan="2">&nbsp;<td></tr>
##</table>
#  
#<table border="1" cellspacing="1" cellpadding="5" align="center">
#  <tr>
#    <td>惑星半径</td>
#    <td>6.371e+6 [m]</td>
#  </tr>
#  <tr>
#    <td>回転角速度</td>
#    <td>7.292e-5 [1/s]</td>
#  </tr>
#  <tr>
#    <td>重力加速度</td>
#    <td> 9.8 [m/s^2]</td>
#  </tr>
#  <tr>
#    <td>大気定圧比熱</td>
#    <td>1004.6  [J K^-1 kg^-1] </td>
#  </tr>
#  <tr>
#    <td>大気気体定数</td>
#    <td>287.04  [J K^-1 kg^-1] </td>
#  </tr>
#  <tr>
#    <td>最大波数に対する e-folding time</td>
#    <td>8640.0  [s] </td>
#  </tr>
#  <tr>
#    <td>超粘性の次数</td>
#    <td>8</td>
#  </tr>
#</table>
}
#message:l,m,rd:{
#== 実行結果に関する付加情報
#}
#message:l,m:{
#<table border="1" cellspacing="1" cellpadding="5" align="center">
#  <tr>
#    <td width="20%">標準出力</td>
#    <td>
#      <a href="info00/dcpam_hs94_T42L20_1200day_stdout.txt">info00/dcpam_hs94_T42L20_1200day_stdout.txt</a>
#    </td>
#  </tr>
#
#  <tr>
#    <td width="20%">標準エラー出力</td>
#    <td>
#      <a href="info00/dcpam_hs94_T42L20_1200day_stderr.txt">info00/dcpam_hs94_T42L20_1200day_stderr.txt</a>
#    </td>
#  </tr>
#  
#  <tr>
#    <td width="20%">プロファイラによる速度計測</td>
#    <td>
#      <a href="info00/ftrace.txt">info00/ftrace.txt</a>
#    </td>
#  </tr>
#</table>
#}
message::<a name="References"></a>
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:{
#<hr>
#初期温度場
#<hr>
#}
#::
#::
#::
#init_Temp_0.97Sig.png:: σ=0.97
#init_Temp_0.52Sig.png:: σ=0.52
#label:m:{
#最上端である σ=0.02 では 250K で一定. 
#}
#::
label:m:{
<hr>
実験結果
<hr>
}
::
::
::
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'))
}