#!/usr/bin/env ruby1.8
#
#= dcmodel thumbnail generate ruby script
#
#  Editor :: deepconv
#  Version:: 2008/02/16 21:07:13
#
#== Overview
#
#This file is generate by following ruby script automatically.
#
#      /GFD_Dennou_Club/ftp/arch/deepconv/arare/sample/2008-02-16_yamasita/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 "/GFD_Dennou_Club/ftp/arch/deepconv/arare/sample/2008-02-16_yamasita/thum-src/dcmodel-thum.rb"

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

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

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

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

#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 = true
                              # 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 = 13
                              # 横にならべるファイル数

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

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

thumb.title     = "deepconv/arare4 : 北守修論再計算(精度向上編)"
                              # 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 での計算結果は arare3 の計算結果とあまり一致していなかった. 
その後, 以下のような計算設定の相違があることが判明した. 
なお, arare3 と arare4 の計算結果を近づけるに当たり, 一方のプログラムのみを書き換えるのではなく, 計算精度が悪いと思われる方のプログラムを改良して両者の歩み寄りを図った. 
* 同じ基本場温度分布を与えているつもりだったが, 持ち上げ凝結高度・乾燥温度減率などに微妙なずれが存在していた. 
  * arare3 では CO2 分子量のデフォルト値を 44.0 としていたが, arare4 では CO2 分子量のデフォルト値としてより真面目な値 44.0095 を使用していた. 
* arare3 では渦粘性係数 Km に対して Assilin フィルタをかけていないが, arare4 ではかけて計算していた. 
* arare3 において Euler 法を用いて第二の初期値を得る init.f90 の温位の計算の所で熱散逸項が抜けていた. 
* arare4 では第二の初期値も main プログラム内で計算しているが, この過程で本来不要である Asselin フィルターがかかってしまっていた. 
  * Asselin フィルターは Leap Frog 法の補正の為にかけるべきものなので, Euler 法に対して用いるのは適切でないと考えられる. 
* arare3 では温位に関する数値粘性・乱流拡散を 4 次精度の中心差分で評価しているが, arare4 では 2 次精度の中心差分で評価していた. 
* arare3 の init.f90 において, 短い時間ステップのループ内での鉛直速度の計算の入力値として ns 値ではなく, nl 値を用いてしまっていた. 


== 計算設定

* 北守修論の設定と同じ. 詳細は[((<こちら|URL:http://www.gfd-dennou.org/arch/deepconv/arare/sample/2006-02-09_kitamo/>))]を参照. 

* 計算領域
  * 水平 20 km
  * 鉛直 20 km

* 基本場の設定
  * 地表気温 165 K, 地表圧力 7 hPa
  * 地表から高度 4 km 付近まで等温位, 高度 4 km 付近から高度 15 km 付近まで湿潤断熱, 高度 15 km 付近から等温(135 K). 

* 空間/時間格子間隔, 境界条件
  * 空間格子間隔
    * 水平 200 m, 鉛直 200 m
    * 格子点数は 100 x 100
  * 時間格子間隔
    * 長い時間ステップ: 2 sec
    * 短い時間ステップ: 0.25 sec
  * 境界条件
      * 水平境界: 周期境界条件
      * 上下境界: 応力なし, w=0 , 熱フラックスなし

* 初期条件
  * x = 10 km を中心として,  ガウス型の正の温位擾乱を与える.
  * 振幅は 2.0 K

* 積分時間
  * 3600 sec



== 結果一覧

計算結果は三段構成で並べている. 
上段が arare3 での計算結果, 
中段が凝結過程移植後の arare4 での計算結果(転置なし), 
中段が凝結過程移植後の arare4 での計算結果(転置あり)である. 



=end

Message









