rrtm_kgb16_00.F90

Go to the documentation of this file.

!OCL SCALAR
SUBROUTINE rrtm_kgb16_00

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

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

#include "tsmbkind.h"

USE yoerrto16, ONLY : kao     ,selfrefo   ,fracrefao
USE yoerrta16, ONLY : strrat

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


IMPLICIT NONE
strrat = 830.411_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.26476e-03_jprb, 1.04962e-03_jprb, 8.71080e-04_jprb, 7.22907e-04_jprb, 5.99939e-04_jprb,&
&4.97888e-04_jprb, 4.13197e-04_jprb, 3.42911e-04_jprb, 2.84581e-04_jprb, 2.36173e-04_jprb/)
selfrefo( :, 2) = (/&
&1.27380e-03_jprb, 1.05867e-03_jprb, 8.79870e-04_jprb, 7.31270e-04_jprb, 6.07767e-04_jprb,&
&5.05121e-04_jprb, 4.19812e-04_jprb, 3.48911e-04_jprb, 2.89983e-04_jprb, 2.41008e-04_jprb/)
selfrefo( :, 3) = (/&
&1.26332e-03_jprb, 1.05867e-03_jprb, 8.87172e-04_jprb, 7.43454e-04_jprb, 6.23018e-04_jprb,&
&5.22092e-04_jprb, 4.37515e-04_jprb, 3.66640e-04_jprb, 3.07246e-04_jprb, 2.57474e-04_jprb/)
selfrefo( :, 4) = (/&
&1.28078e-03_jprb, 1.07802e-03_jprb, 9.07363e-04_jprb, 7.63722e-04_jprb, 6.42820e-04_jprb,&
&5.41057e-04_jprb, 4.55404e-04_jprb, 3.83311e-04_jprb, 3.22630e-04_jprb, 2.71556e-04_jprb/)
selfrefo( :, 5) = (/&
&1.61188e-03_jprb, 1.32375e-03_jprb, 1.08712e-03_jprb, 8.92794e-04_jprb, 7.33203e-04_jprb,&
&6.02140e-04_jprb, 4.94505e-04_jprb, 4.06110e-04_jprb, 3.33516e-04_jprb, 2.73898e-04_jprb/)
selfrefo( :, 6) = (/&
&1.42549e-03_jprb, 1.25964e-03_jprb, 1.11309e-03_jprb, 9.83595e-04_jprb, 8.69162e-04_jprb,&
&7.68043e-04_jprb, 6.78688e-04_jprb, 5.99728e-04_jprb, 5.29955e-04_jprb, 4.68299e-04_jprb/)
selfrefo( :, 7) = (/&
&1.43253e-03_jprb, 1.28535e-03_jprb, 1.15330e-03_jprb, 1.03480e-03_jprb, 9.28488e-04_jprb,&
&8.33094e-04_jprb, 7.47501e-04_jprb, 6.70702e-04_jprb, 6.01794e-04_jprb, 5.39965e-04_jprb/)
selfrefo( :, 8) = (/&
&1.48973e-03_jprb, 1.33348e-03_jprb, 1.19363e-03_jprb, 1.06844e-03_jprb, 9.56378e-04_jprb,&
&8.56072e-04_jprb, 7.66286e-04_jprb, 6.85918e-04_jprb, 6.13978e-04_jprb, 5.49583e-04_jprb/)
selfrefo( :, 9) = (/&
&1.53002e-03_jprb, 1.44607e-03_jprb, 1.36674e-03_jprb, 1.29175e-03_jprb, 1.22088e-03_jprb,&
&1.15390e-03_jprb, 1.09059e-03_jprb, 1.03076e-03_jprb, 9.74209e-04_jprb, 9.20760e-04_jprb/)
selfrefo( :,10) = (/&
&1.78480e-03_jprb, 1.68865e-03_jprb, 1.59768e-03_jprb, 1.51161e-03_jprb, 1.43018e-03_jprb,&
&1.35313e-03_jprb, 1.28023e-03_jprb, 1.21126e-03_jprb, 1.14601e-03_jprb, 1.08427e-03_jprb/)
selfrefo( :,11) = (/&
&1.85121e-03_jprb, 1.73903e-03_jprb, 1.63364e-03_jprb, 1.53464e-03_jprb, 1.44164e-03_jprb,&
&1.35427e-03_jprb, 1.27220e-03_jprb, 1.19510e-03_jprb, 1.12268e-03_jprb, 1.05464e-03_jprb/)
selfrefo( :,12) = (/&
&2.10772e-03_jprb, 1.96994e-03_jprb, 1.84118e-03_jprb, 1.72082e-03_jprb, 1.60834e-03_jprb,&
&1.50321e-03_jprb, 1.40495e-03_jprb, 1.31311e-03_jprb, 1.22728e-03_jprb, 1.14705e-03_jprb/)
selfrefo( :,13) = (/&
&1.99323e-03_jprb, 1.89279e-03_jprb, 1.79741e-03_jprb, 1.70684e-03_jprb, 1.62083e-03_jprb,&
&1.53915e-03_jprb, 1.46159e-03_jprb, 1.38794e-03_jprb, 1.31800e-03_jprb, 1.25159e-03_jprb/)
selfrefo( :,14) = (/&
&2.05835e-03_jprb, 1.93301e-03_jprb, 1.81529e-03_jprb, 1.70475e-03_jprb, 1.60093e-03_jprb,&
&1.50344e-03_jprb, 1.41189e-03_jprb, 1.32591e-03_jprb, 1.24516e-03_jprb, 1.16934e-03_jprb/)
selfrefo( :,15) = (/&
&2.34120e-03_jprb, 2.20012e-03_jprb, 2.06754e-03_jprb, 1.94295e-03_jprb, 1.82587e-03_jprb,&
&1.71585e-03_jprb, 1.61245e-03_jprb, 1.51529e-03_jprb, 1.42398e-03_jprb, 1.33817e-03_jprb/)
selfrefo( :,16) = (/&
&3.03559e-03_jprb, 2.80023e-03_jprb, 2.58312e-03_jprb, 2.38284e-03_jprb, 2.19809e-03_jprb,&
&2.02767e-03_jprb, 1.87045e-03_jprb, 1.72543e-03_jprb, 1.59165e-03_jprb, 1.46825e-03_jprb/)

fracrefao( :, 1) = (/&
!     From P = 862.6 mb.
    &0.17356300_jprb,0.18880001_jprb,0.17704099_jprb,0.13661300_jprb,&
    &0.10691600_jprb,0.08222480_jprb,0.05939860_jprb,0.04230810_jprb,&
    &0.02526330_jprb,0.00244532_jprb,0.00193541_jprb,0.00150415_jprb,&
    &0.00103528_jprb,0.00067068_jprb,0.00024951_jprb,0.00003348_jprb/)
fracrefao( :, 2) = (/&
    &0.17779499_jprb,0.19837400_jprb,0.16557600_jprb,0.13470000_jprb,&
    &0.11013600_jprb,0.08342720_jprb,0.05987030_jprb,0.03938700_jprb,&
    &0.02293650_jprb,0.00238849_jprb,0.00192400_jprb,0.00149921_jprb,&
    &0.00103539_jprb,0.00067150_jprb,0.00024822_jprb,0.00003348_jprb/)
fracrefao( :, 3) = (/&
    &0.18535601_jprb,0.19407199_jprb,0.16053200_jprb,0.13300700_jprb,&
    &0.10779000_jprb,0.08408500_jprb,0.06480450_jprb,0.04070160_jprb,&
    &0.02203590_jprb,0.00227779_jprb,0.00189074_jprb,0.00146888_jprb,&
    &0.00103147_jprb,0.00066770_jprb,0.00024751_jprb,0.00003348_jprb/)
fracrefao( :, 4) = (/&
    &0.19139200_jprb,0.18917400_jprb,0.15748601_jprb,0.13240699_jprb,&
    &0.10557300_jprb,0.08383260_jprb,0.06724060_jprb,0.04364450_jprb,&
    &0.02175820_jprb,0.00225436_jprb,0.00184421_jprb,0.00143153_jprb,&
    &0.00103027_jprb,0.00066066_jprb,0.00024222_jprb,0.00003148_jprb/)
fracrefao( :, 5) = (/&
    &0.19547801_jprb,0.18539500_jprb,0.15442000_jprb,0.13114899_jprb,&
    &0.10515600_jprb,0.08350350_jprb,0.06909780_jprb,0.04671630_jprb,&
    &0.02168820_jprb,0.00224400_jprb,0.00182009_jprb,0.00139098_jprb,&
    &0.00102582_jprb,0.00065367_jprb,0.00023202_jprb,0.00003148_jprb/)
fracrefao( :, 6) = (/&
    &0.19757500_jprb,0.18266800_jprb,0.15208900_jprb,0.12897800_jprb,&
    &0.10637200_jprb,0.08391220_jprb,0.06989830_jprb,0.04964120_jprb,&
    &0.02155800_jprb,0.00224310_jprb,0.00177358_jprb,0.00138184_jprb,&
    &0.00101538_jprb,0.00063370_jprb,0.00023227_jprb,0.00003148_jprb/)
fracrefao( :, 7) = (/&
    &0.20145500_jprb,0.17692900_jprb,0.14940600_jprb,0.12690400_jprb,&
    &0.10828800_jprb,0.08553720_jprb,0.07004940_jprb,0.05153430_jprb,&
    &0.02268740_jprb,0.00216943_jprb,0.00178603_jprb,0.00137754_jprb,&
    &0.00098344_jprb,0.00063165_jprb,0.00023218_jprb,0.00003148_jprb/)
fracrefao( :, 8) = (/&
    &0.20383500_jprb,0.17047501_jprb,0.14570600_jprb,0.12679300_jprb,&
    &0.11043100_jprb,0.08719150_jprb,0.07045440_jprb,0.05345420_jprb,&
    &0.02448340_jprb,0.00215839_jprb,0.00175893_jprb,0.00138296_jprb,&
    &0.00098318_jprb,0.00063188_jprb,0.00023199_jprb,0.00003148_jprb/)
fracrefao( :, 9) = (/&
    &0.18680701_jprb,0.15961801_jprb,0.15092900_jprb,0.13049100_jprb,&
    &0.11418400_jprb,0.09380540_jprb,0.07093450_jprb,0.05664280_jprb,&
    &0.02938410_jprb,0.00217751_jprb,0.00176766_jprb,0.00138275_jprb,&
    &0.00098377_jprb,0.00063181_jprb,0.00023193_jprb,0.00003148_jprb/)

!     -----------------------------------------------------------------
RETURN
END SUBROUTINE rrtm_kgb16_00