惑星大気放射モデル 描画・解析 2 鉛直一次元フラックス分布, 加熱率分布の比較

結果の確認

二つの鉛直一次元フラックス分布・加熱率分布を比べる方法として, 二つのプロファイルの差を二乗平均誤差 (Root Mean Square Error; RMSE) として評価ことを考えるかもしれない. そのためにはこの ruby スクリプト を用いると良い.

スクリプト disp_comp_1Dresults_v1.rb ファイル内の以下の部分を書き換えて使う.

grav           = 9.8     # gravitational acceleration (m -2)
mmw            = 29.0e-3 # mean molecular weight (kg mol-1)
ncfn_t         = "lbl_k-dist_test_TRAPPIST-1/prog01.0_mkprofile_ascii/out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv.nc"
ncfn1_flux     = "lbl_k-dist_test_TRAPPIST-1/prog03.0_calc_rte/out/TRAPPIST-1_planet_Flux.nc"
ncfn1_tendency = "lbl_k-dist_test_TRAPPIST-1/prog03.0_calc_rte/out/TRAPPIST-1_planet_Tendency.nc"
ncfn2_flux     = "lbl_k-dist_test_TRAPPIST-1_dwn1e2/prog03.0_calc_rte/out/TRAPPIST-1_planet_Flux.nc"
ncfn2_tendency = "lbl_k-dist_test_TRAPPIST-1_dwn1e2/prog03.0_calc_rte/out/TRAPPIST-1_planet_Tendency.nc"

変数はそれぞれ以下の意味.

grav           : 重力加速度 (m s-2)
mmw            : 平均分子量 (kg mol-1)
ncfn_t         : 温度分布 (prog01.0_mkprofile_ascii/out 以下の NetCDF ファイルを指定すればよい)
ncfn1_flux     : 比較する一つ目の鉛直一次元フラックス (Flux) 分布のファイル
ncfn1_tendency : 比較する一つ目の鉛直一次元加熱率 (Tendency) 分布のファイル
ncfn2_flux     : 比較する二つ目の鉛直一次元フラックス (Flux) 分布のファイル
ncfn2_tendency : 比較する二つ目の鉛直一次元加熱率 (Tendency) 分布のファイル

上記 (の特にファイル名) を適宜変更し, 下のように使う.

$ 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  : lbl_k-dist_test_TRAPPIST-1/prog03.0_calc_rte/out/TRAPPIST-1_planet_Flux.nc
           : lbl_k-dist_test_TRAPPIST-1/prog03.0_calc_rte/out/TRAPPIST-1_planet_Tendency.nc
 2nd file  : lbl_k-dist_test_TRAPPIST-1_dwn1e2/prog03.0_calc_rte/out/TRAPPIST-1_planet_Flux.nc
           : lbl_k-dist_test_TRAPPIST-1_dwn1e2/prog03.0_calc_rte/out/TRAPPIST-1_planet_Tendency.nc
 variable : Planetary Rad.      Solar Rad.
  upflx   : 1.25820e+00    |    1.32184e-01 W m-2
  dnflx   : 2.41800e+00    |    1.99857e+00 W m-2
  flxcnv  : 5.63158e-04    |    3.27598e-04 W m-3
  tend    : 5.49075e-05    |    1.12680e-04 K s-1

それぞれの意味は以下.

Planetary Rad. : 惑星放射 (惑星表面および惑星大気が射出した放射) の分布とそれによる加熱率の RMSE.
Solar Rad.     : 太陽 (中心星/恒星) 放射 (太陽 (中心星/恒星) が射出した放射) の分布とそれによる加熱率の RMSE.
upflx          : 上向きフラックス (W m-2)
dnflx          : 下向きフラックス (W m-2)
flxcnv         : フラックスの収束 (W m-3)
tend           : 温度変化率 (K s-1)