message:l,m:{
<small>
  [ <a href="/">地球流体電脳倶楽部</a> /
    <a href="/library/dcpam">DCPAM</a> /
    <a href="/library/dcpam/sample">計算結果</a> /
    <a href="SIGEN.htm">SIGEN.htm</a> ] 
</small>
}
message:l,m,rd:{
= 水惑星実験 (H1998con)

#  * 2009/01/20 (森川 靖大) 作成

== 概要

((<AGCM 5.3 による水惑星実験 (H1998con)|URL:/arch/ape/AGCM5/old.2005/#adj>))
と同設定 (調査中) での計算での dcpam5 (20090108-1) の動作試験である. 
}

message:l,m,rd:{
現在動作試験中につき, 一部の画像しか作成していない. 
}

message:l,m,rd:{

計算結果は
#((<初期値および実験結果>))
((<実験結果>))
に記す. 

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

  <tr>
    <td width="20%">依存するライブラリとバージョン</td>
    <td>
      <ul>
	<li> netCDF ver 3.6.1
	<li> gtool5 ver 20081109
	<li> ISPACK ver 0.71
	<li> spml ver 0.4.2
      </ul>
    </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/src">ソースコード</a>
      [<a href="info00/src.tar.gz">TAR.GZ</a>] <br>
      <a href="info00/Config.mk">ビルド時の設定ファイル</a>
    </td>
  </tr>

  <tr>
    <td width="20%">NAMELIST ファイル</td>
    <td>
      <a href="info00/dcpam_ape_T42L16_rst1000.nml">info00/dcpam_ape_T42L16_rst1000.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>20 分</td>
  </tr>
  <tr>
    <td>積分時間</td>
    <td>1000 日</td>
  </tr>
  <tr>
    <td>境界条件</td>
    <td>
      <ul>
        <li>水平(λ,φ) : 周期境界条件<br>
        <li>鉛直(σ)    : dσ/dt = 0 (σ = 0, 1)
      </ul>
    </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>

  <tr>
    <td>空間解像度</td>
    <td>
      東西格子点数 : 128<br>
      南北格子点数 : 63<br>
      鉛直格子点数 : 16
    </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 [s-1] </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> 10800.0  [s] </td>
  </tr>
  <tr>
    <td>超粘性の次数</td>
    <td>8</td>
  </tr>

  <tr><td colspan="2">&nbsp;</td></tr>
  <tr>
    <td>力学過程</td>
    <td>プリミティブ方程式系. 詳しくは
        <a href="/library/dcpam/dcpam5/dcpam5_current/doc/basic_equations/htm/basic_equations.htm">
        dcpam5  支配方程式系とその離散化 
        </a>
        を参照のこと. 
    </td> 
  </tr>
  <tr>
    <td width="20%">物理過程</td>
    <td>
      <ul>
      	<li> 放射 (バンドモデル)
	    <ul>
	      <li> 長波放射 (水蒸気 3 バンド, 乾燥大気 1 バンド)
	      <li> 短波放射吸収なし
	      <!-- <li> 短波放射 (水蒸気 1 バンド) -->
	    </ul>
      	<li> 鉛直拡散 (Mellor and Yamada, 1974, レベル 2)
        <li> 湿潤対流調節 (Manabe et al., 1965)
        <li> 診断型大規模凝結 
        <li> 乾燥対流調節 
      </ul>
    </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/stdout.txt">info00/stdout.txt</a>
    </td>
  </tr>

  <tr>
    <td width="20%">標準エラー出力</td>
    <td>
      <a href="info00/stderr.txt">info00/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.
* Manabe, S., Smagorinsky, J., Strickler, R.F., 1965:
  Simulated climatology of a general circulation model
  with a hydrologic cycle.
  ((|Mon. Weather Rev.|)), ((*93*)), 769--798.
* Mellor, G. L, Yamada, T., 1974:
  A hierarchy of turbulence closure models for planetary boundary layers.
  ((|J. Atmos. Sci.|)), ((*31*)), 1791--1806.
#* 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>
##label:m:{
##<hr>
##初期値
##<hr>
##}
##label:: <hr>
##::
##::
##::
###
##init_Temp_LonLatMean.png:: 温度の基本場の高度分布
##init_QVap_LonLatMean.png:: 比湿の高度分布
##::
##::
##::
##::
###
label:: <hr>
label:m:{
<hr>
ML
<hr>
}
label:: <hr>
#
ml/ml_u.png::ml_u <br /> eastward wind
ml/ml_v.png::ml_v <br /> northward wind
ml/ml_t.png::ml_t <br /> air temperature
#
label::ml_om <br /> omega
label::ml_phi <br /> geopotential height
ml/ml_q.png::ml_t <br />  specific humidity

label:: <hr>
label:m:{
<hr>
TR
<hr>
}
label:: <hr>
#
tr/tr_u250.png::tr_u250 <br /> eastward wind
tr/tr_v250.png::tr_v250 <br /> northward wind
tr/tr_mslp.png::tr_mspl <br /> air pressure at sea level