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-suarez1994.htm">Held and Suarez (1994) テスト</a> /
    <a href="SIGEN.htm">SIGEN.htm</a> ] 
</small>
}
message:l,m,rd:{
= Held and Suarez (1994) テスト

#* 森川 靖大
#  * 2007/12/10 (森川 靖大) 作成

== 概要
DCPAM の力学過程の実装実験として、
Held and Suarez (1994) の力学コアベンチマークテストを行った。
このテストは、大気大循環モデルの力学過程を
評価するためのベンチマーク計算である。
この実験では、温度場にニュートン冷却を与え、σ < 0.7 の速度場に
レイリー摩擦を与え、
ニュートン冷却で大気をある温度分布に近づける。
(詳細は Held and Suarez (1994) を参照のこと)。

環境, およびパラメータは以下の通り. 
結果は下記に記す.

== 実験環境
}
message:l,m:{
<table border="1" cellspacing="1" cellpadding="5" align="center">
  <tr>
    <td width="20%">システム (uname -a)</td>
    <td>
      SUPER-UX unix 17.1  SX-8R
    </td>
  </tr>
  
  <tr>
    <td width="20%">コンパイラ</td>
    <td>
      FORTRAN90/SX
    </td>
  </tr>

  <tr>
    <td width="20%">コンパイラオプション</td>
    <td>
      -w -sx6 -Wf'-pvctl nodivloop noloopchg -O'
      <br>
      (致命的レベルの文法診断メッセージのみを出力.
      SX-6向けの命令列を生成.
      ループ分割によるベクトル化を行わないことを指定.
      ループ入れ換えによるベクトル化を行わないことを指定.
      最適化処理を指示すると共に, 最適化処理に関する個別処理を指定)
    </td>
  </tr>

  <tr>
    <td width="20%">実行時間等に関する情報</td>
    <td>
      <a href="information/sx_comp_info.txt">information/sx_comp_info.txt</a>
      <br>
      <a href="information/dycore_clock.txt">information/dycore_clock.txt</a>
    </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="/library/dcpam/dcpam2/">Version: dcpam2-20050623-4</a>
    </td>
  </tr>
  <tr>
    <td width="20%">初期値生成の NAMELIST</td>
    <td>
      <a href="information/init.nml">information/init.nml</a>
    </td>
  </tr>
  <tr>
    <td width="20%">長期ランの NAMELIST</td>
    <td>
      <a href="information/dycore.nml">information/dycore.nml</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>240.0 [s]</td>
  </tr>
  <tr>
    <td>積分時間</td>
    <td>1200 [day]</td>
  </tr>
  <tr>
    <td>境界条件</td>
    <td>
      水平(λ,φ) : 周期境界条件<br>
      鉛直(σ)    : dσ/dt = 0 (σ = 0, 1)
    </td>
  </tr>
  <tr>
    <td>数値解法</td>
    <td>
      水平(λ,φ) : スペクトル法 (球面調和関数展開・三角形切断)<br>
      鉛直(σ)    : 差分法
                    (<a href="#References">Arakawa and Suarez (1983)</a>)<br>
      時間(t)     : リープフロッグスキーム 
                    (<a href="#References">Asselin (1972)</a> のフィルター係数 : 0.05)
    </td>
  </tr>
  <tr><td colspan="2">&nbsp;<td></tr>

  <tr>
    <td>空間解像度</td>
    <td>
      東西格子点数 : 192<br>
      南北格子点数 : 96<br>
      鉛直格子点数 : 20
    </td>
  </tr>
  <tr>
    <td>最大波数</td>
    <td>63</td>
  </tr>
  <tr><td colspan="2">&nbsp;<td></tr>
  
  <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>
  <tr>
    <td>初期温度</td>
    <td>温度一様(250.00 [K])</td>
  </tr>
  <tr>
    <td>初期気圧</td>
    <td>気圧一様(1000.00 [hPa])</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.
=end
}
message:m,rd:{
== 結果
}
#
# Following parts are Thumbnail comments
#
Ps_GlobalMean_1-1200Day.png::
Temp_LonMean_200-1200DayMean.png::
U_LonMean_200-1200DayMean.png::
V_LonMean_200-1200DayMean.png::