rrtm_kgb5_00.F90

Go to the documentation of this file.

!OCL SCALAR
SUBROUTINE rrtm_kgb5_00

!     Originally by Eli J. Mlawer, Atmospheric & Environmental Research.
!     BAND 5:  700-820 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 yoerrto5 , ONLY : kao     ,kbo     ,selfrefo   ,fracrefao  ,&
           &fracrefbo, ccl4o
USE yoerrta5 , ONLY : strrat1   ,strrat2

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


IMPLICIT NONE
ccl4o( :) = (/&
     &26.1407_jprb,  53.9776_jprb,  63.8085_jprb,  36.1701_jprb,&
     &15.4099_jprb, 10.23116_jprb,  4.82948_jprb,  5.03836_jprb,&
     &1.75558_jprb,  _zero_     ,  _zero_     ,  _zero_     ,&
     &_zero_     ,  _zero_     ,  _zero_     ,  _zero_      /)

strrat1 = 90.4894_jprb
strrat2 = 0.900502_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.27664e-01_jprb, 1.09296e-01_jprb, 9.35703e-02_jprb, 8.01075e-02_jprb, 6.85818e-02_jprb,&
&5.87143e-02_jprb, 5.02666e-02_jprb, 4.30343e-02_jprb, 3.68426e-02_jprb, 3.15417e-02_jprb/)
selfrefo( :, 2) = (/&
&1.39620e-01_jprb, 1.20381e-01_jprb, 1.03793e-01_jprb, 8.94908e-02_jprb, 7.71595e-02_jprb,&
&6.65273e-02_jprb, 5.73601e-02_jprb, 4.94562e-02_jprb, 4.26414e-02_jprb, 3.67656e-02_jprb/)
selfrefo( :, 3) = (/&
&1.42628e-01_jprb, 1.23043e-01_jprb, 1.06148e-01_jprb, 9.15726e-02_jprb, 7.89986e-02_jprb,&
&6.81511e-02_jprb, 5.87931e-02_jprb, 5.07201e-02_jprb, 4.37556e-02_jprb, 3.77474e-02_jprb/)
selfrefo( :, 4) = (/&
&1.53569e-01_jprb, 1.33143e-01_jprb, 1.15435e-01_jprb, 1.00082e-01_jprb, 8.67706e-02_jprb,&
&7.52299e-02_jprb, 6.52241e-02_jprb, 5.65491e-02_jprb, 4.90279e-02_jprb, 4.25070e-02_jprb/)
selfrefo( :, 5) = (/&
&1.70491e-01_jprb, 1.46448e-01_jprb, 1.25796e-01_jprb, 1.08056e-01_jprb, 9.28182e-02_jprb,&
&7.97290e-02_jprb, 6.84856e-02_jprb, 5.88278e-02_jprb, 5.05319e-02_jprb, 4.34059e-02_jprb/)
selfrefo( :, 6) = (/&
&1.76394e-01_jprb, 1.51432e-01_jprb, 1.30003e-01_jprb, 1.11606e-01_jprb, 9.58127e-02_jprb,&
&8.22542e-02_jprb, 7.06144e-02_jprb, 6.06217e-02_jprb, 5.20431e-02_jprb, 4.46784e-02_jprb/)
selfrefo( :, 7) = (/&
&1.85706e-01_jprb, 1.59172e-01_jprb, 1.36429e-01_jprb, 1.16936e-01_jprb, 1.00228e-01_jprb,&
&8.59068e-02_jprb, 7.36322e-02_jprb, 6.31114e-02_jprb, 5.40939e-02_jprb, 4.63648e-02_jprb/)
selfrefo( :, 8) = (/&
&1.88647e-01_jprb, 1.61657e-01_jprb, 1.38529e-01_jprb, 1.18710e-01_jprb, 1.01726e-01_jprb,&
&8.71722e-02_jprb, 7.47005e-02_jprb, 6.40132e-02_jprb, 5.48549e-02_jprb, 4.70068e-02_jprb/)
selfrefo( :, 9) = (/&
&1.90074e-01_jprb, 1.62793e-01_jprb, 1.39427e-01_jprb, 1.19415e-01_jprb, 1.02275e-01_jprb,&
&8.75959e-02_jprb, 7.50233e-02_jprb, 6.42552e-02_jprb, 5.50327e-02_jprb, 4.71338e-02_jprb/)
selfrefo( :,10) = (/&
&1.94769e-01_jprb, 1.66338e-01_jprb, 1.42057e-01_jprb, 1.21320e-01_jprb, 1.03611e-01_jprb,&
&8.84863e-02_jprb, 7.55696e-02_jprb, 6.45384e-02_jprb, 5.51175e-02_jprb, 4.70718e-02_jprb/)
selfrefo( :,11) = (/&
&1.90624e-01_jprb, 1.64229e-01_jprb, 1.41488e-01_jprb, 1.21896e-01_jprb, 1.05017e-01_jprb,&
&9.04757e-02_jprb, 7.79475e-02_jprb, 6.71542e-02_jprb, 5.78554e-02_jprb, 4.98442e-02_jprb/)
selfrefo( :,12) = (/&
&1.90502e-01_jprb, 1.64025e-01_jprb, 1.41228e-01_jprb, 1.21599e-01_jprb, 1.04699e-01_jprb,&
&9.01472e-02_jprb, 7.76181e-02_jprb, 6.68303e-02_jprb, 5.75419e-02_jprb, 4.95444e-02_jprb/)
selfrefo( :,13) = (/&
&1.86786e-01_jprb, 1.61636e-01_jprb, 1.39872e-01_jprb, 1.21039e-01_jprb, 1.04741e-01_jprb,&
&9.06380e-02_jprb, 7.84338e-02_jprb, 6.78729e-02_jprb, 5.87340e-02_jprb, 5.08256e-02_jprb/)
selfrefo( :,14) = (/&
&1.99149e-01_jprb, 1.71475e-01_jprb, 1.47646e-01_jprb, 1.27129e-01_jprb, 1.09462e-01_jprb,&
&9.42512e-02_jprb, 8.11538e-02_jprb, 6.98764e-02_jprb, 6.01662e-02_jprb, 5.18053e-02_jprb/)
selfrefo( :,15) = (/&
&2.02676e-01_jprb, 1.73701e-01_jprb, 1.48869e-01_jprb, 1.27587e-01_jprb, 1.09347e-01_jprb,&
&9.37144e-02_jprb, 8.03170e-02_jprb, 6.88348e-02_jprb, 5.89941e-02_jprb, 5.05603e-02_jprb/)
selfrefo( :,16) = (/&
&1.99865e-01_jprb, 1.72699e-01_jprb, 1.49225e-01_jprb, 1.28942e-01_jprb, 1.11416e-01_jprb,&
&9.62721e-02_jprb, 8.31866e-02_jprb, 7.18797e-02_jprb, 6.21097e-02_jprb, 5.36676e-02_jprb/)      

fracrefao( :, 1) = (/&
!     From P = 387.6 mb.
    &0.13966499_jprb,0.14138900_jprb,0.13763399_jprb,0.13076700_jprb,&
    &0.12299100_jprb,0.10747700_jprb,0.08942000_jprb,0.06769200_jprb,&
    &0.04587610_jprb,0.00501173_jprb,0.00415809_jprb,0.00328398_jprb,&
    &0.00240015_jprb,0.00156222_jprb,0.00059104_jprb,0.00008323_jprb/)
fracrefao( :, 2) = (/&
    &0.13958199_jprb,0.14332899_jprb,0.13785399_jprb,0.13205400_jprb,&
    &0.12199700_jprb,0.10679600_jprb,0.08861080_jprb,0.06712320_jprb,&
    &0.04556030_jprb,0.00500863_jprb,0.00416315_jprb,0.00328629_jprb,&
    &0.00240023_jprb,0.00156220_jprb,0.00059104_jprb,0.00008323_jprb/)
fracrefao( :, 3) = (/&
    &0.13907100_jprb,0.14250501_jprb,0.13889600_jprb,0.13297300_jprb,&
    &0.12218700_jprb,0.10683800_jprb,0.08839260_jprb,0.06677310_jprb,&
    &0.04538570_jprb,0.00495402_jprb,0.00409863_jprb,0.00328219_jprb,&
    &0.00240805_jprb,0.00156266_jprb,0.00059104_jprb,0.00008323_jprb/)
fracrefao( :, 4) = (/&
    &0.13867700_jprb,0.14190100_jprb,0.13932300_jprb,0.13327099_jprb,&
    &0.12280800_jprb,0.10692500_jprb,0.08844510_jprb,0.06658510_jprb,&
    &0.04519340_jprb,0.00492276_jprb,0.00408832_jprb,0.00323856_jprb,&
    &0.00239289_jprb,0.00155698_jprb,0.00059104_jprb,0.00008323_jprb/)
fracrefao( :, 5) = (/&
    &0.13845000_jprb,0.14158800_jprb,0.13929300_jprb,0.13295600_jprb,&
    &0.12348300_jprb,0.10736700_jprb,0.08859480_jprb,0.06650610_jprb,&
    &0.04498230_jprb,0.00491335_jprb,0.00406968_jprb,0.00322901_jprb,&
    &0.00234666_jprb,0.00155235_jprb,0.00058813_jprb,0.00008323_jprb/)
fracrefao( :, 6) = (/&
    &0.13837101_jprb,0.14113200_jprb,0.13930500_jprb,0.13283101_jprb,&
    &0.12349200_jprb,0.10796400_jprb,0.08890490_jprb,0.06646480_jprb,&
    &0.04485990_jprb,0.00489554_jprb,0.00405264_jprb,0.00320313_jprb,&
    &0.00234742_jprb,0.00151159_jprb,0.00058438_jprb,0.00008253_jprb/)
fracrefao( :, 7) = (/&
    &0.13834500_jprb,0.14093500_jprb,0.13896500_jprb,0.13262001_jprb,&
    &0.12326900_jprb,0.10828900_jprb,0.08950050_jprb,0.06674610_jprb,&
    &0.04476560_jprb,0.00489624_jprb,0.00400962_jprb,0.00317423_jprb,&
    &0.00233479_jprb,0.00148249_jprb,0.00058590_jprb,0.00008253_jprb/)
fracrefao( :, 8) = (/&
    &0.13831300_jprb,0.14069000_jprb,0.13871400_jprb,0.13247600_jprb,&
    &0.12251400_jprb,0.10831300_jprb,0.08977090_jprb,0.06776920_jprb,&
    &0.04498390_jprb,0.00484111_jprb,0.00398948_jprb,0.00316069_jprb,&
    &0.00229741_jprb,0.00150104_jprb,0.00058608_jprb,0.00008253_jprb/)
fracrefao( :, 9) = (/&
    &0.14027201_jprb,0.14420401_jprb,0.14215700_jprb,0.13446601_jprb,&
    &0.12303700_jprb,0.10596100_jprb,0.08650370_jprb,0.06409570_jprb,&
    &0.04312310_jprb,0.00471110_jprb,0.00393954_jprb,0.00310850_jprb,&
    &0.00229588_jprb,0.00146366_jprb,0.00058194_jprb,0.00008253_jprb/)

fracrefbo( :, 1) = (/&
!     From P = 1.17 mb.
    &0.14339100_jprb,0.14358699_jprb,0.13935301_jprb,0.13306700_jprb,&
    &0.12135700_jprb,0.10590600_jprb,0.08688240_jprb,0.06553220_jprb,&
    &0.04446740_jprb,0.00483580_jprb,0.00399413_jprb,0.00316225_jprb,&
    &0.00233007_jprb,0.00149135_jprb,0.00056246_jprb,0.00008059_jprb/)
fracrefbo( :, 2) = (/&
    &0.14330500_jprb,0.14430299_jprb,0.14053699_jprb,0.13355300_jprb,&
    &0.12151200_jprb,0.10529100_jprb,0.08627630_jprb,0.06505230_jprb,&
    &0.04385850_jprb,0.00476555_jprb,0.00395010_jprb,0.00313878_jprb,&
    &0.00232273_jprb,0.00149354_jprb,0.00056246_jprb,0.00008059_jprb/)
fracrefbo( :, 3) = (/&
    &0.14328399_jprb,0.14442700_jprb,0.14078601_jprb,0.13390100_jprb,&
    &0.12132600_jprb,0.10510600_jprb,0.08613660_jprb,0.06494630_jprb,&
    &0.04381310_jprb,0.00475378_jprb,0.00394166_jprb,0.00313076_jprb,&
    &0.00231235_jprb,0.00149159_jprb,0.00056301_jprb,0.00008059_jprb/)
fracrefbo( :, 4) = (/&
    &0.14326900_jprb,0.14453100_jprb,0.14114200_jprb,0.13397101_jprb,&
    &0.12127200_jprb,0.10493400_jprb,0.08601380_jprb,0.06483360_jprb,&
    &0.04378900_jprb,0.00474655_jprb,0.00393549_jprb,0.00312583_jprb,&
    &0.00230686_jprb,0.00148433_jprb,0.00056502_jprb,0.00008059_jprb/)
fracrefbo( :, 5) = (/&
    &0.14328900_jprb,0.14532700_jprb,0.14179000_jprb,0.13384600_jprb,&
    &0.12093700_jprb,0.10461500_jprb,0.08573010_jprb,0.06461340_jprb,&
    &0.04366570_jprb,0.00473087_jprb,0.00392539_jprb,0.00311238_jprb,&
    &0.00229865_jprb,0.00147572_jprb,0.00056517_jprb,0.00007939_jprb/)

!     -----------------------------------------------------------------
RETURN
END SUBROUTINE rrtm_kgb5_00