source: LMDZ6/branches/DYNAMICO-conv-GC/libf/phymar/rrtm_kgb9_00.F90 @ 4249

!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