message:l,m,rd:{
= 解像度 T42L20：初期値温度場 300K の場合

* 森川 靖大
  * 2005/07/07 (森川 靖大) 実験環境について追記
#  * 2005/05/01 (森川 靖大) 更新
  * 2005/04/30 (森川 靖大) ページ作成

== 概要
下記の環境, およびパラメータで ((<Held and Suarez (1994) |URL:#References>))
実験を行った結果を記す。

== 実験環境
}
message:l,m:{
<table border="1" cellspacing="1" cellpadding="5" align="center">
  <tr>
    <td width="20%">システム (uname -a)</td>
    <td>
      Linux leaf.ep.sci.hokudai.ac.jp 2.4.20-31.9smp #1 SMP Tue Apr 13 17:40:10 EDT 2004 i686 i686 i386 GNU/Linux
    </td>
  </tr>
  
  <tr>
    <td width="20%">コンパイラ</td>
    <td>
      Fujitsu Fortran Compiler バージョン 5.0
    </td>
  </tr>

  <tr>
    <td width="20%">コンパイラオプション</td>
    <td>
      -Am -Kfast -E fitu  (モジュール読み込み, 自動高速化, ループ化などの表示)
    </td>
  </tr>

  <tr>
    <td width="20%">実行時間等に関する情報</td>
    <td>
      <a href="nc/clock.txt">nc/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="src/">dcpam2 (2005/02/14)</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>1250 [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>数値解法</td>
    <td>
      <dl>
        <dt>水平(λ,φ)：</dt>
	<dd>スペクトル法 (球面調和関数展開・三角形切断)</dd>
	
        <dt>鉛直(σ)：</dt>
	<dd>差分法 (<a href="#References">Arakawa and Suarez (1983)</a>)</dd>
	
	<dt>時間(t)：</dt>
	<dd>リープフロッグスキーム</dd>
        <dd>(<a href="#References">Asselin (1972)</a> のフィルター係数 : 0.05)</dd>
      </dl>
    </td>
  </tr>
-->
  <tr><td colspan="2">&nbsp;<td></tr>

  <tr>
    <td>空間解像度</td>
    <td>
      東西格子点数 : 128<br>
      南北格子点数 : 64<br>
      鉛直格子点数 : 20
    </td>
  </tr>
  <tr>
    <td>最大波数</td>
    <td>42</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>温度一様(300.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
#
label:: <hr> 東西風 <br> <a href="nc/result_vlon.nc">result_vlon.nc</a> <hr>
label:: <hr>
label:: <hr>
VelLon_DayMean200-1250_LonMean.png:m:{
<br>
<div align="left">
  200 日目から 1250 日目までを平均した、帯状平均東西風。
  <a href="#References">Held and Suarez(1994)</a> で示される T63 の
  結果とかなり近い結果が出ている。
</div>
}
VelLon_Day200-1250_LonMean_Sigma0.2.png:m:{
<br>
<div align="left">
  帯状平均東西風の日変化。σ = 0.2。
  東西ジェットの最大値の軸が南北の移動していることが分かる、
  詳細な解析はおこなっていないが、
  平均的すれば南北対象であることは左の図からわかる。
</div>
}
:
label:: <hr> 南北風 <br> <a href="nc/result_vlat.nc">result_vlat.nc</a> <hr>
label:: <hr> 鉛直風 <br> <a href="nc/result_vsigmahalf.nc">result_sigmahalf.nc</a> <hr>
label:: <hr> 温度 <br> <a href="nc/result_temp.nc">result_temp.nc</a> <hr>
VelLat_DayMean200-1250_LonMean.png:m:{
<br>
<div align="left">
  200 日目から 1250 日目までを平均した、帯状平均南北風。
</div>
}
VelSigma_DayMean200-1250_LonMean.png:m:{
<br>
<div align="left">
  200 日目から 1250 日目までを平均した、帯状平均鉛直風。
  σ = 1 が地上、σ = 0 が大気上端なので、
  σ &gt; 0 は下降流、 σ &lt; 0 は上昇流である。
</div>
}
Temp_DayMean200-1250_LonMean.png:m:{
<br>
<div align="left">
  200 日目から 1250 日目までを平均した、帯状平均温度。
</div>
}