惑星大気放射モデル 2 波数解像度依存性実験

概要

• 高い波数解像度で計算する.
  • 各プログラムの設定ファイルで設定されている波数解像度を変更する.
  • 波数解像度が高い時にどの程度の計算時間がかかるかを体験する

設定の変更

前回からの変更箇所は以下. 設定ファイルを書き換えて前回と同様に実行する.

prog02.0_calc_ac

• 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
  /

  を

  &wn_nml
    DelWaveNum        =     1.0d2,  ! wavenumber interval in m-1
    WaveNumS          =     0.5d2,  ! start wavenumber in m-1
    WaveNumE          = 49999.5d2,  ! end wavenumber in m-1
  /

  とすることで波数解像度が変更される.

• calc_ac_file.conf の変更

  &file_nml
    OutputFileName = "./out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_ac.nc",
  /

  を

  &file_nml
    OutputFileName = "./out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_01.0_ac.nc",
  /

  とすることで計算結果のファイルを変更する.

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

  • out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_01.0_ac.nc

prog02.1_calc_optprop_particle

• calc_optprop.conf の変更

  &wn_nml
    DelWaveNum        =    10.0d2,
    WaveNumS          =     5.0d2,
    WaveNumE          = 49995.0d2,
  /

  を

  &wn_nml
    DelWaveNum        =     1.0d2,
    WaveNumS          =     0.5d2,
    WaveNumE          = 49999.5d2,
  /

  とすることで波数解像度を変更する.

• calc_optprop_file.conf の変更

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

  を

  &file_nml
    OutputFileName    = "./out/Earth_particle_opt_prop_01.0.nc",
  /

  とすることで計算結果のファイルを変更する.

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

  • out/Earth_particle_opt_prop_01.0.nc

prog02.5_calc_stellarspectrum

• calc_str_spe.conf の変更

  &wn_nml
    DelWaveNum        =    10.0d2,
    WaveNumS          =     5.0d2,
    WaveNumE          = 49995.0d2,
  /

  を

  &wn_nml
    DelWaveNum        =     1.0d2,
    WaveNumS          =     0.5d2,
    WaveNumE          = 49999.5d2,
  /

  とすることで波数解像度を変更する.

  &spectrum_nml
    OutFileName = "./out/solar_flux_Gueymard2004_1366.1Wm-2.nc",
  /

  を

  &spectrum_nml
    OutFileName = "./out/solar_flux_Gueymard2004_1366.1Wm-2_01.0.nc",
  /

  とすることで計算結果のファイルを変更する.

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

  • out/solar_flux_Gueymard2004_1366.1Wm-2_01.0.nc

prog03.0_calc_rte

• calc_rte.conf の変更

  &rte_nml
    StrSpeNcFN  = "../prog02.5_calc_stellarspectrum/out/solar_flux_Gueymard2004_1366.1Wm-2.nc",
  /

  を

  &rte_nml
    StrSpeNcFN  = "../prog02.5_calc_stellarspectrum/out/solar_flux_Gueymard2004_1366.1Wm-2.nc",
  /

  とすることで入射放射スペクトルのファイルを変更する.

• calc_rte_file.conf の変更

  &rte_file_nml
    PtclOptPropNcFN = "../prog02.1_calc_optprop_particle/out/Earth_particle_opt_prop.nc"
    InOptPropNcFN   = '../prog02.0_calc_ac/out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_ac.nc',
    OutFluxNcFn     = "out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_Flux.nc",
    OutTendencyNcFn = "out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_Tendency.nc"
    OutSpeFluxNcFn  = "out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_Spectral_Flux.nc
  /

  を

  &rte_file_nml
    PtclOptPropNcFN = "../prog02.1_calc_optprop_particle/out/Earth_particle_opt_prop_01.0.nc"
    InOptPropNcFN   = '../prog02.0_calc_ac/out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_01.0_ac.nc',
    OutFluxNcFn     = "out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_01.0_Flux.nc",
    OutTendencyNcFn = "out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_01.0_Tendency.nc"
    OutSpeFluxNcFn  = "out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_01.0_Spectral_Flux.nc
        /

  とすることで, 粒子の光学特性のファイルと計算結果のファイルを変更する.

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

  • out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_01.0_Flux.nc
  • out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_01.0_Spectral_Flux.nc
  • out/Earth_ICRCCM_LW_Case27_MLS_CO2-300ppmv_01.0_Tendency.nc

前回の結果との違いを確認する.