LMDZ
rrtm_kgb3_00.F90
Go to the documentation of this file.
1 !OCL SCALAR
2 SUBROUTINE rrtm_kgb3_00
3 
4 ! Originally by Eli J. Mlawer, Atmospheric & Environmental Research.
5 ! BAND 3: 500-630 cm-1 (low - H2O,CO2; high - H2O,CO2)
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 
13 USE yoerrto3 , ONLY : kao ,kbo ,selfrefo ,fracrefao ,&
14  &fracrefbo ,forrefo ,absn2oao ,absn2obo
15 USE yoerrta3 , ONLY : absn2oa ,absn2ob ,etaref ,h2oref ,&
16  &n2oref ,co2ref ,strrat
17 
18 ! ------------------------------------------------------------------
19 
20 
21 IMPLICIT NONE
22 forrefo( :) = (/&
23  &1.76842e-04_jprb, 1.77913e-04_jprb, 1.25186e-04_jprb, 1.07912e-04_jprb,&
24  &1.05217e-04_jprb, 7.48726e-05_jprb, 1.11701e-04_jprb, 7.68921e-05_jprb,&
25  &9.87242e-05_jprb, 9.85711e-05_jprb, 6.16557e-05_jprb,-1.61291e-05_jprb,&
26  &-1.26794e-04_jprb,-1.19011e-04_jprb,-2.67814e-04_jprb, 6.95005e-05_jprb/)
27 
28 absn2oao( :) = (/&
29  &1.50387e-01_jprb,2.91407e-01_jprb,6.28803e-01_jprb,9.65619e-01_jprb,&
30  &1.15054e-00_jprb,2.23424e-00_jprb,1.83392e-00_jprb,1.39033e-00_jprb,&
31  &4.28457e-01_jprb,2.73502e-01_jprb,1.84307e-01_jprb,1.61325e-01_jprb,&
32  &7.66314e-02_jprb,1.33862e-01_jprb,6.71196e-07_jprb,1.59293e-06_jprb/)
33 
34 absn2obo( :) = (/&
35  &9.37044e-05_jprb,1.23318e-03_jprb,7.91720e-03_jprb,5.33005e-02_jprb,&
36  &1.72343e-01_jprb,4.29571e-01_jprb,1.01288e+00_jprb,3.83863e+00_jprb,&
37  &1.15312e+01_jprb,1.08383e+00_jprb,2.24847e+00_jprb,1.51268e+00_jprb,&
38  &3.33177e-01_jprb,7.82102e-01_jprb,3.44631e-01_jprb,1.61039e-03_jprb/)
39 
40 etaref( :) = (/&
41  &_zero_,0.125_jprb,0.25_jprb,0.375_jprb,_half_,0.625_jprb,0.75_jprb,0.875_jprb,0.9875_jprb,_one_/)
42 
43 h2oref( :) = (/&
44  &1.87599e-02_jprb,1.22233e-02_jprb,5.89086e-03_jprb,2.76753e-03_jprb,1.40651e-03_jprb, &
45  &7.59698e-04_jprb,3.88758e-04_jprb,1.65422e-04_jprb,3.71895e-05_jprb,7.47648e-06_jprb, &
46  &4.30818e-06_jprb,3.33194e-06_jprb,3.20393e-06_jprb,3.16186e-06_jprb,3.25235e-06_jprb, &
47  &3.42258e-06_jprb,3.62884e-06_jprb,3.91482e-06_jprb,4.14875e-06_jprb,4.30810e-06_jprb,&
48  &4.44204e-06_jprb,4.57783e-06_jprb,4.70865e-06_jprb,4.79432e-06_jprb,4.86971e-06_jprb, &
49  &4.92603e-06_jprb,4.96688e-06_jprb,4.99628e-06_jprb,5.05266e-06_jprb,5.12658e-06_jprb, &
50  &5.25028e-06_jprb,5.35708e-06_jprb,5.45085e-06_jprb,5.48304e-06_jprb,5.50000e-06_jprb, &
51  &5.50000e-06_jprb,5.45359e-06_jprb,5.40468e-06_jprb,5.35576e-06_jprb,5.25327e-06_jprb,&
52  &5.14362e-06_jprb,5.03396e-06_jprb,4.87662e-06_jprb,4.69787e-06_jprb,4.51911e-06_jprb, &
53  &4.33600e-06_jprb,4.14416e-06_jprb,3.95232e-06_jprb,3.76048e-06_jprb,3.57217e-06_jprb, &
54  &3.38549e-06_jprb,3.19881e-06_jprb,3.01212e-06_jprb,2.82621e-06_jprb,2.64068e-06_jprb, &
55  &2.45515e-06_jprb,2.26962e-06_jprb,2.08659e-06_jprb,1.93029e-06_jprb/)
56 
57 n2oref( :) = (/&
58  &3.20000e-07_jprb,3.20000e-07_jprb,3.20000e-07_jprb,3.20000e-07_jprb,3.20000e-07_jprb,&
59  &3.19652e-07_jprb,3.15324e-07_jprb,3.03830e-07_jprb,2.94221e-07_jprb,2.84953e-07_jprb,&
60  &2.76714e-07_jprb,2.64709e-07_jprb,2.42847e-07_jprb,2.09547e-07_jprb,1.71945e-07_jprb,&
61  &1.37491e-07_jprb,1.13319e-07_jprb,1.00354e-07_jprb,9.12812e-08_jprb,8.54633e-08_jprb,&
62  &8.03631e-08_jprb,7.33718e-08_jprb,6.59754e-08_jprb,5.60386e-08_jprb,4.70901e-08_jprb,&
63  &3.99774e-08_jprb,3.29786e-08_jprb,2.60642e-08_jprb,2.10663e-08_jprb,1.65918e-08_jprb,&
64  &1.30167e-08_jprb,1.00900e-08_jprb,7.62490e-09_jprb,6.11592e-09_jprb,4.66725e-09_jprb,&
65  &3.28574e-09_jprb,2.84838e-09_jprb,2.46198e-09_jprb,2.07557e-09_jprb,1.85507e-09_jprb,&
66  &1.65675e-09_jprb,1.45843e-09_jprb,1.31948e-09_jprb,1.20716e-09_jprb,1.09485e-09_jprb,&
67  &9.97803e-10_jprb,9.31260e-10_jprb,8.64721e-10_jprb,7.98181e-10_jprb,7.51380e-10_jprb,&
68  &7.13670e-10_jprb,6.75960e-10_jprb,6.38250e-10_jprb,6.09811e-10_jprb,5.85998e-10_jprb,&
69  &5.62185e-10_jprb,5.38371e-10_jprb,5.15183e-10_jprb,4.98660e-10_jprb/)
70 
71 co2ref( :) = (/&
72  &3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,&
73  &3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,&
74  &3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,&
75  &3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,&
76  &3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,&
77  &3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,&
78  &3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,&
79  &3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,&
80  &3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,&
81  &3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,&
82  &3.55000e-04_jprb,3.55000e-04_jprb,3.55000e-04_jprb,&
83  &3.5470873e-04_jprb, 3.5427220e-04_jprb, 3.5383567e-04_jprb,&
84  &3.5339911e-04_jprb, 3.5282588e-04_jprb, 3.5079606e-04_jprb/)
85 
86 strrat = 1.19268_jprb
87 
88 ! ------------------------------------------------------------------
89 
90 ! The array SELFREFO contains the coefficient of the water vapor
91 ! self-continuum (including the energy term). The first index
92 ! refers to temperature in 7.2 degree increments. For instance,
93 ! JT = 1 refers to a temperature of 245.6, JT = 2 refers to 252.8,
94 ! etc. The second index runs over the g-channel (1 to 16).
95 
96 selfrefo( :, 1) = (/&
97 &5.10286e-01_jprb, 4.30424e-01_jprb, 3.63061e-01_jprb, 3.06241e-01_jprb, 2.58313e-01_jprb,&
98 &2.17886e-01_jprb, 1.83786e-01_jprb, 1.55023e-01_jprb, 1.30761e-01_jprb, 1.10297e-01_jprb/)
99 selfrefo( :, 2) = (/&
100 &4.96777e-01_jprb, 4.30758e-01_jprb, 3.73512e-01_jprb, 3.23874e-01_jprb, 2.80833e-01_jprb,&
101 &2.43512e-01_jprb, 2.11150e-01_jprb, 1.83089e-01_jprb, 1.58758e-01_jprb, 1.37659e-01_jprb/)
102 selfrefo( :, 3) = (/&
103 &4.65595e-01_jprb, 4.04778e-01_jprb, 3.51905e-01_jprb, 3.05939e-01_jprb, 2.65976e-01_jprb,&
104 &2.31234e-01_jprb, 2.01030e-01_jprb, 1.74771e-01_jprb, 1.51942e-01_jprb, 1.32095e-01_jprb/)
105 selfrefo( :, 4) = (/&
106 &4.34629e-01_jprb, 3.86709e-01_jprb, 3.44073e-01_jprb, 3.06137e-01_jprb, 2.72384e-01_jprb,&
107 &2.42353e-01_jprb, 2.15633e-01_jprb, 1.91858e-01_jprb, 1.70705e-01_jprb, 1.51884e-01_jprb/)
108 selfrefo( :, 5) = (/&
109 &4.17625e-01_jprb, 3.75442e-01_jprb, 3.37519e-01_jprb, 3.03427e-01_jprb, 2.72778e-01_jprb,&
110 &2.45225e-01_jprb, 2.20456e-01_jprb, 1.98188e-01_jprb, 1.78169e-01_jprb, 1.60173e-01_jprb/)
111 selfrefo( :, 6) = (/&
112 &3.97412e-01_jprb, 3.61249e-01_jprb, 3.28376e-01_jprb, 2.98494e-01_jprb, 2.71332e-01_jprb,&
113 &2.46641e-01_jprb, 2.24198e-01_jprb, 2.03796e-01_jprb, 1.85251e-01_jprb, 1.68394e-01_jprb/)
114 selfrefo( :, 7) = (/&
115 &3.89795e-01_jprb, 3.56178e-01_jprb, 3.25460e-01_jprb, 2.97391e-01_jprb, 2.71743e-01_jprb,&
116 &2.48307e-01_jprb, 2.26892e-01_jprb, 2.07324e-01_jprb, 1.89444e-01_jprb, 1.73105e-01_jprb/)
117 selfrefo( :, 8) = (/&
118 &3.79095e-01_jprb, 3.46828e-01_jprb, 3.17308e-01_jprb, 2.90301e-01_jprb, 2.65592e-01_jprb,&
119 &2.42986e-01_jprb, 2.22305e-01_jprb, 2.03384e-01_jprb, 1.86073e-01_jprb, 1.70235e-01_jprb/)
120 selfrefo( :, 9) = (/&
121 &3.87086e-01_jprb, 3.53988e-01_jprb, 3.23720e-01_jprb, 2.96040e-01_jprb, 2.70727e-01_jprb,&
122 &2.47578e-01_jprb, 2.26409e-01_jprb, 2.07050e-01_jprb, 1.89346e-01_jprb, 1.73156e-01_jprb/)
123 selfrefo( :,10) = (/&
124 &4.07529e-01_jprb, 3.71284e-01_jprb, 3.38263e-01_jprb, 3.08179e-01_jprb, 2.80770e-01_jprb,&
125 &2.55799e-01_jprb, 2.33049e-01_jprb, 2.12322e-01_jprb, 1.93438e-01_jprb, 1.76234e-01_jprb/)
126 selfrefo( :,11) = (/&
127 &4.09315e-01_jprb, 3.74420e-01_jprb, 3.42500e-01_jprb, 3.13302e-01_jprb, 2.86592e-01_jprb,&
128 &2.62160e-01_jprb, 2.39811e-01_jprb, 2.19367e-01_jprb, 2.00665e-01_jprb, 1.83558e-01_jprb/)
129 selfrefo( :,12) = (/&
130 &3.73996e-01_jprb, 3.51462e-01_jprb, 3.30286e-01_jprb, 3.10385e-01_jprb, 2.91684e-01_jprb,&
131 &2.74109e-01_jprb, 2.57593e-01_jprb, 2.42072e-01_jprb, 2.27487e-01_jprb, 2.13780e-01_jprb/)
132 selfrefo( :,13) = (/&
133 &4.87965e-01_jprb, 4.44525e-01_jprb, 4.04952e-01_jprb, 3.68902e-01_jprb, 3.36061e-01_jprb,&
134 &3.06144e-01_jprb, 2.78891e-01_jprb, 2.54063e-01_jprb, 2.31446e-01_jprb, 2.10842e-01_jprb/)
135 selfrefo( :,14) = (/&
136 &5.67999e-01_jprb, 5.03336e-01_jprb, 4.46034e-01_jprb, 3.95256e-01_jprb, 3.50259e-01_jprb,&
137 &3.10384e-01_jprb, 2.75048e-01_jprb, 2.43736e-01_jprb, 2.15988e-01_jprb, 1.91399e-01_jprb/)
138 selfrefo( :,15) = (/&
139 &5.92492e-01_jprb, 5.27072e-01_jprb, 4.68874e-01_jprb, 4.17103e-01_jprb, 3.71048e-01_jprb,&
140 &3.30079e-01_jprb, 2.93633e-01_jprb, 2.61211e-01_jprb, 2.32369e-01_jprb, 2.06712e-01_jprb/)
141 selfrefo( :,16) = (/&
142 &5.04176e-01_jprb, 4.64845e-01_jprb, 4.28582e-01_jprb, 3.95148e-01_jprb, 3.64322e-01_jprb,&
143 &3.35901e-01_jprb, 3.09697e-01_jprb, 2.85538e-01_jprb, 2.63263e-01_jprb, 2.42725e-01_jprb/)
144 
145 fracrefao( :, 1) = (/&
146 ! From P = 1053.6 mb.
147  &0.15116400_jprb,0.14875700_jprb,0.14232300_jprb,0.13234501_jprb,&
148  &0.11881600_jprb,0.10224100_jprb,0.08345580_jprb,0.06267490_jprb,&
149  &0.04250650_jprb,0.00462650_jprb,0.00382259_jprb,0.00302600_jprb,&
150  &0.00222004_jprb,0.00141397_jprb,0.00053379_jprb,0.00007421_jprb/)
151 fracrefao( :, 2) = (/&
152  &0.15266000_jprb,0.14888400_jprb,0.14195900_jprb,0.13179500_jprb,&
153  &0.11842700_jprb,0.10209000_jprb,0.08336130_jprb,0.06264370_jprb,&
154  &0.04247660_jprb,0.00461946_jprb,0.00381536_jprb,0.00302601_jprb,&
155  &0.00222004_jprb,0.00141397_jprb,0.00053302_jprb,0.00007498_jprb/)
156 fracrefao( :, 3) = (/&
157  &0.15282799_jprb,0.14903000_jprb,0.14192399_jprb,0.13174300_jprb,&
158  &0.11835300_jprb,0.10202700_jprb,0.08329830_jprb,0.06264830_jprb,&
159  &0.04246910_jprb,0.00460242_jprb,0.00381904_jprb,0.00301573_jprb,&
160  &0.00222004_jprb,0.00141397_jprb,0.00053379_jprb,0.00007421_jprb/)
161 fracrefao( :, 4) = (/&
162  &0.15298399_jprb,0.14902800_jprb,0.14193401_jprb,0.13173500_jprb,&
163  &0.11833300_jprb,0.10195800_jprb,0.08324730_jprb,0.06264770_jprb,&
164  &0.04246490_jprb,0.00460489_jprb,0.00381123_jprb,0.00301893_jprb,&
165  &0.00221093_jprb,0.00141397_jprb,0.00053379_jprb,0.00007421_jprb/)
166 fracrefao( :, 5) = (/&
167  &0.15307599_jprb,0.14907201_jprb,0.14198899_jprb,0.13169800_jprb,&
168  &0.11827300_jprb,0.10192300_jprb,0.08321600_jprb,0.06263490_jprb,&
169  &0.04245600_jprb,0.00460846_jprb,0.00380836_jprb,0.00301663_jprb,&
170  &0.00221402_jprb,0.00141167_jprb,0.00052807_jprb,0.00007376_jprb/)
171 fracrefao( :, 6) = (/&
172  &0.15311401_jprb,0.14915401_jprb,0.14207301_jprb,0.13167299_jprb,&
173  &0.11819300_jprb,0.10188900_jprb,0.08318760_jprb,0.06261960_jprb,&
174  &0.04243890_jprb,0.00461584_jprb,0.00380929_jprb,0.00300815_jprb,&
175  &0.00221736_jprb,0.00140588_jprb,0.00052776_jprb,0.00007376_jprb/)
176 fracrefao( :, 7) = (/&
177  &0.15316001_jprb,0.14925499_jprb,0.14213000_jprb,0.13170999_jprb,&
178  &0.11807700_jprb,0.10181400_jprb,0.08317400_jprb,0.06260300_jprb,&
179  &0.04242720_jprb,0.00461520_jprb,0.00381381_jprb,0.00301285_jprb,&
180  &0.00220275_jprb,0.00140371_jprb,0.00052776_jprb,0.00007376_jprb/)
181 fracrefao( :, 8) = (/&
182  &0.15321200_jprb,0.14940999_jprb,0.14222500_jprb,0.13164200_jprb,&
183  &0.11798200_jprb,0.10174500_jprb,0.08317500_jprb,0.06253640_jprb,&
184  &0.04243130_jprb,0.00461724_jprb,0.00381534_jprb,0.00300320_jprb,&
185  &0.00220091_jprb,0.00140364_jprb,0.00052852_jprb,0.00007300_jprb/)
186 fracrefao( :, 9) = (/&
187  &0.15312800_jprb,0.14973100_jprb,0.14234400_jprb,0.13168900_jprb,&
188  &0.11795200_jprb,0.10156100_jprb,0.08302990_jprb,0.06252240_jprb,&
189  &0.04240980_jprb,0.00461035_jprb,0.00381381_jprb,0.00300176_jprb,&
190  &0.00220160_jprb,0.00140284_jprb,0.00052774_jprb,0.00007376_jprb/)
191 fracrefao( :,10) = (/&
192  &0.15292500_jprb,0.14978001_jprb,0.14242400_jprb,0.13172600_jprb,&
193  &0.11798800_jprb,0.10156400_jprb,0.08303050_jprb,0.06251670_jprb,&
194  &0.04240970_jprb,0.00461302_jprb,0.00381452_jprb,0.00300250_jprb,&
195  &0.00220126_jprb,0.00140324_jprb,0.00052850_jprb,0.00007300_jprb/)
196 
197 fracrefbo( :, 1) = (/&
198 ! From P = 64.1 mb.
199  &0.16340201_jprb,0.15607700_jprb,0.14601400_jprb,0.13182700_jprb,&
200  &0.11524700_jprb,0.09666570_jprb,0.07825360_jprb,0.05849780_jprb,&
201  &0.03949650_jprb,0.00427980_jprb,0.00353719_jprb,0.00279303_jprb,&
202  &0.00204788_jprb,0.00130139_jprb,0.00049055_jprb,0.00006904_jprb/)
203 fracrefbo( :, 2) = (/&
204  &0.15762900_jprb,0.15494700_jprb,0.14659800_jprb,0.13267800_jprb,&
205  &0.11562700_jprb,0.09838360_jprb,0.07930420_jprb,0.05962700_jprb,&
206  &0.04036360_jprb,0.00438053_jprb,0.00361463_jprb,0.00285723_jprb,&
207  &0.00208345_jprb,0.00132135_jprb,0.00050528_jprb,0.00008003_jprb/)
208 fracrefbo( :, 3) = (/&
209  &0.15641500_jprb,0.15394500_jprb,0.14633600_jprb,0.13180400_jprb,&
210  &0.11617100_jprb,0.09924170_jprb,0.08000510_jprb,0.06021420_jprb,&
211  &0.04082730_jprb,0.00441694_jprb,0.00365364_jprb,0.00287723_jprb,&
212  &0.00210914_jprb,0.00135784_jprb,0.00054651_jprb,0.00008003_jprb/)
213 fracrefbo( :, 4) = (/&
214  &0.15482700_jprb,0.15286300_jprb,0.14392500_jprb,0.13244100_jprb,&
215  &0.11712000_jprb,0.09994920_jprb,0.08119200_jprb,0.06104360_jprb,&
216  &0.04135600_jprb,0.00446685_jprb,0.00368377_jprb,0.00290767_jprb,&
217  &0.00215445_jprb,0.00142865_jprb,0.00056142_jprb,0.00008003_jprb/)
218 fracrefbo( :, 5) = (/&
219  &0.15975100_jprb,0.15653500_jprb,0.14214399_jprb,0.12892200_jprb,&
220  &0.11508400_jprb,0.09906020_jprb,0.08087940_jprb,0.06078190_jprb,&
221  &0.04140530_jprb,0.00452724_jprb,0.00374558_jprb,0.00295328_jprb,&
222  &0.00218509_jprb,0.00138644_jprb,0.00056018_jprb,0.00008003_jprb/)
223 
224 ! -----------------------------------------------------------------
225 RETURN
226 END SUBROUTINE rrtm_kgb3_00
yoerrto3::forrefo
real(kind=jprb), dimension(no3) forrefo
Definition: yoerrto3.F90:18
rrtm_kgb3_00
subroutine rrtm_kgb3_00
Definition: rrtm_kgb3_00.F90:3
yoerrta3::absn2ob
real(kind=jprb), dimension(16) absn2ob
Definition: yoerrta3.F90:20
yoerrto3::fracrefbo
real(kind=jprb), dimension(no3, 5) fracrefbo
Definition: yoerrto3.F90:16
yoerrta3::h2oref
real(kind=jprb), dimension(59) h2oref
Definition: yoerrta3.F90:22
yoerrta3::co2ref
real(kind=jprb), dimension(59) co2ref
Definition: yoerrta3.F90:24
yoerrta3::n2oref
real(kind=jprb), dimension(59) n2oref
Definition: yoerrta3.F90:23
yoerrto3::kao
real(kind=jprb), dimension(10, 5, 13, no3) kao
Definition: yoerrto3.F90:22
yoerrta3::absn2oa
real(kind=jprb), dimension(16) absn2oa
Definition: yoerrta3.F90:19
yoerrto3::absn2obo
real(kind=jprb), dimension(no3) absn2obo
Definition: yoerrto3.F90:20
yoerrta3::etaref
real(kind=jprb), dimension(10) etaref
Definition: yoerrta3.F90:21
yoerrto3::kbo
real(kind=jprb), dimension(5, 5, 13:59, no3) kbo
Definition: yoerrto3.F90:23
yoerrta3::strrat
real(kind=jprb) strrat
Definition: yoerrta3.F90:29
yoerrto3::selfrefo
real(kind=jprb), dimension(10, no3) selfrefo
Definition: yoerrto3.F90:24
yoerrto3::fracrefao
real(kind=jprb), dimension(no3, 10) fracrefao
Definition: yoerrto3.F90:16
yoerrto3::absn2oao
real(kind=jprb), dimension(no3) absn2oao
Definition: yoerrto3.F90:19