#!/usr/bin/env ruby1.8
#
#= dcmodel thumbnail generate ruby script
#
#  Editor :: sugiyama
#  Version:: 2006/06/16 00:32:05
#
#== Overview
#
#This file is generate by following ruby script automatically.
#
#      /GFD_Dennou_Club/ftp/arch/dcmodel/bin/dcmodel-thum.rb
#      version "$Revision: 1.28 $ : $Date: 2006/02/11 19:48:19 $"
#
#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 "/GFD_Dennou_Club/ftp/arch/dcmodel/bin/dcmodel-thum.rb"

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

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

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

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

#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 
                              # が省略可能になる。

#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_height = 150
                              # 画像ファイルサイズ (高さ)

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

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

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

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

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

                              # フッターメッセージ。
                              # "=" ではなく、"<<" で代入することに注意！！
thumb.footer    = Array.new
#thumb.footer    << <<-Footer
# 
#Footer

                              # メッセージ。
                              # "=" ではなく、"<<" で代入することに注意！！
thumb.message   = Array.new
thumb.message   << <<-Message
=begin
[((<地球流体電脳倶楽部|URL:http://www.gfd-dennou.org>))]
[((<dcmodel|URL:http://www.gfd-dennou.org/library/dcmodel>)) |
((<dcmodel-tools|URL:http://www.gfd-dennou.org/library/dcmodel/doc/dcmodel-tools/SIGEN.htm>))]
[((<dcmodel-thum.rb|URL:http://www.gfd-dennou.org/library/dcmodel/doc/dcmodel-tools/dcmodel-thum-rdoc>))]

= deepconv/arare4 : 飽和蒸気圧の確認

== 計算の概要

arare4 にコーディングされた飽和蒸気圧式をプロットする. 

Titens の式はコーディングしていないが, 
地球の雲対流モデルで良く使われている式なので, 
一緒にプロットすることにした. 

== 利用した飽和蒸気圧の式

=== Antoine の式

  log( P[mmHg] ) = A - B / C + T[℃]
  
    H2O(l) の場合: A = 7.9186968, B = 1636.909, C = 224.92
    NH3(s) の場合: A = 9.96382, B = 1617.907, C = 272.55

    出典: 化学便覧 改訂第 4 版 

=== AMP 式

  ln( P [dyn/cm^2 = 0.1 Pa] ) = A / T[K]**2 + B + C ln( T[K] ) + D T[K] + E T[K]**2

    H2O(l) の場合: A = -2313.0338, B = -164.03307, C = 38.053682, D = -0.13844344, E = 7.4465367e-5
    H2O(s) の場合: A = -5631.1206, B = -8.363602, C = 8.2312, D = -0.03861449, E = 2.77494e-5
    NH3(s) の場合: A = -4122.0, B = 41.67871, C = -1.8163, D = 0.0, E = 0.0

    出典: Briggs and Sacket (1989, Icarus)

=== Tetens の式

  P[mbar] = 6.11 * 10 ** ( ( A * T[℃] ) / ( B + T[℃] ) )

    H2O(l) の場合: A = 7.5, B = 237.3
    H2O(s) の場合: A = 9.5, B = 265.3

    出典: 大気科学講座 2 「雲や降水を伴う大気」, 東京大学出版


== 備忘録

  * 一般に Antoine の式の適用範囲は 100 -- 150000 Pa 程度である. 

  * Antoine の式より求めた H2O の飽和蒸気圧は, 低温領域において, 
    非常に小さな値を持つ傾向にある. 

  * Antoine の式と Tetens の式では, 低温領域において, H2O(l) の飽和蒸気圧が H2O(s) の飽和蒸気圧より小さくなる. これは物理的におかしい. 
    * Antoine の式の適用範囲外に外挿しているので仕方ないかな.
    * AMP 式の方が性質のよい関数が良いと言えるかな.

  * モデル中で Antoine の式と AMP 式のどちらを使うか悩みどころ.
    * Antoine の式を用いると, dq/dz が大きくなりすぎ, 数値拡散が効き過ぎる恐れがある

== 結果一覧



=end

Message