鉛直一次元放射計算 1

ひとまず何も変更せずに実行してみる.

プログラムでは, 長波放射モデル相互比較実験 (Ellingson et al., 1991) の Case 5 (Isothermal atmosphere 300 K, CO2 only 300 ppmv) の条件で実行される.

実験の実行

$ ./rad1d
 *** MESSAGE [namelist_util] ***  ----- Initialization Messages -----
 *** MESSAGE [namelist_util] ***    MaxNmlArySize = 256
 *** MESSAGE [namelist_util] ***  -- version = $Name:  $$Id: namelist_util.f90,v 1.4 2013/09/21 14:40:07 yot Exp $
 *** WARNING [constants] ***  NAMELIST group "constants_nml" is not found in "radvert1d.cntl" (iostat=-1).
...
(中略)
...
# Result
# k   Press  RadSUwFlux  RadSDwFlux  RadLUwFlux  RadLDwFlux
          0   100000.00000000000        226.67065909624384        1133.3532954812194        459.16384997273155        83.280785041062714
          1   81873.075307798179        235.42754659018502        1142.1101829786237        459.16384997273155        80.553433599087356
          2   67032.004603563924        242.45292781539044        1149.1355642068588        459.16384997273155        77.893018026728413
          3   54881.163609402633        248.09402498509931        1154.7766613791659        459.16384997273150        75.076977760998318
          4   44932.896411722148        252.62859994690777        1159.3112363431731        459.16384997273155        71.856094015041080
          5   36787.944117144223        256.27831895736216        1162.9609553554699        459.16384997273155        68.185073136662837
          6   30119.421191220201        259.21984258994041        1165.9024789895811        459.16384997273155        64.249431761026884
          7   24659.696394160641        261.59390780145480        1168.2765442023647        459.16384997273155        60.279715216393114
          8   20189.651799465533        263.51266985846729        1170.1953062604239        459.16384997273150        56.318654043420587
          9   16529.888822158646        265.06557934893732        1171.7482157517550        459.16384997273155        52.219202155515966
         10   13533.528323661261        266.32405843121569        1173.0066948347401        459.16384997273155        47.902665412321305
         11   11080.315836233382        267.34521063262844        1174.0278470367327        459.16384997273150        43.515153656759288
         12   9071.7953289412435        268.17476038225539        1174.8573967868344        459.16384997273155        39.312369270531413
         13   7427.3578214333820        268.84938112808652        1175.5320175330546        459.16384997273155        35.405909676054478
         14   6081.0062625217915        269.39853884726347        1176.0811752525501        459.16384997273150        31.611228092659800
         15   4978.7068367863903        269.84595180294326        1176.5285882084902        459.16384997273155        27.633342186295540
         16   4076.2203978366183        270.21074679110535        1176.8933831968654        459.16384997273144        23.372786632377256
         17   3337.3269960326052        270.50837565790135        1177.1910120638361        459.16384997273155        18.953933651820616
         18   2732.3722447292535        270.75134269942930        1177.4339791055068        459.16384997273150        14.440955306367613
         19   2237.0771856165575        270.94978298125903        1177.6324193874534        459.16384997273150        9.3948990224138438
         20   1831.5638888734163        271.11192309260338        1177.7945594988932        459.16384997273155        0.0000000000000000
# Comparison with Case 5 of Ellingson et al. (1991) (E91)
# Case 5 : Isothermal atmosphere 300 K, CO2 only 300 ppmv
# Results of E91 is shown in Table10 of E91.
# Unit : W m-2
# UpLW@Sfc =    459.16384997273155       (E91 Ave=458.85, StdDev=1.05)
# DnLW@Sfc =    83.280785041062714       (E91 Ave= 90.16, StdDev=7.66)
# UpLW@TOA =    459.16384997273155       (E91 Ave=458.90, StdDev=2.08)
  • 上記のうち, 以下の部分が, 計算された放射フラックスの鉛直分布
    • 左から
      • 鉛直面番号 [k]
      • 圧力 (Pa) [Press]
      • 上向き短波フラックス (W m-2) [RadSUwFlux]
      • 下向き短波フラックス (W m-2) [RadSDwFlux]
      • 上向き長波フラックス (W m-2) [RadLUwFlux]
      • 下向き長波フラックス (W m-2) [RadLDwFlux]
# Result
# k   Press  RadSUwFlux  RadSDwFlux  RadLUwFlux  RadLDwFlux
          0   100000.00000000000        226.67065909624384        1133.3532954812194        459.16384997273155        83.280785041062714
          1   81873.075307798179        235.42754659018502        1142.1101829786237        459.16384997273155        80.553433599087356
...
         19   2237.0771856165575        270.94978298125903        1177.6324193874534        459.16384997273150        9.3948990224138438
         20   1831.5638888734163        271.11192309260338        1177.7945594988932        459.16384997273155        0.0000000000000000
  • 以下はサンプル計算結果と他のモデルの結果との比較 (詳細は Ellingson et al. (1991) を参照すること)
# Comparison with Case 5 of Ellingson et al. (1991) (E91)
# Case 5 : Isothermal atmosphere 300 K, CO2 only 300 ppmv
# Results of E91 is shown in Table10 f E91.
# Unit : W m-2
# UpLW@Sfc =    459.16384997273155       (E91 Ave=458.85, StdDev=1.05)
# DnLW@Sfc =    83.280785041062714       (E91 Ave= 90.16, StdDev=7.66)
# UpLW@TOA =    459.16384997273155       (E91 Ave=458.90, StdDev=2.08)