#!/usr/bin/env ruby1.8
#
#= dcmodel thumbnail generate ruby script
#
#  Editor :: momoko
#  Version:: 2013/05/22 13:50:50
#
#== Overview
#
#This file is generate by following ruby script automatically.
#
#      /home/momoko/OneDimEquilibrium/2013-05-01_momoko/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 "/home/momoko/util/dcpam/OneDimEquilibrium/2013-05-22_momoko/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 = 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 = 3
                              # 横にならべるファイル数

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

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

thumb.title     = "1 次元放射対流平衡モデルを用いた射出限界の計算"
                              # html ヘッダのタイトル

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

                              # フッターメッセージ。
                              # "=" ではなく、"<<" で代入することに注意！！
thumb.footer    = Array.new
#thumb.footer    << <<-Footer
# 
#Footer
thumb.message   = Array.new
thumb.message   << <<-Message
=begin
[((<dcmodel|URL:http://www.gfd-dennou.org/library/dcmodel>)) |
[((<dcpam|URL:http://www.gfd-dennou.org/library/dcpam>)) |


= 1 次元放射対流平衡モデルを用いた射出限界の計算

== 概要

dcpam5 を用いて鉛直 1 次元放射対流平衡解を求める計算をおこなった.
地球放射スキームと灰色放射スキームを用いて, それぞれ惑星放射量
の最大値を求めてみた.

灰色放射を用いた場合では, 表面温度を変えても惑星放射量が
およそ 350 W/m2 を越えないことが確認された. 
これは Nakajima et al. (1992, JAS, 49, 2256) と同様の結果である.

地球放射スキームを用いた場合では, 表面温度が 400K を越えた場合に
正しく平衡解が求まらなかった (全層が対流圏になってしまった). 
また, L32 の場合では, 高温状態における放射計算もうまくできない.

== 実験設定

* プログラム: dcpam5-20130302 版に変更を加えたもの.
  * ((<tgz ファイル|URL:./src-code/dcpam5-20130302-1deq.tgz>))
* 軸対称モデル計算
  * 緯度格子点数 : 32
  * 鉛直層数: 32 または 1000
* 更に詳細は ((<発表資料|URL:./presen/pub/index.html>)) を参照してください.
* NAMELIST ファイル
  * ((<灰色放射, L32|URL:./namelists/GrayRadEqL32.conf>))
  * ((<灰色放射, L1000|URL:./namelists/GrayRadEqL1000.conf>))
  * ((<地球放射, L32|URL:./namelists/EarthConstStratL32.conf>))
  * ((<地球放射, L1000|URL:./namelists/EarthConstStratL1000.conf>))
  * ((<地球放射, CO2 混合比 10 倍, L32|URL:./namelists/EarthConstStratMRCO2times0.1L32.conf>))
  * ((<地球放射, CO2 混合比 0.1 倍, L32|URL:./namelists/EarthConstStratMRCO2times10L32.conf>))
  * ((<地球放射, 対流圏〜中間圏 CO2 放射スキームを用いる場合, L32|URL:./namelists/EarthConstStratHighAltL32.conf>))

  * ((<初期値生成用(L32)|URL:./namelists/init_data_1dim_L32.conf>))
  * ((<初期値生成用(L1000)|URL:./namelists/init_data_1dim_L1000.conf>))
  * ((<SST 生成用|URL:./namelists/sst_data_1dim.conf>))

== dcpam5-20130302 版からの変更点
* 平衡モデル用ルーチンを追加.
  * src/equilibrium ディレクトリを作成し, 以下のファイルを置いた.
      Makefile, equilib_modify_ps.f90, equilib_strat_const_temp.f90, equilib_tropopause.f90, equilib_moist_adiabat.f90, equilib_strat_optdep.f90
  * src/Makefile の LIBDIRS = の設定に equilibrium を追加.

* sst ファイルの作成の際に, SST の値を格子点ごとに指定するようにした.
  * src/prepare_data/surface_data.f90 を編集

* 灰色放射ルーチンを追加.
  * src/radiation/rad_gray.f90 を追加
  * src/radiation/Makefile の OBJS = の設定に rad_gray.o を追加

* namelist で読みこむ配列サイズの上限を大きくした.
  * setup/namelist_util.f90 において
        integer, parameter, public:: MaxNmlArySize = 2001
    と変更.

* 1 次元平衡計算ようメインルーチンを作成.
  * src/main/dcpam_main_1deq.f90  を作成
  * src/main/Makefile の OBJS = の設定に dcpam_main_1deq.o dcpam_main_1deq を追加

== 作図に用いたスクリプトたち
* ((<スクリプト置き場|URL:./scripts/>))


== 計算結果


=end

Message