rrtm_kgb1_01.F90

Go to the documentation of this file.

!OCL SCALAR
SUBROUTINE rrtm_kgb1_01

!     Originally by Eli J. Mlawer, Atmospheric & Environmental Research.
!     BAND 1:  10-250 cm-1 (low - H2O; high - H2O)
!     Reformatted for F90 by JJMorcrette, ECMWF

!     ------------------------------------------------------------------

#include "tsmbkind.h"

USE yoerrto1 , ONLY : kao     ,kbo     ,selfrefo   ,fracrefao ,&
           &fracrefbo  ,forrefo

!     ------------------------------------------------------------------


IMPLICIT NONE
fracrefao( :) = (/&
    &0.08452097_jprb,0.17952873_jprb,0.16214369_jprb,0.13602182_jprb,&
    &0.12760490_jprb,0.10302561_jprb,0.08392423_jprb,0.06337652_jprb,&
    &0.04206551_jprb,0.00487497_jprb,0.00410743_jprb,0.00344421_jprb,&
    &0.00285731_jprb,0.00157327_jprb,0.00080648_jprb,0.00012406_jprb/)

fracrefbo( :) = (/&
    &0.15492001_jprb,0.17384727_jprb,0.15165100_jprb,0.12675308_jprb,&
    &0.10986247_jprb,0.09006091_jprb,0.07584465_jprb,0.05990077_jprb,&
    &0.04113461_jprb,0.00438638_jprb,0.00374754_jprb,0.00313924_jprb,&
    &0.00234381_jprb,0.00167167_jprb,0.00062744_jprb,0.00010889_jprb/)

forrefo( :) = (/&
     &-4.50470e-02_jprb,-1.18908e-01_jprb,-7.21730e-02_jprb,-2.83862e-02_jprb,&
     &-3.01961e-02_jprb,-1.56877e-02_jprb,-1.53684e-02_jprb,-1.29135e-02_jprb,&
     &-1.27963e-02_jprb,-1.81742e-03_jprb, 4.40008e-05_jprb, 1.05260e-02_jprb,&
      &2.17290e-02_jprb, 1.65571e-02_jprb, 7.60751e-02_jprb, 1.47405e-01_jprb/)

!     ------------------------------------------------------------------

!     The array KAO contains absorption coefs at the 16 chosen g-values 
!     for a range of pressure levels > ~100mb and temperatures.  The first
!     index in the array, JT, which runs from 1 to 5, corresponds to 
!     different temperatures.  More specifically, JT = 3 means that the 
!     data are for the corresponding TREF for this  pressure level, 
!     JT = 2 refers to the temperatureTREF-15, JT = 1 is for TREF-30, 
!     JT = 4 is for TREF+15, and JT = 5 is for TREF+30.  The second 
!     index, JP, runs from 1 to 13 and refers to the corresponding 
!     pressure level in PREF (e.g. JP = 1 is for a pressure of 1053.63 mb).  
!     The third index, IG, goes from 1 to 16, and tells us which 
!     g-interval the absorption coefficients are for.

kao( :, 1, 1) = (/&
&2.9728e-01_jprb,2.7117e-01_jprb,2.4824e-01_jprb,2.2823e-01_jprb,2.1073e-01_jprb/)
kao( :, 2, 1) = (/&
&2.5383e-01_jprb,2.3094e-01_jprb,2.1122e-01_jprb,1.9387e-01_jprb,1.7870e-01_jprb/)
kao( :, 3, 1) = (/&
&2.2079e-01_jprb,2.0034e-01_jprb,1.8264e-01_jprb,1.6726e-01_jprb,1.5378e-01_jprb/)
kao( :, 4, 1) = (/&
&1.9385e-01_jprb,1.7531e-01_jprb,1.5932e-01_jprb,1.4544e-01_jprb,1.3338e-01_jprb/)
kao( :, 5, 1) = (/&
&1.7059e-01_jprb,1.5374e-01_jprb,1.3930e-01_jprb,1.2682e-01_jprb,1.1600e-01_jprb/)
kao( :, 6, 1) = (/&
&1.5085e-01_jprb,1.3549e-01_jprb,1.2235e-01_jprb,1.1112e-01_jprb,1.0142e-01_jprb/)
kao( :, 7, 1) = (/&
&1.3366e-01_jprb,1.1960e-01_jprb,1.0771e-01_jprb,9.7560e-02_jprb,8.8804e-02_jprb/)
kao( :, 8, 1) = (/&
&1.1833e-01_jprb,1.0551e-01_jprb,9.4748e-02_jprb,8.5642e-02_jprb,7.7808e-02_jprb/)
kao( :, 9, 1) = (/&
&1.0407e-01_jprb,9.2585e-02_jprb,8.3046e-02_jprb,7.5077e-02_jprb,6.8338e-02_jprb/)
kao( :,10, 1) = (/&
&9.6249e-02_jprb,8.6268e-02_jprb,7.7918e-02_jprb,7.0909e-02_jprb,6.5026e-02_jprb/)
kao( :,11, 1) = (/&
&8.5798e-02_jprb,7.7607e-02_jprb,7.0781e-02_jprb,6.5084e-02_jprb,6.0324e-02_jprb/)
kao( :,12, 1) = (/&
&7.5259e-02_jprb,6.8525e-02_jprb,6.2945e-02_jprb,5.8330e-02_jprb,5.4507e-02_jprb/)
kao( :,13, 1) = (/&
&6.6022e-02_jprb,6.0316e-02_jprb,5.5584e-02_jprb,5.1665e-02_jprb,4.8461e-02_jprb/)
kao( :, 1, 2) = (/&
&1.0868e+00_jprb,1.0180e+00_jprb,9.5619e-01_jprb,9.0008e-01_jprb,8.4968e-01_jprb/)
kao( :, 2, 2) = (/&
&9.1280e-01_jprb,8.5427e-01_jprb,8.0215e-01_jprb,7.5470e-01_jprb,7.1212e-01_jprb/)
kao( :, 3, 2) = (/&
&7.7699e-01_jprb,7.2607e-01_jprb,6.8084e-01_jprb,6.3997e-01_jprb,6.0293e-01_jprb/)
kao( :, 4, 2) = (/&
&6.6708e-01_jprb,6.2264e-01_jprb,5.8280e-01_jprb,5.4723e-01_jprb,5.1498e-01_jprb/)
kao( :, 5, 2) = (/&
&5.7399e-01_jprb,5.3509e-01_jprb,5.0018e-01_jprb,4.6892e-01_jprb,4.4080e-01_jprb/)
kao( :, 6, 2) = (/&
&4.9530e-01_jprb,4.6110e-01_jprb,4.3040e-01_jprb,4.0285e-01_jprb,3.7815e-01_jprb/)
kao( :, 7, 2) = (/&
&4.2745e-01_jprb,3.9732e-01_jprb,3.7040e-01_jprb,3.4614e-01_jprb,3.2442e-01_jprb/)
kao( :, 8, 2) = (/&
&3.6869e-01_jprb,3.4268e-01_jprb,3.1900e-01_jprb,2.9776e-01_jprb,2.7880e-01_jprb/)
kao( :, 9, 2) = (/&
&3.1986e-01_jprb,2.9734e-01_jprb,2.7688e-01_jprb,2.5843e-01_jprb,2.4196e-01_jprb/)
kao( :,10, 2) = (/&
&2.8291e-01_jprb,2.6327e-01_jprb,2.4583e-01_jprb,2.3014e-01_jprb,2.1607e-01_jprb/)
kao( :,11, 2) = (/&
&2.4260e-01_jprb,2.2644e-01_jprb,2.1226e-01_jprb,1.9956e-01_jprb,1.8806e-01_jprb/)
kao( :,12, 2) = (/&
&2.0558e-01_jprb,1.9256e-01_jprb,1.8115e-01_jprb,1.7097e-01_jprb,1.6176e-01_jprb/)
kao( :,13, 2) = (/&
&1.7296e-01_jprb,1.6235e-01_jprb,1.5298e-01_jprb,1.4476e-01_jprb,1.3735e-01_jprb/)
kao( :, 1, 3) = (/&
&1.9978e+00_jprb,1.9059e+00_jprb,1.8189e+00_jprb,1.7366e+00_jprb,1.6586e+00_jprb/)
kao( :, 2, 3) = (/&
&1.6606e+00_jprb,1.5839e+00_jprb,1.5107e+00_jprb,1.4431e+00_jprb,1.3779e+00_jprb/)
kao( :, 3, 3) = (/&
&1.3892e+00_jprb,1.3254e+00_jprb,1.2641e+00_jprb,1.2074e+00_jprb,1.1536e+00_jprb/)
kao( :, 4, 3) = (/&
&1.1729e+00_jprb,1.1192e+00_jprb,1.0684e+00_jprb,1.0206e+00_jprb,9.7581e-01_jprb/)
kao( :, 5, 3) = (/&
&9.9394e-01_jprb,9.4858e-01_jprb,9.0600e-01_jprb,8.6583e-01_jprb,8.2779e-01_jprb/)
kao( :, 6, 3) = (/&
&8.4450e-01_jprb,8.0530e-01_jprb,7.6915e-01_jprb,7.3568e-01_jprb,7.0371e-01_jprb/)
kao( :, 7, 3) = (/&
&7.1829e-01_jprb,6.8411e-01_jprb,6.5331e-01_jprb,6.2452e-01_jprb,5.9760e-01_jprb/)
kao( :, 8, 3) = (/&
&6.1224e-01_jprb,5.8194e-01_jprb,5.5498e-01_jprb,5.3049e-01_jprb,5.0749e-01_jprb/)
kao( :, 9, 3) = (/&
&5.2649e-01_jprb,4.9920e-01_jprb,4.7516e-01_jprb,4.5390e-01_jprb,4.3400e-01_jprb/)
kao( :,10, 3) = (/&
&4.6036e-01_jprb,4.3661e-01_jprb,4.1510e-01_jprb,3.9653e-01_jprb,3.7973e-01_jprb/)
kao( :,11, 3) = (/&
&3.9552e-01_jprb,3.7556e-01_jprb,3.5783e-01_jprb,3.4242e-01_jprb,3.2854e-01_jprb/)
kao( :,12, 3) = (/&
&3.3619e-01_jprb,3.1976e-01_jprb,3.0528e-01_jprb,2.9270e-01_jprb,2.8130e-01_jprb/)
kao( :,13, 3) = (/&
&2.8224e-01_jprb,2.6905e-01_jprb,2.5767e-01_jprb,2.4740e-01_jprb,2.3812e-01_jprb/)
kao( :, 1, 4) = (/&
&3.1134e+00_jprb,2.9768e+00_jprb,2.8528e+00_jprb,2.7371e+00_jprb,2.6298e+00_jprb/)
kao( :, 2, 4) = (/&
&2.5980e+00_jprb,2.4825e+00_jprb,2.3758e+00_jprb,2.2768e+00_jprb,2.1885e+00_jprb/)
kao( :, 3, 4) = (/&
&2.1823e+00_jprb,2.0815e+00_jprb,1.9888e+00_jprb,1.9036e+00_jprb,1.8278e+00_jprb/)
kao( :, 4, 4) = (/&
&1.8482e+00_jprb,1.7602e+00_jprb,1.6798e+00_jprb,1.6064e+00_jprb,1.5402e+00_jprb/)
kao( :, 5, 4) = (/&
&1.5691e+00_jprb,1.4937e+00_jprb,1.4240e+00_jprb,1.3603e+00_jprb,1.3038e+00_jprb/)
kao( :, 6, 4) = (/&
&1.3353e+00_jprb,1.2714e+00_jprb,1.2114e+00_jprb,1.1566e+00_jprb,1.1076e+00_jprb/)
kao( :, 7, 4) = (/&
&1.1365e+00_jprb,1.0823e+00_jprb,1.0312e+00_jprb,9.8386e-01_jprb,9.4149e-01_jprb/)
kao( :, 8, 4) = (/&
&9.6785e-01_jprb,9.2216e-01_jprb,8.7848e-01_jprb,8.3819e-01_jprb,8.0171e-01_jprb/)
kao( :, 9, 4) = (/&
&8.2906e-01_jprb,7.9098e-01_jprb,7.5375e-01_jprb,7.1927e-01_jprb,6.8769e-01_jprb/)
kao( :,10, 4) = (/&
&7.2091e-01_jprb,6.8867e-01_jprb,6.5764e-01_jprb,6.2852e-01_jprb,6.0157e-01_jprb/)
kao( :,11, 4) = (/&
&6.1639e-01_jprb,5.8956e-01_jprb,5.6471e-01_jprb,5.4117e-01_jprb,5.1965e-01_jprb/)
kao( :,12, 4) = (/&
&5.2694e-01_jprb,5.0446e-01_jprb,4.8409e-01_jprb,4.6534e-01_jprb,4.4808e-01_jprb/)
kao( :,13, 4) = (/&
&4.4495e-01_jprb,4.2623e-01_jprb,4.0962e-01_jprb,3.9458e-01_jprb,3.8091e-01_jprb/)
kao( :, 1, 5) = (/&
&5.1530e+00_jprb,4.9766e+00_jprb,4.8081e+00_jprb,4.6516e+00_jprb,4.5093e+00_jprb/)
kao( :, 2, 5) = (/&
&4.3491e+00_jprb,4.2045e+00_jprb,4.0708e+00_jprb,3.9418e+00_jprb,3.8163e+00_jprb/)
kao( :, 3, 5) = (/&
&3.6708e+00_jprb,3.5507e+00_jprb,3.4414e+00_jprb,3.3312e+00_jprb,3.2173e+00_jprb/)
kao( :, 4, 5) = (/&
&3.1175e+00_jprb,3.0159e+00_jprb,2.9223e+00_jprb,2.8267e+00_jprb,2.7279e+00_jprb/)
kao( :, 5, 5) = (/&
&2.6518e+00_jprb,2.5625e+00_jprb,2.4801e+00_jprb,2.3975e+00_jprb,2.3146e+00_jprb/)
kao( :, 6, 5) = (/&
&2.2564e+00_jprb,2.1773e+00_jprb,2.1039e+00_jprb,2.0333e+00_jprb,1.9631e+00_jprb/)
kao( :, 7, 5) = (/&
&1.9182e+00_jprb,1.8476e+00_jprb,1.7832e+00_jprb,1.7230e+00_jprb,1.6643e+00_jprb/)
kao( :, 8, 5) = (/&
&1.6337e+00_jprb,1.5705e+00_jprb,1.5139e+00_jprb,1.4617e+00_jprb,1.4128e+00_jprb/)
kao( :, 9, 5) = (/&
&1.4083e+00_jprb,1.3508e+00_jprb,1.3006e+00_jprb,1.2545e+00_jprb,1.2118e+00_jprb/)
kao( :,10, 5) = (/&
&1.2328e+00_jprb,1.1811e+00_jprb,1.1365e+00_jprb,1.0948e+00_jprb,1.0570e+00_jprb/)
kao( :,11, 5) = (/&
&1.0545e+00_jprb,1.0121e+00_jprb,9.7473e-01_jprb,9.3945e-01_jprb,9.0734e-01_jprb/)
kao( :,12, 5) = (/&
&9.0496e-01_jprb,8.7115e-01_jprb,8.4024e-01_jprb,8.1130e-01_jprb,7.8449e-01_jprb/)
kao( :,13, 5) = (/&
&7.7164e-01_jprb,7.4436e-01_jprb,7.1910e-01_jprb,6.9522e-01_jprb,6.7297e-01_jprb/)
kao( :, 1, 6) = (/&
&9.6973e+00_jprb,9.3787e+00_jprb,9.0588e+00_jprb,8.7497e+00_jprb,8.4590e+00_jprb/)
kao( :, 2, 6) = (/&
&8.3422e+00_jprb,8.0414e+00_jprb,7.7516e+00_jprb,7.4834e+00_jprb,7.2408e+00_jprb/)
kao( :, 3, 6) = (/&
&7.1314e+00_jprb,6.8664e+00_jprb,6.6149e+00_jprb,6.3897e+00_jprb,6.1962e+00_jprb/)
kao( :, 4, 6) = (/&
&6.1234e+00_jprb,5.8927e+00_jprb,5.6767e+00_jprb,5.4877e+00_jprb,5.3211e+00_jprb/)
kao( :, 5, 6) = (/&
&5.2539e+00_jprb,5.0548e+00_jprb,4.8732e+00_jprb,4.7110e+00_jprb,4.5668e+00_jprb/)
kao( :, 6, 6) = (/&
&4.5026e+00_jprb,4.3343e+00_jprb,4.1792e+00_jprb,4.0419e+00_jprb,3.9124e+00_jprb/)
kao( :, 7, 6) = (/&
&3.8524e+00_jprb,3.7093e+00_jprb,3.5777e+00_jprb,3.4570e+00_jprb,3.3411e+00_jprb/)
kao( :, 8, 6) = (/&
&3.2864e+00_jprb,3.1665e+00_jprb,3.0553e+00_jprb,2.9496e+00_jprb,2.8463e+00_jprb/)
kao( :, 9, 6) = (/&
&2.8265e+00_jprb,2.7253e+00_jprb,2.6293e+00_jprb,2.5362e+00_jprb,2.4454e+00_jprb/)
kao( :,10, 6) = (/&
&2.4703e+00_jprb,2.3851e+00_jprb,2.3012e+00_jprb,2.2197e+00_jprb,2.1390e+00_jprb/)
kao( :,11, 6) = (/&
&2.1039e+00_jprb,2.0324e+00_jprb,1.9625e+00_jprb,1.8926e+00_jprb,1.8236e+00_jprb/)
kao( :,12, 6) = (/&
&1.7889e+00_jprb,1.7295e+00_jprb,1.6689e+00_jprb,1.6099e+00_jprb,1.5530e+00_jprb/)
kao( :,13, 6) = (/&
&1.5191e+00_jprb,1.4695e+00_jprb,1.4192e+00_jprb,1.3697e+00_jprb,1.3233e+00_jprb/)
kao( :, 1, 7) = (/&
&1.9309e+01_jprb,1.8684e+01_jprb,1.8099e+01_jprb,1.7538e+01_jprb,1.7003e+01_jprb/)
kao( :, 2, 7) = (/&
&1.7226e+01_jprb,1.6672e+01_jprb,1.6118e+01_jprb,1.5594e+01_jprb,1.5067e+01_jprb/)
kao( :, 3, 7) = (/&
&1.5276e+01_jprb,1.4769e+01_jprb,1.4264e+01_jprb,1.3763e+01_jprb,1.3263e+01_jprb/)
kao( :, 4, 7) = (/&
&1.3460e+01_jprb,1.3002e+01_jprb,1.2544e+01_jprb,1.2084e+01_jprb,1.1633e+01_jprb/)
kao( :, 5, 7) = (/&
&1.1768e+01_jprb,1.1365e+01_jprb,1.0948e+01_jprb,1.0535e+01_jprb,1.0140e+01_jprb/)
kao( :, 6, 7) = (/&
&1.0233e+01_jprb,9.8691e+00_jprb,9.5008e+00_jprb,9.1399e+00_jprb,8.8037e+00_jprb/)
kao( :, 7, 7) = (/&
&8.8418e+00_jprb,8.5251e+00_jprb,8.2044e+00_jprb,7.8974e+00_jprb,7.6097e+00_jprb/)
kao( :, 8, 7) = (/&
&7.5932e+00_jprb,7.3228e+00_jprb,7.0494e+00_jprb,6.7890e+00_jprb,6.5466e+00_jprb/)
kao( :, 9, 7) = (/&
&6.5438e+00_jprb,6.3102e+00_jprb,6.0785e+00_jprb,5.8572e+00_jprb,5.6536e+00_jprb/)
kao( :,10, 7) = (/&
&5.7701e+00_jprb,5.5632e+00_jprb,5.3602e+00_jprb,5.1657e+00_jprb,4.9886e+00_jprb/)
kao( :,11, 7) = (/&
&4.9229e+00_jprb,4.7468e+00_jprb,4.5754e+00_jprb,4.4156e+00_jprb,4.2640e+00_jprb/)
kao( :,12, 7) = (/&
&4.1712e+00_jprb,4.0245e+00_jprb,3.8835e+00_jprb,3.7490e+00_jprb,3.6186e+00_jprb/)
kao( :,13, 7) = (/&
&3.5142e+00_jprb,3.3938e+00_jprb,3.2767e+00_jprb,3.1632e+00_jprb,3.0535e+00_jprb/)
kao( :, 1, 8) = (/&
&4.3792e+01_jprb,4.2562e+01_jprb,4.1326e+01_jprb,4.0173e+01_jprb,3.9029e+01_jprb/)
kao( :, 2, 8) = (/&
&4.1083e+01_jprb,3.9961e+01_jprb,3.8875e+01_jprb,3.7761e+01_jprb,3.6732e+01_jprb/)
kao( :, 3, 8) = (/&
&3.8232e+01_jprb,3.7158e+01_jprb,3.6113e+01_jprb,3.5094e+01_jprb,3.4130e+01_jprb/)
kao( :, 4, 8) = (/&
&3.5298e+01_jprb,3.4286e+01_jprb,3.3318e+01_jprb,3.2377e+01_jprb,3.1455e+01_jprb/)
kao( :, 5, 8) = (/&
&3.2307e+01_jprb,3.1366e+01_jprb,3.0496e+01_jprb,2.9612e+01_jprb,2.8733e+01_jprb/)
kao( :, 6, 8) = (/&
&2.9245e+01_jprb,2.8430e+01_jprb,2.7640e+01_jprb,2.6836e+01_jprb,2.6022e+01_jprb/)
kao( :, 7, 8) = (/&
&2.6088e+01_jprb,2.5397e+01_jprb,2.4695e+01_jprb,2.3991e+01_jprb,2.3257e+01_jprb/)
kao( :, 8, 8) = (/&
&2.2981e+01_jprb,2.2412e+01_jprb,2.1823e+01_jprb,2.1203e+01_jprb,2.0533e+01_jprb/)
kao( :, 9, 8) = (/&
&2.0055e+01_jprb,1.9585e+01_jprb,1.9088e+01_jprb,1.8539e+01_jprb,1.7936e+01_jprb/)
kao( :,10, 8) = (/&
&1.7813e+01_jprb,1.7392e+01_jprb,1.6941e+01_jprb,1.6433e+01_jprb,1.5881e+01_jprb/)
kao( :,11, 8) = (/&
&1.5517e+01_jprb,1.5141e+01_jprb,1.4721e+01_jprb,1.4250e+01_jprb,1.3762e+01_jprb/)
kao( :,12, 8) = (/&
&1.3307e+01_jprb,1.2968e+01_jprb,1.2587e+01_jprb,1.2173e+01_jprb,1.1765e+01_jprb/)
kao( :,13, 8) = (/&
&1.1246e+01_jprb,1.0950e+01_jprb,1.0619e+01_jprb,1.0276e+01_jprb,9.9339e+00_jprb/)
kao( :, 1, 9) = (/&
&1.1805e+02_jprb,1.1497e+02_jprb,1.1183e+02_jprb,1.0901e+02_jprb,1.0608e+02_jprb/)
kao( :, 2, 9) = (/&
&1.2044e+02_jprb,1.1710e+02_jprb,1.1402e+02_jprb,1.1089e+02_jprb,1.0784e+02_jprb/)
kao( :, 3, 9) = (/&
&1.2224e+02_jprb,1.1882e+02_jprb,1.1548e+02_jprb,1.1225e+02_jprb,1.0912e+02_jprb/)
kao( :, 4, 9) = (/&
&1.2253e+02_jprb,1.1914e+02_jprb,1.1570e+02_jprb,1.1248e+02_jprb,1.0936e+02_jprb/)
kao( :, 5, 9) = (/&
&1.2110e+02_jprb,1.1772e+02_jprb,1.1430e+02_jprb,1.1120e+02_jprb,1.0822e+02_jprb/)
kao( :, 6, 9) = (/&
&1.1824e+02_jprb,1.1487e+02_jprb,1.1155e+02_jprb,1.0851e+02_jprb,1.0564e+02_jprb/)
kao( :, 7, 9) = (/&
&1.1353e+02_jprb,1.1029e+02_jprb,1.0721e+02_jprb,1.0431e+02_jprb,1.0150e+02_jprb/)
kao( :, 8, 9) = (/&
&1.0728e+02_jprb,1.0430e+02_jprb,1.0145e+02_jprb,9.8717e+01_jprb,9.6110e+01_jprb/)
kao( :, 9, 9) = (/&
&9.9961e+01_jprb,9.7260e+01_jprb,9.4661e+01_jprb,9.2129e+01_jprb,8.9711e+01_jprb/)
kao( :,10, 9) = (/&
&9.2244e+01_jprb,8.9827e+01_jprb,8.7426e+01_jprb,8.5062e+01_jprb,8.2760e+01_jprb/)
kao( :,11, 9) = (/&
&8.3999e+01_jprb,8.1726e+01_jprb,7.9419e+01_jprb,7.7163e+01_jprb,7.4981e+01_jprb/)
kao( :,12, 9) = (/&
&7.5642e+01_jprb,7.3608e+01_jprb,7.1481e+01_jprb,6.9416e+01_jprb,6.7422e+01_jprb/)
kao( :,13, 9) = (/&
&6.6934e+01_jprb,6.5046e+01_jprb,6.3162e+01_jprb,6.1319e+01_jprb,5.9545e+01_jprb/)
kao( :, 1,10) = (/&
&2.5219e+02_jprb,2.4571e+02_jprb,2.3918e+02_jprb,2.3279e+02_jprb,2.2858e+02_jprb/)
kao( :, 2,10) = (/&
&2.6944e+02_jprb,2.6208e+02_jprb,2.5593e+02_jprb,2.5059e+02_jprb,2.4739e+02_jprb/)
kao( :, 3,10) = (/&
&2.8444e+02_jprb,2.7711e+02_jprb,2.7148e+02_jprb,2.6629e+02_jprb,2.6255e+02_jprb/)
kao( :, 4,10) = (/&
&3.0272e+02_jprb,2.9452e+02_jprb,2.8714e+02_jprb,2.7941e+02_jprb,2.7276e+02_jprb/)
kao( :, 5,10) = (/&
&3.1679e+02_jprb,3.0877e+02_jprb,3.0015e+02_jprb,2.9068e+02_jprb,2.8102e+02_jprb/)
kao( :, 6,10) = (/&
&3.2322e+02_jprb,3.1468e+02_jprb,3.0521e+02_jprb,2.9571e+02_jprb,2.8597e+02_jprb/)
kao( :, 7,10) = (/&
&3.2856e+02_jprb,3.1908e+02_jprb,3.0869e+02_jprb,2.9873e+02_jprb,2.8901e+02_jprb/)
kao( :, 8,10) = (/&
&3.2929e+02_jprb,3.1918e+02_jprb,3.0900e+02_jprb,2.9900e+02_jprb,2.8972e+02_jprb/)
kao( :, 9,10) = (/&
&3.2403e+02_jprb,3.1482e+02_jprb,3.0502e+02_jprb,2.9579e+02_jprb,2.8721e+02_jprb/)
kao( :,10,10) = (/&
&3.1039e+02_jprb,3.0245e+02_jprb,2.9340e+02_jprb,2.8545e+02_jprb,2.7823e+02_jprb/)
kao( :,11,10) = (/&
&2.8971e+02_jprb,2.8193e+02_jprb,2.7523e+02_jprb,2.6877e+02_jprb,2.6210e+02_jprb/)
kao( :,12,10) = (/&
&2.7119e+02_jprb,2.6399e+02_jprb,2.5760e+02_jprb,2.5064e+02_jprb,2.4326e+02_jprb/)
kao( :,13,10) = (/&
&2.5228e+02_jprb,2.4612e+02_jprb,2.3937e+02_jprb,2.3200e+02_jprb,2.2456e+02_jprb/)
kao( :, 1,11) = (/&
&3.2436e+02_jprb,3.1705e+02_jprb,3.0936e+02_jprb,3.0225e+02_jprb,2.9528e+02_jprb/)
kao( :, 2,11) = (/&
&3.6387e+02_jprb,3.5499e+02_jprb,3.4471e+02_jprb,3.3477e+02_jprb,3.2544e+02_jprb/)
kao( :, 3,11) = (/&
&3.9089e+02_jprb,3.8297e+02_jprb,3.7268e+02_jprb,3.6169e+02_jprb,3.5104e+02_jprb/)
kao( :, 4,11) = (/&
&4.0614e+02_jprb,3.9591e+02_jprb,3.8584e+02_jprb,3.7702e+02_jprb,3.6787e+02_jprb/)
kao( :, 5,11) = (/&
&4.2111e+02_jprb,4.0953e+02_jprb,3.9791e+02_jprb,3.8868e+02_jprb,3.8088e+02_jprb/)
kao( :, 6,11) = (/&
&4.4478e+02_jprb,4.3170e+02_jprb,4.1878e+02_jprb,4.0704e+02_jprb,3.9607e+02_jprb/)
kao( :, 7,11) = (/&
&4.6528e+02_jprb,4.5296e+02_jprb,4.4095e+02_jprb,4.2830e+02_jprb,4.1414e+02_jprb/)
kao( :, 8,11) = (/&
&4.8103e+02_jprb,4.6875e+02_jprb,4.5586e+02_jprb,4.4138e+02_jprb,4.2642e+02_jprb/)
kao( :, 9,11) = (/&
&4.9211e+02_jprb,4.7875e+02_jprb,4.6419e+02_jprb,4.4833e+02_jprb,4.3292e+02_jprb/)
kao( :,10,11) = (/&
&4.9428e+02_jprb,4.7899e+02_jprb,4.6351e+02_jprb,4.4779e+02_jprb,4.3245e+02_jprb/)
kao( :,11,11) = (/&
&4.7778e+02_jprb,4.6338e+02_jprb,4.4778e+02_jprb,4.3344e+02_jprb,4.2003e+02_jprb/)
kao( :,12,11) = (/&
&4.5095e+02_jprb,4.3803e+02_jprb,4.2442e+02_jprb,4.1242e+02_jprb,4.0150e+02_jprb/)
kao( :,13,11) = (/&
&4.2162e+02_jprb,4.0890e+02_jprb,3.9818e+02_jprb,3.8811e+02_jprb,3.7807e+02_jprb/)
kao( :, 1,12) = (/&
&4.0765e+02_jprb,4.0438e+02_jprb,4.0112e+02_jprb,3.9703e+02_jprb,3.9278e+02_jprb/)
kao( :, 2,12) = (/&
&4.5408e+02_jprb,4.4591e+02_jprb,4.3895e+02_jprb,4.3184e+02_jprb,4.2260e+02_jprb/)
kao( :, 3,12) = (/&
&5.0754e+02_jprb,4.9609e+02_jprb,4.8458e+02_jprb,4.7355e+02_jprb,4.6248e+02_jprb/)
kao( :, 4,12) = (/&
&5.5729e+02_jprb,5.4777e+02_jprb,5.3611e+02_jprb,5.2186e+02_jprb,5.0830e+02_jprb/)
kao( :, 5,12) = (/&
&5.9884e+02_jprb,5.9059e+02_jprb,5.8025e+02_jprb,5.6682e+02_jprb,5.5220e+02_jprb/)
kao( :, 6,12) = (/&
&6.3069e+02_jprb,6.2342e+02_jprb,6.1247e+02_jprb,5.9971e+02_jprb,5.8566e+02_jprb/)
kao( :, 7,12) = (/&
&6.6576e+02_jprb,6.5415e+02_jprb,6.3907e+02_jprb,6.2535e+02_jprb,6.1174e+02_jprb/)
kao( :, 8,12) = (/&
&7.0612e+02_jprb,6.9050e+02_jprb,6.7253e+02_jprb,6.5500e+02_jprb,6.3767e+02_jprb/)
kao( :, 9,12) = (/&
&7.4151e+02_jprb,7.2711e+02_jprb,7.0818e+02_jprb,6.8887e+02_jprb,6.6982e+02_jprb/)
kao( :,10,12) = (/&
&7.6876e+02_jprb,7.5445e+02_jprb,7.3627e+02_jprb,7.1768e+02_jprb,6.9714e+02_jprb/)
kao( :,11,12) = (/&
&7.8345e+02_jprb,7.6652e+02_jprb,7.4553e+02_jprb,7.2384e+02_jprb,7.0114e+02_jprb/)
kao( :,12,12) = (/&
&7.8623e+02_jprb,7.6501e+02_jprb,7.4277e+02_jprb,7.1917e+02_jprb,6.9536e+02_jprb/)
kao( :,13,12) = (/&
&7.6967e+02_jprb,7.4813e+02_jprb,7.2429e+02_jprb,7.0066e+02_jprb,6.7901e+02_jprb/)
kao( :, 1,13) = (/&
&5.1329e+02_jprb,5.0716e+02_jprb,5.0382e+02_jprb,4.9878e+02_jprb,4.9356e+02_jprb/)
kao( :, 2,13) = (/&
&5.7148e+02_jprb,5.6846e+02_jprb,5.6431e+02_jprb,5.6063e+02_jprb,5.5653e+02_jprb/)
kao( :, 3,13) = (/&
&6.5145e+02_jprb,6.4483e+02_jprb,6.3790e+02_jprb,6.3094e+02_jprb,6.2717e+02_jprb/)
kao( :, 4,13) = (/&
&7.5002e+02_jprb,7.3461e+02_jprb,7.2296e+02_jprb,7.1138e+02_jprb,7.0213e+02_jprb/)
kao( :, 5,13) = (/&
&8.5963e+02_jprb,8.3647e+02_jprb,8.1752e+02_jprb,7.9952e+02_jprb,7.8427e+02_jprb/)
kao( :, 6,13) = (/&
&9.6412e+02_jprb,9.4095e+02_jprb,9.1679e+02_jprb,8.9731e+02_jprb,8.7656e+02_jprb/)
kao( :, 7,13) = (/&
&1.0550e+03_jprb,1.0327e+03_jprb,1.0118e+03_jprb,9.9247e+02_jprb,9.6895e+02_jprb/)
kao( :, 8,13) = (/&
&1.1351e+03_jprb,1.1102e+03_jprb,1.0930e+03_jprb,1.0757e+03_jprb,1.0518e+03_jprb/)
kao( :, 9,13) = (/&
&1.2061e+03_jprb,1.1763e+03_jprb,1.1581e+03_jprb,1.1387e+03_jprb,1.1187e+03_jprb/)
kao( :,10,13) = (/&
&1.2746e+03_jprb,1.2431e+03_jprb,1.2211e+03_jprb,1.1969e+03_jprb,1.1743e+03_jprb/)
kao( :,11,13) = (/&
&1.3270e+03_jprb,1.3002e+03_jprb,1.2748e+03_jprb,1.2457e+03_jprb,1.2158e+03_jprb/)
kao( :,12,13) = (/&
&1.3778e+03_jprb,1.3525e+03_jprb,1.3230e+03_jprb,1.2880e+03_jprb,1.2527e+03_jprb/)
kao( :,13,13) = (/&
&1.4203e+03_jprb,1.3956e+03_jprb,1.3627e+03_jprb,1.3265e+03_jprb,1.2880e+03_jprb/)
kao( :, 1,14) = (/&
&7.0659e+02_jprb,6.8649e+02_jprb,6.7585e+02_jprb,6.6280e+02_jprb,6.4777e+02_jprb/)
kao( :, 2,14) = (/&
&8.1314e+02_jprb,7.9545e+02_jprb,7.8563e+02_jprb,7.7042e+02_jprb,7.6159e+02_jprb/)
kao( :, 3,14) = (/&
&9.2212e+02_jprb,9.1802e+02_jprb,9.0672e+02_jprb,8.9703e+02_jprb,8.8878e+02_jprb/)
kao( :, 4,14) = (/&
&1.0533e+03_jprb,1.0471e+03_jprb,1.0407e+03_jprb,1.0252e+03_jprb,1.0216e+03_jprb/)
kao( :, 5,14) = (/&
&1.1902e+03_jprb,1.1874e+03_jprb,1.1827e+03_jprb,1.1705e+03_jprb,1.1629e+03_jprb/)
kao( :, 6,14) = (/&
&1.3538e+03_jprb,1.3468e+03_jprb,1.3305e+03_jprb,1.3230e+03_jprb,1.3152e+03_jprb/)
kao( :, 7,14) = (/&
&1.5575e+03_jprb,1.5374e+03_jprb,1.5135e+03_jprb,1.4950e+03_jprb,1.4827e+03_jprb/)
kao( :, 8,14) = (/&
&1.7978e+03_jprb,1.7650e+03_jprb,1.7335e+03_jprb,1.7038e+03_jprb,1.6772e+03_jprb/)
kao( :, 9,14) = (/&
&2.0608e+03_jprb,2.0242e+03_jprb,1.9804e+03_jprb,1.9383e+03_jprb,1.9020e+03_jprb/)
kao( :,10,14) = (/&
&2.3189e+03_jprb,2.2831e+03_jprb,2.2391e+03_jprb,2.1847e+03_jprb,2.1380e+03_jprb/)
kao( :,11,14) = (/&
&2.5479e+03_jprb,2.4992e+03_jprb,2.4561e+03_jprb,2.4037e+03_jprb,2.3531e+03_jprb/)
kao( :,12,14) = (/&
&2.7503e+03_jprb,2.6955e+03_jprb,2.6487e+03_jprb,2.6000e+03_jprb,2.5509e+03_jprb/)
kao( :,13,14) = (/&
&2.9348e+03_jprb,2.8771e+03_jprb,2.8224e+03_jprb,2.7731e+03_jprb,2.7137e+03_jprb/)
kao( :, 1,15) = (/&
&1.0087e+03_jprb,9.8320e+02_jprb,9.7093e+02_jprb,9.5615e+02_jprb,9.3254e+02_jprb/)
kao( :, 2,15) = (/&
&1.2013e+03_jprb,1.1684e+03_jprb,1.1519e+03_jprb,1.1213e+03_jprb,1.1065e+03_jprb/)
kao( :, 3,15) = (/&
&1.4276e+03_jprb,1.4001e+03_jprb,1.3700e+03_jprb,1.3257e+03_jprb,1.3068e+03_jprb/)
kao( :, 4,15) = (/&
&1.6836e+03_jprb,1.6523e+03_jprb,1.6210e+03_jprb,1.5698e+03_jprb,1.5625e+03_jprb/)
kao( :, 5,15) = (/&
&1.9769e+03_jprb,1.9333e+03_jprb,1.9035e+03_jprb,1.8613e+03_jprb,1.8254e+03_jprb/)
kao( :, 6,15) = (/&
&2.3172e+03_jprb,2.2683e+03_jprb,2.2153e+03_jprb,2.1871e+03_jprb,2.1654e+03_jprb/)
kao( :, 7,15) = (/&
&2.6839e+03_jprb,2.6608e+03_jprb,2.5827e+03_jprb,2.5461e+03_jprb,2.5090e+03_jprb/)
kao( :, 8,15) = (/&
&3.0767e+03_jprb,3.0634e+03_jprb,3.0007e+03_jprb,2.9429e+03_jprb,2.9079e+03_jprb/)
kao( :, 9,15) = (/&
&3.4941e+03_jprb,3.4634e+03_jprb,3.4221e+03_jprb,3.3814e+03_jprb,3.3456e+03_jprb/)
kao( :,10,15) = (/&
&3.9537e+03_jprb,3.9378e+03_jprb,3.9008e+03_jprb,3.8414e+03_jprb,3.8169e+03_jprb/)
kao( :,11,15) = (/&
&4.5046e+03_jprb,4.4637e+03_jprb,4.4210e+03_jprb,4.3583e+03_jprb,4.3332e+03_jprb/)
kao( :,12,15) = (/&
&5.1528e+03_jprb,5.0823e+03_jprb,5.0105e+03_jprb,4.9465e+03_jprb,4.9071e+03_jprb/)
kao( :,13,15) = (/&
&5.8749e+03_jprb,5.7913e+03_jprb,5.6995e+03_jprb,5.6331e+03_jprb,5.5557e+03_jprb/)
kao( :, 1,16) = (/&
&1.2189e+03_jprb,1.2224e+03_jprb,1.2164e+03_jprb,1.2064e+03_jprb,1.1935e+03_jprb/)
kao( :, 2,16) = (/&
&1.4789e+03_jprb,1.4888e+03_jprb,1.4867e+03_jprb,1.4809e+03_jprb,1.4687e+03_jprb/)
kao( :, 3,16) = (/&
&1.8114e+03_jprb,1.8246e+03_jprb,1.8350e+03_jprb,1.8138e+03_jprb,1.8019e+03_jprb/)
kao( :, 4,16) = (/&
&2.1704e+03_jprb,2.2074e+03_jprb,2.2239e+03_jprb,2.2048e+03_jprb,2.2173e+03_jprb/)
kao( :, 5,16) = (/&
&2.5886e+03_jprb,2.6353e+03_jprb,2.6810e+03_jprb,2.6789e+03_jprb,2.6704e+03_jprb/)
kao( :, 6,16) = (/&
&3.0777e+03_jprb,3.1359e+03_jprb,3.1850e+03_jprb,3.2159e+03_jprb,3.2625e+03_jprb/)
kao( :, 7,16) = (/&
&3.6815e+03_jprb,3.7549e+03_jprb,3.7825e+03_jprb,3.8441e+03_jprb,3.9118e+03_jprb/)
kao( :, 8,16) = (/&
&4.4826e+03_jprb,4.4982e+03_jprb,4.5205e+03_jprb,4.5598e+03_jprb,4.6673e+03_jprb/)
kao( :, 9,16) = (/&
&5.4808e+03_jprb,5.4311e+03_jprb,5.4325e+03_jprb,5.4480e+03_jprb,5.4909e+03_jprb/)
kao( :,10,16) = (/&
&6.6128e+03_jprb,6.5719e+03_jprb,6.5761e+03_jprb,6.5449e+03_jprb,6.5647e+03_jprb/)
kao( :,11,16) = (/&
&7.9375e+03_jprb,7.9476e+03_jprb,7.9279e+03_jprb,7.8912e+03_jprb,7.8532e+03_jprb/)
kao( :,12,16) = (/&
&9.5403e+03_jprb,9.5531e+03_jprb,9.5122e+03_jprb,9.4385e+03_jprb,9.3648e+03_jprb/)
kao( :,13,16) = (/&
&1.1398e+04_jprb,1.1400e+04_jprb,1.1334e+04_jprb,1.1241e+04_jprb,1.1075e+04_jprb/)

!     The array SELFREFO contains the coefficient of the water vapor
!     self-continuum (including the energy term).  The first index
!     refers to temperature in 7.2 degree increments.  For instance,
!     JT = 1 refers to a temperature of 245.6, JT = 2 refers to 252.8,
!     etc.  The second index runs over the g-channel (1 to 16).

selfrefo( :, 1) = (/&
 &1.84925e+00_jprb,1.65087e+00_jprb,1.47377e+00_jprb,1.31567e+00_jprb,1.17453e+00_jprb,&
 &1.04853e+00_jprb,9.36042e-01_jprb,8.35626e-01_jprb,7.45982e-01_jprb,6.65955e-01_jprb/)
selfrefo( :, 2) = (/&
 &5.91460e+00_jprb,5.37056e+00_jprb,4.87657e+00_jprb,4.42802e+00_jprb,4.02073e+00_jprb,&
 &3.65090e+00_jprb,3.31508e+00_jprb,3.01016e+00_jprb,2.73328e+00_jprb,2.48187e+00_jprb/)
selfrefo( :, 3) = (/&
 &6.37319e+00_jprb,5.76807e+00_jprb,5.22041e+00_jprb,4.72475e+00_jprb,4.27614e+00_jprb,&
 &3.87014e+00_jprb,3.50268e+00_jprb,3.17011e+00_jprb,2.86911e+00_jprb,2.59670e+00_jprb/)
selfrefo( :, 4) = (/&
 &6.24391e+00_jprb,5.67335e+00_jprb,5.15493e+00_jprb,4.68388e+00_jprb,4.25588e+00_jprb,&
 &3.86699e+00_jprb,3.51363e+00_jprb,3.19256e+00_jprb,2.90083e+00_jprb,2.63576e+00_jprb/)
selfrefo( :, 5) = (/&
 &6.28855e+00_jprb,5.70697e+00_jprb,5.17918e+00_jprb,4.70019e+00_jprb,4.26551e+00_jprb,&
 &3.87103e+00_jprb,3.51303e+00_jprb,3.18813e+00_jprb,2.89329e+00_jprb,2.62571e+00_jprb/)
selfrefo( :, 6) = (/&
 &6.15500e+00_jprb,5.58229e+00_jprb,5.06287e+00_jprb,4.59178e+00_jprb,4.16453e+00_jprb,&
 &3.77703e+00_jprb,3.42558e+00_jprb,3.10684e+00_jprb,2.81776e+00_jprb,2.55557e+00_jprb/)
selfrefo( :, 7) = (/&
 &6.25544e+00_jprb,5.67036e+00_jprb,5.14000e+00_jprb,4.65925e+00_jprb,4.22346e+00_jprb,&
 &3.82844e+00_jprb,3.47036e+00_jprb,3.14577e+00_jprb,2.85154e+00_jprb,2.58483e+00_jprb/)
selfrefo( :, 8) = (/&
 &6.08431e+00_jprb,5.55653e+00_jprb,5.07454e+00_jprb,4.63436e+00_jprb,4.23236e+00_jprb,&
 &3.86523e+00_jprb,3.52994e+00_jprb,3.22374e+00_jprb,2.94410e+00_jprb,2.68872e+00_jprb/)
selfrefo( :, 9) = (/&
 &6.22772e+00_jprb,5.71381e+00_jprb,5.24231e+00_jprb,4.80972e+00_jprb,4.41282e+00_jprb,&
 &4.04868e+00_jprb,3.71458e+00_jprb,3.40806e+00_jprb,3.12683e+00_jprb,2.86880e+00_jprb/)
selfrefo( :,10) = (/&
 &6.63397e+00_jprb,6.10398e+00_jprb,5.61633e+00_jprb,5.16764e+00_jprb,4.75480e+00_jprb,&
 &4.37494e+00_jprb,4.02543e+00_jprb,3.70384e+00_jprb,3.40794e+00_jprb,3.13568e+00_jprb/)
selfrefo( :,11) = (/&
 &7.05802e+00_jprb,6.40896e+00_jprb,5.81959e+00_jprb,5.28442e+00_jprb,4.79846e+00_jprb,&
 &4.35719e+00_jprb,3.95651e+00_jprb,3.59266e+00_jprb,3.26228e+00_jprb,2.96228e+00_jprb/)
selfrefo( :,12) = (/&
 &7.26302e+00_jprb,6.57623e+00_jprb,5.95438e+00_jprb,5.39134e+00_jprb,4.88153e+00_jprb,&
 &4.41994e+00_jprb,4.00199e+00_jprb,3.62356e+00_jprb,3.28092e+00_jprb,2.97067e+00_jprb/)
selfrefo( :,13) = (/&
 &6.76462e+00_jprb,6.23813e+00_jprb,5.75262e+00_jprb,5.30490e+00_jprb,4.89202e+00_jprb,&
 &4.51128e+00_jprb,4.16017e+00_jprb,3.83638e+00_jprb,3.53780e+00_jprb,3.26245e+00_jprb/)
selfrefo( :,14) = (/&
 &7.35977e+00_jprb,6.71848e+00_jprb,6.13306e+00_jprb,5.59866e+00_jprb,5.11082e+00_jprb,&
 &4.66549e+00_jprb,4.25896e+00_jprb,3.88786e+00_jprb,3.54909e+00_jprb,3.23984e+00_jprb/)
selfrefo( :,15) = (/&
 &7.51978e+00_jprb,6.83155e+00_jprb,6.20632e+00_jprb,5.63830e+00_jprb,5.12227e+00_jprb,&
 &4.65347e+00_jprb,4.22758e+00_jprb,3.84066e+00_jprb,3.48915e+00_jprb,3.16982e+00_jprb/)
selfrefo( :,16) = (/&
 &6.95561e+00_jprb,6.36239e+00_jprb,5.81977e+00_jprb,5.32342e+00_jprb,4.86941e+00_jprb,&
 &4.45411e+00_jprb,4.07424e+00_jprb,3.72676e+00_jprb,3.40892e+00_jprb,3.11819e+00_jprb/)

!     -----------------------------------------------------------------
RETURN
END SUBROUTINE rrtm_kgb1_01