#!/usr/bin/env ruby1.8
# -*- coding: euc-jp -*-
#
#= 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/2010-07-27_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 = 37
                              # 横にならべるファイル数

#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 : 主成分凝結を考慮した湿潤対流実験(中層一様冷却・下層一様加熱, 臨界飽和比 1.0, 地表面温度 165 K の場合) "
                              # 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 : 主成分凝結を考慮した湿潤対流実験(中層一様冷却・下層一様加熱, 臨界飽和比 1.0, 地表面温度 165 K の場合)

雲粒落下・荷重効果なしで, 計算領域の鉛直距離を 50 km にして計算を行なったところ, 140 秒程度の時間積分を行なうと計算が破綻してしまった. 
これまでに行った試みと結果を整理する. 

== 計算設定

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

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

* 空間/時間格子間隔, 境界条件
  * 空間格子間隔
    * 水平 200 m, 鉛直 200 m
    * 格子点数は 250 x 250
  * 時間格子間隔
    * 長い時間ステップ 1.0 sec, 短い時間ステップ 0.125 sec
  * 境界条件
      * 水平境界: 周期境界条件
      * 上下境界: 応力なし, w=0 , 熱フラックスなし
  * 放射加熱・冷却
    * 高度 1 -- 15 km で -5.0 K/day の冷却を与える. 
    * 高度 0 -- 1 km で放射冷却と釣り合うだけの加熱を与える
  * スポンジ層
    * 厚さ : 15 km
    * e-folding time : 300 sec

* 初期条件
  * 最下層で振幅 1.0 K のランダムな温位擾乱を与える. 

* 積分時間
  * 3600 sec

== まとめ

* 症状
  * 鉛直距離を 50 km にすると, 上部境界付近で two-grid wave が発生, 
    計算が破綻してしまう
* 試してみたこと
  * スポンジ層の e-folding time を 300 sec -> 3000, 100, 30, 1 sec と変更
    * いずれも 140 sec 弱で破綻
  * 長いタイムステップを 1.0 sec -> 0.5 sec, 短い時間ステップを 0.125 sec 
    -> 0.1 sec に変更
    * 130 sec 程度で破綻
  * 数値粘性係数を 5 倍(alpha を 1.0e-4 -> 5.0e-4 )に
    * 140 sec 弱で破綻
  * スポンジ層の厚さを 15 km -> 18 km に変更
    * 140 sec 弱で破綻
  * 鉛直距離を 40 km, スポンジ層の厚さを 12 km に変更
    * 10 日間時間積分をすることが出来た
    * 上部境界でギザギザも発生していない模様

== 結果一覧(絵は後程更新予定)



=end

Message




