地球大気お試し放射計算

地球大気放射場のラインバイライン計算を行う方法を説明する.

ラインバイライン放射フラックス・加熱率計算 (prog03.0)

prog02.X で計算された結果を用いて放射フラックスと加熱率を計算する.

$ cd prog03.0_calc_rte

設定ファイルの準備

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

• calc_rte.conf : 計算条件の設定 [変更なし]
• calc_rte_file.conf : 計算に用いるファイル等の設定 [変更なし]
• planetconst.conf : 惑星パラメータの設定 [変更なし]

コンパイル

$ make

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

• calc_rte

実行

$ ./calc_rte

Finish と表示されれば無事実行終了.

計算結果のファイルは以下.

• out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_Flux.nc
• out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_Spectral_Flux.nc
• out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_Tendency.nc

それぞれのファイルに保存されている変数は下のような意味である.

out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_Flux.nc
  Press       : (層境界)圧力 (Pa)
  RadUwFlux   : 上向き(全)放射フラックス (W m-2)
  RadDwFlux   : 下向き(全)放射フラックス (W m-2)
  RadUwFluxLW : 上向き惑星放射フラックス (W m-2)
  RadDwFluxLW : 下向き惑星放射フラックス (W m-2)
  RadUwFluxSW : 上向き太陽放射フラックス (W m-2)
  RadDwFluxSW : 下向き太陽放射フラックス (W m-2)
out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_Spectral_Flux.nc
  Press          : (層境界)圧力 (Pa)
  PressLC        : (層中心)圧力 (Pa)
  WaveNum        : 波数 (m-1)
  RadUwFlux      : 上向き(全)放射フラックス密度 (W m-2 (m-1)-1)
  RadDwFlux      : 下向き(全)放射フラックス密度 (W m-2 (m-1)-1)
  RadUwFluxLW    : 上向き惑星放射フラックス密度 (W m-2 (m-1)-1)
  RadDwFluxLW    : 下向き惑星放射フラックス密度 (W m-2 (m-1)-1)
  RadUwFluxSW    : 上向き太陽放射フラックス密度 (W m-2 (m-1)-1)
  RadDwFluxSW    : 下向き太陽放射フラックス密度 (W m-2 (m-1)-1)
  OptDep         : 光学的厚さ (1)
  DelOptDep      : 各層の光学的厚さ
  DRadFluxDPress : (全)放射フラックスの圧力微分 (W m-2 Pa-1 (m-1)-1)
out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_Tendency.nc
  Press       : (層境界)圧力 (Pa)
  DTempDt     : (全)放射による温度変化率 (K s-1)
  DTempDtLW   : 惑星放射による温度変化率 (K s-1)
  DTempDtSW   : 太陽放射による温度変化率 (K s-1)

描画に関わる細かい事

結果を描画する際に単位を変更したり, 波数軸を波長軸にの変更するにはこちらを参照すると良い.

親ディレクトリに戻る

$ cd ..