惑星大気放射モデル 1

地球大気計算を行う方法を説明する.

[OLD VERSION] 大気プロファイルの準備 (prog01.?)

今回は既に準備されている

prog01.0_mkprofile/out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv.nc

を用いる. これは, 地球中緯度夏大気のプロファイル.

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

$ 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 ..

Line-by-line 吸収係数計算 (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,  ! wavenumber interval in m-1
  WaveNumS          =     5.0d2,  ! start wavenumber in m-1
  WaveNumE          = 49995.0d2,  ! end wavenumber in m-1
/
&gas_nml
  FlagNetCDF        = .true.,  ! format of gas parameter
  ! HITFileName(:) : HITRAN file name for each specie
  HITFileName(1)    = "01_hit12.nc",
  HITFileName(2)    = "02_hit12.nc",
  HITFileName(3)    = "03_hit12.nc",
  !
  ! FlagCalcLine(:) : flag for calculation of line absorption for each specie
  ! FlagCalcCont(:) : flag for calculation of continuum absorption for each specie
  FlagCalcLine(1:3) = 3*.true.,
  FlagCalcCont(1:3) = 3*.true.,
  !
  ! AirBroadHalfWidFac(:) : factor for air broadened half width for each specie
  AirBroadHalfWidFac(1:3) = 3*1.0d0,
  !
  ! CIAFileName(:) : CIA filename
  CIAFileName(1) = ""
  !
  ! XSCFileName(:) : absorption coefficient filename
  XSCFileName(1) = ""
  !
  ! LineCutOffWaveNum(:) : line cut off wavenumber in m-1 for each specie
  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 ..

Line-by-line 粒子光学パラメータ計算 (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 ..

Line-by-line 入射放射フラックス計算 (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 ..

Line-by-line 放射フラックス・加熱率計算 (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 ..