#!/usr/bin/env ruby1.8
# -*- coding: euc-jp -*-
#
#= dcmodel thumbnail generate ruby script
#
#  Editor :: yukai
#  Version:: 2011/08/24 16:34:43
#
#== Overview
#
#This file is generate by following ruby script automatically.
#
#      /GFD_Dennou_Club/ftp/arch/dcpam/sample/2011-08-24_yukai/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/dcpam/sample/2011-08-24_yukai/thum-src/dcmodel-thum.rb"

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

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

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

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

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

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

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

thumb.title     = "サブグリッドスケールの過程に関するテスト計算"
                              # html ヘッダのタイトル

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

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

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

=begin サブグリッドスケールの過程に関するテスト計算

= サブグリッドスケールの過程に関するテスト計算

* ((<概要>))
* ((<実験環境>))
* ((<実験設定>))
* ((<考察メモ>))
* ((<結果>))

=== 更新履歴

* 2011.08.24 / ページ作成

== 概要
 
鉛直乱流混合の評価に用いるリチャードソン数の計算と, 積雲対流調節のスキーム, 大規模凝結のスキームを変更し, ((<雲の寿命と移流の有無による放射の変化|URL:../2011-06-21_yukai_1>)) と同様の設定で実験を行った.

今回の変更前/後の計算方法は以下の通りである.

* 変更前
  * 鉛直乱流混合
    * 温位を用いてリチャードソン数を計算
  * 積雲対流調節
    * Manabe (1965) の湿潤対流調節
  * 大規模凝結
    * Manabe (1965)

* 変更後
  * 鉛直乱流混合
    * 仮温位を用いてリチャードソン数を計算
      * 湿潤大気における空気と水蒸気の分子量の違いを考慮
  * 積雲対流調節
    * Relaxed Arakawa-Schubert scheme (Moorthi and Suarez, 1992)
  * 大規模凝結
    * Le Treut and Li (1991)
  

== 実験環境

下記の 2 台で計算した. コンパイラ, 依存するライブラリは同じ (同じバージョン) である.

=== 計算機名 : joho09-itpass

  * CPU : Intel(R) Core(TM) i7-2600 CPU @ 3.40GHz
  * OS  : GNU Debian Linux 5.0 (lenny)
  * uname -a : Linux joho09-itpass 2.6.26-2-amd64 #1 SMP Mon Jun 13 16:29:33 UTC 2011 x86_64 GNU/Linux


=== 計算機名 : joho12-itpass

  * CPU : Intel(R) Core(TM) i7 CPU 860  @ 2.80GHz
  * OS  : GNU Debian Linux 6.0 (sueeze)
  * uname -a : Linux joho12-itpass 2.6.32-5-amd64 #1 SMP Mon Mar 7 21:35:22 UTC 2011 x86_6
4 GNU/Linux


=== コンパイラ

* Fujitsu Fortran 6.0

=== 依存するライブラリ (バージョン)

  * netCDF (3.6.2)
  * gtool5 (20101228-1)
  * ISPACK (0.93)
  * spmodel (0.6.1)

== 実験設定

* ソースファイル
  * dcpam5 ver.20110407 を改変したもの [((<tar.gz|URL:./src/dcpam5-20110714-mod.tar.gz>))]
  * ((<ソースツリー|URL:./src/dcpam5-20110714-mod>))

* NAMELIST ファイル
  * 移流なし
    * 雲の寿命: ((<900 秒|URL:./nml/dcpam_E_T21L22_CLT0900-AdvF_RAS_LL91_Ri.nml.txt>)),\
      ((<1200 秒|URL:./nml/dcpam_E_T21L22_CLT1200-AdvF_RAS_LL91_Ri.nml.txt>)),\
      ((<1500 秒|URL:./nml/dcpam_E_T21L22_CLT1500-AdvF_RAS_LL91_Ri.nml.txt>)), ((<1800 秒|URL:./nml/dcpam_E_T21L22_CLT1800-AdvF_RAS_LL91_Ri.nml.txt>))

* 空間解像度 : T21L22
  * 経度格子点数 : 64 (格子間隔 5.6 度)
  * 緯度格子点数 : 32
  * 鉛直層数 : 22

* タイムステップ : 24 分
* 積分時間 : 20 年

* 初期値
  * 温度場 : 280 K + 小擾乱

* 地表面気圧 : 10^5 Pa
* 東西風速 : 0 m/s
#  * 比湿: ??

#* 境界条件

* 惑星半径 : 6371 km
* 重力加速度 : 9.8 m/s

== 考察メモ

* 放射の全球平均値
  * スキーム変更前よりも絶対値が小さい
    * 10 W/m^2 ほど小さい
  * Threnberth et al. (2007) と合わない
  * 収支は雲の寿命が 900 から 1800 秒のどの場合でも釣り合っている
    * OLR と OSR の和がどの場合も O(1)

* OLR の東西平均分布
  * スキーム変更前と比べて, 低緯度の分布が NCEP 再解析データの分布に近付いた
    * 凹むようになった
  * 北緯 60 度から南緯 60 度における値が NCEP よりも小さくなっ.



== 結果

=== OLR と OSR の全球平均値 (W/m^2)

  雲の寿命(秒)        900.0d0           1200.0d0          1500.0d0          1800.0d0
  --------------------------------------------------------------------------------------------------------
  移流なし   OLR   230.880880004771  227.943034895198  225.779645542562  223.814797857556
             OSR  -236.561760401199 -228.743918216485 -222.717111427596 -217.317611223577


=== 図

* fig.1 から fig.4
  * 積分時間の 11 年目から 20 年目の値を全球平均した OLR, OSR
  * 雲の寿命が 900, 1200, 1500, 1800 秒で移流なしの場合

* fig.5 から fig.8
  * 10 年平均・東西平均した OLR
    * 実線: 雲の寿命が 900, 1200, 1500, 1800 秒で移流なしの場合の計算結果
      * 積分時間の 11 年目から 20 年目の値を平均
    * 破線: NCEP の再解析データ
      * 1988 年から 2007 年の値の平均
      * 再解析データの NetCDF ( ((<上向き長波放射|URL:./data/ulwrftop.mon.mean.nc>)), ((<上向き短波放射|URL:./data/dswrftop.mon.mean.nc>)), ((<入射短波放射|URL:./data/dswrftop.mon.mean.nc>)) )
        * ((<取得元: 京都大学生存圏研究所 の NCEP 再解析データミラーサイト|URL:http://database.rish.kyoto-u.ac.jp/arch/ncep/>))

=end


Message

