rrtm_kgb2_01.F90

Go to the documentation of this file.

!OCL SCALAR
SUBROUTINE rrtm_kgb2_01

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

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

#include "tsmbkind.h"

USE yoerrto2 , ONLY : kao     ,kbo     ,selfrefo   ,fracrefao  ,&
           &fracrefbo  ,forrefo
USE yoerrta2 , ONLY : refparam

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


IMPLICIT NONE
fracrefbo( :) = (/&
!     From P = 0.432 mb.
    &0.17444289_jprb,0.16467269_jprb,0.15021490_jprb,0.12460902_jprb,&
    &0.10400643_jprb,0.09481928_jprb,0.07590704_jprb,0.05752856_jprb,&
    &0.03931715_jprb,0.00428572_jprb,0.00349352_jprb,0.00278938_jprb,&
    &0.00203448_jprb,0.00130037_jprb,0.00051560_jprb,0.00006255_jprb/)

forrefo( :) = (/&
     &-2.34550e-03_jprb,-8.42698e-03_jprb,-2.01816e-02_jprb,-5.66701e-02_jprb,&
     &-8.93189e-02_jprb,-6.37487e-02_jprb,-4.56455e-02_jprb,-4.41417e-02_jprb,&
     &-4.48605e-02_jprb,-4.74696e-02_jprb,-5.16648e-02_jprb,-5.63099e-02_jprb,&
     &-4.74781e-02_jprb,-3.84704e-02_jprb,-2.49905e-02_jprb, 2.02114e-03_jprb/)

!     The following are parameters related to the reference water vapor
!     mixing ratios by REFPARAM(I) = REFH2O(I) / (.002+REFH2O(I)).
!     These parameters are used for the Planck function interpolation.
refparam( :) = (/&
 &0.903661_jprb   , 0.859386_jprb   , 0.746542_jprb   , 0.580496_jprb   , 0.412889_jprb   ,&
 &0.275283_jprb   , 0.162745_jprb   , 7.63929e-02_jprb, 1.82553e-02_jprb, 3.72432e-03_jprb, &
 &2.14946e-03_jprb, 1.66320e-03_jprb, 1.59940e-03_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) = (/&
&1.4400e-02_jprb,1.5554e-02_jprb,1.6739e-02_jprb,1.7961e-02_jprb,1.9379e-02_jprb/)
kao( :, 2, 1) = (/&
&1.1452e-02_jprb,1.2288e-02_jprb,1.3232e-02_jprb,1.4190e-02_jprb,1.5277e-02_jprb/)
kao( :, 3, 1) = (/&
&8.9989e-03_jprb,9.5960e-03_jprb,1.0327e-02_jprb,1.1068e-02_jprb,1.1889e-02_jprb/)
kao( :, 4, 1) = (/&
&7.1078e-03_jprb,7.5424e-03_jprb,8.0710e-03_jprb,8.6639e-03_jprb,9.2920e-03_jprb/)
kao( :, 5, 1) = (/&
&5.6632e-03_jprb,5.9616e-03_jprb,6.3421e-03_jprb,6.8076e-03_jprb,7.2949e-03_jprb/)
kao( :, 6, 1) = (/&
&4.5348e-03_jprb,4.7267e-03_jprb,5.0015e-03_jprb,5.3391e-03_jprb,5.7317e-03_jprb/)
kao( :, 7, 1) = (/&
&3.6484e-03_jprb,3.7737e-03_jprb,3.9583e-03_jprb,4.2050e-03_jprb,4.5069e-03_jprb/)
kao( :, 8, 1) = (/&
&2.9485e-03_jprb,3.0305e-03_jprb,3.1509e-03_jprb,3.3270e-03_jprb,3.5470e-03_jprb/)
kao( :, 9, 1) = (/&
&2.4028e-03_jprb,2.4476e-03_jprb,2.5272e-03_jprb,2.6456e-03_jprb,2.8089e-03_jprb/)
kao( :,10, 1) = (/&
&1.9899e-03_jprb,2.0047e-03_jprb,2.0583e-03_jprb,2.1403e-03_jprb,2.2643e-03_jprb/)
kao( :,11, 1) = (/&
&1.6478e-03_jprb,1.6593e-03_jprb,1.7036e-03_jprb,1.7737e-03_jprb,1.8823e-03_jprb/)
kao( :,12, 1) = (/&
&1.3599e-03_jprb,1.3694e-03_jprb,1.4067e-03_jprb,1.4672e-03_jprb,1.5617e-03_jprb/)
kao( :,13, 1) = (/&
&1.1157e-03_jprb,1.1241e-03_jprb,1.1555e-03_jprb,1.2073e-03_jprb,1.2864e-03_jprb/)
kao( :, 1, 2) = (/&
&3.8905e-02_jprb,4.2429e-02_jprb,4.6513e-02_jprb,5.0588e-02_jprb,5.5021e-02_jprb/)
kao( :, 2, 2) = (/&
&3.0803e-02_jprb,3.3461e-02_jprb,3.6546e-02_jprb,3.9922e-02_jprb,4.3440e-02_jprb/)
kao( :, 3, 2) = (/&
&2.4144e-02_jprb,2.6019e-02_jprb,2.8287e-02_jprb,3.0953e-02_jprb,3.3739e-02_jprb/)
kao( :, 4, 2) = (/&
&1.9082e-02_jprb,2.0371e-02_jprb,2.2034e-02_jprb,2.4014e-02_jprb,2.6250e-02_jprb/)
kao( :, 5, 2) = (/&
&1.5203e-02_jprb,1.6092e-02_jprb,1.7271e-02_jprb,1.8721e-02_jprb,2.0471e-02_jprb/)
kao( :, 6, 2) = (/&
&1.2164e-02_jprb,1.2769e-02_jprb,1.3578e-02_jprb,1.4643e-02_jprb,1.5928e-02_jprb/)
kao( :, 7, 2) = (/&
&9.7815e-03_jprb,1.0186e-02_jprb,1.0740e-02_jprb,1.1487e-02_jprb,1.2435e-02_jprb/)
kao( :, 8, 2) = (/&
&7.9245e-03_jprb,8.1687e-03_jprb,8.5497e-03_jprb,9.0644e-03_jprb,9.7618e-03_jprb/)
kao( :, 9, 2) = (/&
&6.4727e-03_jprb,6.5836e-03_jprb,6.8391e-03_jprb,7.1947e-03_jprb,7.6826e-03_jprb/)
kao( :,10, 2) = (/&
&5.3240e-03_jprb,5.3575e-03_jprb,5.5193e-03_jprb,5.7706e-03_jprb,6.1166e-03_jprb/)
kao( :,11, 2) = (/&
&4.3783e-03_jprb,4.4065e-03_jprb,4.5391e-03_jprb,4.7453e-03_jprb,5.0313e-03_jprb/)
kao( :,12, 2) = (/&
&3.5968e-03_jprb,3.6203e-03_jprb,3.7298e-03_jprb,3.8998e-03_jprb,4.1378e-03_jprb/)
kao( :,13, 2) = (/&
&2.9507e-03_jprb,2.9704e-03_jprb,3.0611e-03_jprb,3.2015e-03_jprb,3.4012e-03_jprb/)
kao( :, 1, 3) = (/&
&1.0636e-01_jprb,1.1680e-01_jprb,1.2824e-01_jprb,1.4017e-01_jprb,1.5266e-01_jprb/)
kao( :, 2, 3) = (/&
&8.4443e-02_jprb,9.2480e-02_jprb,1.0148e-01_jprb,1.1092e-01_jprb,1.2100e-01_jprb/)
kao( :, 3, 3) = (/&
&6.5959e-02_jprb,7.1869e-02_jprb,7.8770e-02_jprb,8.6250e-02_jprb,9.4058e-02_jprb/)
kao( :, 4, 3) = (/&
&5.1681e-02_jprb,5.6029e-02_jprb,6.1207e-02_jprb,6.7088e-02_jprb,7.3195e-02_jprb/)
kao( :, 5, 3) = (/&
&4.0734e-02_jprb,4.3938e-02_jprb,4.7777e-02_jprb,5.2288e-02_jprb,5.7243e-02_jprb/)
kao( :, 6, 3) = (/&
&3.2248e-02_jprb,3.4534e-02_jprb,3.7344e-02_jprb,4.0738e-02_jprb,4.4673e-02_jprb/)
kao( :, 7, 3) = (/&
&2.5731e-02_jprb,2.7244e-02_jprb,2.9333e-02_jprb,3.1849e-02_jprb,3.4861e-02_jprb/)
kao( :, 8, 3) = (/&
&2.0693e-02_jprb,2.1657e-02_jprb,2.3154e-02_jprb,2.5012e-02_jprb,2.7268e-02_jprb/)
kao( :, 9, 3) = (/&
&1.6730e-02_jprb,1.7346e-02_jprb,1.8327e-02_jprb,1.9716e-02_jprb,2.1381e-02_jprb/)
kao( :,10, 3) = (/&
&1.3588e-02_jprb,1.4015e-02_jprb,1.4658e-02_jprb,1.5680e-02_jprb,1.6936e-02_jprb/)
kao( :,11, 3) = (/&
&1.1180e-02_jprb,1.1528e-02_jprb,1.2063e-02_jprb,1.2917e-02_jprb,1.3959e-02_jprb/)
kao( :,12, 3) = (/&
&9.1920e-03_jprb,9.4766e-03_jprb,9.9244e-03_jprb,1.0636e-02_jprb,1.1502e-02_jprb/)
kao( :,13, 3) = (/&
&7.5517e-03_jprb,7.7861e-03_jprb,8.1602e-03_jprb,8.7504e-03_jprb,9.4682e-03_jprb/)
kao( :, 1, 4) = (/&
&2.6559e-01_jprb,2.7778e-01_jprb,2.8983e-01_jprb,3.0112e-01_jprb,3.1151e-01_jprb/)
kao( :, 2, 4) = (/&
&2.1415e-01_jprb,2.2376e-01_jprb,2.3372e-01_jprb,2.4291e-01_jprb,2.5148e-01_jprb/)
kao( :, 3, 4) = (/&
&1.7113e-01_jprb,1.7864e-01_jprb,1.8679e-01_jprb,1.9418e-01_jprb,2.0158e-01_jprb/)
kao( :, 4, 4) = (/&
&1.3689e-01_jprb,1.4305e-01_jprb,1.4933e-01_jprb,1.5577e-01_jprb,1.6212e-01_jprb/)
kao( :, 5, 4) = (/&
&1.0951e-01_jprb,1.1484e-01_jprb,1.1993e-01_jprb,1.2527e-01_jprb,1.3048e-01_jprb/)
kao( :, 6, 4) = (/&
&8.7678e-02_jprb,9.2013e-02_jprb,9.6292e-02_jprb,1.0052e-01_jprb,1.0492e-01_jprb/)
kao( :, 7, 4) = (/&
&7.0271e-02_jprb,7.3646e-02_jprb,7.7252e-02_jprb,8.0679e-02_jprb,8.4295e-02_jprb/)
kao( :, 8, 4) = (/&
&5.6292e-02_jprb,5.9027e-02_jprb,6.1915e-02_jprb,6.4784e-02_jprb,6.7662e-02_jprb/)
kao( :, 9, 4) = (/&
&4.5045e-02_jprb,4.7308e-02_jprb,4.9593e-02_jprb,5.1992e-02_jprb,5.4350e-02_jprb/)
kao( :,10, 4) = (/&
&3.6083e-02_jprb,3.8037e-02_jprb,3.9928e-02_jprb,4.1868e-02_jprb,4.3827e-02_jprb/)
kao( :,11, 4) = (/&
&2.9684e-02_jprb,3.1305e-02_jprb,3.2869e-02_jprb,3.4458e-02_jprb,3.6069e-02_jprb/)
kao( :,12, 4) = (/&
&2.4428e-02_jprb,2.5760e-02_jprb,2.7039e-02_jprb,2.8345e-02_jprb,2.9671e-02_jprb/)
kao( :,13, 4) = (/&
&2.0094e-02_jprb,2.1179e-02_jprb,2.2224e-02_jprb,2.3304e-02_jprb,2.4394e-02_jprb/)
kao( :, 1, 5) = (/&
&4.7290e-01_jprb,4.8928e-01_jprb,5.0887e-01_jprb,5.3328e-01_jprb,5.6216e-01_jprb/)
kao( :, 2, 5) = (/&
&3.8475e-01_jprb,3.9762e-01_jprb,4.1280e-01_jprb,4.3219e-01_jprb,4.5473e-01_jprb/)
kao( :, 3, 5) = (/&
&3.0979e-01_jprb,3.1996e-01_jprb,3.3105e-01_jprb,3.4530e-01_jprb,3.6216e-01_jprb/)
kao( :, 4, 5) = (/&
&2.5002e-01_jprb,2.5800e-01_jprb,2.6674e-01_jprb,2.7709e-01_jprb,2.8946e-01_jprb/)
kao( :, 5, 5) = (/&
&2.0218e-01_jprb,2.0847e-01_jprb,2.1546e-01_jprb,2.2311e-01_jprb,2.3239e-01_jprb/)
kao( :, 6, 5) = (/&
&1.6286e-01_jprb,1.6840e-01_jprb,1.7397e-01_jprb,1.7996e-01_jprb,1.8683e-01_jprb/)
kao( :, 7, 5) = (/&
&1.3055e-01_jprb,1.3602e-01_jprb,1.4050e-01_jprb,1.4524e-01_jprb,1.5047e-01_jprb/)
kao( :, 8, 5) = (/&
&1.0442e-01_jprb,1.0954e-01_jprb,1.1353e-01_jprb,1.1737e-01_jprb,1.2143e-01_jprb/)
kao( :, 9, 5) = (/&
&8.3388e-02_jprb,8.7841e-02_jprb,9.1682e-02_jprb,9.4816e-02_jprb,9.8050e-02_jprb/)
kao( :,10, 5) = (/&
&6.7010e-02_jprb,7.0598e-02_jprb,7.4015e-02_jprb,7.6769e-02_jprb,7.9372e-02_jprb/)
kao( :,11, 5) = (/&
&5.5118e-02_jprb,5.8051e-02_jprb,6.0869e-02_jprb,6.3114e-02_jprb,6.5294e-02_jprb/)
kao( :,12, 5) = (/&
&4.5300e-02_jprb,4.7721e-02_jprb,5.0049e-02_jprb,5.1890e-02_jprb,5.3682e-02_jprb/)
kao( :,13, 5) = (/&
&3.7225e-02_jprb,3.9236e-02_jprb,4.1133e-02_jprb,4.2631e-02_jprb,4.4126e-02_jprb/)
kao( :, 1, 6) = (/&
&8.4224e-01_jprb,9.2264e-01_jprb,1.0073e+00_jprb,1.0930e+00_jprb,1.1793e+00_jprb/)
kao( :, 2, 6) = (/&
&6.7595e-01_jprb,7.4232e-01_jprb,8.1248e-01_jprb,8.8302e-01_jprb,9.5616e-01_jprb/)
kao( :, 3, 6) = (/&
&5.2904e-01_jprb,5.8048e-01_jprb,6.3751e-01_jprb,6.9642e-01_jprb,7.5721e-01_jprb/)
kao( :, 4, 6) = (/&
&4.1315e-01_jprb,4.5407e-01_jprb,4.9971e-01_jprb,5.4886e-01_jprb,5.9964e-01_jprb/)
kao( :, 5, 6) = (/&
&3.2467e-01_jprb,3.5531e-01_jprb,3.9159e-01_jprb,4.3161e-01_jprb,4.7423e-01_jprb/)
kao( :, 6, 6) = (/&
&2.5574e-01_jprb,2.7784e-01_jprb,3.0579e-01_jprb,3.3817e-01_jprb,3.7278e-01_jprb/)
kao( :, 7, 6) = (/&
&2.0313e-01_jprb,2.1792e-01_jprb,2.3876e-01_jprb,2.6386e-01_jprb,2.9205e-01_jprb/)
kao( :, 8, 6) = (/&
&1.6285e-01_jprb,1.7201e-01_jprb,1.8675e-01_jprb,2.0598e-01_jprb,2.2843e-01_jprb/)
kao( :, 9, 6) = (/&
&1.3178e-01_jprb,1.3687e-01_jprb,1.4670e-01_jprb,1.6087e-01_jprb,1.7827e-01_jprb/)
kao( :,10, 6) = (/&
&1.0760e-01_jprb,1.1022e-01_jprb,1.1661e-01_jprb,1.2690e-01_jprb,1.4033e-01_jprb/)
kao( :,11, 6) = (/&
&8.8471e-02_jprb,9.0742e-02_jprb,9.6088e-02_jprb,1.0477e-01_jprb,1.1601e-01_jprb/)
kao( :,12, 6) = (/&
&7.2724e-02_jprb,7.4645e-02_jprb,7.9165e-02_jprb,8.6464e-02_jprb,9.5844e-02_jprb/)
kao( :,13, 6) = (/&
&5.9719e-02_jprb,6.1350e-02_jprb,6.5177e-02_jprb,7.1301e-02_jprb,7.9088e-02_jprb/)
kao( :, 1, 7) = (/&
&1.8416e+00_jprb,2.0377e+00_jprb,2.2351e+00_jprb,2.4351e+00_jprb,2.6357e+00_jprb/)
kao( :, 2, 7) = (/&
&1.5057e+00_jprb,1.6691e+00_jprb,1.8342e+00_jprb,2.0031e+00_jprb,2.1705e+00_jprb/)
kao( :, 3, 7) = (/&
&1.1980e+00_jprb,1.3329e+00_jprb,1.4704e+00_jprb,1.6099e+00_jprb,1.7488e+00_jprb/)
kao( :, 4, 7) = (/&
&9.4408e-01_jprb,1.0566e+00_jprb,1.1708e+00_jprb,1.2874e+00_jprb,1.4014e+00_jprb/)
kao( :, 5, 7) = (/&
&7.4052e-01_jprb,8.3332e-01_jprb,9.3007e-01_jprb,1.0263e+00_jprb,1.1221e+00_jprb/)
kao( :, 6, 7) = (/&
&5.7502e-01_jprb,6.5145e-01_jprb,7.3102e-01_jprb,8.1228e-01_jprb,8.9283e-01_jprb/)
kao( :, 7, 7) = (/&
&4.4413e-01_jprb,5.0565e-01_jprb,5.7113e-01_jprb,6.3832e-01_jprb,7.0636e-01_jprb/)
kao( :, 8, 7) = (/&
&3.4215e-01_jprb,3.9125e-01_jprb,4.4414e-01_jprb,4.9963e-01_jprb,5.5600e-01_jprb/)
kao( :, 9, 7) = (/&
&2.6278e-01_jprb,3.0155e-01_jprb,3.4409e-01_jprb,3.8940e-01_jprb,4.3634e-01_jprb/)
kao( :,10, 7) = (/&
&2.0399e-01_jprb,2.3450e-01_jprb,2.6863e-01_jprb,3.0559e-01_jprb,3.4410e-01_jprb/)
kao( :,11, 7) = (/&
&1.6899e-01_jprb,1.9452e-01_jprb,2.2306e-01_jprb,2.5367e-01_jprb,2.8600e-01_jprb/)
kao( :,12, 7) = (/&
&1.4000e-01_jprb,1.6132e-01_jprb,1.8504e-01_jprb,2.1049e-01_jprb,2.3736e-01_jprb/)
kao( :,13, 7) = (/&
&1.1587e-01_jprb,1.3364e-01_jprb,1.5332e-01_jprb,1.7448e-01_jprb,1.9695e-01_jprb/)
kao( :, 1, 8) = (/&
&4.7751e+00_jprb,5.2791e+00_jprb,5.7933e+00_jprb,6.3083e+00_jprb,6.8326e+00_jprb/)
kao( :, 2, 8) = (/&
&4.0048e+00_jprb,4.4443e+00_jprb,4.8940e+00_jprb,5.3580e+00_jprb,5.8267e+00_jprb/)
kao( :, 3, 8) = (/&
&3.2400e+00_jprb,3.6253e+00_jprb,4.0207e+00_jprb,4.4230e+00_jprb,4.8210e+00_jprb/)
kao( :, 4, 8) = (/&
&2.5883e+00_jprb,2.9210e+00_jprb,3.2656e+00_jprb,3.6154e+00_jprb,3.9641e+00_jprb/)
kao( :, 5, 8) = (/&
&2.0550e+00_jprb,2.3390e+00_jprb,2.6351e+00_jprb,2.9336e+00_jprb,3.2370e+00_jprb/)
kao( :, 6, 8) = (/&
&1.6130e+00_jprb,1.8538e+00_jprb,2.1027e+00_jprb,2.3582e+00_jprb,2.6184e+00_jprb/)
kao( :, 7, 8) = (/&
&1.2548e+00_jprb,1.4567e+00_jprb,1.6664e+00_jprb,1.8820e+00_jprb,2.1036e+00_jprb/)
kao( :, 8, 8) = (/&
&9.6687e-01_jprb,1.1348e+00_jprb,1.3110e+00_jprb,1.4930e+00_jprb,1.6805e+00_jprb/)
kao( :, 9, 8) = (/&
&7.3847e-01_jprb,8.7593e-01_jprb,1.0221e+00_jprb,1.1741e+00_jprb,1.3319e+00_jprb/)
kao( :,10, 8) = (/&
&5.6776e-01_jprb,6.8017e-01_jprb,8.0059e-01_jprb,9.2750e-01_jprb,1.0593e+00_jprb/)
kao( :,11, 8) = (/&
&4.7511e-01_jprb,5.6972e-01_jprb,6.7118e-01_jprb,7.7840e-01_jprb,8.8949e-01_jprb/)
kao( :,12, 8) = (/&
&3.9658e-01_jprb,4.7593e-01_jprb,5.6128e-01_jprb,6.5124e-01_jprb,7.4417e-01_jprb/)
kao( :,13, 8) = (/&
&3.3005e-01_jprb,3.9650e-01_jprb,4.6807e-01_jprb,5.4357e-01_jprb,6.2155e-01_jprb/)
kao( :, 1, 9) = (/&
&1.7831e+01_jprb,1.9823e+01_jprb,2.1884e+01_jprb,2.3937e+01_jprb,2.5904e+01_jprb/)
kao( :, 2, 9) = (/&
&1.5983e+01_jprb,1.7908e+01_jprb,1.9844e+01_jprb,2.1750e+01_jprb,2.3639e+01_jprb/)
kao( :, 3, 9) = (/&
&1.3838e+01_jprb,1.5652e+01_jprb,1.7428e+01_jprb,1.9195e+01_jprb,2.0946e+01_jprb/)
kao( :, 4, 9) = (/&
&1.1707e+01_jprb,1.3372e+01_jprb,1.5018e+01_jprb,1.6675e+01_jprb,1.8330e+01_jprb/)
kao( :, 5, 9) = (/&
&9.7619e+00_jprb,1.1244e+01_jprb,1.2752e+01_jprb,1.4279e+01_jprb,1.5837e+01_jprb/)
kao( :, 6, 9) = (/&
&7.9817e+00_jprb,9.2988e+00_jprb,1.0650e+01_jprb,1.2029e+01_jprb,1.3449e+01_jprb/)
kao( :, 7, 9) = (/&
&6.4478e+00_jprb,7.5864e+00_jprb,8.7717e+00_jprb,1.0008e+01_jprb,1.1289e+01_jprb/)
kao( :, 8, 9) = (/&
&5.1550e+00_jprb,6.1381e+00_jprb,7.1667e+00_jprb,8.2567e+00_jprb,9.3759e+00_jprb/)
kao( :, 9, 9) = (/&
&4.0462e+00_jprb,4.8924e+00_jprb,5.7944e+00_jprb,6.7466e+00_jprb,7.7100e+00_jprb/)
kao( :,10, 9) = (/&
&3.1920e+00_jprb,3.9133e+00_jprb,4.6912e+00_jprb,5.5091e+00_jprb,6.3432e+00_jprb/)
kao( :,11, 9) = (/&
&2.7717e+00_jprb,3.4039e+00_jprb,4.0800e+00_jprb,4.7818e+00_jprb,5.5081e+00_jprb/)
kao( :,12, 9) = (/&
&2.3978e+00_jprb,2.9436e+00_jprb,3.5235e+00_jprb,4.1293e+00_jprb,4.7600e+00_jprb/)
kao( :,13, 9) = (/&
&2.0575e+00_jprb,2.5263e+00_jprb,3.0244e+00_jprb,3.5477e+00_jprb,4.0896e+00_jprb/)
kao( :, 1,10) = (/&
&4.9748e+01_jprb,5.5112e+01_jprb,5.9938e+01_jprb,6.5023e+01_jprb,7.0458e+01_jprb/)
kao( :, 2,10) = (/&
&4.6588e+01_jprb,5.1661e+01_jprb,5.6901e+01_jprb,6.2463e+01_jprb,6.8487e+01_jprb/)
kao( :, 3,10) = (/&
&4.1628e+01_jprb,4.6699e+01_jprb,5.2314e+01_jprb,5.8353e+01_jprb,6.4458e+01_jprb/)
kao( :, 4,10) = (/&
&3.7456e+01_jprb,4.2469e+01_jprb,4.7867e+01_jprb,5.3384e+01_jprb,5.9111e+01_jprb/)
kao( :, 5,10) = (/&
&3.3436e+01_jprb,3.8793e+01_jprb,4.3992e+01_jprb,4.9077e+01_jprb,5.4046e+01_jprb/)
kao( :, 6,10) = (/&
&2.8948e+01_jprb,3.4066e+01_jprb,3.9093e+01_jprb,4.3943e+01_jprb,4.8577e+01_jprb/)
kao( :, 7,10) = (/&
&2.4411e+01_jprb,2.9169e+01_jprb,3.3772e+01_jprb,3.8328e+01_jprb,4.2674e+01_jprb/)
kao( :, 8,10) = (/&
&2.0213e+01_jprb,2.4379e+01_jprb,2.8547e+01_jprb,3.2632e+01_jprb,3.6717e+01_jprb/)
kao( :, 9,10) = (/&
&1.6770e+01_jprb,2.0275e+01_jprb,2.3726e+01_jprb,2.7289e+01_jprb,3.1149e+01_jprb/)
kao( :,10,10) = (/&
&1.3740e+01_jprb,1.6814e+01_jprb,1.9874e+01_jprb,2.3039e+01_jprb,2.6529e+01_jprb/)
kao( :,11,10) = (/&
&1.2161e+01_jprb,1.4820e+01_jprb,1.7667e+01_jprb,2.0767e+01_jprb,2.3948e+01_jprb/)
kao( :,12,10) = (/&
&1.0666e+01_jprb,1.3093e+01_jprb,1.5727e+01_jprb,1.8473e+01_jprb,2.1195e+01_jprb/)
kao( :,13,10) = (/&
&9.3915e+00_jprb,1.1557e+01_jprb,1.3847e+01_jprb,1.6172e+01_jprb,1.8558e+01_jprb/)
kao( :, 1,11) = (/&
&7.2842e+01_jprb,8.2127e+01_jprb,9.1751e+01_jprb,1.0081e+02_jprb,1.0868e+02_jprb/)
kao( :, 2,11) = (/&
&7.0912e+01_jprb,8.0211e+01_jprb,8.9133e+01_jprb,9.7767e+01_jprb,1.0537e+02_jprb/)
kao( :, 3,11) = (/&
&6.7169e+01_jprb,7.5895e+01_jprb,8.3957e+01_jprb,9.1467e+01_jprb,9.8740e+01_jprb/)
kao( :, 4,11) = (/&
&6.1046e+01_jprb,6.9086e+01_jprb,7.7071e+01_jprb,8.4625e+01_jprb,9.2325e+01_jprb/)
kao( :, 5,11) = (/&
&5.4291e+01_jprb,6.1624e+01_jprb,6.9146e+01_jprb,7.6916e+01_jprb,8.5050e+01_jprb/)
kao( :, 6,11) = (/&
&4.7940e+01_jprb,5.4754e+01_jprb,6.1943e+01_jprb,6.9739e+01_jprb,7.7765e+01_jprb/)
kao( :, 7,11) = (/&
&4.1424e+01_jprb,4.8323e+01_jprb,5.5884e+01_jprb,6.3434e+01_jprb,7.1064e+01_jprb/)
kao( :, 8,11) = (/&
&3.4808e+01_jprb,4.1718e+01_jprb,4.9240e+01_jprb,5.6728e+01_jprb,6.4032e+01_jprb/)
kao( :, 9,11) = (/&
&2.8724e+01_jprb,3.5333e+01_jprb,4.2395e+01_jprb,4.9372e+01_jprb,5.6186e+01_jprb/)
kao( :,10,11) = (/&
&2.4060e+01_jprb,2.9886e+01_jprb,3.6171e+01_jprb,4.2511e+01_jprb,4.8855e+01_jprb/)
kao( :,11,11) = (/&
&2.2466e+01_jprb,2.7842e+01_jprb,3.3276e+01_jprb,3.8707e+01_jprb,4.4411e+01_jprb/)
kao( :,12,11) = (/&
&2.0635e+01_jprb,2.5402e+01_jprb,3.0130e+01_jprb,3.5065e+01_jprb,4.0456e+01_jprb/)
kao( :,13,11) = (/&
&1.8375e+01_jprb,2.2528e+01_jprb,2.6944e+01_jprb,3.1721e+01_jprb,3.6642e+01_jprb/)
kao( :, 1,12) = (/&
&1.2814e+02_jprb,1.3938e+02_jprb,1.4994e+02_jprb,1.6150e+02_jprb,1.7273e+02_jprb/)
kao( :, 2,12) = (/&
&1.2028e+02_jprb,1.3171e+02_jprb,1.4381e+02_jprb,1.5628e+02_jprb,1.6936e+02_jprb/)
kao( :, 3,12) = (/&
&1.0908e+02_jprb,1.2221e+02_jprb,1.3662e+02_jprb,1.5108e+02_jprb,1.6591e+02_jprb/)
kao( :, 4,12) = (/&
&9.9689e+01_jprb,1.1393e+02_jprb,1.2899e+02_jprb,1.4446e+02_jprb,1.5932e+02_jprb/)
kao( :, 5,12) = (/&
&9.1574e+01_jprb,1.0647e+02_jprb,1.2179e+02_jprb,1.3681e+02_jprb,1.5090e+02_jprb/)
kao( :, 6,12) = (/&
&8.3159e+01_jprb,9.8372e+01_jprb,1.1370e+02_jprb,1.2787e+02_jprb,1.4184e+02_jprb/)
kao( :, 7,12) = (/&
&7.5130e+01_jprb,8.9279e+01_jprb,1.0335e+02_jprb,1.1696e+02_jprb,1.3088e+02_jprb/)
kao( :, 8,12) = (/&
&6.6942e+01_jprb,7.9719e+01_jprb,9.2547e+01_jprb,1.0530e+02_jprb,1.1887e+02_jprb/)
kao( :, 9,12) = (/&
&5.7722e+01_jprb,6.9837e+01_jprb,8.2082e+01_jprb,9.4593e+01_jprb,1.0770e+02_jprb/)
kao( :,10,12) = (/&
&4.9040e+01_jprb,6.0620e+01_jprb,7.2577e+01_jprb,8.5187e+01_jprb,9.7844e+01_jprb/)
kao( :,11,12) = (/&
&4.5433e+01_jprb,5.6352e+01_jprb,6.8223e+01_jprb,8.0499e+01_jprb,9.2512e+01_jprb/)
kao( :,12,12) = (/&
&4.1971e+01_jprb,5.2467e+01_jprb,6.3533e+01_jprb,7.4543e+01_jprb,8.5084e+01_jprb/)
kao( :,13,12) = (/&
&3.9013e+01_jprb,4.8568e+01_jprb,5.8272e+01_jprb,6.7741e+01_jprb,7.7194e+01_jprb/)
kao( :, 1,13) = (/&
&2.2210e+02_jprb,2.4767e+02_jprb,2.7177e+02_jprb,2.9506e+02_jprb,3.1428e+02_jprb/)
kao( :, 2,13) = (/&
&2.2387e+02_jprb,2.5269e+02_jprb,2.7872e+02_jprb,2.9937e+02_jprb,3.1679e+02_jprb/)
kao( :, 3,13) = (/&
&2.1746e+02_jprb,2.4620e+02_jprb,2.7168e+02_jprb,2.9352e+02_jprb,3.1488e+02_jprb/)
kao( :, 4,13) = (/&
&2.0527e+02_jprb,2.3513e+02_jprb,2.5998e+02_jprb,2.8303e+02_jprb,3.0506e+02_jprb/)
kao( :, 5,13) = (/&
&1.8964e+02_jprb,2.1890e+02_jprb,2.4311e+02_jprb,2.6810e+02_jprb,2.9248e+02_jprb/)
kao( :, 6,13) = (/&
&1.7090e+02_jprb,1.9758e+02_jprb,2.2284e+02_jprb,2.4839e+02_jprb,2.7463e+02_jprb/)
kao( :, 7,13) = (/&
&1.5109e+02_jprb,1.7710e+02_jprb,2.0295e+02_jprb,2.2932e+02_jprb,2.5634e+02_jprb/)
kao( :, 8,13) = (/&
&1.3270e+02_jprb,1.5826e+02_jprb,1.8450e+02_jprb,2.1121e+02_jprb,2.3873e+02_jprb/)
kao( :, 9,13) = (/&
&1.1570e+02_jprb,1.4021e+02_jprb,1.6653e+02_jprb,1.9389e+02_jprb,2.2105e+02_jprb/)
kao( :,10,13) = (/&
&1.0259e+02_jprb,1.2584e+02_jprb,1.5091e+02_jprb,1.7718e+02_jprb,2.0281e+02_jprb/)
kao( :,11,13) = (/&
&1.0076e+02_jprb,1.2370e+02_jprb,1.4756e+02_jprb,1.7147e+02_jprb,1.9575e+02_jprb/)
kao( :,12,13) = (/&
&9.7325e+01_jprb,1.1939e+02_jprb,1.4174e+02_jprb,1.6449e+02_jprb,1.8836e+02_jprb/)
kao( :,13,13) = (/&
&9.2070e+01_jprb,1.1299e+02_jprb,1.3489e+02_jprb,1.5734e+02_jprb,1.8036e+02_jprb/)
kao( :, 1,14) = (/&
&4.2241e+02_jprb,4.5820e+02_jprb,4.8815e+02_jprb,5.1730e+02_jprb,5.4765e+02_jprb/)
kao( :, 2,14) = (/&
&4.2461e+02_jprb,4.6681e+02_jprb,4.9982e+02_jprb,5.3726e+02_jprb,5.6640e+02_jprb/)
kao( :, 3,14) = (/&
&4.1742e+02_jprb,4.5924e+02_jprb,4.9930e+02_jprb,5.3630e+02_jprb,5.8090e+02_jprb/)
kao( :, 4,14) = (/&
&3.9739e+02_jprb,4.4336e+02_jprb,4.8832e+02_jprb,5.2886e+02_jprb,5.8035e+02_jprb/)
kao( :, 5,14) = (/&
&3.7398e+02_jprb,4.2048e+02_jprb,4.7438e+02_jprb,5.2486e+02_jprb,5.7134e+02_jprb/)
kao( :, 6,14) = (/&
&3.4725e+02_jprb,3.9940e+02_jprb,4.5648e+02_jprb,5.1230e+02_jprb,5.6617e+02_jprb/)
kao( :, 7,14) = (/&
&3.2072e+02_jprb,3.7728e+02_jprb,4.3646e+02_jprb,4.9362e+02_jprb,5.5200e+02_jprb/)
kao( :, 8,14) = (/&
&2.9404e+02_jprb,3.5487e+02_jprb,4.1802e+02_jprb,4.7825e+02_jprb,5.4076e+02_jprb/)
kao( :, 9,14) = (/&
&2.7034e+02_jprb,3.3232e+02_jprb,3.9722e+02_jprb,4.5997e+02_jprb,5.2290e+02_jprb/)
kao( :,10,14) = (/&
&2.5045e+02_jprb,3.1309e+02_jprb,3.7582e+02_jprb,4.4096e+02_jprb,5.0588e+02_jprb/)
kao( :,11,14) = (/&
&2.5624e+02_jprb,3.1911e+02_jprb,3.8269e+02_jprb,4.5034e+02_jprb,5.1713e+02_jprb/)
kao( :,12,14) = (/&
&2.6046e+02_jprb,3.2208e+02_jprb,3.8851e+02_jprb,4.5644e+02_jprb,5.2468e+02_jprb/)
kao( :,13,14) = (/&
&2.6376e+02_jprb,3.2552e+02_jprb,3.9097e+02_jprb,4.5804e+02_jprb,5.2614e+02_jprb/)
kao( :, 1,15) = (/&
&8.5706e+02_jprb,9.2142e+02_jprb,9.6878e+02_jprb,1.0281e+03_jprb,1.0795e+03_jprb/)
kao( :, 2,15) = (/&
&9.5189e+02_jprb,1.0335e+03_jprb,1.0913e+03_jprb,1.1632e+03_jprb,1.2011e+03_jprb/)
kao( :, 3,15) = (/&
&1.0232e+03_jprb,1.0861e+03_jprb,1.1646e+03_jprb,1.2285e+03_jprb,1.3057e+03_jprb/)
kao( :, 4,15) = (/&
&1.0392e+03_jprb,1.1356e+03_jprb,1.2544e+03_jprb,1.3115e+03_jprb,1.4338e+03_jprb/)
kao( :, 5,15) = (/&
&1.0616e+03_jprb,1.1713e+03_jprb,1.3147e+03_jprb,1.4146e+03_jprb,1.4989e+03_jprb/)
kao( :, 6,15) = (/&
&1.0788e+03_jprb,1.1886e+03_jprb,1.3302e+03_jprb,1.4680e+03_jprb,1.5638e+03_jprb/)
kao( :, 7,15) = (/&
&1.0593e+03_jprb,1.1909e+03_jprb,1.3423e+03_jprb,1.4840e+03_jprb,1.6000e+03_jprb/)
kao( :, 8,15) = (/&
&9.9873e+02_jprb,1.1655e+03_jprb,1.3439e+03_jprb,1.4887e+03_jprb,1.6318e+03_jprb/)
kao( :, 9,15) = (/&
&9.4071e+02_jprb,1.1173e+03_jprb,1.3103e+03_jprb,1.4801e+03_jprb,1.6507e+03_jprb/)
kao( :,10,15) = (/&
&8.8647e+02_jprb,1.0727e+03_jprb,1.2703e+03_jprb,1.4721e+03_jprb,1.6711e+03_jprb/)
kao( :,11,15) = (/&
&9.2131e+02_jprb,1.1303e+03_jprb,1.3379e+03_jprb,1.5675e+03_jprb,1.7706e+03_jprb/)
kao( :,12,15) = (/&
&9.6245e+02_jprb,1.1976e+03_jprb,1.4205e+03_jprb,1.6577e+03_jprb,1.8724e+03_jprb/)
kao( :,13,15) = (/&
&1.0149e+03_jprb,1.2514e+03_jprb,1.4928e+03_jprb,1.7376e+03_jprb,1.9688e+03_jprb/)
kao( :, 1,16) = (/&
&1.1893e+03_jprb,1.2688e+03_jprb,1.3351e+03_jprb,1.3959e+03_jprb,1.4433e+03_jprb/)
kao( :, 2,16) = (/&
&1.4043e+03_jprb,1.5015e+03_jprb,1.5933e+03_jprb,1.6702e+03_jprb,1.7343e+03_jprb/)
kao( :, 3,16) = (/&
&1.6203e+03_jprb,1.7272e+03_jprb,1.8338e+03_jprb,1.9414e+03_jprb,2.0288e+03_jprb/)
kao( :, 4,16) = (/&
&1.7997e+03_jprb,1.9426e+03_jprb,2.1186e+03_jprb,2.2171e+03_jprb,2.3429e+03_jprb/)
kao( :, 5,16) = (/&
&1.9840e+03_jprb,2.1906e+03_jprb,2.3875e+03_jprb,2.5360e+03_jprb,2.6578e+03_jprb/)
kao( :, 6,16) = (/&
&2.1706e+03_jprb,2.4031e+03_jprb,2.6185e+03_jprb,2.8410e+03_jprb,2.9858e+03_jprb/)
kao( :, 7,16) = (/&
&2.3367e+03_jprb,2.6044e+03_jprb,2.8764e+03_jprb,3.1365e+03_jprb,3.3254e+03_jprb/)
kao( :, 8,16) = (/&
&2.4353e+03_jprb,2.7922e+03_jprb,3.1601e+03_jprb,3.4455e+03_jprb,3.6918e+03_jprb/)
kao( :, 9,16) = (/&
&2.5375e+03_jprb,2.9571e+03_jprb,3.4043e+03_jprb,3.7567e+03_jprb,4.0876e+03_jprb/)
kao( :,10,16) = (/&
&2.6751e+03_jprb,3.1240e+03_jprb,3.6015e+03_jprb,4.0631e+03_jprb,4.4658e+03_jprb/)
kao( :,11,16) = (/&
&3.0220e+03_jprb,3.5297e+03_jprb,4.0306e+03_jprb,4.6185e+03_jprb,5.0516e+03_jprb/)
kao( :,12,16) = (/&
&3.3138e+03_jprb,4.0094e+03_jprb,4.5368e+03_jprb,5.2030e+03_jprb,5.7257e+03_jprb/)
kao( :,13,16) = (/&
&3.6770e+03_jprb,4.3534e+03_jprb,5.0748e+03_jprb,5.8067e+03_jprb,6.4888e+03_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) = (/&
&7.54232e-01_jprb, 6.87231e-01_jprb, 6.26181e-01_jprb, 5.70555e-01_jprb, 5.19870e-01_jprb,&
&4.73688e-01_jprb, 4.31608e-01_jprb, 3.93266e-01_jprb, 3.58331e-01_jprb, 3.26499e-01_jprb/)
selfrefo( :, 2) = (/&
&1.01377e+00_jprb, 9.30266e-01_jprb, 8.53637e-01_jprb, 7.83320e-01_jprb, 7.18796e-01_jprb,&
&6.59587e-01_jprb, 6.05254e-01_jprb, 5.55398e-01_jprb, 5.09648e-01_jprb, 4.67667e-01_jprb/)
selfrefo( :, 3) = (/&
&1.31954e+00_jprb, 1.22523e+00_jprb, 1.13766e+00_jprb, 1.05635e+00_jprb, 9.80853e-01_jprb,&
&9.10749e-01_jprb, 8.45656e-01_jprb, 7.85215e-01_jprb, 7.29095e-01_jprb, 6.76985e-01_jprb/)
selfrefo( :, 4) = (/&
&2.18726e+00_jprb, 2.05960e+00_jprb, 1.93940e+00_jprb, 1.82621e+00_jprb, 1.71962e+00_jprb,&
&1.61926e+00_jprb, 1.52476e+00_jprb, 1.43577e+00_jprb, 1.35197e+00_jprb, 1.27306e+00_jprb/)
selfrefo( :, 5) = (/&
&3.17440e+00_jprb, 2.92794e+00_jprb, 2.70061e+00_jprb, 2.49094e+00_jprb, 2.29754e+00_jprb,&
&2.11916e+00_jprb, 1.95463e+00_jprb, 1.80287e+00_jprb, 1.66290e+00_jprb, 1.53379e+00_jprb/)
selfrefo( :, 6) = (/&
&2.69705e+00_jprb, 2.48684e+00_jprb, 2.29302e+00_jprb, 2.11431e+00_jprb, 1.94952e+00_jprb,&
&1.79758e+00_jprb, 1.65748e+00_jprb, 1.52830e+00_jprb, 1.40919e+00_jprb, 1.29936e+00_jprb/)
selfrefo( :, 7) = (/&
&2.91893e+00_jprb, 2.68371e+00_jprb, 2.46745e+00_jprb, 2.26862e+00_jprb, 2.08581e+00_jprb,&
&1.91773e+00_jprb, 1.76319e+00_jprb, 1.62111e+00_jprb, 1.49048e+00_jprb, 1.37037e+00_jprb/)
selfrefo( :, 8) = (/&
&2.91847e+00_jprb, 2.69138e+00_jprb, 2.48195e+00_jprb, 2.28883e+00_jprb, 2.11073e+00_jprb,&
&1.94648e+00_jprb, 1.79502e+00_jprb, 1.65535e+00_jprb, 1.52654e+00_jprb, 1.40776e+00_jprb/)
selfrefo( :, 9) = (/&
&2.73994e+00_jprb, 2.54834e+00_jprb, 2.37014e+00_jprb, 2.20440e+00_jprb, 2.05026e+00_jprb,&
&1.90689e+00_jprb, 1.77354e+00_jprb, 1.64952e+00_jprb, 1.53418e+00_jprb, 1.42690e+00_jprb/)
selfrefo( :,10) = (/&
&2.81902e+00_jprb, 2.66010e+00_jprb, 2.51014e+00_jprb, 2.36863e+00_jprb, 2.23510e+00_jprb,&
&2.10909e+00_jprb, 1.99020e+00_jprb, 1.87800e+00_jprb, 1.77213e+00_jprb, 1.67222e+00_jprb/)
selfrefo( :,11) = (/&
&3.19430e+00_jprb, 2.97671e+00_jprb, 2.77395e+00_jprb, 2.58500e+00_jprb, 2.40891e+00_jprb,&
&2.24483e+00_jprb, 2.09192e+00_jprb, 1.94942e+00_jprb, 1.81663e+00_jprb, 1.69289e+00_jprb/)
selfrefo( :,12) = (/&
&3.24778e+00_jprb, 3.04535e+00_jprb, 2.85553e+00_jprb, 2.67755e+00_jprb, 2.51066e+00_jprb,&
&2.35417e+00_jprb, 2.20744e+00_jprb, 2.06985e+00_jprb, 1.94084e+00_jprb, 1.81987e+00_jprb/)
selfrefo( :,13) = (/&
&3.67142e+00_jprb, 3.39875e+00_jprb, 3.14634e+00_jprb, 2.91267e+00_jprb, 2.69636e+00_jprb,&
&2.49611e+00_jprb, 2.31073e+00_jprb, 2.13913e+00_jprb, 1.98026e+00_jprb, 1.83319e+00_jprb/)
selfrefo( :,14) = (/&
&3.73679e+00_jprb, 3.49113e+00_jprb, 3.26162e+00_jprb, 3.04719e+00_jprb, 2.84686e+00_jprb,&
&2.65970e+00_jprb, 2.48485e+00_jprb, 2.32149e+00_jprb, 2.16887e+00_jprb, 2.02628e+00_jprb/)
selfrefo( :,15) = (/&
&4.25684e+00_jprb, 3.95885e+00_jprb, 3.68171e+00_jprb, 3.42397e+00_jprb, 3.18428e+00_jprb,&
&2.96137e+00_jprb, 2.75406e+00_jprb, 2.56126e+00_jprb, 2.38196e+00_jprb, 2.21521e+00_jprb/)
selfrefo( :,16) = (/&
&4.97048e+00_jprb, 4.59255e+00_jprb, 4.24336e+00_jprb, 3.92072e+00_jprb, 3.62261e+00_jprb,&
&3.34717e+00_jprb, 3.09267e+00_jprb, 2.85752e+00_jprb, 2.64025e+00_jprb, 2.43950e+00_jprb/)

!     The ith set of reference fractions are from the ith reference
!     pressure level.
fracrefao( :, 1) = (/&
    &0.18068060_jprb,0.16803175_jprb,0.15140158_jprb,0.12221480_jprb,&
    &0.10240850_jprb,0.09330297_jprb,0.07518960_jprb,0.05611294_jprb,&
    &0.03781487_jprb,0.00387192_jprb,0.00321285_jprb,0.00244440_jprb,&
    &0.00179546_jprb,0.00107704_jprb,0.00038798_jprb,0.00005060_jprb/)
fracrefao( :, 2) = (/&
    &0.17927621_jprb,0.16731168_jprb,0.15129538_jprb,0.12328085_jprb,&
    &0.10243484_jprb,0.09354796_jprb,0.07538418_jprb,0.05633071_jprb,&
    &0.03810832_jprb,0.00398347_jprb,0.00320262_jprb,0.00250029_jprb,&
    &0.00178666_jprb,0.00111127_jprb,0.00039438_jprb,0.00005169_jprb/)
fracrefao( :, 3) = (/&
    &0.17762886_jprb,0.16638555_jprb,0.15115446_jprb,0.12470623_jprb,&
    &0.10253213_jprb,0.09383459_jprb,0.07560240_jprb,0.05646568_jprb,&
    &0.03844077_jprb,0.00409142_jprb,0.00322521_jprb,0.00254918_jprb,&
    &0.00179296_jprb,0.00113652_jprb,0.00040169_jprb,0.00005259_jprb/)
fracrefao( :, 4) = (/&
    &0.17566043_jprb,0.16539773_jprb,0.15092199_jprb,0.12571971_jprb,&
    &0.10340609_jprb,0.09426189_jprb,0.07559051_jprb,0.05678188_jprb,&
    &0.03881499_jprb,0.00414102_jprb,0.00328551_jprb,0.00258795_jprb,&
    &0.00181648_jprb,0.00115145_jprb,0.00040969_jprb,0.00005357_jprb/)
fracrefao( :, 5) = (/&
    &0.17335825_jprb,0.16442548_jprb,0.15070701_jprb,0.12667464_jprb,&
    &0.10452303_jprb,0.09450833_jprb,0.07599410_jprb,0.05706393_jprb,&
    &0.03910370_jprb,0.00417880_jprb,0.00335256_jprb,0.00261708_jprb,&
    &0.00185491_jprb,0.00116627_jprb,0.00041759_jprb,0.00005464_jprb/)
fracrefao( :, 6) = (/&
    &0.17082544_jprb,0.16321516_jprb,0.15044247_jprb,0.12797612_jprb,&
    &0.10574646_jprb,0.09470057_jprb,0.07647423_jprb,0.05738756_jprb,&
    &0.03935621_jprb,0.00423789_jprb,0.00342651_jprb,0.00264549_jprb,&
    &0.00190188_jprb,0.00118281_jprb,0.00042592_jprb,0.00005583_jprb/)
fracrefao( :, 7) = (/&
    &0.16809277_jprb,0.16193336_jprb,0.15013184_jprb,0.12937409_jprb,&
    &0.10720784_jprb,0.09485368_jprb,0.07692636_jprb,0.05771774_jprb,&
    &0.03966988_jprb,0.00427754_jprb,0.00349696_jprb,0.00268946_jprb,&
    &0.00193536_jprb,0.00120222_jprb,0.00043462_jprb,0.00005712_jprb/)
fracrefao( :, 8) = (/&
    &0.16517997_jprb,0.16059248_jprb,0.14984852_jprb,0.13079269_jprb,&
    &0.10865030_jprb,0.09492947_jprb,0.07759736_jprb,0.05812201_jprb,&
    &0.03997169_jprb,0.00432356_jprb,0.00355308_jprb,0.00274031_jprb,&
    &0.00197243_jprb,0.00122401_jprb,0.00044359_jprb,0.00005849_jprb/)
fracrefao( :, 9) = (/&
    &0.16209179_jprb,0.15912023_jprb,0.14938223_jprb,0.13198245_jprb,&
    &0.11077233_jprb,0.09487948_jprb,0.07831636_jprb,0.05863440_jprb,&
    &0.04028239_jprb,0.00436804_jprb,0.00360407_jprb,0.00279885_jprb,&
    &0.00200364_jprb,0.00124861_jprb,0.00045521_jprb,0.00005996_jprb/)
fracrefao( :,10) = (/&
    &0.15962425_jprb,0.15789343_jprb,0.14898103_jprb,0.13275230_jprb,&
    &0.11253940_jprb,0.09503502_jprb,0.07884382_jprb,0.05908009_jprb,&
    &0.04053524_jprb,0.00439971_jprb,0.00364269_jprb,0.00284965_jprb,&
    &0.00202758_jprb,0.00127076_jprb,0.00046408_jprb,0.00006114_jprb/)
fracrefao( :,11) = (/&
    &0.15926200_jprb,0.15770932_jprb,0.14891729_jprb,0.13283882_jprb,&
    &0.11276010_jprb,0.09507311_jprb,0.07892222_jprb,0.05919230_jprb,&
    &0.04054824_jprb,0.00440833_jprb,0.00365575_jprb,0.00286459_jprb,&
    &0.00203786_jprb,0.00128405_jprb,0.00046504_jprb,0.00006146_jprb/)
fracrefao( :,12) = (/&
    &0.15926351_jprb,0.15770483_jprb,0.14891177_jprb,0.13279966_jprb,&
    &0.11268171_jprb,0.09515216_jprb,0.07890341_jprb,0.05924807_jprb,&
    &0.04052851_jprb,0.00440870_jprb,0.00365425_jprb,0.00286878_jprb,&
    &0.00205747_jprb,0.00128916_jprb,0.00046589_jprb,0.00006221_jprb/)
fracrefao( :,13) = (/&
    &0.15937765_jprb,0.15775780_jprb,0.14892603_jprb,0.13273248_jprb,&
    &0.11252731_jprb,0.09521657_jprb,0.07885858_jprb,0.05927679_jprb,&
    &0.04050184_jprb,0.00440285_jprb,0.00365748_jprb,0.00286791_jprb,&
    &0.00207507_jprb,0.00129193_jprb,0.00046679_jprb,0.00006308_jprb/)

!     -----------------------------------------------------------------
RETURN
END SUBROUTINE rrtm_kgb2_01