rrtm_kgb4_00.F90

Go to the documentation of this file.

!OCL SCALAR
SUBROUTINE rrtm_kgb4_00

!     Originally by Eli J. Mlawer, Atmospheric & Environmental Research.
!     BAND 4:  630-700 cm-1 (low - H2O,CO2; high - O3,CO2)
!     Reformatted for F90 by JJMorcrette, ECMWF
!     Reformatted for NEC by H.Gallée   , LGGE  (splitting)

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

#include "tsmbkind.h"

USE yoerrto4 , ONLY : kao     ,kbo     ,selfrefo   ,fracrefao  ,fracrefbo
USE yoerrta4 , ONLY : strrat1   ,strrat2

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


IMPLICIT NONE
strrat1 = 850.577_jprb
strrat2 = 35.7416_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) = (/&
&2.62628e-01_jprb, 2.29008e-01_jprb, 1.99692e-01_jprb, 1.74129e-01_jprb, 1.51838e-01_jprb,&
&1.32400e-01_jprb, 1.15451e-01_jprb, 1.00672e-01_jprb, 8.77845e-02_jprb, 7.65469e-02_jprb/)
selfrefo( :, 2) = (/&
&2.45051e-01_jprb, 2.12961e-01_jprb, 1.85073e-01_jprb, 1.60837e-01_jprb, 1.39775e-01_jprb,&
&1.21471e-01_jprb, 1.05564e-01_jprb, 9.17397e-02_jprb, 7.97260e-02_jprb, 6.92856e-02_jprb/)
selfrefo( :, 3) = (/&
&2.42194e-01_jprb, 2.09976e-01_jprb, 1.82044e-01_jprb, 1.57827e-01_jprb, 1.36832e-01_jprb,&
&1.18630e-01_jprb, 1.02849e-01_jprb, 8.91673e-02_jprb, 7.73057e-02_jprb, 6.70221e-02_jprb/)
selfrefo( :, 4) = (/&
&2.44485e-01_jprb, 2.11926e-01_jprb, 1.83702e-01_jprb, 1.59237e-01_jprb, 1.38030e-01_jprb,&
&1.19648e-01_jprb, 1.03714e-01_jprb, 8.99014e-02_jprb, 7.79286e-02_jprb, 6.75503e-02_jprb/)
selfrefo( :, 5) = (/&
&2.43120e-01_jprb, 2.10743e-01_jprb, 1.82679e-01_jprb, 1.58351e-01_jprb, 1.37263e-01_jprb,&
&1.18984e-01_jprb, 1.03139e-01_jprb, 8.94038e-02_jprb, 7.74978e-02_jprb, 6.71774e-02_jprb/)
selfrefo( :, 6) = (/&
&2.40558e-01_jprb, 2.08922e-01_jprb, 1.81446e-01_jprb, 1.57583e-01_jprb, 1.36859e-01_jprb,&
&1.18860e-01_jprb, 1.03229e-01_jprb, 8.96529e-02_jprb, 7.78624e-02_jprb, 6.76225e-02_jprb/)
selfrefo( :, 7) = (/&
&2.42496e-01_jprb, 2.10386e-01_jprb, 1.82528e-01_jprb, 1.58359e-01_jprb, 1.37390e-01_jprb,&
&1.19198e-01_jprb, 1.03415e-01_jprb, 8.97211e-02_jprb, 7.78409e-02_jprb, 6.75337e-02_jprb/)
selfrefo( :, 8) = (/&
&2.39781e-01_jprb, 2.08227e-01_jprb, 1.80825e-01_jprb, 1.57029e-01_jprb, 1.36365e-01_jprb,&
&1.18419e-01_jprb, 1.02836e-01_jprb, 8.93030e-02_jprb, 7.75510e-02_jprb, 6.73456e-02_jprb/)
selfrefo( :, 9) = (/&
&2.38707e-01_jprb, 2.07058e-01_jprb, 1.79605e-01_jprb, 1.55792e-01_jprb, 1.35136e-01_jprb,&
&1.17219e-01_jprb, 1.01677e-01_jprb, 8.81962e-02_jprb, 7.65026e-02_jprb, 6.63594e-02_jprb/)
selfrefo( :,10) = (/&
&2.29942e-01_jprb, 2.00668e-01_jprb, 1.75121e-01_jprb, 1.52826e-01_jprb, 1.33370e-01_jprb,&
&1.16390e-01_jprb, 1.01572e-01_jprb, 8.86410e-02_jprb, 7.73560e-02_jprb, 6.75077e-02_jprb/)
selfrefo( :,11) = (/&
&2.39870e-01_jprb, 2.08120e-01_jprb, 1.80573e-01_jprb, 1.56671e-01_jprb, 1.35934e-01_jprb,&
&1.17941e-01_jprb, 1.02330e-01_jprb, 8.87854e-02_jprb, 7.70335e-02_jprb, 6.68371e-02_jprb/)
selfrefo( :,12) = (/&
&2.40196e-01_jprb, 2.08400e-01_jprb, 1.80812e-01_jprb, 1.56877e-01_jprb, 1.36110e-01_jprb,&
&1.18092e-01_jprb, 1.02460e-01_jprb, 8.88962e-02_jprb, 7.71284e-02_jprb, 6.69184e-02_jprb/)
selfrefo( :,13) = (/&
&2.40426e-01_jprb, 2.08603e-01_jprb, 1.80991e-01_jprb, 1.57035e-01_jprb, 1.36249e-01_jprb,&
&1.18214e-01_jprb, 1.02567e-01_jprb, 8.89909e-02_jprb, 7.72117e-02_jprb, 6.69917e-02_jprb/)
selfrefo( :,14) = (/&
&2.40590e-01_jprb, 2.08742e-01_jprb, 1.81110e-01_jprb, 1.57135e-01_jprb, 1.36334e-01_jprb,&
&1.18287e-01_jprb, 1.02628e-01_jprb, 8.90428e-02_jprb, 7.72556e-02_jprb, 6.70288e-02_jprb/)
selfrefo( :,15) = (/&
&2.40634e-01_jprb, 2.08779e-01_jprb, 1.81141e-01_jprb, 1.57162e-01_jprb, 1.36357e-01_jprb,&
&1.18306e-01_jprb, 1.02645e-01_jprb, 8.90565e-02_jprb, 7.72673e-02_jprb, 6.70387e-02_jprb/)
selfrefo( :,16) = (/&
&2.40652e-01_jprb, 2.08793e-01_jprb, 1.81151e-01_jprb, 1.57169e-01_jprb, 1.36362e-01_jprb,&
&1.18309e-01_jprb, 1.02647e-01_jprb, 8.90576e-02_jprb, 7.72675e-02_jprb, 6.70383e-02_jprb/)

fracrefao( :, 1) = (/&
!     From P = 
    &0.15579100_jprb,0.14918099_jprb,0.14113800_jprb,0.13127001_jprb,&
    &0.11796300_jprb,0.10174300_jprb,0.08282370_jprb,0.06238150_jprb,&
    &0.04213440_jprb,0.00458968_jprb,0.00377949_jprb,0.00298736_jprb,&
    &0.00220743_jprb,0.00140644_jprb,0.00053024_jprb,0.00007459_jprb/)
fracrefao( :, 2) = (/&
    &0.15292799_jprb,0.15004000_jprb,0.14211500_jprb,0.13176700_jprb,&
    &0.11821100_jprb,0.10186300_jprb,0.08288040_jprb,0.06241390_jprb,&
    &0.04220720_jprb,0.00459006_jprb,0.00377919_jprb,0.00298743_jprb,&
    &0.00220743_jprb,0.00140644_jprb,0.00053024_jprb,0.00007459_jprb/)
fracrefao( :, 3) = (/&
    &0.14386199_jprb,0.15125300_jprb,0.14650001_jprb,0.13377000_jprb,&
    &0.11895900_jprb,0.10229400_jprb,0.08312110_jprb,0.06239520_jprb,&
    &0.04225560_jprb,0.00459428_jprb,0.00378865_jprb,0.00298860_jprb,&
    &0.00220743_jprb,0.00140644_jprb,0.00053024_jprb,0.00007459_jprb/)
fracrefao( :, 4) = (/&
    &0.14359100_jprb,0.14561599_jprb,0.14479300_jprb,0.13740200_jprb,&
    &0.12150100_jprb,0.10315400_jprb,0.08355480_jprb,0.06247240_jprb,&
    &0.04230980_jprb,0.00459916_jprb,0.00378373_jprb,0.00300063_jprb,&
    &0.00221111_jprb,0.00140644_jprb,0.00053024_jprb,0.00007459_jprb/)
fracrefao( :, 5) = (/&
    &0.14337599_jprb,0.14451601_jprb,0.14238000_jprb,0.13520500_jprb,&
    &0.12354200_jprb,0.10581200_jprb,0.08451810_jprb,0.06262440_jprb,&
    &0.04239590_jprb,0.00460297_jprb,0.00378701_jprb,0.00300466_jprb,&
    &0.00221899_jprb,0.00141020_jprb,0.00053024_jprb,0.00007459_jprb/)
fracrefao( :, 6) = (/&
    &0.14322001_jprb,0.14397401_jprb,0.14117201_jprb,0.13401900_jprb,&
    &0.12255500_jprb,0.10774100_jprb,0.08617650_jprb,0.06296420_jprb,&
    &0.04249590_jprb,0.00463406_jprb,0.00378241_jprb,0.00302037_jprb,&
    &0.00221583_jprb,0.00141103_jprb,0.00053814_jprb,0.00007991_jprb/)
fracrefao( :, 7) = (/&
    &0.14309500_jprb,0.14364301_jprb,0.14043900_jprb,0.13348100_jprb,&
    &0.12211600_jprb,0.10684700_jprb,0.08820590_jprb,0.06374610_jprb,&
    &0.04264730_jprb,0.00464231_jprb,0.00384022_jprb,0.00303427_jprb,&
    &0.00221825_jprb,0.00140943_jprb,0.00055564_jprb,0.00007991_jprb/)
fracrefao( :, 8) = (/&
    &0.15579100_jprb,0.14918099_jprb,0.14113800_jprb,0.13127001_jprb,&
    &0.11796300_jprb,0.10174300_jprb,0.08282370_jprb,0.06238150_jprb,&
    &0.04213440_jprb,0.00458968_jprb,0.00377949_jprb,0.00298736_jprb,&
    &0.00220743_jprb,0.00140644_jprb,0.00053024_jprb,0.00007459_jprb/)
fracrefao( :, 9) = (/&
    &0.15937001_jprb,0.15159500_jprb,0.14242800_jprb,0.13078900_jprb,&
    &0.11671300_jprb,0.10035700_jprb,0.08143450_jprb,0.06093850_jprb,&
    &0.04105320_jprb,0.00446233_jprb,0.00369844_jprb,0.00293784_jprb,&
    &0.00216425_jprb,0.00143403_jprb,0.00054571_jprb,0.00007991_jprb/)

fracrefbo( :, 1) = (/&
!     From P = 1.17 mb.
    &0.15558299_jprb,0.14930600_jprb,0.14104301_jprb,0.13124099_jprb,&
    &0.11792900_jprb,0.10159200_jprb,0.08314130_jprb,0.06240450_jprb,&
    &0.04217020_jprb,0.00459313_jprb,0.00379798_jprb,0.00299835_jprb,&
    &0.00218950_jprb,0.00140615_jprb,0.00053010_jprb,0.00007457_jprb/)
fracrefbo( :, 2) = (/&
    &0.15592700_jprb,0.14918999_jprb,0.14095700_jprb,0.13115700_jprb,&
    &0.11788900_jprb,0.10158000_jprb,0.08313780_jprb,0.06240240_jprb,&
    &0.04217000_jprb,0.00459313_jprb,0.00379798_jprb,0.00299835_jprb,&
    &0.00218950_jprb,0.00140615_jprb,0.00053010_jprb,0.00007457_jprb/)
fracrefbo( :, 3) = (/&
    &0.15949000_jprb,0.15014900_jprb,0.14162201_jprb,0.13080800_jprb,&
    &0.11713500_jprb,0.10057100_jprb,0.08170080_jprb,0.06128110_jprb,&
    &0.04165600_jprb,0.00459202_jprb,0.00379835_jprb,0.00299717_jprb,&
    &0.00218958_jprb,0.00140616_jprb,0.00053010_jprb,0.00007457_jprb/)
fracrefbo( :, 4) = (/&
    &0.15967900_jprb,0.15038200_jprb,0.14196999_jprb,0.13074800_jprb,&
    &0.11701700_jprb,0.10053000_jprb,0.08160790_jprb,0.06122690_jprb,&
    &0.04128310_jprb,0.00456598_jprb,0.00379486_jprb,0.00299457_jprb,&
    &0.00219016_jprb,0.00140619_jprb,0.00053011_jprb,0.00007456_jprb/)
fracrefbo( :, 5) = (/&
    &0.15989800_jprb,0.15057300_jprb,0.14207700_jprb,0.13068600_jprb,&
    &0.11682900_jprb,0.10053900_jprb,0.08163610_jprb,0.06121870_jprb,&
    &0.04121690_jprb,0.00449061_jprb,0.00371235_jprb,0.00294207_jprb,&
    &0.00217778_jprb,0.00139877_jprb,0.00053011_jprb,0.00007455_jprb/)
fracrefbo( :, 6) = (/&
    &0.15950100_jprb,0.15112500_jprb,0.14199100_jprb,0.13071300_jprb,&
    &0.11680800_jprb,0.10054600_jprb,0.08179050_jprb,0.06120910_jprb,&
    &0.04126050_jprb,0.00444324_jprb,0.00366843_jprb,0.00289369_jprb,&
    &0.00211550_jprb,0.00134746_jprb,0.00050874_jprb,0.00007863_jprb/)

!     -----------------------------------------------------------------
RETURN
END SUBROUTINE rrtm_kgb4_00