rrtm_kgb13_00.F90

Go to the documentation of this file.

!OCL SCALAR
SUBROUTINE rrtm_kgb13_00

!     Originally by Eli J. Mlawer, Atmospheric & Environmental Research.
!     BAND 13:  2080-2250 cm-1 (low - H2O,N2O; high - nothing)
!     Reformatted for F90 by JJMorcrette, ECMWF
!     Reformatted for NEC by H.Gallée   , LGGE  (splitting)

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

#include "tsmbkind.h"

USE yoerrto13, ONLY : kao     ,selfrefo   ,fracrefao
USE yoerrta13, ONLY : strrat

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


IMPLICIT NONE
strrat = 16658.87_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) = (/&
&9.83745e-03_jprb, 8.45687e-03_jprb, 7.27004e-03_jprb, 6.24977e-03_jprb, 5.37269e-03_jprb,&
&4.61869e-03_jprb, 3.97051e-03_jprb, 3.41329e-03_jprb, 2.93427e-03_jprb, 2.52248e-03_jprb/)
selfrefo( :, 2) = (/&
&9.79014e-03_jprb, 8.48557e-03_jprb, 7.35483e-03_jprb, 6.37477e-03_jprb, 5.52531e-03_jprb,&
&4.78905e-03_jprb, 4.15089e-03_jprb, 3.59777e-03_jprb, 3.11835e-03_jprb, 2.70282e-03_jprb/)
selfrefo( :, 3) = (/&
&9.36931e-03_jprb, 8.10880e-03_jprb, 7.01787e-03_jprb, 6.07371e-03_jprb, 5.25658e-03_jprb,&
&4.54938e-03_jprb, 3.93732e-03_jprb, 3.40761e-03_jprb, 2.94916e-03_jprb, 2.55239e-03_jprb/)
selfrefo( :, 4) = (/&
&8.44896e-03_jprb, 7.26401e-03_jprb, 6.24525e-03_jprb, 5.36937e-03_jprb, 4.61632e-03_jprb,&
&3.96889e-03_jprb, 3.41226e-03_jprb, 2.93370e-03_jprb, 2.52226e-03_jprb, 2.16851e-03_jprb/)
selfrefo( :, 5) = (/&
&6.69018e-03_jprb, 5.93911e-03_jprb, 5.27237e-03_jprb, 4.68047e-03_jprb, 4.15502e-03_jprb,&
&3.68856e-03_jprb, 3.27447e-03_jprb, 2.90686e-03_jprb, 2.58053e-03_jprb, 2.29083e-03_jprb/)
selfrefo( :, 6) = (/&
&6.26590e-03_jprb, 5.54076e-03_jprb, 4.89954e-03_jprb, 4.33252e-03_jprb, 3.83113e-03_jprb,&
&3.38776e-03_jprb, 2.99570e-03_jprb, 2.64901e-03_jprb, 2.34244e-03_jprb, 2.07136e-03_jprb/)
selfrefo( :, 7) = (/&
&6.01943e-03_jprb, 5.31006e-03_jprb, 4.68429e-03_jprb, 4.13226e-03_jprb, 3.64529e-03_jprb,&
&3.21570e-03_jprb, 2.83674e-03_jprb, 2.50244e-03_jprb, 2.20754e-03_jprb, 1.94739e-03_jprb/)
selfrefo( :, 8) = (/&
&5.80205e-03_jprb, 5.07696e-03_jprb, 4.44248e-03_jprb, 3.88730e-03_jprb, 3.40149e-03_jprb,&
&2.97640e-03_jprb, 2.60444e-03_jprb, 2.27895e-03_jprb, 1.99415e-03_jprb, 1.74494e-03_jprb/)
selfrefo( :, 9) = (/&
&4.90962e-03_jprb, 4.54613e-03_jprb, 4.20956e-03_jprb, 3.89790e-03_jprb, 3.60931e-03_jprb,&
&3.34210e-03_jprb, 3.09466e-03_jprb, 2.86555e-03_jprb, 2.65339e-03_jprb, 2.45695e-03_jprb/)
selfrefo( :,10) = (/&
&4.51541e-03_jprb, 4.34841e-03_jprb, 4.18758e-03_jprb, 4.03271e-03_jprb, 3.88356e-03_jprb,&
&3.73992e-03_jprb, 3.60160e-03_jprb, 3.46840e-03_jprb, 3.34012e-03_jprb, 3.21658e-03_jprb/)
selfrefo( :,11) = (/&
&4.41859e-03_jprb, 4.25299e-03_jprb, 4.09360e-03_jprb, 3.94018e-03_jprb, 3.79252e-03_jprb,&
&3.65038e-03_jprb, 3.51357e-03_jprb, 3.38189e-03_jprb, 3.25515e-03_jprb, 3.13316e-03_jprb/)
selfrefo( :,12) = (/&
&4.41867e-03_jprb, 4.22140e-03_jprb, 4.03294e-03_jprb, 3.85289e-03_jprb, 3.68088e-03_jprb,&
&3.51655e-03_jprb, 3.35956e-03_jprb, 3.20957e-03_jprb, 3.06628e-03_jprb, 2.92939e-03_jprb/)
selfrefo( :,13) = (/&
&4.42776e-03_jprb, 4.34604e-03_jprb, 4.26583e-03_jprb, 4.18710e-03_jprb, 4.10983e-03_jprb,&
&4.03398e-03_jprb, 3.95953e-03_jprb, 3.88645e-03_jprb, 3.81472e-03_jprb, 3.74432e-03_jprb/)
selfrefo( :,14) = (/&
&4.46410e-03_jprb, 4.46210e-03_jprb, 4.46011e-03_jprb, 4.45812e-03_jprb, 4.45613e-03_jprb,&
&4.45414e-03_jprb, 4.45216e-03_jprb, 4.45017e-03_jprb, 4.44818e-03_jprb, 4.44620e-03_jprb/)
selfrefo( :,15) = (/&
&3.21270e-03_jprb, 3.42098e-03_jprb, 3.64276e-03_jprb, 3.87892e-03_jprb, 4.13040e-03_jprb,&
&4.39817e-03_jprb, 4.68330e-03_jprb, 4.98692e-03_jprb, 5.31023e-03_jprb, 5.65449e-03_jprb/)
selfrefo( :,16) = (/&
&3.11411e-03_jprb, 3.35874e-03_jprb, 3.62259e-03_jprb, 3.90717e-03_jprb, 4.21410e-03_jprb,&
&4.54514e-03_jprb, 4.90219e-03_jprb, 5.28729e-03_jprb, 5.70264e-03_jprb, 6.15061e-03_jprb/)

fracrefao( : , 1) = (/&
!   From P = 706.3 mb.
    &0.17683899_jprb,0.17319500_jprb,0.15712699_jprb,0.13604601_jprb,&
    &0.10776200_jprb,0.08750010_jprb,0.06808820_jprb,0.04905150_jprb,&
    &0.03280360_jprb,0.00350836_jprb,0.00281864_jprb,0.00219862_jprb,&
    &0.00160943_jprb,0.00101885_jprb,0.00038147_jprb,0.00005348_jprb/)
fracrefao( : , 2) = (/&
    &0.17535400_jprb,0.16999300_jprb,0.15610200_jprb,0.13589200_jprb,&
    &0.10842100_jprb,0.08988550_jprb,0.06943920_jprb,0.04974900_jprb,&
    &0.03323400_jprb,0.00352752_jprb,0.00289402_jprb,0.00231003_jprb,&
    &0.00174659_jprb,0.00101884_jprb,0.00038147_jprb,0.00005348_jprb/)
fracrefao( : , 3) = (/&
    &0.17409500_jprb,0.16846400_jprb,0.15641899_jprb,0.13503000_jprb,&
    &0.10838600_jprb,0.08985800_jprb,0.07092720_jprb,0.05075710_jprb,&
    &0.03364180_jprb,0.00354241_jprb,0.00303507_jprb,0.00243391_jprb,&
    &0.00177502_jprb,0.00114638_jprb,0.00043585_jprb,0.00005348_jprb/)
fracrefao( : , 4) = (/&
    &0.17248300_jprb,0.16778600_jprb,0.15543500_jprb,0.13496999_jprb,&
    &0.10826300_jprb,0.09028740_jprb,0.07156720_jprb,0.05187120_jprb,&
    &0.03424890_jprb,0.00363933_jprb,0.00324715_jprb,0.00255030_jprb,&
    &0.00187380_jprb,0.00116978_jprb,0.00051229_jprb,0.00009768_jprb/)
fracrefao( : , 5) = (/&
    &0.17061099_jprb,0.16715799_jprb,0.15405200_jprb,0.13471501_jprb,&
    &0.10896400_jprb,0.09069460_jprb,0.07229760_jprb,0.05218280_jprb,&
    &0.03555340_jprb,0.00379576_jprb,0.00330240_jprb,0.00274693_jprb,&
    &0.00201587_jprb,0.00119598_jprb,0.00061885_jprb,0.00009768_jprb/)
fracrefao( : , 6) = (/&
    &0.16789700_jprb,0.16629100_jprb,0.15270300_jprb,0.13360199_jprb,&
    &0.11047200_jprb,0.09151080_jprb,0.07325000_jprb,0.05261450_jprb,&
    &0.03657990_jprb,0.00450092_jprb,0.00349537_jprb,0.00283321_jprb,&
    &0.00208396_jprb,0.00140354_jprb,0.00066587_jprb,0.00009768_jprb/)
fracrefao( : , 7) = (/&
    &0.16412200_jprb,0.16387400_jprb,0.15211500_jprb,0.13062200_jprb,&
    &0.11325100_jprb,0.09348130_jprb,0.07381380_jprb,0.05434740_jprb,&
    &0.03803160_jprb,0.00481346_jprb,0.00393592_jprb,0.00296633_jprb,&
    &0.00222532_jprb,0.00163762_jprb,0.00066648_jprb,0.00009768_jprb/)
fracrefao( : , 8) = (/&
    &0.15513401_jprb,0.15768200_jprb,0.14850400_jprb,0.13330200_jprb,&
    &0.11446500_jprb,0.09868230_jprb,0.07642050_jprb,0.05624170_jprb,&
    &0.04197810_jprb,0.00502288_jprb,0.00429452_jprb,0.00315347_jprb,&
    &0.00263559_jprb,0.00171772_jprb,0.00066860_jprb,0.00009768_jprb/)
fracrefao( : , 9) = (/&
    &0.15732600_jprb,0.15223300_jprb,0.14271900_jprb,0.13563600_jprb,&
    &0.11859600_jprb,0.10274200_jprb,0.07934560_jprb,0.05763410_jprb,&
    &0.03921740_jprb,0.00437741_jprb,0.00337921_jprb,0.00280212_jprb,&
    &0.00200156_jprb,0.00124812_jprb,0.00064664_jprb,0.00009768_jprb/)

!     ------------------------------------------------------------------
RETURN
END SUBROUTINE rrtm_kgb13_00