#!/usr/bin/env ruby1.8
#
#= dcmodel thumbnail generate ruby script
#
#  Editor :: sugiyama
#  Version:: 2007/05/17 21:27:32
#
#== Overview
#
#This file is generate by following ruby script automatically.
#
#      /usr/local/dcmodel-tools/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 "/usr/local/dcmodel-tools/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 
                              # が省略可能になる。

thumb.blankfig = true
                              # 画像ファイル名を表示しない

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

#thumb.thumbnaildir    = "../figdir-thumb"
                              # サムネイル画像を置くディレクトリの名前

#thumb.thumbnailtail    = "_thumb.png"
                              # サムネイル画像の拡張子名

#thumb.thumbnailpage = "1"
                              # サムネイル画像化する元画像のページ番号

#thumb.convert_cmd = "convert -depth 8 -geometry"
                              # convert コマンドおよびオプション。この
                              # 後ろに画像サイズが 200x150 のような形
                              # 式で指定される。

#thumb.convert_overwrite = false
                              # convert コマンドが実行される際、既に出
                              # 力先のファイルがある場合に上書きするた
                              # めのフラグ。これが false の場合、出力
                              # 先のファイルが存在し、且つ元のファイル
                              # よりも新しい場合は変換を行わない。

#thumb.css       = "/usr/local/dcmodeltools/dcmodel.css"
thumb.css       = "../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  = 400
                              # 画像ファイルサイズ (幅)

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

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

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

#thumb.html_author  = "sugiyama"
                              # 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>))]

= 統計量の出力

== はじめに

宿題となっていた, skewness (歪み度), krutosis (尖り度) をプロットする. 

=== 歪み度

分布の対称･非対称性をあらわすパラメタとして用いられる. 対称な分布については値がゼロとなる. 

γ1=Σ(Xi-μ)3/(Nσ3)
(
σ : 標準偏差 , 
μ : 平均
)
=== 尖り度

正規分布からの形状のずれを正負で表わすようにしている. 
平均付近の形状が尖っているほどクルトシスは大きくなる. 

γ2=Σ(Xi-μ)4/(Nσ4)

=== 参考

 * http://www.wwq.jp/stacalcul2/skewness.htm
 * http://www.wwq.jp/stacalcul2/kurtosis.htm
 * http://zisin.geophys.tohoku.ac.jp/~naka/timeseries.pdf

== 図の解釈

* 歪み度
  * 活発期(t = 840, 890, 920 day) において, H2O 凝結高度より下で値が負. 
    => 下降流が卓越していることを意味する
  * 静穏期において, 雲対流層での値は基本的に負. 
    => 下降流が卓越. 上で冷やされて下降する系となっている. 
  * 活発期の雲対流領域で瞬間的に値が正
    => 活発期は H2O 凝結高度から対流が発生

* 尖り度
  * 自乗平均で値が大きい領域(= 対流運動が見られる領域)に対応して, 
    尖り度も大きくなっている. 
    => その領域での速度は平均値からあまりずれていない. 

* 結果として....
  * 歪み度より, 下降流によって駆動される系である, といえる. 

== 結果一覧

=end

Message