rrtm_kgb15_00.F90

Go to the documentation of this file.

!OCL SCALAR
SUBROUTINE rrtm_kgb15_00

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

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

#include "tsmbkind.h"

USE yoerrto15, ONLY : kao     ,selfrefo   ,fracrefao
USE yoerrta15, ONLY : strrat

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


IMPLICIT NONE
strrat = 0.2883201_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.72146e-03_jprb, 1.41213e-03_jprb, 1.15839e-03_jprb, 9.50236e-04_jprb, 7.79489e-04_jprb,&
&6.39423e-04_jprb, 5.24525e-04_jprb, 4.30274e-04_jprb, 3.52958e-04_jprb, 2.89535e-04_jprb/)
selfrefo( :, 2) = (/&
&1.83092e-03_jprb, 1.50082e-03_jprb, 1.23023e-03_jprb, 1.00843e-03_jprb, 8.26611e-04_jprb,&
&6.77578e-04_jprb, 5.55414e-04_jprb, 4.55276e-04_jprb, 3.73192e-04_jprb, 3.05908e-04_jprb/)
selfrefo( :, 3) = (/&
&1.92620e-03_jprb, 1.57871e-03_jprb, 1.29391e-03_jprb, 1.06048e-03_jprb, 8.69169e-04_jprb,&
&7.12370e-04_jprb, 5.83857e-04_jprb, 4.78528e-04_jprb, 3.92201e-04_jprb, 3.21447e-04_jprb/)
selfrefo( :, 4) = (/&
&2.11224e-03_jprb, 1.72966e-03_jprb, 1.41638e-03_jprb, 1.15984e-03_jprb, 9.49760e-04_jprb,&
&7.77734e-04_jprb, 6.36867e-04_jprb, 5.21514e-04_jprb, 4.27055e-04_jprb, 3.49704e-04_jprb/)
selfrefo( :, 5) = (/&
&2.31367e-03_jprb, 1.88988e-03_jprb, 1.54371e-03_jprb, 1.26095e-03_jprb, 1.02999e-03_jprb,&
&8.41325e-04_jprb, 6.87221e-04_jprb, 5.61344e-04_jprb, 4.58523e-04_jprb, 3.74536e-04_jprb/)
selfrefo( :, 6) = (/&
&2.47614e-03_jprb, 2.03107e-03_jprb, 1.66599e-03_jprb, 1.36654e-03_jprb, 1.12091e-03_jprb,&
&9.19433e-04_jprb, 7.54170e-04_jprb, 6.18612e-04_jprb, 5.07419e-04_jprb, 4.16213e-04_jprb/)
selfrefo( :, 7) = (/&
&2.55137e-03_jprb, 2.10353e-03_jprb, 1.73430e-03_jprb, 1.42987e-03_jprb, 1.17889e-03_jprb,&
&9.71958e-04_jprb, 8.01350e-04_jprb, 6.60688e-04_jprb, 5.44717e-04_jprb, 4.49103e-04_jprb/)
selfrefo( :, 8) = (/&
&2.59317e-03_jprb, 2.13957e-03_jprb, 1.76532e-03_jprb, 1.45653e-03_jprb, 1.20175e-03_jprb,&
&9.91539e-04_jprb, 8.18098e-04_jprb, 6.74996e-04_jprb, 5.56926e-04_jprb, 4.59508e-04_jprb/)
selfrefo( :, 9) = (/&
&2.83495e-03_jprb, 2.34704e-03_jprb, 1.94309e-03_jprb, 1.60867e-03_jprb, 1.33181e-03_jprb,&
&1.10259e-03_jprb, 9.12826e-04_jprb, 7.55721e-04_jprb, 6.25656e-04_jprb, 5.17975e-04_jprb/)
selfrefo( :,10) = (/&
&2.98196e-03_jprb, 2.46746e-03_jprb, 2.04173e-03_jprb, 1.68946e-03_jprb, 1.39796e-03_jprb,&
&1.15676e-03_jprb, 9.57177e-04_jprb, 7.92028e-04_jprb, 6.55373e-04_jprb, 5.42297e-04_jprb/)
selfrefo( :,11) = (/&
&2.98241e-03_jprb, 2.46776e-03_jprb, 2.04191e-03_jprb, 1.68955e-03_jprb, 1.39800e-03_jprb,&
&1.15675e-03_jprb, 9.57141e-04_jprb, 7.91974e-04_jprb, 6.55308e-04_jprb, 5.42226e-04_jprb/)
selfrefo( :,12) = (/&
&2.97837e-03_jprb, 2.46464e-03_jprb, 2.03952e-03_jprb, 1.68773e-03_jprb, 1.39662e-03_jprb,&
&1.15572e-03_jprb, 9.56375e-04_jprb, 7.91413e-04_jprb, 6.54904e-04_jprb, 5.41942e-04_jprb/)
selfrefo( :,13) = (/&
&2.97777e-03_jprb, 2.46461e-03_jprb, 2.03989e-03_jprb, 1.68835e-03_jprb, 1.39740e-03_jprb,&
&1.15659e-03_jprb, 9.57271e-04_jprb, 7.92305e-04_jprb, 6.55767e-04_jprb, 5.42759e-04_jprb/)
selfrefo( :,14) = (/&
&2.98331e-03_jprb, 2.46947e-03_jprb, 2.04413e-03_jprb, 1.69205e-03_jprb, 1.40061e-03_jprb,&
&1.15937e-03_jprb, 9.59683e-04_jprb, 7.94388e-04_jprb, 6.57564e-04_jprb, 5.44306e-04_jprb/)
selfrefo( :,15) = (/&
&2.99414e-03_jprb, 2.47830e-03_jprb, 2.05133e-03_jprb, 1.69792e-03_jprb, 1.40540e-03_jprb,&
&1.16328e-03_jprb, 9.62863e-04_jprb, 7.96978e-04_jprb, 6.59673e-04_jprb, 5.46023e-04_jprb/)
selfrefo( :,16) = (/&
&3.00006e-03_jprb, 2.48297e-03_jprb, 2.05501e-03_jprb, 1.70081e-03_jprb, 1.40766e-03_jprb,&
&1.16503e-03_jprb, 9.64227e-04_jprb, 7.98033e-04_jprb, 6.60484e-04_jprb, 5.46643e-04_jprb/)

fracrefao( :, 1) = (/&
!     From P = 1053.6 mb.
    &0.11287100_jprb,0.12070200_jprb,0.12729000_jprb,0.12858100_jprb,&
    &0.12743001_jprb,0.11961800_jprb,0.10290400_jprb,0.07888980_jprb,&
    &0.05900120_jprb,0.00667979_jprb,0.00552926_jprb,0.00436993_jprb,&
    &0.00320611_jprb,0.00204765_jprb,0.00077371_jprb,0.00010894_jprb/)
fracrefao( :, 2) = (/&
    &0.13918801_jprb,0.16353001_jprb,0.16155800_jprb,0.14090499_jprb,&
    &0.11322300_jprb,0.08757720_jprb,0.07225720_jprb,0.05173390_jprb,&
    &0.04731360_jprb,0.00667979_jprb,0.00552926_jprb,0.00436993_jprb,&
    &0.00320611_jprb,0.00204765_jprb,0.00077371_jprb,0.00010894_jprb/)
fracrefao( :, 3) = (/&
    &0.14687300_jprb,0.17853101_jprb,0.15664500_jprb,0.13351700_jprb,&
    &0.10791200_jprb,0.08684320_jprb,0.07158090_jprb,0.05198410_jprb,&
    &0.04340110_jprb,0.00667979_jprb,0.00552926_jprb,0.00436993_jprb,&
    &0.00320611_jprb,0.00204765_jprb,0.00077371_jprb,0.00010894_jprb/)
fracrefao( :, 4) = (/&
    &0.15760700_jprb,0.17759100_jprb,0.15158001_jprb,0.13193300_jprb,&
    &0.10742800_jprb,0.08693760_jprb,0.07159490_jprb,0.05196250_jprb,&
    &0.04065270_jprb,0.00667979_jprb,0.00552926_jprb,0.00436993_jprb,&
    &0.00320611_jprb,0.00204765_jprb,0.00077371_jprb,0.00010894_jprb/)
fracrefao( :, 5) = (/&
    &0.16646700_jprb,0.17299300_jprb,0.15018500_jprb,0.13138700_jprb,&
    &0.10735900_jprb,0.08713110_jprb,0.07130330_jprb,0.05279420_jprb,&
    &0.03766730_jprb,0.00667979_jprb,0.00552926_jprb,0.00436993_jprb,&
    &0.00320611_jprb,0.00204765_jprb,0.00077371_jprb,0.00010894_jprb/)
fracrefao( :, 6) = (/&
    &0.17546000_jprb,0.16666500_jprb,0.14969499_jprb,0.13105400_jprb,&
    &0.10782500_jprb,0.08718610_jprb,0.07156770_jprb,0.05308320_jprb,&
    &0.03753960_jprb,0.00432465_jprb,0.00509623_jprb,0.00436993_jprb,&
    &0.00320611_jprb,0.00204765_jprb,0.00077371_jprb,0.00010894_jprb/)
fracrefao( :, 7) = (/&
    &0.18378501_jprb,0.16064601_jprb,0.14940400_jprb,0.13146400_jprb,&
    &0.10810300_jprb,0.08775740_jprb,0.07115360_jprb,0.05400040_jprb,&
    &0.03689970_jprb,0.00388333_jprb,0.00323610_jprb,0.00353414_jprb,&
    &0.00320611_jprb,0.00204765_jprb,0.00077371_jprb,0.00010894_jprb/)
fracrefao( :, 8) = (/&
    &0.18966800_jprb,0.15744300_jprb,0.14993000_jprb,0.13152599_jprb,&
    &0.10899200_jprb,0.08858690_jprb,0.07142920_jprb,0.05399600_jprb,&
    &0.03433460_jprb,0.00374886_jprb,0.00302066_jprb,0.00240653_jprb,&
    &0.00199205_jprb,0.00204765_jprb,0.00077371_jprb,0.00010894_jprb/)
fracrefao( :, 9) = (/&
    &0.11887100_jprb,0.12479600_jprb,0.12569501_jprb,0.12839900_jprb,&
    &0.12473500_jprb,0.12012800_jprb,0.11086700_jprb,0.08493590_jprb,&
    &0.05063770_jprb,0.00328723_jprb,0.00266849_jprb,0.00210232_jprb,&
    &0.00152114_jprb,0.00095635_jprb,0.00035374_jprb,0.00004980_jprb/)

!     -----------------------------------------------------------------
RETURN
END SUBROUTINE rrtm_kgb15_00