rrtm_kgb12_00.F90

Go to the documentation of this file.

!OCL SCALAR
SUBROUTINE rrtm_kgb12_00

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

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

#include "tsmbkind.h"

USE yoerrto12, ONLY : kao     ,selfrefo   ,fracrefao
USE yoerrta12, ONLY : strrat

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


IMPLICIT NONE
strrat = 0.009736757_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.34973e-02_jprb, 2.07877e-02_jprb, 1.83905e-02_jprb, 1.62698e-02_jprb, 1.43936e-02_jprb,&
&1.27338e-02_jprb, 1.12654e-02_jprb, 9.96632e-03_jprb, 8.81704e-03_jprb, 7.80030e-03_jprb/)
selfrefo( :, 2) = (/&
&3.10090e-02_jprb, 2.82626e-02_jprb, 2.57594e-02_jprb, 2.34779e-02_jprb, 2.13985e-02_jprb,&
&1.95033e-02_jprb, 1.77759e-02_jprb, 1.62015e-02_jprb, 1.47665e-02_jprb, 1.34587e-02_jprb/)
selfrefo( :, 3) = (/&
&5.36676e-02_jprb, 4.95621e-02_jprb, 4.57706e-02_jprb, 4.22692e-02_jprb, 3.90356e-02_jprb,&
&3.60494e-02_jprb, 3.32917e-02_jprb, 3.07449e-02_jprb, 2.83929e-02_jprb, 2.62209e-02_jprb/)
selfrefo( :, 4) = (/&
&9.05844e-02_jprb, 8.33349e-02_jprb, 7.66656e-02_jprb, 7.05300e-02_jprb, 6.48855e-02_jprb,&
&5.96927e-02_jprb, 5.49155e-02_jprb, 5.05206e-02_jprb, 4.64774e-02_jprb, 4.27578e-02_jprb/)
selfrefo( :, 5) = (/&
&1.10531e-01_jprb, 1.03430e-01_jprb, 9.67851e-02_jprb, 9.05670e-02_jprb, 8.47484e-02_jprb,&
&7.93037e-02_jprb, 7.42087e-02_jprb, 6.94411e-02_jprb, 6.49797e-02_jprb, 6.08050e-02_jprb/)
selfrefo( :, 6) = (/&
&1.24308e-01_jprb, 1.16014e-01_jprb, 1.08274e-01_jprb, 1.01050e-01_jprb, 9.43078e-02_jprb,&
&8.80155e-02_jprb, 8.21431e-02_jprb, 7.66624e-02_jprb, 7.15475e-02_jprb, 6.67738e-02_jprb/)
selfrefo( :, 7) = (/&
&1.34594e-01_jprb, 1.25054e-01_jprb, 1.16190e-01_jprb, 1.07954e-01_jprb, 1.00302e-01_jprb,&
&9.31920e-02_jprb, 8.65863e-02_jprb, 8.04488e-02_jprb, 7.47464e-02_jprb, 6.94481e-02_jprb/)
selfrefo( :, 8) = (/&
&1.26966e-01_jprb, 1.18270e-01_jprb, 1.10169e-01_jprb, 1.02624e-01_jprb, 9.55948e-02_jprb,&
&8.90474e-02_jprb, 8.29484e-02_jprb, 7.72672e-02_jprb, 7.19750e-02_jprb, 6.70453e-02_jprb/)
selfrefo( :, 9) = (/&
&1.29806e-01_jprb, 1.21296e-01_jprb, 1.13344e-01_jprb, 1.05913e-01_jprb, 9.89690e-02_jprb,&
&9.24805e-02_jprb, 8.64175e-02_jprb, 8.07519e-02_jprb, 7.54578e-02_jprb, 7.05108e-02_jprb/)
selfrefo( :,10) = (/&
&1.37781e-01_jprb, 1.26830e-01_jprb, 1.16751e-01_jprb, 1.07472e-01_jprb, 9.89303e-02_jprb,&
&9.10678e-02_jprb, 8.38301e-02_jprb, 7.71677e-02_jprb, 7.10347e-02_jprb, 6.53892e-02_jprb/)
selfrefo( :,11) = (/&
&1.30675e-01_jprb, 1.24254e-01_jprb, 1.18149e-01_jprb, 1.12344e-01_jprb, 1.06824e-01_jprb,&
&1.01575e-01_jprb, 9.65839e-02_jprb, 9.18383e-02_jprb, 8.73258e-02_jprb, 8.30351e-02_jprb/)
selfrefo( :,12) = (/&
&1.55998e-01_jprb, 1.45155e-01_jprb, 1.35067e-01_jprb, 1.25679e-01_jprb, 1.16944e-01_jprb,&
&1.08816e-01_jprb, 1.01253e-01_jprb, 9.42160e-02_jprb, 8.76677e-02_jprb, 8.15746e-02_jprb/)
selfrefo( :,13) = (/&
&1.56508e-01_jprb, 1.46390e-01_jprb, 1.36925e-01_jprb, 1.28072e-01_jprb, 1.19792e-01_jprb,&
&1.12047e-01_jprb, 1.04803e-01_jprb, 9.80266e-02_jprb, 9.16888e-02_jprb, 8.57608e-02_jprb/)
selfrefo( :,14) = (/&
&1.55206e-01_jprb, 1.43469e-01_jprb, 1.32620e-01_jprb, 1.22591e-01_jprb, 1.13321e-01_jprb,&
&1.04752e-01_jprb, 9.68305e-02_jprb, 8.95082e-02_jprb, 8.27396e-02_jprb, 7.64829e-02_jprb/)
selfrefo( :,15) = (/&
&1.46694e-01_jprb, 1.35458e-01_jprb, 1.25082e-01_jprb, 1.15501e-01_jprb, 1.06654e-01_jprb,&
&9.84840e-02_jprb, 9.09403e-02_jprb, 8.39744e-02_jprb, 7.75420e-02_jprb, 7.16024e-02_jprb/)
selfrefo( :,16) = (/&
&1.65126e-01_jprb, 1.53281e-01_jprb, 1.42287e-01_jprb, 1.32081e-01_jprb, 1.22607e-01_jprb,&
&1.13812e-01_jprb, 1.05649e-01_jprb, 9.80707e-02_jprb, 9.10362e-02_jprb, 8.45063e-02_jprb/)

fracrefao( : ,1) = (/&
!     From P = 706.3 mb.
    &0.21245100_jprb,0.15164700_jprb,0.14486700_jprb,0.13075501_jprb,&
    &0.11629600_jprb,0.09266050_jprb,0.06579930_jprb,0.04524000_jprb,&
    &0.03072870_jprb,0.00284297_jprb,0.00234660_jprb,0.00185208_jprb,&
    &0.00133978_jprb,0.00082214_jprb,0.00031016_jprb,0.00004363_jprb/)
fracrefao( : ,2) = (/&
    &0.14703900_jprb,0.16937999_jprb,0.15605700_jprb,0.14159000_jprb,&
    &0.12088500_jprb,0.10058500_jprb,0.06809110_jprb,0.05131470_jprb,&
    &0.03487040_jprb,0.00327281_jprb,0.00250183_jprb,0.00190024_jprb,&
    &0.00133978_jprb,0.00082214_jprb,0.00031016_jprb,0.00004363_jprb/)
fracrefao( : ,3) = (/&
    &0.13689300_jprb,0.16610400_jprb,0.15723500_jprb,0.14299500_jprb,&
    &0.12399400_jprb,0.09907820_jprb,0.07169690_jprb,0.05367370_jprb,&
    &0.03671630_jprb,0.00378148_jprb,0.00290510_jprb,0.00221076_jprb,&
    &0.00142810_jprb,0.00093527_jprb,0.00031016_jprb,0.00004363_jprb/)
fracrefao( : ,4) = (/&
    &0.13054299_jprb,0.16273800_jprb,0.15874299_jprb,0.14279599_jprb,&
    &0.12674300_jprb,0.09664900_jprb,0.07462200_jprb,0.05620080_jprb,&
    &0.03789090_jprb,0.00411690_jprb,0.00322920_jprb,0.00245036_jprb,&
    &0.00178303_jprb,0.00098595_jprb,0.00040802_jprb,0.00010150_jprb/)
fracrefao( : ,5) = (/&
    &0.12828299_jprb,0.15824600_jprb,0.15688400_jprb,0.14449100_jprb,&
    &0.12787800_jprb,0.09517830_jprb,0.07679350_jprb,0.05890820_jprb,&
    &0.03883570_jprb,0.00442304_jprb,0.00346796_jprb,0.00255333_jprb,&
    &0.00212519_jprb,0.00116168_jprb,0.00067065_jprb,0.00010150_jprb/)
fracrefao( : ,6) = (/&
    &0.12649800_jprb,0.15195100_jprb,0.15646499_jprb,0.14569700_jprb,&
    &0.12669300_jprb,0.09653520_jprb,0.07887920_jprb,0.06106920_jprb,&
    &0.04043910_jprb,0.00430390_jprb,0.00364453_jprb,0.00314360_jprb,&
    &0.00203206_jprb,0.00187787_jprb,0.00067075_jprb,0.00010150_jprb/)
fracrefao( : ,7) = (/&
    &0.12500300_jprb,0.14460599_jprb,0.15672199_jprb,0.14724600_jprb,&
    &0.11978900_jprb,0.10190200_jprb,0.08196710_jprb,0.06315770_jprb,&
    &0.04240100_jprb,0.00433645_jprb,0.00404097_jprb,0.00329466_jprb,&
    &0.00288491_jprb,0.00187803_jprb,0.00067093_jprb,0.00010150_jprb/)
fracrefao( : ,8) = (/&
    &0.12317200_jprb,0.14118700_jprb,0.15242000_jprb,0.13794300_jprb,&
    &0.12119200_jprb,0.10655400_jprb,0.08808350_jprb,0.06521370_jprb,&
    &0.04505680_jprb,0.00485949_jprb,0.00477105_jprb,0.00401468_jprb,&
    &0.00288491_jprb,0.00187786_jprb,0.00067110_jprb,0.00010150_jprb/)
fracrefao( : ,9) = (/&
    &0.10193600_jprb,0.11693000_jprb,0.13236099_jprb,0.14053200_jprb,&
    &0.13749801_jprb,0.12193100_jprb,0.10221000_jprb,0.07448910_jprb,&
    &0.05205320_jprb,0.00572312_jprb,0.00476882_jprb,0.00403380_jprb,&
    &0.00288871_jprb,0.00187396_jprb,0.00067218_jprb,0.00010150_jprb/)

!     -----------------------------------------------------------------
RETURN
END SUBROUTINE rrtm_kgb12_00