rrtm_kgb6.F90

Go to the documentation of this file.

!OCL SCALAR
SUBROUTINE rrtm_kgb6

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

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

#include "tsmbkind.h"

USE yoerrto6 , ONLY : kao     ,selfrefo   ,fracrefao  ,cfc11adjo ,&
           &cfc12o   , absco2o 

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


IMPLICIT NONE
fracrefao( :) = (/&
!     From P = 706 mb.
    &0.13739009_jprb,0.14259538_jprb,0.14033118_jprb,0.13547136_jprb,&
    &0.12569460_jprb,0.11028396_jprb,0.08626066_jprb,0.06245148_jprb,&
    &0.04309394_jprb,0.00473551_jprb,0.00403920_jprb,0.00321695_jprb,&
    &0.00232470_jprb,0.00147662_jprb,0.00056095_jprb,0.00007373_jprb/)

!     CFC11 is multiplied by 1.385 to account for the 1060-1107 cm-1 band.
cfc11adjo( :) = (/&
          &_zero_,      _zero_, 36.7627_jprb, 150.757_jprb,    &
     &81.4109_jprb, 74.9112_jprb, 56.9325_jprb, 49.3226_jprb,  &
     &57.1074_jprb, 66.1202_jprb, 109.557_jprb, 89.0562_jprb,  &
     &149.865_jprb, 196.140_jprb, 258.393_jprb, 80.9923_jprb/)

cfc12o( :) = (/&
     &62.8368_jprb, 43.2626_jprb, 26.7549_jprb, 22.2487_jprb,&
     &23.5029_jprb, 34.8323_jprb, 26.2335_jprb, 23.2306_jprb,&
     &18.4062_jprb, 13.9534_jprb, 22.6268_jprb, 24.2604_jprb,&
     &30.0088_jprb, 26.3634_jprb, 15.8237_jprb, 57.5050_jprb/)

absco2o( :) = (/&
     &7.44852e-05_jprb, 6.29208e-05_jprb, 7.34031e-05_jprb, 6.65218e-05_jprb,&
     &7.87511e-05_jprb, 1.22489e-04_jprb, 3.39785e-04_jprb, 9.33040e-04_jprb,&
     &1.54323e-03_jprb, 4.07220e-04_jprb, 4.34332e-04_jprb, 8.76418e-05_jprb,&
     &9.80381e-05_jprb, 3.51680e-05_jprb, 5.31766e-05_jprb, 1.01542e-05_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.2568e-05_jprb,2.5521e-05_jprb,2.4392e-05_jprb,2.3686e-05_jprb,2.3850e-05_jprb/)
kao( :, 2, 1) = (/&
&1.8872e-05_jprb,1.8763e-05_jprb,1.8106e-05_jprb,1.7702e-05_jprb,1.7681e-05_jprb/)
kao( :, 3, 1) = (/&
&1.2996e-05_jprb,1.3086e-05_jprb,1.3378e-05_jprb,1.3260e-05_jprb,1.1771e-05_jprb/)
kao( :, 4, 1) = (/&
&9.6276e-06_jprb,9.9774e-06_jprb,9.6694e-06_jprb,8.4861e-06_jprb,8.9160e-06_jprb/)
kao( :, 5, 1) = (/&
&7.1690e-06_jprb,7.2984e-06_jprb,5.8013e-06_jprb,6.2225e-06_jprb,6.6628e-06_jprb/)
kao( :, 6, 1) = (/&
&5.6038e-06_jprb,4.4095e-06_jprb,4.4722e-06_jprb,4.4395e-06_jprb,5.4630e-06_jprb/)
kao( :, 7, 1) = (/&
&3.4775e-06_jprb,3.0437e-06_jprb,3.3876e-06_jprb,3.8636e-06_jprb,4.9257e-06_jprb/)
kao( :, 8, 1) = (/&
&2.7310e-06_jprb,3.2224e-06_jprb,3.6552e-06_jprb,4.3975e-06_jprb,5.5359e-06_jprb/)
kao( :, 9, 1) = (/&
&6.6325e-06_jprb,7.7076e-06_jprb,8.9070e-06_jprb,1.0179e-05_jprb,1.1909e-05_jprb/)
kao( :,10, 1) = (/&
&2.4079e-05_jprb,2.6989e-05_jprb,3.0753e-05_jprb,3.4229e-05_jprb,3.8060e-05_jprb/)
kao( :,11, 1) = (/&
&3.3970e-05_jprb,3.7965e-05_jprb,4.2993e-05_jprb,4.7630e-05_jprb,5.2629e-05_jprb/)
kao( :,12, 1) = (/&
&3.5907e-05_jprb,4.0095e-05_jprb,4.5267e-05_jprb,5.0045e-05_jprb,5.5178e-05_jprb/)
kao( :,13, 1) = (/&
&3.0600e-05_jprb,3.4162e-05_jprb,3.8530e-05_jprb,4.2579e-05_jprb,4.6937e-05_jprb/)
kao( :, 1, 2) = (/&
&2.9639e-05_jprb,3.3221e-05_jprb,3.2806e-05_jprb,2.8964e-05_jprb,2.6726e-05_jprb/)
kao( :, 2, 2) = (/&
&2.3653e-05_jprb,2.4355e-05_jprb,2.4028e-05_jprb,2.0887e-05_jprb,2.2925e-05_jprb/)
kao( :, 3, 2) = (/&
&1.6734e-05_jprb,1.6216e-05_jprb,1.4570e-05_jprb,1.6223e-05_jprb,1.7762e-05_jprb/)
kao( :, 4, 2) = (/&
&1.0615e-05_jprb,1.0115e-05_jprb,1.1508e-05_jprb,1.1885e-05_jprb,1.3929e-05_jprb/)
kao( :, 5, 2) = (/&
&7.5406e-06_jprb,7.3085e-06_jprb,8.6006e-06_jprb,9.3888e-06_jprb,1.1719e-05_jprb/)
kao( :, 6, 2) = (/&
&4.8289e-06_jprb,6.2547e-06_jprb,5.8818e-06_jprb,7.9263e-06_jprb,1.0542e-05_jprb/)
kao( :, 7, 2) = (/&
&4.3622e-06_jprb,4.7605e-06_jprb,5.4656e-06_jprb,7.2224e-06_jprb,1.0212e-05_jprb/)
kao( :, 8, 2) = (/&
&4.2923e-06_jprb,4.6584e-06_jprb,6.2268e-06_jprb,8.3462e-06_jprb,1.1797e-05_jprb/)
kao( :, 9, 2) = (/&
&9.7259e-06_jprb,1.1243e-05_jprb,1.3832e-05_jprb,1.8290e-05_jprb,2.4869e-05_jprb/)
kao( :,10, 2) = (/&
&3.3596e-05_jprb,3.8056e-05_jprb,4.3895e-05_jprb,5.4699e-05_jprb,7.1889e-05_jprb/)
kao( :,11, 2) = (/&
&4.7132e-05_jprb,5.3142e-05_jprb,6.1177e-05_jprb,7.5188e-05_jprb,9.8046e-05_jprb/)
kao( :,12, 2) = (/&
&4.9755e-05_jprb,5.6067e-05_jprb,6.4461e-05_jprb,7.8769e-05_jprb,1.0238e-04_jprb/)
kao( :,13, 2) = (/&
&4.2448e-05_jprb,4.7895e-05_jprb,5.4990e-05_jprb,6.7099e-05_jprb,8.7152e-05_jprb/)
kao( :, 1, 3) = (/&
&2.4391e-05_jprb,2.6179e-05_jprb,3.0228e-05_jprb,3.9877e-05_jprb,4.1318e-05_jprb/)
kao( :, 2, 3) = (/&
&2.0506e-05_jprb,2.0015e-05_jprb,2.6885e-05_jprb,2.9284e-05_jprb,2.9837e-05_jprb/)
kao( :, 3, 3) = (/&
&1.3542e-05_jprb,1.5721e-05_jprb,1.9501e-05_jprb,1.8947e-05_jprb,2.4397e-05_jprb/)
kao( :, 4, 3) = (/&
&9.0261e-06_jprb,1.3633e-05_jprb,1.2485e-05_jprb,1.7278e-05_jprb,2.0007e-05_jprb/)
kao( :, 5, 3) = (/&
&8.7945e-06_jprb,9.7631e-06_jprb,1.1634e-05_jprb,1.3600e-05_jprb,1.9745e-05_jprb/)
kao( :, 6, 3) = (/&
&7.5275e-06_jprb,7.1532e-06_jprb,9.4521e-06_jprb,1.3168e-05_jprb,1.9311e-05_jprb/)
kao( :, 7, 3) = (/&
&5.6021e-06_jprb,6.6005e-06_jprb,8.8148e-06_jprb,1.3299e-05_jprb,1.9798e-05_jprb/)
kao( :, 8, 3) = (/&
&5.6076e-06_jprb,7.6899e-06_jprb,1.0237e-05_jprb,1.5421e-05_jprb,2.3499e-05_jprb/)
kao( :, 9, 3) = (/&
&1.3625e-05_jprb,1.7122e-05_jprb,2.2283e-05_jprb,3.1506e-05_jprb,4.8193e-05_jprb/)
kao( :,10, 3) = (/&
&4.6234e-05_jprb,5.4855e-05_jprb,6.8590e-05_jprb,8.9964e-05_jprb,1.3093e-04_jprb/)
kao( :,11, 3) = (/&
&6.5003e-05_jprb,7.6686e-05_jprb,9.5408e-05_jprb,1.2507e-04_jprb,1.8079e-04_jprb/)
kao( :,12, 3) = (/&
&6.8759e-05_jprb,8.1014e-05_jprb,1.0079e-04_jprb,1.3196e-04_jprb,1.9017e-04_jprb/)
kao( :,13, 3) = (/&
&5.8741e-05_jprb,6.9297e-05_jprb,8.6397e-05_jprb,1.1292e-04_jprb,1.6288e-04_jprb/)
kao( :, 1, 4) = (/&
&3.6952e-05_jprb,3.9167e-05_jprb,4.5764e-05_jprb,4.7701e-05_jprb,5.5871e-05_jprb/)
kao( :, 2, 4) = (/&
&2.6008e-05_jprb,3.6209e-05_jprb,2.8658e-05_jprb,3.8449e-05_jprb,3.8668e-05_jprb/)
kao( :, 3, 4) = (/&
&2.1596e-05_jprb,2.1370e-05_jprb,2.2422e-05_jprb,2.7123e-05_jprb,3.0425e-05_jprb/)
kao( :, 4, 4) = (/&
&1.4408e-05_jprb,1.3094e-05_jprb,1.9460e-05_jprb,2.1634e-05_jprb,3.3208e-05_jprb/)
kao( :, 5, 4) = (/&
&8.9120e-06_jprb,1.2874e-05_jprb,1.5169e-05_jprb,2.3591e-05_jprb,3.5032e-05_jprb/)
kao( :, 6, 4) = (/&
&7.0427e-06_jprb,1.1114e-05_jprb,1.4987e-05_jprb,2.2856e-05_jprb,3.4461e-05_jprb/)
kao( :, 7, 4) = (/&
&7.4919e-06_jprb,9.5389e-06_jprb,1.5038e-05_jprb,2.3474e-05_jprb,3.5967e-05_jprb/)
kao( :, 8, 4) = (/&
&8.6186e-06_jprb,1.1413e-05_jprb,1.7846e-05_jprb,2.8336e-05_jprb,4.4637e-05_jprb/)
kao( :, 9, 4) = (/&
&1.9661e-05_jprb,2.5521e-05_jprb,3.7417e-05_jprb,5.9323e-05_jprb,9.6439e-05_jprb/)
kao( :,10, 4) = (/&
&6.3994e-05_jprb,8.1146e-05_jprb,1.0850e-04_jprb,1.6659e-04_jprb,2.7179e-04_jprb/)
kao( :,11, 4) = (/&
&9.0362e-05_jprb,1.1392e-04_jprb,1.5220e-04_jprb,2.3565e-04_jprb,3.8761e-04_jprb/)
kao( :,12, 4) = (/&
&9.5893e-05_jprb,1.2069e-04_jprb,1.6199e-04_jprb,2.5340e-04_jprb,4.1839e-04_jprb/)
kao( :,13, 4) = (/&
&8.2144e-05_jprb,1.0344e-04_jprb,1.3998e-04_jprb,2.2132e-04_jprb,3.6483e-04_jprb/)
kao( :, 1, 5) = (/&
&3.3532e-05_jprb,4.1940e-05_jprb,5.7696e-05_jprb,6.0186e-05_jprb,6.3029e-05_jprb/)
kao( :, 2, 5) = (/&
&2.3830e-05_jprb,3.5587e-05_jprb,4.6056e-05_jprb,4.1505e-05_jprb,5.7960e-05_jprb/)
kao( :, 3, 5) = (/&
&2.1812e-05_jprb,2.5706e-05_jprb,2.8233e-05_jprb,3.5229e-05_jprb,5.8555e-05_jprb/)
kao( :, 4, 5) = (/&
&1.7315e-05_jprb,1.9695e-05_jprb,2.2530e-05_jprb,3.8597e-05_jprb,5.8973e-05_jprb/)
kao( :, 5, 5) = (/&
&1.1890e-05_jprb,1.3867e-05_jprb,2.3935e-05_jprb,3.9940e-05_jprb,5.8993e-05_jprb/)
kao( :, 6, 5) = (/&
&1.0677e-05_jprb,1.4192e-05_jprb,2.4839e-05_jprb,4.0149e-05_jprb,5.9735e-05_jprb/)
kao( :, 7, 5) = (/&
&9.2784e-06_jprb,1.5528e-05_jprb,2.6050e-05_jprb,4.1864e-05_jprb,6.4186e-05_jprb/)
kao( :, 8, 5) = (/&
&1.1359e-05_jprb,1.8672e-05_jprb,3.1354e-05_jprb,5.1149e-05_jprb,8.0918e-05_jprb/)
kao( :, 9, 5) = (/&
&2.6841e-05_jprb,4.0944e-05_jprb,6.8221e-05_jprb,1.1515e-04_jprb,1.8821e-04_jprb/)
kao( :,10, 5) = (/&
&8.7415e-05_jprb,1.2053e-04_jprb,1.9441e-04_jprb,3.3166e-04_jprb,5.5740e-04_jprb/)
kao( :,11, 5) = (/&
&1.2386e-04_jprb,1.7143e-04_jprb,2.7887e-04_jprb,4.7716e-04_jprb,8.0477e-04_jprb/)
kao( :,12, 5) = (/&
&1.3204e-04_jprb,1.8364e-04_jprb,3.0111e-04_jprb,5.1685e-04_jprb,8.7520e-04_jprb/)
kao( :,13, 5) = (/&
&1.1358e-04_jprb,1.5915e-04_jprb,2.6323e-04_jprb,4.5343e-04_jprb,7.7073e-04_jprb/)
kao( :, 1, 6) = (/&
&5.6626e-05_jprb,6.2952e-05_jprb,5.6997e-05_jprb,7.8387e-05_jprb,1.4023e-04_jprb/)
kao( :, 2, 6) = (/&
&4.2639e-05_jprb,4.0472e-05_jprb,5.1521e-05_jprb,7.2367e-05_jprb,1.2444e-04_jprb/)
kao( :, 3, 6) = (/&
&3.3191e-05_jprb,2.9821e-05_jprb,4.1035e-05_jprb,7.5835e-05_jprb,1.1292e-04_jprb/)
kao( :, 4, 6) = (/&
&1.7563e-05_jprb,2.2621e-05_jprb,4.3898e-05_jprb,7.3541e-05_jprb,1.0514e-04_jprb/)
kao( :, 5, 6) = (/&
&1.4086e-05_jprb,2.3545e-05_jprb,4.5696e-05_jprb,7.1090e-05_jprb,1.0366e-04_jprb/)
kao( :, 6, 6) = (/&
&1.2590e-05_jprb,2.4048e-05_jprb,4.4731e-05_jprb,6.9646e-05_jprb,1.0367e-04_jprb/)
kao( :, 7, 6) = (/&
&1.3429e-05_jprb,2.6924e-05_jprb,4.6412e-05_jprb,7.4051e-05_jprb,1.1343e-04_jprb/)
kao( :, 8, 6) = (/&
&1.8185e-05_jprb,3.3168e-05_jprb,5.7618e-05_jprb,9.5901e-05_jprb,1.5240e-04_jprb/)
kao( :, 9, 6) = (/&
&4.0915e-05_jprb,7.1523e-05_jprb,1.2899e-04_jprb,2.2341e-04_jprb,3.6849e-04_jprb/)
kao( :,10, 6) = (/&
&1.2433e-04_jprb,2.0393e-04_jprb,3.7307e-04_jprb,6.6855e-04_jprb,1.1425e-03_jprb/)
kao( :,11, 6) = (/&
&1.8043e-04_jprb,3.0126e-04_jprb,5.5504e-04_jprb,9.9966e-04_jprb,1.7159e-03_jprb/)
kao( :,12, 6) = (/&
&1.9679e-04_jprb,3.3338e-04_jprb,6.1765e-04_jprb,1.1153e-03_jprb,1.9175e-03_jprb/)
kao( :,13, 6) = (/&
&1.7363e-04_jprb,2.9831e-04_jprb,5.5435e-04_jprb,1.0025e-03_jprb,1.7267e-03_jprb/)
kao( :, 1, 7) = (/&
&5.4728e-05_jprb,6.3216e-05_jprb,9.7227e-05_jprb,1.9830e-04_jprb,3.0115e-04_jprb/)
kao( :, 2, 7) = (/&
&3.7424e-05_jprb,4.6891e-05_jprb,9.3922e-05_jprb,1.8575e-04_jprb,2.6393e-04_jprb/)
kao( :, 3, 7) = (/&
&3.2885e-05_jprb,5.4240e-05_jprb,1.1210e-04_jprb,1.7073e-04_jprb,2.3633e-04_jprb/)
kao( :, 4, 7) = (/&
&2.6560e-05_jprb,6.2341e-05_jprb,1.0718e-04_jprb,1.5691e-04_jprb,2.1887e-04_jprb/)
kao( :, 5, 7) = (/&
&2.8005e-05_jprb,6.1158e-05_jprb,9.7995e-05_jprb,1.4461e-04_jprb,2.0623e-04_jprb/)
kao( :, 6, 7) = (/&
&2.6304e-05_jprb,5.4858e-05_jprb,8.8237e-05_jprb,1.3524e-04_jprb,2.0068e-04_jprb/)
kao( :, 7, 7) = (/&
&2.8057e-05_jprb,5.2849e-05_jprb,8.8492e-05_jprb,1.4199e-04_jprb,2.1851e-04_jprb/)
kao( :, 8, 7) = (/&
&3.4390e-05_jprb,6.3412e-05_jprb,1.1163e-04_jprb,1.8592e-04_jprb,2.9376e-04_jprb/)
kao( :, 9, 7) = (/&
&7.2790e-05_jprb,1.4039e-04_jprb,2.5926e-04_jprb,4.5257e-04_jprb,7.4795e-04_jprb/)
kao( :,10, 7) = (/&
&2.1021e-04_jprb,4.1798e-04_jprb,8.0790e-04_jprb,1.4658e-03_jprb,2.5048e-03_jprb/)
kao( :,11, 7) = (/&
&3.2720e-04_jprb,6.6181e-04_jprb,1.2928e-03_jprb,2.3683e-03_jprb,4.0754e-03_jprb/)
kao( :,12, 7) = (/&
&3.7693e-04_jprb,7.7171e-04_jprb,1.5176e-03_jprb,2.7982e-03_jprb,4.8311e-03_jprb/)
kao( :,13, 7) = (/&
&3.5194e-04_jprb,7.2710e-04_jprb,1.4369e-03_jprb,2.6518e-03_jprb,4.5779e-03_jprb/)
kao( :, 1, 8) = (/&
&8.8622e-05_jprb,2.2200e-04_jprb,4.1760e-04_jprb,5.9916e-04_jprb,8.1306e-04_jprb/)
kao( :, 2, 8) = (/&
&7.7484e-05_jprb,1.9083e-04_jprb,3.2988e-04_jprb,4.7824e-04_jprb,6.5145e-04_jprb/)
kao( :, 3, 8) = (/&
&6.8249e-05_jprb,1.8888e-04_jprb,2.9496e-04_jprb,4.2064e-04_jprb,5.7727e-04_jprb/)
kao( :, 4, 8) = (/&
&1.0307e-04_jprb,1.8157e-04_jprb,2.7450e-04_jprb,3.9411e-04_jprb,5.4734e-04_jprb/)
kao( :, 5, 8) = (/&
&1.0013e-04_jprb,1.6578e-04_jprb,2.5412e-04_jprb,3.7165e-04_jprb,5.2629e-04_jprb/)
kao( :, 6, 8) = (/&
&8.7661e-05_jprb,1.4688e-04_jprb,2.3121e-04_jprb,3.4816e-04_jprb,5.0439e-04_jprb/)
kao( :, 7, 8) = (/&
&7.9171e-05_jprb,1.3672e-04_jprb,2.2306e-04_jprb,3.4738e-04_jprb,5.2076e-04_jprb/)
kao( :, 8, 8) = (/&
&8.3889e-05_jprb,1.5157e-04_jprb,2.5895e-04_jprb,4.2268e-04_jprb,6.6361e-04_jprb/)
kao( :, 9, 8) = (/&
&1.6488e-04_jprb,3.1976e-04_jprb,5.8317e-04_jprb,1.0082e-03_jprb,1.6627e-03_jprb/)
kao( :,10, 8) = (/&
&4.7725e-04_jprb,9.7517e-04_jprb,1.8746e-03_jprb,3.3981e-03_jprb,5.7914e-03_jprb/)
kao( :,11, 8) = (/&
&7.9688e-04_jprb,1.6560e-03_jprb,3.2130e-03_jprb,5.8470e-03_jprb,9.9511e-03_jprb/)
kao( :,12, 8) = (/&
&9.7776e-04_jprb,2.0671e-03_jprb,4.0593e-03_jprb,7.3885e-03_jprb,1.2564e-02_jprb/)
kao( :,13, 8) = (/&
&9.8124e-04_jprb,2.1024e-03_jprb,4.1582e-03_jprb,7.5870e-03_jprb,1.2921e-02_jprb/)
kao( :, 1, 9) = (/&
&7.4912e-04_jprb,1.1099e-03_jprb,1.5644e-03_jprb,2.1464e-03_jprb,2.8897e-03_jprb/)
kao( :, 2, 9) = (/&
&6.2597e-04_jprb,9.2912e-04_jprb,1.3142e-03_jprb,1.8181e-03_jprb,2.4773e-03_jprb/)
kao( :, 3, 9) = (/&
&4.7892e-04_jprb,7.1171e-04_jprb,1.0243e-03_jprb,1.4397e-03_jprb,1.9820e-03_jprb/)
kao( :, 4, 9) = (/&
&3.8453e-04_jprb,5.9140e-04_jprb,8.7342e-04_jprb,1.2581e-03_jprb,1.7584e-03_jprb/)
kao( :, 5, 9) = (/&
&3.6098e-04_jprb,5.7696e-04_jprb,8.8197e-04_jprb,1.2954e-03_jprb,1.8431e-03_jprb/)
kao( :, 6, 9) = (/&
&3.4542e-04_jprb,5.7060e-04_jprb,8.9328e-04_jprb,1.3402e-03_jprb,1.9357e-03_jprb/)
kao( :, 7, 9) = (/&
&3.3888e-04_jprb,5.8373e-04_jprb,9.4422e-04_jprb,1.4513e-03_jprb,2.1408e-03_jprb/)
kao( :, 8, 9) = (/&
&3.7707e-04_jprb,6.8408e-04_jprb,1.1538e-03_jprb,1.8357e-03_jprb,2.7831e-03_jprb/)
kao( :, 9, 9) = (/&
&7.0739e-04_jprb,1.3735e-03_jprb,2.4588e-03_jprb,4.1112e-03_jprb,6.4788e-03_jprb/)
kao( :,10, 9) = (/&
&1.8575e-03_jprb,3.8603e-03_jprb,7.2941e-03_jprb,1.2720e-02_jprb,2.0839e-02_jprb/)
kao( :,11, 9) = (/&
&2.9958e-03_jprb,6.2197e-03_jprb,1.1775e-02_jprb,2.0597e-02_jprb,3.3906e-02_jprb/)
kao( :,12, 9) = (/&
&3.6668e-03_jprb,7.6016e-03_jprb,1.4363e-02_jprb,2.5216e-02_jprb,4.1719e-02_jprb/)
kao( :,13, 9) = (/&
&3.7360e-03_jprb,7.7341e-03_jprb,1.4600e-02_jprb,2.5684e-02_jprb,4.2617e-02_jprb/)
kao( :, 1,10) = (/&
&2.2204e-03_jprb,3.1215e-03_jprb,4.1592e-03_jprb,5.5251e-03_jprb,7.2462e-03_jprb/)
kao( :, 2,10) = (/&
&1.8794e-03_jprb,2.6515e-03_jprb,3.6524e-03_jprb,4.9115e-03_jprb,6.4747e-03_jprb/)
kao( :, 3,10) = (/&
&1.4689e-03_jprb,2.0942e-03_jprb,2.9433e-03_jprb,4.0441e-03_jprb,5.4540e-03_jprb/)
kao( :, 4,10) = (/&
&1.1490e-03_jprb,1.6495e-03_jprb,2.3294e-03_jprb,3.2354e-03_jprb,4.4282e-03_jprb/)
kao( :, 5,10) = (/&
&8.6570e-04_jprb,1.3179e-03_jprb,1.9604e-03_jprb,2.8379e-03_jprb,4.0441e-03_jprb/)
kao( :, 6,10) = (/&
&8.5652e-04_jprb,1.4272e-03_jprb,2.2416e-03_jprb,3.3643e-03_jprb,4.8824e-03_jprb/)
kao( :, 7,10) = (/&
&9.3347e-04_jprb,1.6326e-03_jprb,2.6435e-03_jprb,4.0816e-03_jprb,6.0203e-03_jprb/)
kao( :, 8,10) = (/&
&1.1389e-03_jprb,2.0806e-03_jprb,3.5187e-03_jprb,5.6514e-03_jprb,8.6161e-03_jprb/)
kao( :, 9,10) = (/&
&2.3318e-03_jprb,4.5800e-03_jprb,8.2437e-03_jprb,1.3780e-02_jprb,2.1744e-02_jprb/)
kao( :,10,10) = (/&
&6.2326e-03_jprb,1.3165e-02_jprb,2.5232e-02_jprb,4.4615e-02_jprb,7.3492e-02_jprb/)
kao( :,11,10) = (/&
&9.8704e-03_jprb,2.0942e-02_jprb,4.0187e-02_jprb,7.0789e-02_jprb,1.1609e-01_jprb/)
kao( :,12,10) = (/&
&1.1999e-02_jprb,2.5297e-02_jprb,4.8337e-02_jprb,8.4705e-02_jprb,1.3826e-01_jprb/)
kao( :,13,10) = (/&
&1.1789e-02_jprb,2.4642e-02_jprb,4.6747e-02_jprb,8.1645e-02_jprb,1.3284e-01_jprb/)
kao( :, 1,11) = (/&
&3.1261e-03_jprb,4.2813e-03_jprb,5.7811e-03_jprb,7.6623e-03_jprb,1.0021e-02_jprb/)
kao( :, 2,11) = (/&
&2.7730e-03_jprb,3.9256e-03_jprb,5.2937e-03_jprb,7.0385e-03_jprb,9.1561e-03_jprb/)
kao( :, 3,11) = (/&
&2.2211e-03_jprb,3.2332e-03_jprb,4.4680e-03_jprb,6.0516e-03_jprb,7.9945e-03_jprb/)
kao( :, 4,11) = (/&
&1.7260e-03_jprb,2.5271e-03_jprb,3.5974e-03_jprb,5.0208e-03_jprb,6.7957e-03_jprb/)
kao( :, 5,11) = (/&
&1.3640e-03_jprb,2.0174e-03_jprb,2.8988e-03_jprb,4.1037e-03_jprb,5.6676e-03_jprb/)
kao( :, 6,11) = (/&
&1.1158e-03_jprb,1.8312e-03_jprb,2.8771e-03_jprb,4.3425e-03_jprb,6.2982e-03_jprb/)
kao( :, 7,11) = (/&
&1.2321e-03_jprb,2.1507e-03_jprb,3.5466e-03_jprb,5.5748e-03_jprb,8.3259e-03_jprb/)
kao( :, 8,11) = (/&
&1.5691e-03_jprb,2.9089e-03_jprb,5.0376e-03_jprb,8.1424e-03_jprb,1.2407e-02_jprb/)
kao( :, 9,11) = (/&
&3.5241e-03_jprb,7.0356e-03_jprb,1.2741e-02_jprb,2.1461e-02_jprb,3.3839e-02_jprb/)
kao( :,10,11) = (/&
&1.0170e-02_jprb,2.1914e-02_jprb,4.2506e-02_jprb,7.5273e-02_jprb,1.2386e-01_jprb/)
kao( :,11,11) = (/&
&1.6549e-02_jprb,3.5671e-02_jprb,6.9003e-02_jprb,1.2202e-01_jprb,2.0089e-01_jprb/)
kao( :,12,11) = (/&
&2.0073e-02_jprb,4.2914e-02_jprb,8.2517e-02_jprb,1.4566e-01_jprb,2.3923e-01_jprb/)
kao( :,13,11) = (/&
&1.9754e-02_jprb,4.1795e-02_jprb,8.0135e-02_jprb,1.4074e-01_jprb,2.2980e-01_jprb/)
kao( :, 1,12) = (/&
&4.1564e-03_jprb,5.9925e-03_jprb,8.3120e-03_jprb,1.1209e-02_jprb,1.4597e-02_jprb/)
kao( :, 2,12) = (/&
&3.9567e-03_jprb,5.4909e-03_jprb,7.6176e-03_jprb,1.0371e-02_jprb,1.3856e-02_jprb/)
kao( :, 3,12) = (/&
&3.3974e-03_jprb,4.7965e-03_jprb,6.6252e-03_jprb,8.9930e-03_jprb,1.2050e-02_jprb/)
kao( :, 4,12) = (/&
&2.7271e-03_jprb,4.0009e-03_jprb,5.6537e-03_jprb,7.7402e-03_jprb,1.0333e-02_jprb/)
kao( :, 5,12) = (/&
&2.1241e-03_jprb,3.2132e-03_jprb,4.6676e-03_jprb,6.6066e-03_jprb,9.0309e-03_jprb/)
kao( :, 6,12) = (/&
&1.6495e-03_jprb,2.5482e-03_jprb,3.8834e-03_jprb,5.7278e-03_jprb,8.1979e-03_jprb/)
kao( :, 7,12) = (/&
&1.5776e-03_jprb,2.7768e-03_jprb,4.6314e-03_jprb,7.2948e-03_jprb,1.0930e-02_jprb/)
kao( :, 8,12) = (/&
&2.1394e-03_jprb,4.0088e-03_jprb,6.9630e-03_jprb,1.1391e-02_jprb,1.7592e-02_jprb/)
kao( :, 9,12) = (/&
&5.1508e-03_jprb,1.0419e-02_jprb,1.9252e-02_jprb,3.2952e-02_jprb,5.2981e-02_jprb/)
kao( :,10,12) = (/&
&1.6525e-02_jprb,3.6418e-02_jprb,7.1731e-02_jprb,1.2850e-01_jprb,2.1301e-01_jprb/)
kao( :,11,12) = (/&
&2.8843e-02_jprb,6.3252e-02_jprb,1.2384e-01_jprb,2.2076e-01_jprb,3.6471e-01_jprb/)
kao( :,12,12) = (/&
&3.6629e-02_jprb,7.9908e-02_jprb,1.5563e-01_jprb,2.7617e-01_jprb,4.5458e-01_jprb/)
kao( :,13,12) = (/&
&3.6750e-02_jprb,7.9425e-02_jprb,1.5346e-01_jprb,2.7137e-01_jprb,4.4636e-01_jprb/)
kao( :, 1,13) = (/&
&6.2842e-03_jprb,8.5901e-03_jprb,1.1770e-02_jprb,1.5848e-02_jprb,2.1231e-02_jprb/)
kao( :, 2,13) = (/&
&5.8832e-03_jprb,8.3961e-03_jprb,1.1688e-02_jprb,1.5821e-02_jprb,2.1081e-02_jprb/)
kao( :, 3,13) = (/&
&5.0693e-03_jprb,7.3487e-03_jprb,1.0502e-02_jprb,1.4654e-02_jprb,1.9671e-02_jprb/)
kao( :, 4,13) = (/&
&4.3064e-03_jprb,6.2164e-03_jprb,8.9638e-03_jprb,1.2692e-02_jprb,1.7434e-02_jprb/)
kao( :, 5,13) = (/&
&3.5860e-03_jprb,5.3245e-03_jprb,7.6791e-03_jprb,1.0860e-02_jprb,1.5079e-02_jprb/)
kao( :, 6,13) = (/&
&2.8314e-03_jprb,4.3212e-03_jprb,6.4657e-03_jprb,9.2062e-03_jprb,1.2862e-02_jprb/)
kao( :, 7,13) = (/&
&2.2377e-03_jprb,3.6343e-03_jprb,5.9331e-03_jprb,9.1864e-03_jprb,1.3701e-02_jprb/)
kao( :, 8,13) = (/&
&2.8263e-03_jprb,5.3372e-03_jprb,9.3113e-03_jprb,1.5184e-02_jprb,2.3750e-02_jprb/)
kao( :, 9,13) = (/&
&7.1632e-03_jprb,1.4853e-02_jprb,2.8075e-02_jprb,4.9370e-02_jprb,8.0768e-02_jprb/)
kao( :,10,13) = (/&
&2.5324e-02_jprb,5.6755e-02_jprb,1.1356e-01_jprb,2.0699e-01_jprb,3.4818e-01_jprb/)
kao( :,11,13) = (/&
&4.7551e-02_jprb,1.0638e-01_jprb,2.1201e-01_jprb,3.8439e-01_jprb,6.4430e-01_jprb/)
kao( :,12,13) = (/&
&6.5536e-02_jprb,1.4596e-01_jprb,2.8971e-01_jprb,5.2186e-01_jprb,8.7097e-01_jprb/)
kao( :,13,13) = (/&
&7.1040e-02_jprb,1.5739e-01_jprb,3.0996e-01_jprb,5.5635e-01_jprb,9.2532e-01_jprb/)
kao( :, 1,14) = (/&
&8.9073e-03_jprb,1.2577e-02_jprb,1.7104e-02_jprb,2.3298e-02_jprb,3.0760e-02_jprb/)
kao( :, 2,14) = (/&
&9.1449e-03_jprb,1.2712e-02_jprb,1.7315e-02_jprb,2.3477e-02_jprb,3.1268e-02_jprb/)
kao( :, 3,14) = (/&
&8.6458e-03_jprb,1.1907e-02_jprb,1.6470e-02_jprb,2.2100e-02_jprb,2.9699e-02_jprb/)
kao( :, 4,14) = (/&
&7.6642e-03_jprb,1.0832e-02_jprb,1.5035e-02_jprb,2.0201e-02_jprb,2.7549e-02_jprb/)
kao( :, 5,14) = (/&
&6.6532e-03_jprb,9.4973e-03_jprb,1.3436e-02_jprb,1.8486e-02_jprb,2.5410e-02_jprb/)
kao( :, 6,14) = (/&
&5.6659e-03_jprb,8.2188e-03_jprb,1.1655e-02_jprb,1.6380e-02_jprb,2.2833e-02_jprb/)
kao( :, 7,14) = (/&
&4.7388e-03_jprb,7.0553e-03_jprb,1.0065e-02_jprb,1.4371e-02_jprb,2.0476e-02_jprb/)
kao( :, 8,14) = (/&
&4.0074e-03_jprb,6.9317e-03_jprb,1.1650e-02_jprb,1.9270e-02_jprb,3.0109e-02_jprb/)
kao( :, 9,14) = (/&
&9.4290e-03_jprb,1.9762e-02_jprb,3.7876e-02_jprb,6.6863e-02_jprb,1.1024e-01_jprb/)
kao( :,10,14) = (/&
&3.5639e-02_jprb,8.0302e-02_jprb,1.6119e-01_jprb,2.9567e-01_jprb,5.0285e-01_jprb/)
kao( :,11,14) = (/&
&7.0994e-02_jprb,1.5991e-01_jprb,3.2149e-01_jprb,5.8889e-01_jprb,1.0009e+00_jprb/)
kao( :,12,14) = (/&
&1.0472e-01_jprb,2.3582e-01_jprb,4.7386e-01_jprb,8.6654e-01_jprb,1.4714e+00_jprb/)
kao( :,13,14) = (/&
&1.2287e-01_jprb,2.7626e-01_jprb,5.5413e-01_jprb,1.0143e+00_jprb,1.7166e+00_jprb/)
kao( :, 1,15) = (/&
&1.0915e-02_jprb,1.5406e-02_jprb,2.1787e-02_jprb,3.0490e-02_jprb,4.2541e-02_jprb/)
kao( :, 2,15) = (/&
&1.1315e-02_jprb,1.6166e-02_jprb,2.2798e-02_jprb,3.1923e-02_jprb,4.3992e-02_jprb/)
kao( :, 3,15) = (/&
&1.0941e-02_jprb,1.5701e-02_jprb,2.1967e-02_jprb,3.0980e-02_jprb,4.2137e-02_jprb/)
kao( :, 4,15) = (/&
&1.0488e-02_jprb,1.4945e-02_jprb,2.1284e-02_jprb,2.9684e-02_jprb,4.1028e-02_jprb/)
kao( :, 5,15) = (/&
&9.8738e-03_jprb,1.4407e-02_jprb,2.0449e-02_jprb,2.8357e-02_jprb,3.9999e-02_jprb/)
kao( :, 6,15) = (/&
&8.9244e-03_jprb,1.3460e-02_jprb,1.9534e-02_jprb,2.7000e-02_jprb,3.8027e-02_jprb/)
kao( :, 7,15) = (/&
&8.0056e-03_jprb,1.2241e-02_jprb,1.8433e-02_jprb,2.6084e-02_jprb,3.5910e-02_jprb/)
kao( :, 8,15) = (/&
&7.2557e-03_jprb,1.0997e-02_jprb,1.6847e-02_jprb,2.4941e-02_jprb,3.5979e-02_jprb/)
kao( :, 9,15) = (/&
&1.1206e-02_jprb,2.3492e-02_jprb,4.5074e-02_jprb,7.9749e-02_jprb,1.3179e-01_jprb/)
kao( :,10,15) = (/&
&4.3714e-02_jprb,9.8800e-02_jprb,1.9893e-01_jprb,3.6488e-01_jprb,6.2367e-01_jprb/)
kao( :,11,15) = (/&
&9.0564e-02_jprb,2.0463e-01_jprb,4.1227e-01_jprb,7.5790e-01_jprb,1.2930e+00_jprb/)
kao( :,12,15) = (/&
&1.4027e-01_jprb,3.1640e-01_jprb,6.3737e-01_jprb,1.1714e+00_jprb,1.9962e+00_jprb/)
kao( :,13,15) = (/&
&1.7360e-01_jprb,3.9198e-01_jprb,7.8920e-01_jprb,1.4509e+00_jprb,2.4685e+00_jprb/)
kao( :, 1,16) = (/&
&1.1824e-02_jprb,1.7312e-02_jprb,2.6044e-02_jprb,4.2271e-02_jprb,6.3915e-02_jprb/)
kao( :, 2,16) = (/&
&1.2333e-02_jprb,1.8642e-02_jprb,2.7212e-02_jprb,4.2155e-02_jprb,7.0740e-02_jprb/)
kao( :, 3,16) = (/&
&1.1914e-02_jprb,1.8039e-02_jprb,2.6483e-02_jprb,3.8087e-02_jprb,6.0410e-02_jprb/)
kao( :, 4,16) = (/&
&1.1617e-02_jprb,1.6736e-02_jprb,2.5551e-02_jprb,3.7243e-02_jprb,5.5742e-02_jprb/)
kao( :, 5,16) = (/&
&1.2207e-02_jprb,1.6058e-02_jprb,2.3899e-02_jprb,3.5249e-02_jprb,5.2521e-02_jprb/)
kao( :, 6,16) = (/&
&1.2552e-02_jprb,1.6221e-02_jprb,2.2453e-02_jprb,3.3451e-02_jprb,4.8959e-02_jprb/)
kao( :, 7,16) = (/&
&1.2682e-02_jprb,1.6583e-02_jprb,2.1720e-02_jprb,3.1647e-02_jprb,4.6528e-02_jprb/)
kao( :, 8,16) = (/&
&1.2800e-02_jprb,1.6916e-02_jprb,2.2044e-02_jprb,3.0239e-02_jprb,4.5143e-02_jprb/)
kao( :, 9,16) = (/&
&1.2606e-02_jprb,2.5371e-02_jprb,4.8585e-02_jprb,8.6042e-02_jprb,1.4214e-01_jprb/)
kao( :,10,16) = (/&
&4.7940e-02_jprb,1.0811e-01_jprb,2.1696e-01_jprb,3.9856e-01_jprb,6.7922e-01_jprb/)
kao( :,11,16) = (/&
&1.0038e-01_jprb,2.2692e-01_jprb,4.5634e-01_jprb,8.3917e-01_jprb,1.4277e+00_jprb/)
kao( :,12,16) = (/&
&1.5780e-01_jprb,3.5540e-01_jprb,7.1875e-01_jprb,1.3169e+00_jprb,2.2414e+00_jprb/)
kao( :,13,16) = (/&
&1.9907e-01_jprb,4.4847e-01_jprb,9.0510e-01_jprb,1.6605e+00_jprb,2.8238e+00_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.71391e-02_jprb, 6.43175e-02_jprb, 5.36271e-02_jprb, 4.47136e-02_jprb, 3.72816e-02_jprb,&
&3.10849e-02_jprb, 2.59182e-02_jprb, 2.16102e-02_jprb, 1.80183e-02_jprb, 1.50234e-02_jprb/)
selfrefo( :, 2) = (/&
&8.39275e-02_jprb, 7.04019e-02_jprb, 5.90561e-02_jprb, 4.95388e-02_jprb, 4.15552e-02_jprb,&
&3.48583e-02_jprb, 2.92406e-02_jprb, 2.45283e-02_jprb, 2.05754e-02_jprb, 1.72595e-02_jprb/)
selfrefo( :, 3) = (/&
&8.77747e-02_jprb, 7.41411e-02_jprb, 6.26251e-02_jprb, 5.28978e-02_jprb, 4.46814e-02_jprb,&
&3.77413e-02_jprb, 3.18791e-02_jprb, 2.69275e-02_jprb, 2.27449e-02_jprb, 1.92121e-02_jprb/)
selfrefo( :, 4) = (/&
&9.35855e-02_jprb, 7.94116e-02_jprb, 6.73843e-02_jprb, 5.71787e-02_jprb, 4.85187e-02_jprb,&
&4.11704e-02_jprb, 3.49349e-02_jprb, 2.96439e-02_jprb, 2.51542e-02_jprb, 2.13445e-02_jprb/)
selfrefo( :, 5) = (/&
&9.31058e-02_jprb, 7.97629e-02_jprb, 6.83321e-02_jprb, 5.85395e-02_jprb, 5.01503e-02_jprb,&
&4.29633e-02_jprb, 3.68062e-02_jprb, 3.15316e-02_jprb, 2.70128e-02_jprb, 2.31416e-02_jprb/)
selfrefo( :, 6) = (/&
&7.92454e-02_jprb, 7.03070e-02_jprb, 6.23769e-02_jprb, 5.53412e-02_jprb, 4.90991e-02_jprb,&
&4.35611e-02_jprb, 3.86477e-02_jprb, 3.42885e-02_jprb, 3.04210e-02_jprb, 2.69898e-02_jprb/)
selfrefo( :, 7) = (/&
&7.29322e-02_jprb, 6.57740e-02_jprb, 5.93183e-02_jprb, 5.34963e-02_jprb, 4.82457e-02_jprb,&
&4.35104e-02_jprb, 3.92399e-02_jprb, 3.53886e-02_jprb, 3.19152e-02_jprb, 2.87828e-02_jprb/)
selfrefo( :, 8) = (/&
&7.53680e-02_jprb, 6.63817e-02_jprb, 5.84669e-02_jprb, 5.14958e-02_jprb, 4.53558e-02_jprb,&
&3.99480e-02_jprb, 3.51849e-02_jprb, 3.09897e-02_jprb, 2.72948e-02_jprb, 2.40404e-02_jprb/)
selfrefo( :, 9) = (/&
&7.17087e-02_jprb, 6.34545e-02_jprb, 5.61505e-02_jprb, 4.96871e-02_jprb, 4.39678e-02_jprb,&
&3.89068e-02_jprb, 3.44283e-02_jprb, 3.04654e-02_jprb, 2.69586e-02_jprb, 2.38554e-02_jprb/)
selfrefo( :,10) = (/&
&7.20376e-02_jprb, 6.42071e-02_jprb, 5.72277e-02_jprb, 5.10070e-02_jprb, 4.54625e-02_jprb,&
&4.05207e-02_jprb, 3.61161e-02_jprb, 3.21903e-02_jprb, 2.86912e-02_jprb, 2.55724e-02_jprb/)
selfrefo( :,11) = (/&
&8.56200e-02_jprb, 7.44248e-02_jprb, 6.46934e-02_jprb, 5.62344e-02_jprb, 4.88815e-02_jprb,&
&4.24900e-02_jprb, 3.69342e-02_jprb, 3.21049e-02_jprb, 2.79070e-02_jprb, 2.42581e-02_jprb/)
selfrefo( :,12) = (/&
&9.30364e-02_jprb, 8.00633e-02_jprb, 6.88992e-02_jprb, 5.92918e-02_jprb, 5.10240e-02_jprb,&
&4.39092e-02_jprb, 3.77864e-02_jprb, 3.25174e-02_jprb, 2.79832e-02_jprb, 2.40811e-02_jprb/)
selfrefo( :,13) = (/&
&9.58838e-02_jprb, 8.17150e-02_jprb, 6.96399e-02_jprb, 5.93492e-02_jprb, 5.05792e-02_jprb,&
&4.31051e-02_jprb, 3.67354e-02_jprb, 3.13070e-02_jprb, 2.66808e-02_jprb, 2.27381e-02_jprb/)
selfrefo( :,14) = (/&
&8.92791e-02_jprb, 7.69435e-02_jprb, 6.63123e-02_jprb, 5.71500e-02_jprb, 4.92536e-02_jprb,&
&4.24483e-02_jprb, 3.65833e-02_jprb, 3.15286e-02_jprb, 2.71723e-02_jprb, 2.34179e-02_jprb/)
selfrefo( :,15) = (/&
&9.09947e-02_jprb, 7.81455e-02_jprb, 6.71107e-02_jprb, 5.76340e-02_jprb, 4.94956e-02_jprb,&
&4.25064e-02_jprb, 3.65041e-02_jprb, 3.13494e-02_jprb, 2.69226e-02_jprb, 2.31209e-02_jprb/)
selfrefo( :,16) = (/&
&7.99990e-02_jprb, 7.02289e-02_jprb, 6.16519e-02_jprb, 5.41224e-02_jprb, 4.75125e-02_jprb,&
&4.17099e-02_jprb, 3.66159e-02_jprb, 3.21440e-02_jprb, 2.82183e-02_jprb, 2.47720e-02_jprb/)

!     -----------------------------------------------------------------
RETURN
END SUBROUTINE rrtm_kgb6