木星大気計算

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

3. 粒子光学パラメータの計算

について説明する.

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

$ cd prog02.1_calc_optprop_particle

設定ファイルの準備

$ cp sample/None.calc_optprop.conf calc_optprop.conf
$ cp sample/None.calc_optprop_file.conf calc_optprop_file.conf

• calc_optprop.conf : 計算条件の設定 (! 以降はコメント)

  &wn_nml
    DelWaveNum        =     1.0d2,
    WaveNumS          =    10.5d2,
    WaveNumE          = 49999.5d2,
  /
  &ptcl_nml
    NMolNum     = 3,
    MolNum(1:3) = 1, 2, 3,
    ! Jupiter
  !  NMolNum     = 5,
  !  MolNum(1:5) = 1, 6, 11, 45, 91,
    ! Titan
  !  NMolNum     = 3,
  !  MolNum(1:3) = 6, 22, 45,
  !
    FlagReadQExtSSAAFFromFile = .true.,
    NPtcl      = 0,
    PtclName   = "",
    PtclDataFN = "",
    PtclNEffRadius = 0,
    PtclEffRadiusS = 1d-6,
    PtclEffRadiusE = 1d-6,
  /

  変数の説明と変数の変更.

  • wn_nml
    • prog02.0_calc_ac/calc_ac.conf の wn_nml と同じ値を指定する.

      • つまり

        DelWaveNum        =    10.0d2,
        WaveNumS          =    15.0d2,
        WaveNumE          = 49995.0d2,

  • ptcl_nml
    • NMolNum : 大気プロファイルのファイルに指定した気体種の数

      • 下のように指定.

        NMolNum     = 6,

    • MolNum(X) : 大気プロファイルのファイルに指定した気体の分子番号

      • 下のように指定.

        MolNum(1:6) = 1, 6, 11, 31, 45, 91,

    • FlagReadQExtSSAAFFromFile
      • 変更不要
    • NPtcl
      • 変更不要
    • PtclName
      • 変更不要
    • PtclDataFN
      • 変更不要
    • PtclNEffRadius
      • 変更不要
    • PtclEffRadiusS
      • 変更不要
    • PtclEffRadiusE
      • 変更不要

• calc_optprop_file.conf : 計算に用いるファイル等の設定 (! 以降はコメント)

  &file_nml
    OutputFileName    = "./out/Jupiter_particle_opt_prop.nc",
  /

  変数の説明と変数の変更.

  • file_nml
    • OutputFileName : 出力ファイルのパス
      • 変更不要

コンパイル

$ make

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

• calc_optprop
• merge

実行

$ ./calc_optprop

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

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

• out/Jupiter_particle_opt_prop.nc

注意:

• ただし, このプログラムでは, H2S のレイリー散乱係数を空気のレイリー散乱係数と等しいとしている.

親ディレクトリに戻る

$ cd ..