rrtm_kgb7_00.F90

Go to the documentation of this file.

!OCL SCALAR
SUBROUTINE rrtm_kgb7_00

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

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

#include "tsmbkind.h"

USE yoerrto7 , ONLY : kao     ,kbo     ,selfrefo   ,fracrefao  ,&
           &fracrefbo  ,absco2o
USE yoerrta7 , ONLY : strrat

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


IMPLICIT NONE
fracrefbo( :) = (/&
    &0.15355594_jprb,0.15310939_jprb,0.14274909_jprb,0.13129812_jprb,&
    &0.11736792_jprb,0.10118213_jprb,0.08215259_jprb,0.06165591_jprb,&
    &0.04164486_jprb,0.00451141_jprb,0.00372837_jprb,0.00294095_jprb,&
    &0.00215259_jprb,0.00136792_jprb,0.00051233_jprb,0.00007075_jprb/)

absco2o( :) = (/&
     &9.30038e-05_jprb, 1.74061e-04_jprb, 2.09293e-04_jprb, 2.52360e-04_jprb,&
     &3.13404e-04_jprb, 4.16619e-04_jprb, 6.27394e-04_jprb, 1.29386e-03_jprb,&
     &4.05192e-03_jprb, 3.97050e-03_jprb, 7.00634e-04_jprb, 6.06617e-04_jprb,&
     &7.66978e-04_jprb, 6.70661e-04_jprb, 7.89971e-04_jprb, 7.55709e-04_jprb/)

strrat = 8.21104e+04_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) = (/&
&5.19175e-02_jprb, 4.28815e-02_jprb, 3.54182e-02_jprb, 2.92539e-02_jprb, 2.41624e-02_jprb,&
&1.99571e-02_jprb, 1.64836e-02_jprb, 1.36148e-02_jprb, 1.12452e-02_jprb, 9.28801e-03_jprb/)
selfrefo( :, 2) = (/&
&4.36985e-02_jprb, 3.78503e-02_jprb, 3.27848e-02_jprb, 2.83973e-02_jprb, 2.45969e-02_jprb,&
&2.13051e-02_jprb, 1.84538e-02_jprb, 1.59842e-02_jprb, 1.38450e-02_jprb, 1.19921e-02_jprb/)
selfrefo( :, 3) = (/&
&4.25575e-02_jprb, 3.71214e-02_jprb, 3.23797e-02_jprb, 2.82437e-02_jprb, 2.46360e-02_jprb,&
&2.14891e-02_jprb, 1.87442e-02_jprb, 1.63499e-02_jprb, 1.42615e-02_jprb, 1.24398e-02_jprb/)
selfrefo( :, 4) = (/&
&4.03854e-02_jprb, 3.54677e-02_jprb, 3.11488e-02_jprb, 2.73558e-02_jprb, 2.40247e-02_jprb,&
&2.10992e-02_jprb, 1.85300e-02_jprb, 1.62736e-02_jprb, 1.42920e-02_jprb, 1.25517e-02_jprb/)
selfrefo( :, 5) = (/&
&3.95070e-02_jprb, 3.46726e-02_jprb, 3.04297e-02_jprb, 2.67060e-02_jprb, 2.34380e-02_jprb,&
&2.05699e-02_jprb, 1.80528e-02_jprb, 1.58437e-02_jprb, 1.39049e-02_jprb, 1.22034e-02_jprb/)
selfrefo( :, 6) = (/&
&3.89748e-02_jprb, 3.40186e-02_jprb, 2.96926e-02_jprb, 2.59168e-02_jprb, 2.26211e-02_jprb,&
&1.97445e-02_jprb, 1.72337e-02_jprb, 1.50422e-02_jprb, 1.31293e-02_jprb, 1.14597e-02_jprb/)
selfrefo( :, 7) = (/&
&3.85527e-02_jprb, 3.36561e-02_jprb, 2.93815e-02_jprb, 2.56497e-02_jprb, 2.23920e-02_jprb,&
&1.95480e-02_jprb, 1.70652e-02_jprb, 1.48977e-02_jprb, 1.30056e-02_jprb, 1.13537e-02_jprb/)
selfrefo( :, 8) = (/&
&3.78212e-02_jprb, 3.30614e-02_jprb, 2.89005e-02_jprb, 2.52633e-02_jprb, 2.20839e-02_jprb,&
&1.93046e-02_jprb, 1.68751e-02_jprb, 1.47513e-02_jprb, 1.28948e-02_jprb, 1.12720e-02_jprb/)
selfrefo( :, 9) = (/&
&3.68624e-02_jprb, 3.22676e-02_jprb, 2.82456e-02_jprb, 2.47248e-02_jprb, 2.16429e-02_jprb,&
&1.89452e-02_jprb, 1.65838e-02_jprb, 1.45166e-02_jprb, 1.27072e-02_jprb, 1.11233e-02_jprb/)
selfrefo( :,10) = (/&
&3.63135e-02_jprb, 3.19332e-02_jprb, 2.80813e-02_jprb, 2.46940e-02_jprb, 2.17154e-02_jprb,&
&1.90960e-02_jprb, 1.67926e-02_jprb, 1.47670e-02_jprb, 1.29857e-02_jprb, 1.14193e-02_jprb/)
selfrefo( :,11) = (/&
&3.60018e-02_jprb, 3.17294e-02_jprb, 2.79640e-02_jprb, 2.46454e-02_jprb, 2.17207e-02_jprb,&
&1.91430e-02_jprb, 1.68713e-02_jprb, 1.48691e-02_jprb, 1.31045e-02_jprb, 1.15494e-02_jprb/)
selfrefo( :,12) = (/&
&3.65147e-02_jprb, 3.19158e-02_jprb, 2.78962e-02_jprb, 2.43827e-02_jprb, 2.13118e-02_jprb,&
&1.86277e-02_jprb, 1.62816e-02_jprb, 1.42310e-02_jprb, 1.24386e-02_jprb, 1.08720e-02_jprb/)
selfrefo( :,13) = (/&
&3.72684e-02_jprb, 3.23732e-02_jprb, 2.81210e-02_jprb, 2.44274e-02_jprb, 2.12189e-02_jprb,&
&1.84318e-02_jprb, 1.60108e-02_jprb, 1.39078e-02_jprb, 1.20810e-02_jprb, 1.04942e-02_jprb/)
selfrefo( :,14) = (/&
&3.87604e-02_jprb, 3.34099e-02_jprb, 2.87979e-02_jprb, 2.48226e-02_jprb, 2.13960e-02_jprb,&
&1.84425e-02_jprb, 1.58967e-02_jprb, 1.37023e-02_jprb, 1.18108e-02_jprb, 1.01804e-02_jprb/)
selfrefo( :,15) = (/&
&3.51227e-02_jprb, 3.05052e-02_jprb, 2.64947e-02_jprb, 2.30115e-02_jprb, 1.99862e-02_jprb,&
&1.73586e-02_jprb, 1.50765e-02_jprb, 1.30944e-02_jprb, 1.13729e-02_jprb, 9.87775e-03_jprb/)
selfrefo( :,16) = (/&
&3.51227e-02_jprb, 3.05052e-02_jprb, 2.64947e-02_jprb, 2.30115e-02_jprb, 1.99862e-02_jprb,&
&1.73586e-02_jprb, 1.50765e-02_jprb, 1.30944e-02_jprb, 1.13729e-02_jprb, 9.87775e-03_jprb/)

fracrefao( :, 1) = (/&
  &0.16461779_jprb, 0.14889984_jprb, 0.14233345_jprb, 0.13156526_jprb,&
  &0.11679733_jprb, 0.09988949_jprb, 0.08078653_jprb, 0.06006384_jprb,&
  &0.04028391_jprb, 0.00435899_jprb, 0.00359173_jprb, 0.00281707_jprb,&
  &0.00206767_jprb, 0.00135012_jprb, 0.00050720_jprb, 0.00007146_jprb/)
fracrefao( :, 2) = (/&
  &0.16442357_jprb, 0.14944240_jprb, 0.14245804_jprb, 0.13111183_jprb,&
  &0.11688625_jprb, 0.09983791_jprb, 0.08085148_jprb, 0.05993948_jprb,&
  &0.04028057_jprb, 0.00435939_jprb, 0.00358708_jprb, 0.00284036_jprb,&
  &0.00208869_jprb, 0.00133256_jprb, 0.00049260_jprb, 0.00006931_jprb/)
fracrefao( :, 3) = (/&
  &0.16368519_jprb, 0.15018989_jprb, 0.14262174_jprb, 0.13084342_jprb,&
  &0.11682195_jprb, 0.09996257_jprb, 0.08074036_jprb, 0.05985692_jprb,&
  &0.04045362_jprb, 0.00436208_jprb, 0.00358257_jprb, 0.00287122_jprb,&
  &0.00211004_jprb, 0.00133804_jprb, 0.00049260_jprb, 0.00006931_jprb/)
fracrefao( :, 4) = (/&
  &0.16274056_jprb, 0.15133780_jprb, 0.14228874_jprb, 0.13081114_jprb,&
  &0.11688486_jprb, 0.09979610_jprb, 0.08073687_jprb, 0.05996741_jprb,&
  &0.04040616_jprb, 0.00439869_jprb, 0.00368910_jprb, 0.00293041_jprb,&
  &0.00211604_jprb, 0.00133536_jprb, 0.00049260_jprb, 0.00006931_jprb/)
fracrefao( :, 5) = (/&
  &0.16176532_jprb, 0.15207882_jprb, 0.14226955_jprb, 0.13079646_jprb,&
  &0.11688191_jprb, 0.09966998_jprb, 0.08066384_jprb, 0.06020275_jprb,&
  &0.04047901_jprb, 0.00446696_jprb, 0.00377456_jprb, 0.00294410_jprb,&
  &0.00211082_jprb, 0.00133536_jprb, 0.00049260_jprb, 0.00006931_jprb/)
fracrefao( :, 6) = (/&
  &0.15993737_jprb, 0.15305527_jprb, 0.14259829_jprb, 0.13078023_jprb,&
  &0.11686983_jprb, 0.09980131_jprb, 0.08058286_jprb, 0.06031430_jprb,&
  &0.04082833_jprb, 0.00450509_jprb, 0.00377574_jprb, 0.00294823_jprb,&
  &0.00210977_jprb, 0.00133302_jprb, 0.00049260_jprb, 0.00006931_jprb/)
fracrefao( :, 7) = (/&
  &0.15371189_jprb, 0.15592396_jprb, 0.14430280_jprb, 0.13076764_jprb,&
  &0.11720382_jprb, 0.10023471_jprb, 0.08066396_jprb, 0.06073554_jprb,&
  &0.04121581_jprb, 0.00451202_jprb, 0.00377832_jprb, 0.00294609_jprb,&
  &0.00210943_jprb, 0.00133336_jprb, 0.00049260_jprb, 0.00006931_jprb/)
fracrefao( :, 8) = (/&
  &0.14262275_jprb, 0.14572631_jprb, 0.14560597_jprb, 0.13736825_jprb,&
  &0.12271351_jprb, 0.10419556_jprb, 0.08294533_jprb, 0.06199794_jprb,&
  &0.04157615_jprb, 0.00452842_jprb, 0.00377704_jprb, 0.00293852_jprb,&
  &0.00211034_jprb, 0.00133278_jprb, 0.00049259_jprb, 0.00006931_jprb/)
fracrefao( :, 9) = (/&
  &0.14500433_jprb, 0.14590444_jprb, 0.14430299_jprb, 0.13770708_jprb,&
  &0.12288283_jprb, 0.10350952_jprb, 0.08269450_jprb, 0.06130579_jprb,&
  &0.04144571_jprb, 0.00452096_jprb, 0.00377382_jprb, 0.00294532_jprb,&
  &0.00210943_jprb, 0.00133228_jprb, 0.00049260_jprb, 0.00006931_jprb/)

!     -----------------------------------------------------------------
RETURN
END SUBROUTINE rrtm_kgb7_00