#!/usr/bin/env ruby1.8 # #= dcmodel thumbnail generate ruby script # # Editor :: murashin # Version:: 2008/01/27 19:33:32 # #== Overview # #This file is generate by following ruby script automatically. # # /home/murashin/test2/thum-src/dcmodel-thum.rb # version "$Revision: 1.31 $ : $Date: 2007/05/22 05:13:14 $" # #Please edit this file according to your purpose. # #== Usage # #Please check following sample page and reference manual. # # http://www.gfd-dennou.org/library/dcmodel/doc/dcmodel-tools/dcmodel-thum-sample # http://www.gfd-dennou.org/library/dcmodel/doc/dcmodel-tools/dcmodel-thum-rdoc # ################################################## require "/home/murashin/test2/thum-src/dcmodel-thum.rb" ###################################################### # DCModelThumbnail のインスタンスを作成 (必須) thumb = DCModelThumbnail.new # サムネイル実行コマンド (必須)。最後に実行すること END{ thumb.create } #thumb.copyright = "Shin'ya Murakami" # コピーライト thumb.index = "../thumb" # 作成されるサムネイルのファイル名 (拡張子除く) #thumb.index_ext = ".htm" # 作成されるサムネイルの拡張子名 thumb.infofile = "#{File.basename(thumb.index)}.txt" # 「情報ファイル (infofile)」のファイル名 #thumb.ext_list.push("bmp") # 画像ファイルとして認識させたい拡張子の # 追加。デフォルトでは # gif, png, jpg, jpeg, # GIF, PNG, JPG, JPEG # が認識される。 #thumb.headlimit = "headlimit_" # headlimit に当てはまらないものは画像と # しては無視され、当てはまるものは「情報 # ファイル」において、接頭部分 headlimit # を省略して指定することが可能となる。 thumb.blankfig = true # 画像ファイル名を表示しない #thumb.figdir = "../figdir" # 絵のあるディレクトリの名前。 # 必ず存在していなければならない。 #thumb.thumbnaildir = "../thumbdir" # サムネイル画像を置くディレクトリの名前 #thumb.thumbnailtail = "_thumb.png" # サムネイル画像の拡張子名 #thumb.thumbnailpage = "1" # サムネイル画像化する元画像のページ番号 #thumb.convert_cmd = "convert -depth 8 -geometry" # convert コマンドおよびオプション。この # 後ろに画像サイズが 200x150 のような形 # 式で指定される。 #thumb.convert_overwrite = false # convert コマンドが実行される際、既に出 # 力先のファイルがある場合に上書きするた # めのフラグ。これが false の場合、出力 # 先のファイルが存在し、且つ元のファイル # よりも新しい場合は変換を行わない。 #thumb.css = "/GFD_Dennou_Club/ftp/arch/dcmodel/htmltools/dcmodel.css" # スタイルシートファイル #thumb.rd2_path = "/usr/bin/rd2" # rd2 コマンドへのパス #thumb.rd2htmlextlib = "/GFD_Dennou_Club/ftp/arch/dcmodel/lib/ruby/1.8" # rd2-ext-lib へのライブラリへのパス thumb.mksigen = false # SIGEN ファイルを作らない場合は false # にセットする thumb.norobots = true # HTML のロボット検索を禁止する場合には # true にセットする #thumb.img_width = 200 thumb.img_width = 160 # 画像ファイルサイズ (幅) #thumb.img_height = 150 thumb.img_height = 120 # 画像ファイルサイズ (高さ) thumb.figtable_num = 4 # 横にならべるファイル数 #thumb.style = <<-STYLE # a:link { color:#269900; } # a:visited { color:#269900; } # a:hover { color:#99FF33; } #STYLE # スタイルシートを直接設定するためのもの # (css ファイルに書き込む情報を直接指定 # できる) #thumb.html_author = "murashin" # html の作成者情報 (デフォルトはユーザ # アカウント名が自動取得される) thumb.title = "simulation results on 2-d beta-plane turbulence" # html ヘッダのタイトル # 本体に書き出すメッセージ。サムネイルの # 部分よりも上に出力される。この変数自体 # は Array オブジェクトで、その内部に # String オブジェクトが格納される。 # フッターメッセージ。 # "=" ではなく、"<<" で代入することに注意！！ thumb.footer = Array.new #thumb.footer << <<-Footer # #Footer # メッセージ。 # "=" ではなく、"<<" で代入することに注意！！ thumb.message = Array.new thumb.message << <<-Message =begin == 実験設定 ベータ面上の順圧渦度方程式をベータを変化させて数値計算した. どの実験も全波数が 98以上102以下の波数成分をマルコフ強制した. 高波数と低波数にカスケードしたエネルギーとエンストロフィーはそれぞれ超粘性とhypofrictionによって散逸されるように設定した. beta = 0, 5, 10, 15, 20, 25, 30 でそれぞれ統計的平衡状態を実現し, 統計的平衡状態における物理量の振舞を観察した. * 空間格子点数 512^2, 切断波数=170, 空間 [2π×2π] * 時間刻幅 dt = 1.0*10^(-3) * 数値計算法: * 空間: 変換法を使って非線形項を評価するスペクトル法 * 時間: 3次 Adams-Bashforth * integrating factorを用いて硬い微分方程式を変形して解いた. * マルコフ強制のパラメタ * 振幅 * 相関時間 * 散逸項パラメタ * 超粘性係数: 20/(170^(8)) * 超粘性次数: 4 * hypofriction 係数: 100 * hypofriction 次数: 2 * 実験名とベータと統計的平衡状態として用いた期間の表 * exp_1_1: beta=0, t=100からt=600 * exp_1_2: beta=5, t=100からt=600 * exp_1_3: beta=10, t=100からt=600 * exp_1_4: beta=15, t=200からt=800 * exp_1_5: beta=20, t=300からt=1000 * exp_1_6: beta=25, t=300からt=1000 * exp_1_7: beta=30, t=300からt=1300 == 粗視化の詳細 === 波数空間を同心円に囲まれた領域に分割 波数空間を全波数 k について k_0=1, λ=2^(1/2)として k_0*λ^i <= k < k_0*λ^(i+1), i=0,1,2,..,n を満たす k = (k_x^2 + k_y^2)^(1/2) の集合を k_i で表す. k_i の集合を粗視化領域 i と呼ぶ. 今の場合は具体的には, * 領域 0: 1 <= k < 1.41 * 領域 1: 1.41 <= k < 2 * 領域 2: 2 <= k < 2.83 * 領域 3: 2.83 <= k < 4 * 領域 4: 4 <= k < 5.66 * 領域 5: 5.66 <= k < 8 * 領域 6: 8 <= k < 11.3 * 領域 7: 11.3 <= k < 16 * 領域 8: 16 <= k < 22.6 * 領域 9: 22.6 <= k < 32 * 領域 10: 32 <= k < 45.3 * 領域 11: 45.3 <= k < 64 * 領域 12: 64 <= k < 90.5 * 領域 13: 90.5 <= k < 128 * 領域 14: 128 <= k < 181 * 領域 15: 181 <= k < 256 * 領域 15: 256 <= k < 362 となる. === 波数空間を線形Rossby波の振動数の大きさの範囲で分割 線形Rossby波の振動数を omega で表す. |omega/beta| の最大値は (k_x, k_y)= (1,0) で 1 になる. 最小値は k_x = 0 で 0 になる. そこで, omega(0, k_y)なる波数ベクトルの集合を omega_0 と定義する. また, λ = 2^(-1/2) として, λ^(i+1) < |omega(k_x, k_y)/beta| <= λ^i, i=1,2,..,n を満たす (k_x, k_y) の集合を領域 i と定義する. i=0を除いて, iが増えるに従って振動数が小さい領域を指すことになる. 具体的には * 領域 0: |omega/beta| == 0 * 領域 1: 0.707 < |omega/beta| <= 1 * 領域 2: 0.5 < |omega/beta| <= 0.707 * 領域 3: 0.354 < |omega/beta| <= 0.5 などとなる. == 結果 * ((<渦度場>)) * ((<帯状平均東西速度 (統計的平衡状態, 時間平均)>)) * ((<帯状平均東西速度 (統計的平衡状態, スナップショット)>)) * ((<帯状平均渦度 (統計的平衡状態, 時間平均)>)) * ((<帯状平均渦度 (統計的平衡状態, スナップショット)>)) * ((<帯状平均渦度の南北微分 (統計的平衡状態, スナップショット)>)) * ((<帯状平均渦度の南北微分 (統計的平衡状態, スナップショット)>)) * ((<エネルギーの時系列, E(t)>)) * ((<エンストロフィーの時系列, Z(t)>)) * ((<パリンストロフィーの時系列, P(t)>)) * ((<帯状平均エネルギーの時系列, E_z(t)>)) * ((<帯状平均エンストロフィーの時系列, Z_z(t)>)) * ((<帯状平均パリンストロフィーの時系列, P_z(t)>)) * ((<エネルギー散逸率の時系列>)) * ((<エンストロフィー散逸率の時系列>)) * ((<パリンストロフィー散逸率の時系列>)) * ((<エネルギースペクトル (統計的平衡状態, 時間平均)>)) * ((<エンストロフィースペクトル (統計的平衡状態, 時間平均)>)) * ((<パリンストロフィースペクトル (統計的平衡状態, 時間平均)>)) * ((<1次元エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<1次元エンストロフィー輸送 (統計的平衡状態, 時間平均)>)) * ((<1次元パリンストロフィー輸送 (統計的平衡状態, 時間平均)>)) * ((<エネルギーフラックス (統計的平衡状態, 時間平均)>)) * ((<エンストロフィーフラックス (統計的平衡状態, 時間平均)>)) * ((<パリンストロフィーフラックス (統計的平衡状態, 時間平均)>)) * ((<2次元エネルギースペクトル (統計的平衡状態, スナップショット)>)) * ((<2次元エネルギースペクトル (統計的平衡状態, スナップショット, 拡大その1)>)) * ((<2次元エネルギースペクトル (統計的平衡状態, スナップショット, 拡大その2)>)) * ((<2次元エネルギースペクトル (統計的平衡状態, スナップショット, 拡大その3)>)) * ((<2次元エンストロフィースペクトル (統計的平衡状態, スナップショット)>)) * ((<2次元エンストロフィースペクトル (統計的平衡状態, スナップショット, 拡大その1)>)) * ((<2次元エンストロフィースペクトル (統計的平衡状態, スナップショット, 拡大その2)>)) * ((<2次元エンストロフィースペクトル (統計的平衡状態, スナップショット, 拡大その3)>)) * ((<2次元パリンストロフィースペクトル (統計的平衡状態, スナップショット)>)) * ((<2次元パリンストロフィースペクトル (統計的平衡状態, スナップショット, 拡大その1)>)) * ((<2次元パリンストロフィースペクトル (統計的平衡状態, スナップショット, 拡大その2)>)) * ((<2次元パリンストロフィースペクトル (統計的平衡状態, スナップショット, 拡大その3)>)) * ((<2次元エネルギースペクトル (統計的平衡状態, 時間平均)>)) * ((<2次元エネルギースペクトル (統計的平衡状態, 時間平均, 拡大その1)>)) * ((<2次元エネルギースペクトル (統計的平衡状態, 時間平均, 拡大その2)>)) * ((<2次元エネルギースペクトル (統計的平衡状態, 時間平均, 拡大その3)>)) * ((<2次元エンストロフィースペクトル (統計的平衡状態, 時間平均)>)) * ((<2次元エンストロフィースペクトル (統計的平衡状態, 時間平均, 拡大その1)>)) * ((<2次元エンストロフィースペクトル (統計的平衡状態, 時間平均, 拡大その2)>)) * ((<2次元エンストロフィースペクトル (統計的平衡状態, 時間平均, 拡大その3)>)) * ((<2次元パリンストロフィースペクトル (統計的平衡状態, 時間平均)>)) * ((<2次元パリンストロフィースペクトル (統計的平衡状態, 時間平均, 拡大その1)>)) * ((<2次元パリンストロフィースペクトル (統計的平衡状態, 時間平均, 拡大その2)>)) * ((<2次元パリンストロフィースペクトル (統計的平衡状態, 時間平均, 拡大その3)>)) * ((<2次元エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<2次元エネルギー輸送 (統計的平衡状態, 時間平均, 拡大その1)>)) * ((<2次元エネルギー輸送 (統計的平衡状態, 時間平均, 拡大その2)>)) * ((<2次元エネルギー輸送 (統計的平衡状態, 時間平均, 拡大その3)>)) * ((<2次元エンストロフィー輸送 (統計的平衡状態, 時間平均)>)) * ((<2次元エンストロフィー輸送 (統計的平衡状態, 時間平均, 拡大その1)>)) * ((<2次元エンストロフィー輸送 (統計的平衡状態, 時間平均, 拡大その2)>)) * ((<2次元エンストロフィー輸送 (統計的平衡状態, 時間平均, 拡大その3)>)) * ((<粗視化領域 0 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 1 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 2 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 3 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 4 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 5 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 6 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 7 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 8 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 9 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 10 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 11 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 12 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 13 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 14 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 15 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<粗視化領域 16 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 0 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 1 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 2 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 3 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 4 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 5 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 6 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 7 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 8 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 9 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 10 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 11 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 12 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 13 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 14 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 15 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 16 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 17 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) * ((<線形Rossby波の振動数で分類した, 粗視化領域 18 への粗視化エネルギー輸送 (統計的平衡状態, 時間平均)>)) =end Message