惑星大気放射モデル 1

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

概要

以下の手順で行う.

1. 大気プロファイルの準備
2. 気体吸収係数の計算
3. 粒子光学パラメータの計算
4. 入射放射フラックスの計算
5. 放射フラックスの計算

大気プロファイルの準備 (prog01.X)

$ cd prog01.0_mkprofile_ascii

設定ファイルの準備

$ cp sample/Earth.mkprofile.conf mkprofile.conf
  • mkprofile.conf : 計算条件の設定

コンパイル

$ make

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

  • mkprofile

実行

$ ./mkprofile

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

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

  • out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv.nc

親ディレクトリに戻る

$ cd ..

ラインバイライン吸収係数計算 (prog02.0)

$ cd prog02.0_calc_ac

設定ファイルの準備

$ cp sample/Earth.calc_ac.conf calc_ac.conf
$ cp sample/Earth.calc_ac_file.conf calc_ac_file.conf

  • calc_ac.conf : 計算条件の設定

    &wn_nml
  DelWaveNum        =    10.0d2,
  WaveNumS          =     5.0d2,
  WaveNumE          = 49995.0d2,
/
&gas_nml
  FlagNetCDF        = .true.,
  HITFileName(1) = "../../../radmodels-data/HITRAN2012/NC_By-Molecule/01_hit12.nc",
  HITFileName(2) = "../../../radmodels-data/HITRAN2012/NC_By-Molecule/02_hit12.nc",
  HITFileName(3) = "../../../radmodels-data/HITRAN2012/NC_By-Molecule/03_hit12.nc",
!
  FlagCalcLine(1:3) = 3*.true.,
  FlagCalcCont(1:3) = 3*.true.,
!
  AirBroadHalfWidFac(1:3) = 3*1.0d0,
!
  CIAFileName(1) = ""
!
  XSCFileName(1) = ""
!
  LineCutOffWaveNum(1:3) = 3*25.0d2,
/
  • calc_ac_file.conf : 計算に用いるファイル等の設定

コンパイル

$ make

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

  • calc_ac
  • merge

実行

$ ./calc_ac

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

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

  • out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_ac.nc

親ディレクトリに戻る

$ cd ..

ラインバイライン粒子光学パラメータ計算 (prog02.1)

$ cd prog02.1_calc_optprop_particle

設定ファイルの準備

$ cp sample/Earth.calc_optprop.conf calc_optprop.conf
$ cp sample/Earth.calc_optprop_file.conf calc_optprop_file.conf
  • calc_optprop.conf : 計算条件の設定
  • calc_optprop_file.conf : 計算に用いるファイル等の設定

コンパイル

$ make

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

  • calc_optprop
  • merge

実行

$ ./calc_optprop

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

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

  • out/Earth_particle_opt_prop.nc

親ディレクトリに戻る

$ cd ..

ラインバイライン入射放射フラックス計算 (prog02.5)

$ cd prog02.5_calc_stellarspectrum

設定ファイルの準備

$ cp sample/Earth.calc_str_spe.conf calc_str_spe.conf
  • calc_str_spe.conf : 計算の条件および計算に用いるファイル等の設定

コンパイル

$ make

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

  • calc_str_spe

実行

$ ./calc_str_spe

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

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

  • out/solar_flux_Gueymard2004_1366.1Wm-2.nc

親ディレクトリに戻る

$ cd ..

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

$ 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

親ディレクトリに戻る

$ cd ..