#!/usr/bin/env ruby
#
# 表題: 長波放射フラックスのお勉強
#
# 履歴: 2003/09/23 やまだ由
#       2003/10/01 やまだ由
#
# < 変数メモ(って途中までしかかいてないけど) >
#
# pres : 気圧 [半]
# tauq : 光学的厚さ, 水   [半] 
# taua : 光学的厚さ, 大気 [半]
# q    : 比湿 [整]
# pib  : \pi B [整]
# pib_g: \pi B (地面)
# rflux: 長波フラックス [半]
# trans: 透過係数
#
#
# ----------------------------------------------------------
require "numru/ggraph"
require "mradl-common"
require "mradl-physics"
require "mradl-dcl"
include NumRu
include NMath

# ----------------------------------------------------------
# main()

# 物理定数, 基本メゾットの指定
mradl_common
# オプション解析
mradl_optdcl

p $tncname

END{

  # 放射フラックス, 放射加熱計算
  mradl_physics

  # 可視化部
  mradl_dcl

}


# ----------------------------------------------------------
# NAMELIST <NMRADL> の指定
$nclng = 5
#rwgtl = [ 0.2  , 0.1  , 0.1  , 0.6  ]
#absal = [ 0.0  , 0.0  , 0.0  , 5.0E-5 ]
#absql = [ 8.0  , 1.0  , 0.1  , 0.0  ]

$rwgtl = [ 0.2  , 0.1  , 0.1  , 0.6, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
$absal = [ 0.0  , 0.0  , 0.0  , 5.0E-5 , 0.0 ]
$absql = [ 8.0  , 1.0  , 0.1  , 0.0  , 0.0]

#absal = [ 0.0  , 0.0  , 0.0  , 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]
#absql = [ 0.0  , 0.0  , 0.0  , 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]


# ----------------------------------------------------------
# 計算に使うデータ指定

# 地表気圧 adj700con1000, x=0, y=32
$ps = 1002.0*100

# <温度> とりあえず, 乾燥断熱減率 10 [K/km]. 
$Gamma   = 10.0/1000
$temp = $Temp_g * $sigma[1..-1]**($R*$Gamma/$Grav)

# <比湿> とりあえず乾燥大気
#q  = 0.0             

# <温度> adj700con1000, x=0, y=32
#tnc = NArray.sfloat(16).fill!(Temp_g) 
$tnc = GPhys::NetCDF_IO.open("#{$tncname}", 't').cut('lat'=>0,'time'=>0).average(0).val

# <比湿> adj700con1000, x=0, y=32
#qnc = NArray.sfloat(16).fill!(1.0) 
$qnc = GPhys::NetCDF_IO.open("#{$qncname}", 'q').cut('lat'=>0,'time'=>0).average(0).val

# <長波放射フラックス> adj700con1000, x=0, y=32
# $lradnc = GPhys::NetCDF_IO.open('lrad.nc', 'lrad').cut('lat'=>0,'time'=>0).average(0).val

# <長波放射加熱> adj700con1000, x=0, y=32
$qradlnc = GPhys::NetCDF_IO.open("#{$qradlncname}", 'qradl').cut('lat'=>0,'time'=>0).average(0).val