k 分布作成 : k 分布計算の確認

k 分布法を用いて放射場を計算し, 作成した k 分布テーブルを確認する方法を説明する.

k 分布計算の実施 (prog08.0v2)

$ cd prog08.0v2_calc_rte_check
$ rm -f out/*

設定ファイルの準備

地球用サンプル設定ファイル

$ cp sample/Venus.calc_rte.conf calc_rte.conf
$ cp sample/Venus.calc_rte_file.conf calc_rte_file.conf
$ cp sample/Venus.planetconst.conf planetconst.conf

calc_rte.conf : 放射伝達計算の設定 [変更なし]
calc_rte_file.conf : 入出力ファイルの設定 [変更なし]
planetconst.conf : 惑星パラメータの設定 [変更なし]

コンパイル

$ make

コンパイルが成功すると, 下の実行ファイルができる.

calc_rte

実行

$ ./calc_rte

標準出力に Finish と表示されれば無事実行終了. (場合によってはスタックが不足してエラーが出るので, "ulimit -s unlimited" とする.) 計算結果は

out/Venus_ProfB_Flux.nc
out/Venus_ProfB_Tendency.nc

結果の確認

k 分布計算とラインバイライン計算の結果の差を Root Mean Square Error (RMSE) として評価する.

ruby/GPhys スクリプトで確認する. この ruby スクリプトを用いる.

disp_comp_1Dresults_v1.rb の以下を編集して入力ファイルを指定.

grav =  8.9    # gravitational acceleration (m -2)
mmw  = 44.0e-3 # mean molecular weight (kg mol-1)
ncfn_t         = '../prog01.0_mkprofile_ascii/out/Venus_ProfB.nc'
#
ncfn1_flux     = 'out/Venus_ProfB_Flux.nc'
ncfn1_tendency = 'out/Venus_ProfB_Tendency.nc'
ncfn2_flux     = '../prog03.0_calc_rte/out/Venus_ProfB_Flux.nc'
ncfn2_tendency = '../prog03.0_calc_rte/out/Venus_ProfB_Tendency.nc'

スクリプトを実行.

$ ruby disp_comp_1Dresults_v1.rb 1
 *** MESSAGE (SWDOPN) ***  GRPH1 : STARTED / IWS =  1.
 *** MESSAGE (UXPTTL) ***  CHARACTER LENGTH IS ZERO.
 ...
 *** MESSAGE (SWPCLS) ***  GRPH1 : PAGE =   1 COMPLETED.
 *** MESSAGE (SWDCLS) ***  GRPH1 : TERMINATED.
RMSE summary
 1st file  : out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_Flux.nc
           : out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_Tendency.nc
 2nd file  : ../prog03.0_calc_rte/out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_Flux.nc
           : ../prog03.0_calc_rte/out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_Tendency.nc
 variable : Planetary Rad.      Solar Rad.
  upflx   : 3.50564e-01    |    5.41147e-01 W m-2
  dnflx   : 2.88914e-01    |    5.02385e-01 W m-2
  flxcnv  : 1.16026e-04    |    6.97129e-05 W m-3
  tend    : 1.32899e-05    |    3.74033e-05 K s-1

雲なしの場合.

RMSE summary
 1st file  : out/Venus_ProfB_nocld_Flux.nc
           : out/Venus_ProfB_nocld_Tendency.nc
 2nd file  : ../prog03.0_calc_rte/out/Venus_ProfB_nocld_Flux.nc
           : ../prog03.0_calc_rte/out/Venus_ProfB_nocld_Tendency.nc
 variable : Planetary Rad.      Solar Rad.
  upflx   : 2.53996e-01    |    3.67580e-01 W m-2
  dnflx   : 3.02244e-01    |    8.35107e-01 W m-2
  flxcnv  : 1.36286e-04    |    1.94423e-04 W m-3
  tend    : 1.32591e-05    |    3.73006e-05 K s-1

親ディレクトリに戻る

$ cd ..