#!/usr/bin/env ruby1.8
#
#= dcmodel thumnail generate ruby script
#
#  Editor :: kitamo
#  Version:: 2005/09/19 15:30:20
#
#== Overview
#
#This file is generate by following ruby script automatically.
#
#      /home/kitamo/work/zatsumu/dcmodel-tools/dcmodel-thum.rb
#      version "$Revision: 1.14 $ : $Date: 2005/05/25 15:54:32 $"
#
#Please edit this file according to your purpose.
#
#== Usage
#
#Please check following sample page and reference manual.
#
#  http://www.gfd-dennou.org/arch/dcmodel/doc/dcmodel-tools/dcmodel-thum-sample
#  http://www.gfd-dennou.org/arch/dcmodel/doc/dcmodel-tools/dcmodel-thum-rdoc
#
##################################################

require "/home/kitamo/util/dcmodel-tools/dcmodel-thum.rb"

######################################################
# DCModelThumbnail のインスタンスを作成 (必須)
thumb = DCModelThumbnail.new

# サムネイル実行コマンド (必須)。最後に実行すること
END{
  thumb.create
}

# コピーライト
#thumb.copyright = "GFD Dennou Club"

# 作成されるサムネイルのファイル名 (拡張子除く)
thumb.index     = "../odaka1998-20051008"

# 作成されるサムネイルの拡張子名
#thumb.index_ext = ".htm"

# 「情報ファイル (infofile)」のファイル名
thumb.infofile  = "#{File.basename(thumb.index)}.txt"

# 画像ファイルとして認識させたい拡張子の追加。
# デフォルトでは gif, png, jpg, jpeg, GIF, PNG, JPG, JPEG が認識される。
#thumb.ext_list.push("bmp")

# headlimit に当てはまらないものは画像としては無視され、
# 当てはまるものは「情報ファイル」において、接頭部分 headlimit が
# 省略可能になる。
#thumb.headlimit = "headlimit_"

# 絵のあるディレクトリの名前。必ず存在していなければならない。
thumb.figdir    = "../fig"

# スタイルシートファイル
#thumb.css       = "/GFD_Dennou_Club/ftp/arch/dcmodel/htmltools/dcmodel.css"

# rd2 コマンドへのパス
#thumb.rd2_path  = "/usr/bin/rd2"

# rd2-ext-lib へのライブラリへのパス
#thumb.rd2htmlextlib = "/GFD_Dennou_Club/ftp/arch/dcmodel/lib"

# SIGEN ファイルを作らない場合は false にセットする
#thumb.mksigen = false

# HTML のロボット検索を禁止する場合には true にセットする
#thumb.norobots = true

# 画像ファイルサイズ (幅)
thumb.img_width  = 400

# 画像ファイルサイズ (高さ)
thumb.img_height = 300

# 横にならべるファイル数
thumb.figtable_num = 2

# スタイルシートを直接設定するためのもの
# (css ファイルに書き込む情報を直接指定できる)
#thumb.style = <<-STYLE
#  a:link { color:#269900; }
#  a:visited { color:#269900; }
#  a:hover { color:#99FF33; }
#STYLE

# html の作成者情報  (デフォルトはユーザアカウント名が自動取得される)
#thumb.html_author  = "kitamo"

# html ヘッダのタイトル
thumb.title     = "dcmodel-thum : 実験結果サムネイル簡易作成スクリプト"

# 本体に書き出すメッセージ。サムネイルの部分よりも上に出力される。
# この変数自体は Array オブジェクトで、その内部に String オブジェクトが
# 格納される。
thumb.message   = Array.new

# "=" ではなく、"<<" で代入することに注意！！
thumb.message   << <<-Message
=begin

= Odaka et al. (1998) の計算

arare3m を用いて Odaka et al. (1998) の計算を行った. 
静止した安定な大気層の最下層に温度擾乱を与えた時の流れ場の様子を 6 モデル時間計算している. 
物理過程にはオリジナルと同じものを用いている. 
すなわち, 放射過程として高度 5 km 以下の領域に一様冷却を与え, 
地表面過程として地表面温度およびバルク係数一定のバルク法を用いた. 

今までのプログラムと比べて異なる点は
乱流拡散係数の計算部分で, 浮力に伴う乱流拡散係数の変化を計算する際, 
基本場の温位の勾配を考慮するようにしたことである. 
この修正に伴い, ((<2005-09-17|URL:../2005-09-17_kitamo/odaka1998-20050917.htm>))など
で示されている, 基本場の温位の拡散を考慮すると温位分布が
階段状になるという問題が解決された. 

計算は Dual Pentium III 1 Ghz を搭載した PC で行った. 
計算に要した実行時間はおよそ 1 時間半である. 

== 計算設定

=end
Message

thumb.message   << <<-Message

    <table BORDER=\"0\" cellspacing=\"10\" align=\"center\">
    <tr>
    <td>ソースプログラム</td>
    <td><a href=\"./src/\">arare3m-20051005</a>
    </tr>
    <tr>
    <td>基本場の温度</td>
    <td>高度 5 km 以下は温位 245 K, 高度 5 km 以上では温度 220 K</td>
    </tr>
    <tr>
    <td>基本場の圧力</td>
    <td>静水圧平衡 (地表面で 700 Pa)</td>
    </tr>
    <tr>
    <td>温度の初期擾乱</td>
    <td>最下層に最大振幅 2 K のランダムな擾乱を与える</td>
    </tr>
    <tr>
    <td>圧力の初期擾乱</td>
    <td>なし</td>
    </tr>
    <tr>
    <td>風速の初期値</td>
    <td>u = 0 m/s , w = 0 m/s</td>
    </tr>
    <tr>
    <td>放射加熱</td>
    <td>高度 5 km 以下で -50 K /day の一様冷却</td>
    </tr>
    <tr>
    <td>地表面フラックス</td>
    <td>バルク法(バルク係数は熱, 運動量ともに 0.01 で一定, 地表面温度 270 K で固定)</td>
    </tr>
    <tr>
    <td>重力加速度</td>
    <td>3.72 m/s<sup>2</sup></td>
    </tr>
    <td>粘性</td>
    <td>Klemp and Wilhelmson (1978). Cm = 0.2</td>
    </tr>
    <tr>
    <td>時間フィルタ</td>
    <td>あり (係数は 0.1)</td>
    </tr>
    <tr>
    <td>数値粘性</td>
    <td>あり (5.0 x 10<sup>-4</sup>)</td>
    </tr>
    <tr>
    <td>計算領域</td>
    <td>水平 25.6 km x 鉛直 10 km</td>
    </tr>
    <tr>
    <td>格子点間隔</td>
    <td>水平 200 m x 鉛直 200 m</td>
    </tr>
    <tr>
    <td>積分時間</td>
    <td>6 時間</td>
    </tr>
    <tr>
    <td>時間ステップ</td>
    <td>長い: 3 秒, 短い: 0.3 秒</td>
    </tr>
    <tr>
    <td>境界条件</td>
    <td>水平: 周期, 鉛直: free-slip</td>
    </tr>
    <tr>
    <td>数値解法</td>
    <td>水平: 陽解法, 鉛直: 陰解法</td>
    </tr>
    </table>

Message

thumb.message << <<-Message

=begin
==計算結果

  * 出力ファイル
    * ((<初期値|URL:./nc/arare-odaka1998-dx200-20051008-init.nc>))
    * ((<時間発展|URL:./nc/arare-odaka1998-dx200-20051008.nc>))
    * ((<リスタート|URL:./nc/arare-odaka1998-dx200-20051008-restart.nc>))
    * ((<温位変化の寄与|URL:./nc/arare_diag-odaka1998-dx200-20051008.nc>))

=end
Message

thumb.create