source: LMDZ5/trunk/libf/phymar/rrtm_kgb12_00.F90 @ 2089

Last change on this file since 2089 was 2089, checked in by Laurent Fairhead, 10 years ago

Inclusion de la physique de MAR


Integration of MAR physics

File size: 7.0 KB
Line 
1!OCL SCALAR
2SUBROUTINE RRTM_KGB12_00
3
4!     Originally by Eli J. Mlawer, Atmospheric & Environmental Research.
5!     BAND 12:  1800-2080 cm-1 (low - H2O,CO2; high - nothing)
6!     Reformatted for F90 by JJMorcrette, ECMWF
7!     Reformatted for NEC by H.Gallée   , LGGE  (splitting)
8
9!     ------------------------------------------------------------------
10
11#include "tsmbkind.h"
12
13USE YOERRTO12, ONLY : KAO     ,SELFREFO   ,FRACREFAO
14USE YOERRTA12, ONLY : STRRAT
15
16!     ------------------------------------------------------------------
17
18
19IMPLICIT NONE
20STRRAT = 0.009736757_JPRB
21
22!     ------------------------------------------------------------------
23
24!     The array SELFREFO contains the coefficient of the water vapor
25!     self-continuum (including the energy term).  The first index
26!     refers to temperature in 7.2 degree increments.  For instance,
27!     JT = 1 refers to a temperature of 245.6, JT = 2 refers to 252.8,
28!     etc.  The second index runs over the g-channel (1 to 16).
29
30SELFREFO( :, 1) = (/&
31&2.34973E-02_JPRB, 2.07877E-02_JPRB, 1.83905E-02_JPRB, 1.62698E-02_JPRB, 1.43936E-02_JPRB,&
32&1.27338E-02_JPRB, 1.12654E-02_JPRB, 9.96632E-03_JPRB, 8.81704E-03_JPRB, 7.80030E-03_JPRB/)
33SELFREFO( :, 2) = (/&
34&3.10090E-02_JPRB, 2.82626E-02_JPRB, 2.57594E-02_JPRB, 2.34779E-02_JPRB, 2.13985E-02_JPRB,&
35&1.95033E-02_JPRB, 1.77759E-02_JPRB, 1.62015E-02_JPRB, 1.47665E-02_JPRB, 1.34587E-02_JPRB/)
36SELFREFO( :, 3) = (/&
37&5.36676E-02_JPRB, 4.95621E-02_JPRB, 4.57706E-02_JPRB, 4.22692E-02_JPRB, 3.90356E-02_JPRB,&
38&3.60494E-02_JPRB, 3.32917E-02_JPRB, 3.07449E-02_JPRB, 2.83929E-02_JPRB, 2.62209E-02_JPRB/)
39SELFREFO( :, 4) = (/&
40&9.05844E-02_JPRB, 8.33349E-02_JPRB, 7.66656E-02_JPRB, 7.05300E-02_JPRB, 6.48855E-02_JPRB,&
41&5.96927E-02_JPRB, 5.49155E-02_JPRB, 5.05206E-02_JPRB, 4.64774E-02_JPRB, 4.27578E-02_JPRB/)
42SELFREFO( :, 5) = (/&
43&1.10531E-01_JPRB, 1.03430E-01_JPRB, 9.67851E-02_JPRB, 9.05670E-02_JPRB, 8.47484E-02_JPRB,&
44&7.93037E-02_JPRB, 7.42087E-02_JPRB, 6.94411E-02_JPRB, 6.49797E-02_JPRB, 6.08050E-02_JPRB/)
45SELFREFO( :, 6) = (/&
46&1.24308E-01_JPRB, 1.16014E-01_JPRB, 1.08274E-01_JPRB, 1.01050E-01_JPRB, 9.43078E-02_JPRB,&
47&8.80155E-02_JPRB, 8.21431E-02_JPRB, 7.66624E-02_JPRB, 7.15475E-02_JPRB, 6.67738E-02_JPRB/)
48SELFREFO( :, 7) = (/&
49&1.34594E-01_JPRB, 1.25054E-01_JPRB, 1.16190E-01_JPRB, 1.07954E-01_JPRB, 1.00302E-01_JPRB,&
50&9.31920E-02_JPRB, 8.65863E-02_JPRB, 8.04488E-02_JPRB, 7.47464E-02_JPRB, 6.94481E-02_JPRB/)
51SELFREFO( :, 8) = (/&
52&1.26966E-01_JPRB, 1.18270E-01_JPRB, 1.10169E-01_JPRB, 1.02624E-01_JPRB, 9.55948E-02_JPRB,&
53&8.90474E-02_JPRB, 8.29484E-02_JPRB, 7.72672E-02_JPRB, 7.19750E-02_JPRB, 6.70453E-02_JPRB/)
54SELFREFO( :, 9) = (/&
55&1.29806E-01_JPRB, 1.21296E-01_JPRB, 1.13344E-01_JPRB, 1.05913E-01_JPRB, 9.89690E-02_JPRB,&
56&9.24805E-02_JPRB, 8.64175E-02_JPRB, 8.07519E-02_JPRB, 7.54578E-02_JPRB, 7.05108E-02_JPRB/)
57SELFREFO( :,10) = (/&
58&1.37781E-01_JPRB, 1.26830E-01_JPRB, 1.16751E-01_JPRB, 1.07472E-01_JPRB, 9.89303E-02_JPRB,&
59&9.10678E-02_JPRB, 8.38301E-02_JPRB, 7.71677E-02_JPRB, 7.10347E-02_JPRB, 6.53892E-02_JPRB/)
60SELFREFO( :,11) = (/&
61&1.30675E-01_JPRB, 1.24254E-01_JPRB, 1.18149E-01_JPRB, 1.12344E-01_JPRB, 1.06824E-01_JPRB,&
62&1.01575E-01_JPRB, 9.65839E-02_JPRB, 9.18383E-02_JPRB, 8.73258E-02_JPRB, 8.30351E-02_JPRB/)
63SELFREFO( :,12) = (/&
64&1.55998E-01_JPRB, 1.45155E-01_JPRB, 1.35067E-01_JPRB, 1.25679E-01_JPRB, 1.16944E-01_JPRB,&
65&1.08816E-01_JPRB, 1.01253E-01_JPRB, 9.42160E-02_JPRB, 8.76677E-02_JPRB, 8.15746E-02_JPRB/)
66SELFREFO( :,13) = (/&
67&1.56508E-01_JPRB, 1.46390E-01_JPRB, 1.36925E-01_JPRB, 1.28072E-01_JPRB, 1.19792E-01_JPRB,&
68&1.12047E-01_JPRB, 1.04803E-01_JPRB, 9.80266E-02_JPRB, 9.16888E-02_JPRB, 8.57608E-02_JPRB/)
69SELFREFO( :,14) = (/&
70&1.55206E-01_JPRB, 1.43469E-01_JPRB, 1.32620E-01_JPRB, 1.22591E-01_JPRB, 1.13321E-01_JPRB,&
71&1.04752E-01_JPRB, 9.68305E-02_JPRB, 8.95082E-02_JPRB, 8.27396E-02_JPRB, 7.64829E-02_JPRB/)
72SELFREFO( :,15) = (/&
73&1.46694E-01_JPRB, 1.35458E-01_JPRB, 1.25082E-01_JPRB, 1.15501E-01_JPRB, 1.06654E-01_JPRB,&
74&9.84840E-02_JPRB, 9.09403E-02_JPRB, 8.39744E-02_JPRB, 7.75420E-02_JPRB, 7.16024E-02_JPRB/)
75SELFREFO( :,16) = (/&
76&1.65126E-01_JPRB, 1.53281E-01_JPRB, 1.42287E-01_JPRB, 1.32081E-01_JPRB, 1.22607E-01_JPRB,&
77&1.13812E-01_JPRB, 1.05649E-01_JPRB, 9.80707E-02_JPRB, 9.10362E-02_JPRB, 8.45063E-02_JPRB/)
78
79FRACREFAO( : ,1) = (/&
80!     From P = 706.3 mb.
81    &0.21245100_JPRB,0.15164700_JPRB,0.14486700_JPRB,0.13075501_JPRB,&
82    &0.11629600_JPRB,0.09266050_JPRB,0.06579930_JPRB,0.04524000_JPRB,&
83    &0.03072870_JPRB,0.00284297_JPRB,0.00234660_JPRB,0.00185208_JPRB,&
84    &0.00133978_JPRB,0.00082214_JPRB,0.00031016_JPRB,0.00004363_JPRB/)
85FRACREFAO( : ,2) = (/&
86    &0.14703900_JPRB,0.16937999_JPRB,0.15605700_JPRB,0.14159000_JPRB,&
87    &0.12088500_JPRB,0.10058500_JPRB,0.06809110_JPRB,0.05131470_JPRB,&
88    &0.03487040_JPRB,0.00327281_JPRB,0.00250183_JPRB,0.00190024_JPRB,&
89    &0.00133978_JPRB,0.00082214_JPRB,0.00031016_JPRB,0.00004363_JPRB/)
90FRACREFAO( : ,3) = (/&
91    &0.13689300_JPRB,0.16610400_JPRB,0.15723500_JPRB,0.14299500_JPRB,&
92    &0.12399400_JPRB,0.09907820_JPRB,0.07169690_JPRB,0.05367370_JPRB,&
93    &0.03671630_JPRB,0.00378148_JPRB,0.00290510_JPRB,0.00221076_JPRB,&
94    &0.00142810_JPRB,0.00093527_JPRB,0.00031016_JPRB,0.00004363_JPRB/)
95FRACREFAO( : ,4) = (/&
96    &0.13054299_JPRB,0.16273800_JPRB,0.15874299_JPRB,0.14279599_JPRB,&
97    &0.12674300_JPRB,0.09664900_JPRB,0.07462200_JPRB,0.05620080_JPRB,&
98    &0.03789090_JPRB,0.00411690_JPRB,0.00322920_JPRB,0.00245036_JPRB,&
99    &0.00178303_JPRB,0.00098595_JPRB,0.00040802_JPRB,0.00010150_JPRB/)
100FRACREFAO( : ,5) = (/&
101    &0.12828299_JPRB,0.15824600_JPRB,0.15688400_JPRB,0.14449100_JPRB,&
102    &0.12787800_JPRB,0.09517830_JPRB,0.07679350_JPRB,0.05890820_JPRB,&
103    &0.03883570_JPRB,0.00442304_JPRB,0.00346796_JPRB,0.00255333_JPRB,&
104    &0.00212519_JPRB,0.00116168_JPRB,0.00067065_JPRB,0.00010150_JPRB/)
105FRACREFAO( : ,6) = (/&
106    &0.12649800_JPRB,0.15195100_JPRB,0.15646499_JPRB,0.14569700_JPRB,&
107    &0.12669300_JPRB,0.09653520_JPRB,0.07887920_JPRB,0.06106920_JPRB,&
108    &0.04043910_JPRB,0.00430390_JPRB,0.00364453_JPRB,0.00314360_JPRB,&
109    &0.00203206_JPRB,0.00187787_JPRB,0.00067075_JPRB,0.00010150_JPRB/)
110FRACREFAO( : ,7) = (/&
111    &0.12500300_JPRB,0.14460599_JPRB,0.15672199_JPRB,0.14724600_JPRB,&
112    &0.11978900_JPRB,0.10190200_JPRB,0.08196710_JPRB,0.06315770_JPRB,&
113    &0.04240100_JPRB,0.00433645_JPRB,0.00404097_JPRB,0.00329466_JPRB,&
114    &0.00288491_JPRB,0.00187803_JPRB,0.00067093_JPRB,0.00010150_JPRB/)
115FRACREFAO( : ,8) = (/&
116    &0.12317200_JPRB,0.14118700_JPRB,0.15242000_JPRB,0.13794300_JPRB,&
117    &0.12119200_JPRB,0.10655400_JPRB,0.08808350_JPRB,0.06521370_JPRB,&
118    &0.04505680_JPRB,0.00485949_JPRB,0.00477105_JPRB,0.00401468_JPRB,&
119    &0.00288491_JPRB,0.00187786_JPRB,0.00067110_JPRB,0.00010150_JPRB/)
120FRACREFAO( : ,9) = (/&
121    &0.10193600_JPRB,0.11693000_JPRB,0.13236099_JPRB,0.14053200_JPRB,&
122    &0.13749801_JPRB,0.12193100_JPRB,0.10221000_JPRB,0.07448910_JPRB,&
123    &0.05205320_JPRB,0.00572312_JPRB,0.00476882_JPRB,0.00403380_JPRB,&
124    &0.00288871_JPRB,0.00187396_JPRB,0.00067218_JPRB,0.00010150_JPRB/)
125
126!     -----------------------------------------------------------------
127RETURN
128END SUBROUTINE RRTM_KGB12_00
Note: See TracBrowser for help on using the repository browser.