rrtm_kgb9_00.F90

Go to the documentation of this file.

!OCL SCALAR
SUBROUTINE rrtm_kgb9_00

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

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

#include "tsmbkind.h"

USE yoerrto9 , ONLY : kao     ,kbo     ,selfrefo   ,fracrefao  ,&
           &fracrefbo, absn2oo
USE yoerrta9 , ONLY : n2oref  ,h2oref  ,ch4ref     ,etaref     ,strrat

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


IMPLICIT NONE
fracrefbo( :) = (/&
!     From P = 0.071 mb.
    &0.20148601_jprb,0.15252700_jprb,0.13376500_jprb,0.12184600_jprb,&
    &0.10767800_jprb,0.09307410_jprb,0.07674570_jprb,0.05876940_jprb,&
    &0.04001480_jprb,0.00424612_jprb,0.00346896_jprb,0.00269954_jprb,&
    &0.00196864_jprb,0.00122562_jprb,0.00043628_jprb,0.00004892_jprb/)

n2oref( :) = (/&
     &3.20000e-07_jprb,3.20000e-07_jprb,3.20000e-07_jprb,3.20000e-07_jprb,3.20000e-07_jprb,&
     &3.19652e-07_jprb,3.15324e-07_jprb,3.03830e-07_jprb,2.94221e-07_jprb,2.84953e-07_jprb,&
     &2.76714e-07_jprb,2.64709e-07_jprb,2.42847e-07_jprb/)

h2oref( :) = (/&
     &1.8759999e-02_jprb, 1.2223309e-02_jprb, 5.8908667e-03_jprb, 2.7675382e-03_jprb,  &
     &1.4065107e-03_jprb, 7.5969833e-04_jprb, 3.8875898e-04_jprb, 1.6542293e-04_jprb,  &
     &3.7189537e-05_jprb, 7.4764857e-06_jprb, 4.3081886e-06_jprb, 3.3319423e-06_jprb,  &
     &3.2039343e-06_jprb/)

ch4ref( :) = (/&
     &1.7000001e-06_jprb, 1.7000001e-06_jprb, 1.6998713e-06_jprb, 1.6904165e-06_jprb,  &
     &1.6671424e-06_jprb, 1.6350652e-06_jprb, 1.6097551e-06_jprb, 1.5590465e-06_jprb,  &
     &1.5119849e-06_jprb, 1.4741138e-06_jprb, 1.4384609e-06_jprb, 1.4002215e-06_jprb,  &
     &1.3573376e-06_jprb/)

etaref( :) = (/&
     &_zero_,0.125_jprb,0.25_jprb,0.375_jprb,_half_,0.625_jprb,0.75_jprb,0.875_jprb,0.96_jprb,0.99_jprb,_one_/)

absn2oo( :) = (/&
!     From P = 952.
     &3.26267e-01_jprb,2.42869e-00_jprb,1.15455e+01_jprb,7.39478e-00_jprb,&
     &5.16550e-00_jprb,2.54474e-00_jprb,3.53082e-00_jprb,3.82278e-00_jprb,&
     &1.81297e-00_jprb,6.65313e-01_jprb,1.23652e-01_jprb,1.83895e-03_jprb,&
     &1.70592e-03_jprb,2.68434e-09_jprb,_zero_,_zero_,&
!     From P = 620.
     &2.08632e-01_jprb,1.11865e+00_jprb,4.95975e+00_jprb,8.10907e+00_jprb,&
     &1.10408e+01_jprb,5.45460e+00_jprb,4.18611e+00_jprb,3.53422e+00_jprb,&
     &2.54164e+00_jprb,3.65093e-01_jprb,5.84480e-01_jprb,2.26918e-01_jprb,&
     &1.36230e-03_jprb,5.54400e-10_jprb,6.83703e-10_jprb,_zero_,&
!     From P=313.
     &6.20022e-02_jprb,2.69521e-01_jprb,9.81928e-01_jprb,1.65004e-00_jprb,&
     &3.08089e-00_jprb,5.38696e-00_jprb,1.14600e+01_jprb,2.41211e+01_jprb,&
     &1.69655e+01_jprb,1.37556e-00_jprb,5.43254e-01_jprb,3.52079e-01_jprb,&
     &4.31888e-01_jprb,4.82523e-06_jprb,5.74747e-11_jprb,_zero_/)

strrat = 21.6282_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.83413e-02_jprb, 2.51443e-02_jprb, 2.23078e-02_jprb, 1.97914e-02_jprb, 1.75588e-02_jprb,&
&1.55780e-02_jprb, 1.38207e-02_jprb, 1.22617e-02_jprb, 1.08785e-02_jprb, 9.65132e-03_jprb/)
selfrefo( :, 2) = (/&
&3.01506e-02_jprb, 2.69506e-02_jprb, 2.40901e-02_jprb, 2.15333e-02_jprb, 1.92478e-02_jprb,&
&1.72049e-02_jprb, 1.53789e-02_jprb, 1.37466e-02_jprb, 1.22876e-02_jprb, 1.09834e-02_jprb/)
selfrefo( :, 3) = (/&
&4.24899e-02_jprb, 3.77138e-02_jprb, 3.34745e-02_jprb, 2.97118e-02_jprb, 2.63720e-02_jprb,&
&2.34077e-02_jprb, 2.07765e-02_jprb, 1.84411e-02_jprb, 1.63682e-02_jprb, 1.45283e-02_jprb/)
selfrefo( :, 4) = (/&
&5.78648e-02_jprb, 5.16988e-02_jprb, 4.61899e-02_jprb, 4.12679e-02_jprb, 3.68705e-02_jprb,&
&3.29416e-02_jprb, 2.94314e-02_jprb, 2.62952e-02_jprb, 2.34932e-02_jprb, 2.09898e-02_jprb/)
selfrefo( :, 5) = (/&
&6.96367e-02_jprb, 6.36716e-02_jprb, 5.82175e-02_jprb, 5.32306e-02_jprb, 4.86709e-02_jprb,&
&4.45018e-02_jprb, 4.06898e-02_jprb, 3.72043e-02_jprb, 3.40174e-02_jprb, 3.11034e-02_jprb/)
selfrefo( :, 6) = (/&
&6.04073e-02_jprb, 5.90856e-02_jprb, 5.77929e-02_jprb, 5.65284e-02_jprb, 5.52916e-02_jprb,&
&5.40819e-02_jprb, 5.28986e-02_jprb, 5.17413e-02_jprb, 5.06092e-02_jprb, 4.95020e-02_jprb/)
selfrefo( :, 7) = (/&
&6.19460e-02_jprb, 6.02812e-02_jprb, 5.86612e-02_jprb, 5.70847e-02_jprb, 5.55506e-02_jprb,&
&5.40577e-02_jprb, 5.26049e-02_jprb, 5.11912e-02_jprb, 4.98155e-02_jprb, 4.84767e-02_jprb/)
selfrefo( :, 8) = (/&
&6.86939e-02_jprb, 6.58429e-02_jprb, 6.31102e-02_jprb, 6.04909e-02_jprb, 5.79804e-02_jprb,&
&5.55740e-02_jprb, 5.32675e-02_jprb, 5.10567e-02_jprb, 4.89377e-02_jprb, 4.69067e-02_jprb/)
selfrefo( :, 9) = (/&
&8.18828e-02_jprb, 7.77172e-02_jprb, 7.37635e-02_jprb, 7.00110e-02_jprb, 6.64494e-02_jprb,&
&6.30689e-02_jprb, 5.98604e-02_jprb, 5.68152e-02_jprb, 5.39248e-02_jprb, 5.11815e-02_jprb/)
selfrefo( :,10) = (/&
&8.41953e-02_jprb, 8.10610e-02_jprb, 7.80434e-02_jprb, 7.51381e-02_jprb, 7.23409e-02_jprb,&
&6.96479e-02_jprb, 6.70551e-02_jprb, 6.45589e-02_jprb, 6.21556e-02_jprb, 5.98417e-02_jprb/)
selfrefo( :,11) = (/&
&9.41304e-02_jprb, 8.97275e-02_jprb, 8.55305e-02_jprb, 8.15299e-02_jprb, 7.77164e-02_jprb,&
&7.40812e-02_jprb, 7.06161e-02_jprb, 6.73131e-02_jprb, 6.41646e-02_jprb, 6.11633e-02_jprb/)
selfrefo( :,12) = (/&
&9.55404e-02_jprb, 9.02311e-02_jprb, 8.52169e-02_jprb, 8.04813e-02_jprb, 7.60089e-02_jprb,&
&7.17850e-02_jprb, 6.77958e-02_jprb, 6.40284e-02_jprb, 6.04703e-02_jprb, 5.71099e-02_jprb/)
selfrefo( :,13) = (/&
&9.32377e-02_jprb, 8.99082e-02_jprb, 8.66977e-02_jprb, 8.36018e-02_jprb, 8.06164e-02_jprb,&
&7.77376e-02_jprb, 7.49617e-02_jprb, 7.22848e-02_jprb, 6.97036e-02_jprb, 6.72145e-02_jprb/)
selfrefo( :,14) = (/&
&1.14746e-01_jprb, 1.08112e-01_jprb, 1.01861e-01_jprb, 9.59719e-02_jprb, 9.04232e-02_jprb,&
&8.51952e-02_jprb, 8.02695e-02_jprb, 7.56285e-02_jprb, 7.12559e-02_jprb, 6.71361e-02_jprb/)
selfrefo( :,15) = (/&
&1.23517e-01_jprb, 1.19527e-01_jprb, 1.15666e-01_jprb, 1.11929e-01_jprb, 1.08313e-01_jprb,&
&1.04814e-01_jprb, 1.01428e-01_jprb, 9.81512e-02_jprb, 9.49803e-02_jprb, 9.19119e-02_jprb/)
selfrefo( :,16) = (/&
&1.44086e-01_jprb, 1.36396e-01_jprb, 1.29117e-01_jprb, 1.22226e-01_jprb, 1.15703e-01_jprb,&
&1.09528e-01_jprb, 1.03683e-01_jprb, 9.81497e-02_jprb, 9.29116e-02_jprb, 8.79530e-02_jprb/)

fracrefao( :, 1) = (/&
!     From P = 1053.6 mb.
    &0.16898900_jprb,0.15898301_jprb,0.13575301_jprb,0.12600900_jprb,&
    &0.11545800_jprb,0.09879170_jprb,0.08106830_jprb,0.06063440_jprb,&
    &0.03988780_jprb,0.00421760_jprb,0.00346635_jprb,0.00278779_jprb,&
    &0.00206225_jprb,0.00132324_jprb,0.00050033_jprb,0.00007038_jprb/)
fracrefao( :, 2) = (/&
    &0.18209399_jprb,0.15315101_jprb,0.13571000_jprb,0.12504999_jprb,&
    &0.11379100_jprb,0.09680810_jprb,0.08008570_jprb,0.05970280_jprb,&
    &0.03942860_jprb,0.00413383_jprb,0.00343186_jprb,0.00275558_jprb,&
    &0.00204657_jprb,0.00130219_jprb,0.00045454_jprb,0.00005664_jprb/)
fracrefao( :, 3) = (/&
    &0.18459500_jprb,0.15512000_jprb,0.13395500_jprb,0.12576801_jprb,&
    &0.11276800_jprb,0.09645190_jprb,0.07956650_jprb,0.05903340_jprb,&
    &0.03887050_jprb,0.00412226_jprb,0.00339453_jprb,0.00273518_jprb,&
    &0.00196922_jprb,0.00119411_jprb,0.00040263_jprb,0.00005664_jprb/)
fracrefao( :, 4) = (/&
    &0.18458800_jprb,0.15859900_jprb,0.13278100_jprb,0.12589300_jprb,&
    &0.11272700_jprb,0.09599660_jprb,0.07903030_jprb,0.05843600_jprb,&
    &0.03843400_jprb,0.00405181_jprb,0.00337980_jprb,0.00263818_jprb,&
    &0.00186869_jprb,0.00111807_jprb,0.00040263_jprb,0.00005664_jprb/)
fracrefao( :, 5) = (/&
    &0.18459301_jprb,0.16176100_jprb,0.13235000_jprb,0.12528200_jprb,&
    &0.11237100_jprb,0.09618840_jprb,0.07833760_jprb,0.05800770_jprb,&
    &0.03787610_jprb,0.00408253_jprb,0.00330363_jprb,0.00250445_jprb,&
    &0.00176725_jprb,0.00111753_jprb,0.00040263_jprb,0.00005664_jprb/)
fracrefao( :, 6) = (/&
    &0.18454400_jprb,0.16505300_jprb,0.13221300_jprb,0.12476600_jprb,&
    &0.11158300_jprb,0.09618120_jprb,0.07797340_jprb,0.05740380_jprb,&
    &0.03742820_jprb,0.00392691_jprb,0.00312208_jprb,0.00246306_jprb,&
    &0.00176735_jprb,0.00111721_jprb,0.00040263_jprb,0.00005664_jprb/)
fracrefao( :, 7) = (/&
    &0.18452001_jprb,0.16697501_jprb,0.13445500_jprb,0.12391300_jprb,&
    &0.11059100_jprb,0.09596890_jprb,0.07761050_jprb,0.05643200_jprb,&
    &0.03686520_jprb,0.00377086_jprb,0.00309351_jprb,0.00246297_jprb,&
    &0.00176765_jprb,0.00111700_jprb,0.00040263_jprb,0.00005664_jprb/)
fracrefao( :, 8) = (/&
    &0.18460999_jprb,0.16854499_jprb,0.13922299_jprb,0.12266400_jprb,&
    &0.10962200_jprb,0.09452030_jprb,0.07653800_jprb,0.05551340_jprb,&
    &0.03609660_jprb,0.00377043_jprb,0.00309367_jprb,0.00246304_jprb,&
    &0.00176749_jprb,0.00111689_jprb,0.00040263_jprb,0.00005664_jprb/)
fracrefao( :, 9) = (/&
    &0.18312500_jprb,0.16787501_jprb,0.14720701_jprb,0.12766500_jprb,&
    &0.10890900_jprb,0.08935530_jprb,0.07310870_jprb,0.05443140_jprb,&
    &0.03566380_jprb,0.00376446_jprb,0.00309521_jprb,0.00246510_jprb,&
    &0.00176139_jprb,0.00111543_jprb,0.00040263_jprb,0.00005664_jprb/)

!     -----------------------------------------------------------------
RETURN
END SUBROUTINE rrtm_kgb9_00