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

* 森川 靖大
  * 2005/07/07 (森川 靖大) 作成

== 概要
下記の環境, およびパラメータで ((<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>
      SUPER-UX unix 12.1  SX-6
    </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>
      東西格子点数 : 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:m:{
<hr> 東西風 <br>
<a href="nc/result_vlon.nc">result_vlon.nc</a>
(<a href="figdir/HS94_T42L20_TempIni300_SX6_VelLon_DayMean200-1200_LonMean.ps">PS file</a>)
<hr>
}
label:m:{
<hr> &nbsp; <br>
(<a href="figdir/HS94_T42L20_TempIni300_SX6_VelLon_DayMean200-1200_LonMean.ps">PS file</a>)
<hr>
}
label:m:{
<hr> 南北風 <br>
<a href="nc/result_vlat.nc">result_vlat.nc</a> 
(<a href="figdir/HS94_T42L20_TempIni300_SX6_VelLat_DayMean200-1200_LonMean.ps">PS file</a>)
<hr>
}
VelLon_DayMean200-1200_LonMean.png:m:{
<br>
<div align="left">
  200 日目から 1250 日目までを平均した、帯状平均東西風。
  <a href="#References">Held and Suarez(1994)</a> で示される T63 の
  結果とかなり近い結果が出ている。
</div>
}
VelLon_Day0-1200_LonMean_Sigma0.2.png:m:{
<br>
<div align="left">
  帯状平均東西風の日変化。σ = 0.2。
  東西ジェットの最大値の軸が南北の移動していることが分かる、
  詳細な解析はおこなっていないが、
  平均的すれば南北対象であることは左の図からわかる。
</div>
}
VelLat_DayMean200-1200_LonMean.png:m:{
<br>
<div align="left">
  200 日目から 1200 日目までを平均した、帯状平均南北風。
</div>
}
##################################################
label:m:{
<hr> 鉛直風 <br>
<a href="nc/result_vsigmahalf.nc">result_vsigmahalf.nc</a>
(<a href="figdir/HS94_T42L20_TempIni300_SX6_VelSigma_DayMean200-1200_LonMean.ps">PS file</a>)
<hr>
}
label:m:{
<hr> 温度 <br>
<a href="nc/result_temp.nc">result_temp.nc</a>
(<a href="figdir/HS94_T42L20_TempIni300_SX6_Temp_DayMean200-1200_LonMean.ps">PS file</a>)
<hr>
}
label:m:{
<hr> 大気の全質量 <br>
<a href="nc/result_tm.nc">result_tm.nc</a>
(<a href="figdir/HS94_T42L20_TempIni300_SX6_TotalMass_Day0-1200.ps">PS file</a>)
<hr>
}
VelSigma_DayMean200-1200_LonMean.png:m:{
<br>
<div align="left">
  200 日目から 1200 日目までを平均した、帯状平均鉛直風。
  σ = 1 が地上、σ = 0 が大気上端なので、
  σ &gt; 0 は下降流、 σ &lt; 0 は上昇流である。
</div>
}
Temp_DayMean200-1200_LonMean.png:m:{
<br>
<div align="left">
  200 日目から 1200 日目までを平均した、帯状平均温度。
</div>
}
TotalMass_Day0-1200.png:m:{
<br>
<div align="left">
  大気の全質量の時間変化。わずかに増加している。
</div>
}