木星大気計算

木星大気放射場のラインバイライン計算を行う際の

4. 入射放射フラックスの計算

について説明する.

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

$ cd prog02.5_calc_stellarspectrum

設定ファイルの準備

$ cp sample/Earth.calc_str_spe.conf calc_str_spe.conf

calc_str_spe.conf : 計算の条件および計算に用いるファイル等の設定 (! 以降はコメント)
```
&wn_nml
  DelWaveNum        =    10.0d2,
  WaveNumS          =    15.0d2,
  WaveNumE          = 49995.0d2,
/
&spectrum_nml
  StrTemp    = -1.0d0,
  InFileName = "data/solar_irradiance_Gueymard2004.nc",
  StrFluxTOA = 1366.1d0
  OutFileName = "./out/solar_flux_Gueymard2004_1366.1Wm-2.nc",
  Comment     = 'Solar spectral irradiance of Gueymard (2004), 1366.1 Wm-2'
!
  WaveNumName = "wn",
  SpeName     = "irradiance"
/
```
変数の説明と変数の変更.
- wn_nml
  - prog02.0_calc_ac/calc_ac.conf の wn_nml と同じ値を指定する.
    - つまり
```
DelWaveNum        =    10.0d2,
WaveNumS          =    15.0d2,
WaveNumE          = 49995.0d2,
```
- spectrum_nml
  - StrTemp : 中心星 (太陽) 表面温度.
    - 正の場合はプランク関数が計算される.
    - 負の場合は無効.
  - InFileName : 中心星 (太陽) スペクトルの入力ファイルのパス
    - 変更不要
  - StrFluxTOA : 中心星 (太陽) 放射の大気上端におけるフラックス
    - 下のように指定
```
StrFluxTOA = 50.5d0
```
    - この値は, 1366.1 ÷ 5.2**2
  - OutFileName : 出力ファイルのパス.
    - 例えば下のように指定.
```
OutFileName = "./out/solar_flux_50.5Wm-2.nc",
```
  - Comment : コメント
    - 例えば下のように指定.
```
Comment = 'Solar spectral irradiance of Gueymard (2004), 50.5 Wm-2'
```
  - WaveNumName
    - 変更不要
  - SpeName
    - 変更不要

コンパイル

$ make

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

calc_str_spe

実行

$ ./calc_str_spe

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

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

out/solar_flux_50.5Wm-2.nc

親ディレクトリに戻る

$ cd ..