message:l,m:{ [ 地球流体電脳倶楽部 / DCPAM / 計算結果 / SIGEN.htm ] } message:l,m,rd:{ = 木星大気を想定した仮想惑星実験のための予備的実験 #* 森川 靖大 # * 2008/03/13 (森川 靖大) 作成 == 概要 ※ 書きかけ ((:
:))((:
:)) 木星大気の 2 次元計算をおこなった Sugiyama et al. (2008) を参考に, 木星大気を想定した仮想惑星実験のための予備的実験を行った. } message:l,m,rd:{ 環境, およびパラメータは以下の通り. 結果は下記に記す. == 実験環境 } message:l,m:{
システム (uname -a) SUPER-UX unix 14.1 SX-6
コンパイラ FORTRAN90/SX
依存するライブラリとバージョン
  • netCDF ver 3.6.1
  • gt4f90io ver 20080220
  • ISPACK ver 0.71
  • spml ver 0.4.2
} message:l,m,rd:{ == 実験設定 } message:l,m:{
 
計算領域 水平(λ,φ) : 全球 (λ=[0, 360], φ=[-90, 90])
鉛直(σ=p/p0) : 大気下端〜上端 (σ = [0, 1]), 大気下端は 3.0e+6 Pa
時間間隔 10 分
積分時間 300 日
境界条件 水平(λ,φ) : 周期境界条件
鉛直(σ) : dσ/dt = 0 (σ = 0, 1)
数値解法
  • 水平(λ,φ) : スペクトル法 (球面調和関数展開・三角形切断)
  • 鉛直(σ) : 差分法 (Arakawa and Suarez, 1983)
  • 時間(t) : リープフロッグスキーム.
    • ただし重力波の項にはセミインプリシット法を用いる. 時間フィルター (Asselin, 1972). を用いる. フィルター係数は 0.05
 
空間解像度 東西格子点数 : 64
南北格子点数 : 32
鉛直格子点数 : 24
最大波数 21
 
惑星半径 6.371e+6 [m] (地球の値)
回転角速度 7.292e-5 [s-1] (地球の値)
重力加速度 23.1 [m s-2] (Sugiyama et al (2008) より)
大気定圧比熱 11900.9264 [J K-1 kg-1] (Sugiyama et al (2008) より)
大気気体定数 3611.44466 [J K^-1 kg^-1] (Sugiyama et al (2008) より)
最大波数に対する e-folding time 10800.0 [s]
超粘性の次数 8
} message:l,m,rd:{ == 実行結果に関する付加情報 } #message:l,m:{ # # # # # # # # # # # # # # # #
標準出力 # info00/agcm5.t42l20.hs94_stdout.txt #
標準エラー出力 # info00/agcm5.t42l20.hs94_stderr.txt #
プロファイラによる速度計測 # info00/ftrace.txt #
#} message:: #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 s-1] label::南北風速 [m s-1] label::温度 [K] label::比湿 [kg kg-1] label::鉛直速度 [s-1] label:m:{
σ=1〜0 での
東西平均,
200〜300 日
時間平均値
} U_LonMean_200-300DayMean_1-0Sig.png:: V_LonMean_200-300DayMean_1-0Sig.png:: Temp_LonMean_200-300DayMean_1-0Sig.png:: QVap_LonMean_200-300DayMean_1-0Sig.png:: SigmaDot_LonMean_200-300DayMean_1-0Sig.png:: label:m:{
σ=1〜0.003 での
東西平均,
200〜300 日
時間平均値
} U_LonMean_200-300DayMean_1-0.003Sig.png:: V_LonMean_200-300DayMean_1-0.003Sig.png:: Temp_LonMean_200-300DayMean_1-0.003Sig.png:: QVap_LonMean_200-300DayMean_1-0.003Sig.png:: SigmaDot_LonMean_200-300DayMean_1-0.003Sig.png:: label:m:{
σ=0.003 での
200〜300 日
時間平均値
} U_0.003Sig_200-300DayMean.png:: V_0.003Sig_200-300DayMean.png:: Temp_0.003Sig_200-300DayMean.png:: QVap_0.003Sig_200-300DayMean.png:: SigmaDot_0.003Sig_200-300DayMean.png:: label:m:{
σ=0.003 での
東西平均, 
200〜300 日
時間平均値
} U_0.003Sig_LonMean_200-300DayMean.png:: V_0.003Sig_LonMean_200-300DayMean.png:: Temp_0.003Sig_LonMean_200-300DayMean.png:: QVap_0.003Sig_LonMean_200-300DayMean.png:: SigmaDot_0.003Sig_LonMean_200-300DayMean.png:: label:m:{
σ=0.1 での
200〜300 日
時間平均値
} U_0.1Sig_200-300DayMean.png:: V_0.1Sig_200-300DayMean.png:: Temp_0.1Sig_200-300DayMean.png:: SigmaDot_0.1Sig_200-300DayMean.png:: QVap_0.1Sig_200-300DayMean.png:: label:m:{
σ=0.2 での
200〜300 日
時間平均値
} U_0.2Sig_200-300DayMean.png:: V_0.2Sig_200-300DayMean.png:: Temp_0.2Sig_200-300DayMean.png:: SigmaDot_0.2Sig_200-300DayMean.png:: QVap_0.2Sig_200-300DayMean.png:: label::
label::
label::
label::
label::
label::
label:: Rain_GlobalMean_0-300Day.png:m:{ 降水量の全球平均値の時間変化. } label:: Temp_LonLatMean_0-300Day.png:m:{ 温度の緯度経度平均値の時間変化. } QVap_LonLatMean_0-300Day.png:m:{ 比湿の緯度経度平均値の時間変化. } :: :: :: :: PotTemp_LonLatMean_1-0.01Sig_0-300Day.png:m:{ 温位の緯度経度平均値の時間変化. } :: ::