LMDZ
suolwn.F90
Go to the documentation of this file.
1 SUBROUTINE suolw
2 !
3 !**** *SULW* - INITIALIZE COMMON YOELW
4 !
5 ! PURPOSE.
6 ! --------
7 ! INITIALIZE YOELW, THE COMMON THAT CONTAINS COEFFICIENTS
8 ! NEEDED TO RUN THE ECMWF LONGWAVE RADIATION ROUTINES
9 !
10 !** INTERFACE.
11 ! ----------
12 ! *CALL* *SULW
13 !
14 ! EXPLICIT ARGUMENTS :
15 ! --------------------
16 ! NONE
17 !
18 ! IMPLICIT ARGUMENTS :
19 ! --------------------
20 ! COMMON YOELW
21 !
22 ! METHOD.
23 ! -------
24 ! SEE DOCUMENTATION
25 !
26 ! EXTERNALS.
27 ! ----------
28 !
29 ! REFERENCE.
30 ! ----------
31 ! ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS
32 !
33 ! AUTHOR.
34 ! -------
35 ! JEAN-JACQUES MORCRETTE *ECMWF*
36 !
37 ! MODIFICATIONS.
38 ! --------------
39 ! ORIGINAL : 88-12-15
40 ! J.-J. Morcrette ECMWF 951215: Pseudo-Voigt Line Profiles
41 ! ------------------------------------------------------------------
42 ! IMPLICIT LOGICAL (L)
43 !
44 !#include "yomlun.h"
45 !#include "yoerad.h"
46 !#include "yoelw.h"
47 
48 #include "tsmbkind.h"
49 
50 USE yoeolw , ONLY : mxixt ,nisp ,nipd ,nipd2 ,&
51  &ntr ,ntra ,nua ,ng1 ,ng1p1 ,&
52  &at ,bt ,ga ,gb ,&
53  &o1h ,o2h ,oct ,rpialf0 ,&
54  &rt1 ,tintp ,tref ,tstand ,tstp ,&
55  &wg1 ,xp ,rvgco2 ,rvgh2o ,rvgo3
56 
57 
58 IMPLICIT NONE
59 
60 
61 ! ----------------------------------------------------------------
62 !
63 !* 1. SET VALUES.
64 ! -----------
65 !
66 ng1=2
67 ng1p1=ng1+1
68 nisp=6
69 nipd=8
70 nipd2=2*nipd
71 ntr=11
72 ntra=15
73 nua=24
74 !
75 !
76 !-- COEFFICIENTS FOR THE PSEUDO-VOIGT LINE PROFILE (PA)
77 !
78 rvgco2= 60.
79 rvgh2o= 30.
80 rvgo3 =400.
81 !
82 !
83 !-- REFERENCE TEMPERATURES FOR THE PLANCK FUNCTIONS ----- TSTAND
84 ! TEMPERATURE DEPENDENCE OF THE ABSORPTION ----- TREF
85 tstand= 250.00
86 tref = 250.00
87 !
88 
89 !-- ROOTS AND WEIGHTS FOR THE 2-POINT GAUSSIAN QUADRATURE
90 rt1 = (/ -0.577350269, +0.577350269 /)
91 wg1 = (/ 1. , 1. /)
92 
93 !-- TEMPERATURE STEP AND MAXIMUM INDEX FOR THE TRANSMISSIVITIES
94 tstp = 12.5
95 mxixt = 10
96 
97 !
98 !-- REFERENCE TEMPERATURE FOR THE TRANSMISSIVITIES
99 
100 tintp = (/ 187.5, 200., 212.5, 225., 237.5, &
101  &250.,&
102  &262.5, 275., 287.5, 300., 312.5 /)
103 !
104 !-- COEFFICIENTS OF THE POLYNOMIALS GIVING THE PLANCK FUNCTIONS
105 !
106 xp( :, 1) = (/&
107  &0.46430621e+02, 0.12928299e+03, 0.20732648e+03,&
108  &0.31398411e+03, 0.18373177e+03,-0.11412303e+03/)
109 xp( :, 2) = (/&
110  &0.73604774e+02, 0.27887914e+03, 0.27076947e+03,&
111  &-.57322111e+02,-0.64742459e+02, 0.87238280e+02/)
112 xp( :, 3) = (/&
113  &0.37050866e+02, 0.20498759e+03, 0.37558029e+03,&
114  &0.17401171e+03,-0.13350302e+03,-0.37651795e+02/)
115 xp( :, 4) = (/&
116  &0.14930141e+02, 0.89161160e+02, 0.17793062e+03,&
117  &0.93433860e+02,-0.70646020e+02,-0.26373150e+02/)
118 xp( :, 5) = (/&
119  &0.40386780e+02, 0.10855270e+03, 0.50755010e+02,&
120 &-0.31496190e+02, 0.12791300e+00, 0.18017770e+01/)
121 xp( :, 6) = (/&
122  &0.90811926e+01, 0.75073923e+02, 0.24654438e+03,&
123  &0.39332612e+03, 0.29385281e+03, 0.89107921e+02/)
124 
125 !
126 !-- COEFFICIENTS OF THE POLYNOMIALS GIVING THE TEMPERATURE DEPENDENCE OF
127 ! THE ABSORPTION
128 !
129 !- from Zhong & Haigh
130 !
131 !-- H2O -- 0 - 350
132 at(1,:)= (/ 0.298199e-02,-.394023e-03,0.319566e-04 /)
133 bt(1,:)= (/-0.106432e-04,0.660324e-06,0.174356e-06 /)
134 !
135 !-- H2O -- 500 - 800
136 at(2,:)= (/ 0.143676e-01,0.366501e-02,-.160822e-02 /)
137 bt(2,:)= (/-0.553979e-04,-.101701e-04,0.920868e-05 /)
138 !
139 !-- H2O -- 800 - 970 + 1110 - 1250
140 at(3,:)= (/ 0.197861e-01,0.315541e-02,-.174547e-02 /)
141 bt(3,:)= (/-0.877012e-04,0.513302e-04,0.523138e-06 /)
142 !
143 !-- H2O -- 970 - 1110
144 at(4,:)= (/ 0.289560e-01,-.208807e-02,-.121943e-02 /)
145 bt(4,:)= (/-0.165960e-03,0.157704e-03,-.146427e-04 /)
146 !
147 !-- H2O -- 350 - 500
148 at(5,:)= (/ 0.103800e-01,0.436296e-02,-.161431e-02 /)
149 bt(5,:)= (/ -.276744e-04,-.327381e-04,0.127646e-04 /)
150 !
151 !-- H2O -- 1250 - 1450 + 1880 - 2820
152 at(6,:)= (/ 0.868859e-02,-.972752e-03,0.000000e-00 /)
153 bt(6,:)= (/ -.278412e-04,-.713940e-06,0.117469e-05 /)
154 !
155 
156 !-------CO2
157 
158 at(7,:)= (/&
159  & 0.250073e-03,0.455875e-03,0.109242e-03 /)
160 bt(7,:)= (/&
161  & 0.199846e-05,-.216313e-05,0.175991e-06 /)
162 at(8,:)= (/&
163  & 0.307423e-01,0.110879e-02,-.322172e-03 /)
164 bt(8,:)= (/&
165  &-0.108482e-03,0.258096e-05,-.814575e-06 /)
166 !
167 !-- COEFFICIENTS FOR OZONE ABSORPTION
168 !
169 o1h = 2230.
170 o2h = 100.
171 rpialf0= 2.
172 oct(1)= -.326e-03
173 oct(2)= -.102e-05
174 oct(3)= .137e-02
175 oct(4)= -.535e-05
176 !
177 !-- COEFFICIENTS FOR H2O E-TYPE CONTINUUM ABSORPTION
178 !
179 ! DATA RNTNU / 70.0, 26.77, 8.45, 5.78 /
180 !
181 !-----------------------------------------------------------------------
182 ! PADE APPROXIMANTS FOR THE TRANSMISSIVITIES
183 ! ==========================================
184 !-----------------------------------------------------------------------
185 !
186 !
187 !
188 !
189 !-----------------------------------------------------------------------
190 ! PADE APPROXIMANTS FOR THE TRANSMISSIVITIES
191 !
192 !- H2O 0-350 cm-1
193 !
194 !-- R.D. -- G = - 0.2 SLA
195 !
196 !
197 !----- INTERVAL = 1 ----- T = 187.5
198 !
199 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
200 ga( 1, 1, :) = (/&
201  & 0.63499072e-02,-0.99506586e-03, 0.00000000e+00/)
202 gb( 1, 1, :) = (/&
203  & 0.63499072e-02, 0.97222852e-01, 0.10000000e+01/)
204 ga( 1, 2, :) = (/&
205  & 0.77266491e-02,-0.11661515e-02, 0.00000000e+00/)
206 gb( 1, 2, :) = (/&
207  & 0.77266491e-02, 0.10681591e+00, 0.10000000e+01/)
208 !
209 !----- INTERVAL = 1 ----- T = 200.0
210 !
211 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
212 ga( 2, 1, :) = (/&
213  & 0.65566348e-02,-0.10184169e-02, 0.00000000e+00/)
214 gb( 2, 1, :) = (/&
215  & 0.65566348e-02, 0.98862238e-01, 0.10000000e+01/)
216 ga( 2, 2, :) = (/&
217  & 0.81323287e-02,-0.11886130e-02, 0.00000000e+00/)
218 gb( 2, 2, :) = (/&
219  & 0.81323287e-02, 0.10921298e+00, 0.10000000e+01/)
220 !
221 !----- INTERVAL = 1 ----- T = 212.5
222 !
223 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
224 ga( 3, 1, :) = (/&
225  & 0.67849730e-02,-0.10404730e-02, 0.00000000e+00/)
226 gb( 3, 1, :) = (/&
227  & 0.67849730e-02, 0.10061504e+00, 0.10000000e+01/)
228 ga( 3, 2, :) = (/&
229  & 0.86507620e-02,-0.12139929e-02, 0.00000000e+00/)
230 gb( 3, 2, :) = (/&
231  & 0.86507620e-02, 0.11198225e+00, 0.10000000e+01/)
232 !
233 !----- INTERVAL = 1 ----- T = 225.0
234 !
235 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
236 ga( 4, 1, :) = (/&
237  & 0.70481947e-02,-0.10621792e-02, 0.00000000e+00/)
238 gb( 4, 1, :) = (/&
239  & 0.70481947e-02, 0.10256222e+00, 0.10000000e+01/)
240 ga( 4, 2, :) = (/&
241  & 0.92776391e-02,-0.12445811e-02, 0.00000000e+00/)
242 gb( 4, 2, :) = (/&
243  & 0.92776391e-02, 0.11487826e+00, 0.10000000e+01/)
244 !
245 !----- INTERVAL = 1 ----- T = 237.5
246 !
247 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
248 ga( 5, 1, :) = (/&
249  & 0.73585943e-02,-0.10847662e-02, 0.00000000e+00/)
250 gb( 5, 1, :) = (/&
251  & 0.73585943e-02, 0.10475952e+00, 0.10000000e+01/)
252 ga( 5, 2, :) = (/&
253  & 0.99806312e-02,-0.12807672e-02, 0.00000000e+00/)
254 gb( 5, 2, :) = (/&
255  & 0.99806312e-02, 0.11751113e+00, 0.10000000e+01/)
256 !
257 !----- INTERVAL = 1 ----- T = 250.0
258 !
259 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
260 ga( 6, 1, :) = (/&
261  & 0.77242818e-02,-0.11094726e-02, 0.00000000e+00/)
262 gb( 6, 1, :) = (/&
263  & 0.77242818e-02, 0.10720986e+00, 0.10000000e+01/)
264 ga( 6, 2, :) = (/&
265  & 0.10709803e-01,-0.13208251e-02, 0.00000000e+00/)
266 gb( 6, 2, :) = (/&
267  & 0.10709803e-01, 0.11951535e+00, 0.10000000e+01/)
268 !
269 !----- INTERVAL = 1 ----- T = 262.5
270 !
271 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
272 ga( 7, 1, :) = (/&
273  & 0.81472693e-02,-0.11372949e-02, 0.00000000e+00/)
274 gb( 7, 1, :) = (/&
275  & 0.81472693e-02, 0.10985370e+00, 0.10000000e+01/)
276 ga( 7, 2, :) = (/&
277  & 0.11414739e-01,-0.13619034e-02, 0.00000000e+00/)
278 gb( 7, 2, :) = (/&
279  & 0.11414739e-01, 0.12069945e+00, 0.10000000e+01/)
280 !
281 !----- INTERVAL = 1 ----- T = 275.0
282 !
283 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
284 ga( 8, 1, :) = (/&
285  & 0.86227527e-02,-0.11687683e-02, 0.00000000e+00/)
286 gb( 8, 1, :) = (/&
287  & 0.86227527e-02, 0.11257633e+00, 0.10000000e+01/)
288 ga( 8, 2, :) = (/&
289  & 0.12058772e-01,-0.14014165e-02, 0.00000000e+00/)
290 gb( 8, 2, :) = (/&
291  & 0.12058772e-01, 0.12108524e+00, 0.10000000e+01/)
292 !
293 !----- INTERVAL = 1 ----- T = 287.5
294 !
295 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
296 ga( 9, 1, :) = (/&
297  & 0.91396814e-02,-0.12038314e-02, 0.00000000e+00/)
298 gb( 9, 1, :) = (/&
299  & 0.91396814e-02, 0.11522980e+00, 0.10000000e+01/)
300 ga( 9, 2, :) = (/&
301  & 0.12623992e-01,-0.14378639e-02, 0.00000000e+00/)
302 gb( 9, 2, :) = (/&
303  & 0.12623992e-01, 0.12084229e+00, 0.10000000e+01/)
304 !
305 !----- INTERVAL = 1 ----- T = 300.0
306 !
307 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
308 ga(10, 1, :) = (/&
309  & 0.96825438e-02,-0.12418367e-02, 0.00000000e+00/)
310 gb(10, 1, :) = (/&
311  & 0.96825438e-02, 0.11766343e+00, 0.10000000e+01/)
312 ga(10, 2, :) = (/&
313  & 0.13108146e-01,-0.14708488e-02, 0.00000000e+00/)
314 gb(10, 2, :) = (/&
315  & 0.13108146e-01, 0.12019005e+00, 0.10000000e+01/)
316 !
317 !----- INTERVAL = 1 ----- T = 312.5
318 !
319 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
320 ga(11, 1, :) = (/&
321  & 0.10233955e-01,-0.12817135e-02, 0.00000000e+00/)
322 gb(11, 1, :) = (/&
323  & 0.10233955e-01, 0.11975320e+00, 0.10000000e+01/)
324 ga(11, 2, :) = (/&
325  & 0.13518390e-01,-0.15006791e-02, 0.00000000e+00/)
326 gb(11, 2, :) = (/&
327  & 0.13518390e-01, 0.11932684e+00, 0.10000000e+01/)
328 !
329 !!
330 !- H2O 500-800 cm-1
331 !
332 !--- R.D. --- G = 0.02 + 0.50 / ( 1 + 4.5 U )
333 !
334 !
335 !----- INTERVAL = 2 ----- T = 187.5
336 !
337 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
338 ga( 1, 3, :) = (/&
339  & 0.11644593e+01, 0.41243390e+00, 0.00000000e+00/)
340 gb( 1, 3, :) = (/&
341  & 0.11644593e+01, 0.10346097e+01, 0.10000000e+01/)
342 ga( 1, 4, :) = (/&
343  & 0.12006968e+01, 0.48318936e+00, 0.00000000e+00/)
344 gb( 1, 4, :) = (/&
345  & 0.12006968e+01, 0.10626130e+01, 0.10000000e+01/)
346 !
347 !----- INTERVAL = 2 ----- T = 200.0
348 !
349 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
350 ga( 2, 3, :) = (/&
351  & 0.11747203e+01, 0.43407282e+00, 0.00000000e+00/)
352 gb( 2, 3, :) = (/&
353  & 0.11747203e+01, 0.10433655e+01, 0.10000000e+01/)
354 ga( 2, 4, :) = (/&
355  & 0.12108196e+01, 0.50501827e+00, 0.00000000e+00/)
356 gb( 2, 4, :) = (/&
357  & 0.12108196e+01, 0.10716026e+01, 0.10000000e+01/)
358 !
359 !----- INTERVAL = 2 ----- T = 212.5
360 !
361 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
362 ga( 3, 3, :) = (/&
363  & 0.11837872e+01, 0.45331413e+00, 0.00000000e+00/)
364 gb( 3, 3, :) = (/&
365  & 0.11837872e+01, 0.10511933e+01, 0.10000000e+01/)
366 ga( 3, 4, :) = (/&
367  & 0.12196717e+01, 0.52409502e+00, 0.00000000e+00/)
368 gb( 3, 4, :) = (/&
369  & 0.12196717e+01, 0.10795108e+01, 0.10000000e+01/)
370 !
371 !----- INTERVAL = 2 ----- T = 225.0
372 !
373 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
374 ga( 4, 3, :) = (/&
375  & 0.11918561e+01, 0.47048604e+00, 0.00000000e+00/)
376 gb( 4, 3, :) = (/&
377  & 0.11918561e+01, 0.10582150e+01, 0.10000000e+01/)
378 ga( 4, 4, :) = (/&
379  & 0.12274493e+01, 0.54085277e+00, 0.00000000e+00/)
380 gb( 4, 4, :) = (/&
381  & 0.12274493e+01, 0.10865006e+01, 0.10000000e+01/)
382 !
383 !----- INTERVAL = 2 ----- T = 237.5
384 !
385 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
386 ga( 5, 3, :) = (/&
387  & 0.11990757e+01, 0.48586286e+00, 0.00000000e+00/)
388 gb( 5, 3, :) = (/&
389  & 0.11990757e+01, 0.10645317e+01, 0.10000000e+01/)
390 ga( 5, 4, :) = (/&
391  & 0.12343189e+01, 0.55565422e+00, 0.00000000e+00/)
392 gb( 5, 4, :) = (/&
393  & 0.12343189e+01, 0.10927103e+01, 0.10000000e+01/)
394 !
395 !----- INTERVAL = 2 ----- T = 250.0
396 !
397 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
398 ga( 6, 3, :) = (/&
399  & 0.12055643e+01, 0.49968044e+00, 0.00000000e+00/)
400 gb( 6, 3, :) = (/&
401  & 0.12055643e+01, 0.10702313e+01, 0.10000000e+01/)
402 ga( 6, 4, :) = (/&
403  & 0.12404147e+01, 0.56878618e+00, 0.00000000e+00/)
404 gb( 6, 4, :) = (/&
405  & 0.12404147e+01, 0.10982489e+01, 0.10000000e+01/)
406 !
407 !----- INTERVAL = 2 ----- T = 262.5
408 !
409 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
410 ga( 7, 3, :) = (/&
411  & 0.12114186e+01, 0.51214132e+00, 0.00000000e+00/)
412 gb( 7, 3, :) = (/&
413  & 0.12114186e+01, 0.10753907e+01, 0.10000000e+01/)
414 ga( 7, 4, :) = (/&
415  & 0.12458431e+01, 0.58047395e+00, 0.00000000e+00/)
416 gb( 7, 4, :) = (/&
417  & 0.12458431e+01, 0.11032019e+01, 0.10000000e+01/)
418 !
419 !----- INTERVAL = 2 ----- T = 275.0
420 !
421 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
422 ga( 8, 3, :) = (/&
423  & 0.12167192e+01, 0.52341830e+00, 0.00000000e+00/)
424 gb( 8, 3, :) = (/&
425  & 0.12167192e+01, 0.10800762e+01, 0.10000000e+01/)
426 ga( 8, 4, :) = (/&
427  & 0.12506907e+01, 0.59089894e+00, 0.00000000e+00/)
428 gb( 8, 4, :) = (/&
429  & 0.12506907e+01, 0.11076379e+01, 0.10000000e+01/)
430 !
431 !----- INTERVAL = 2 ----- T = 287.5
432 !
433 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
434 ga( 9, 3, :) = (/&
435  & 0.12215344e+01, 0.53365803e+00, 0.00000000e+00/)
436 gb( 9, 3, :) = (/&
437  & 0.12215344e+01, 0.10843446e+01, 0.10000000e+01/)
438 ga( 9, 4, :) = (/&
439  & 0.12550299e+01, 0.60021475e+00, 0.00000000e+00/)
440 gb( 9, 4, :) = (/&
441  & 0.12550299e+01, 0.11116160e+01, 0.10000000e+01/)
442 !
443 !----- INTERVAL = 2 ----- T = 300.0
444 !
445 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
446 ga(10, 3, :) = (/&
447  & 0.12259226e+01, 0.54298448e+00, 0.00000000e+00/)
448 gb(10, 3, :) = (/&
449  & 0.12259226e+01, 0.10882439e+01, 0.10000000e+01/)
450 ga(10, 4, :) = (/&
451  & 0.12589256e+01, 0.60856112e+00, 0.00000000e+00/)
452 gb(10, 4, :) = (/&
453  & 0.12589256e+01, 0.11151910e+01, 0.10000000e+01/)
454 !
455 !----- INTERVAL = 2 ----- T = 312.5
456 !
457 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
458 ga(11, 3, :) = (/&
459  & 0.12299344e+01, 0.55150227e+00, 0.00000000e+00/)
460 gb(11, 3, :) = (/&
461  & 0.12299344e+01, 0.10918144e+01, 0.10000000e+01/)
462 ga(11, 4, :) = (/&
463  & 0.12624402e+01, 0.61607594e+00, 0.00000000e+00/)
464 gb(11, 4, :) = (/&
465  & 0.12624402e+01, 0.11184188e+01, 0.10000000e+01/)
466 !
467 !
468 !
469 !
470 !- H2O 350-500 cm-1
471 !
472 !-- G = - 0.2*SLA, 0.0 +0.5/(1+0.5U)
473 !
474 !----- INTERVAL = 5 ----- T = 187.5
475 !
476 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
477 ga( 1, 5, :) = (/&
478  & 0.15750172e+00,-0.22159303e-01, 0.00000000e+00/)
479 gb( 1, 5, :) = (/&
480  & 0.15750172e+00, 0.38103212e+00, 0.10000000e+01/)
481 ga( 1, 6, :) = (/&
482  & 0.17770551e+00,-0.24972399e-01, 0.00000000e+00/)
483 gb( 1, 6, :) = (/&
484  & 0.17770551e+00, 0.41646579e+00, 0.10000000e+01/)
485 !
486 !----- INTERVAL = 5 ----- T = 200.0
487 !
488 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
489 ga( 2, 5, :) = (/&
490  & 0.16174076e+00,-0.22748917e-01, 0.00000000e+00/)
491 gb( 2, 5, :) = (/&
492  & 0.16174076e+00, 0.38913800e+00, 0.10000000e+01/)
493 ga( 2, 6, :) = (/&
494  & 0.18176757e+00,-0.25537247e-01, 0.00000000e+00/)
495 gb( 2, 6, :) = (/&
496  & 0.18176757e+00, 0.42345095e+00, 0.10000000e+01/)
497 !
498 !----- INTERVAL = 5 ----- T = 212.5
499 !
500 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
501 ga( 3, 5, :) = (/&
502  & 0.16548628e+00,-0.23269898e-01, 0.00000000e+00/)
503 gb( 3, 5, :) = (/&
504  & 0.16548628e+00, 0.39613651e+00, 0.10000000e+01/)
505 ga( 3, 6, :) = (/&
506  & 0.18527967e+00,-0.26025624e-01, 0.00000000e+00/)
507 gb( 3, 6, :) = (/&
508  & 0.18527967e+00, 0.42937476e+00, 0.10000000e+01/)
509 !
510 !----- INTERVAL = 5 ----- T = 225.0
511 !
512 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
513 ga( 4, 5, :) = (/&
514  & 0.16881124e+00,-0.23732392e-01, 0.00000000e+00/)
515 gb( 4, 5, :) = (/&
516  & 0.16881124e+00, 0.40222421e+00, 0.10000000e+01/)
517 ga( 4, 6, :) = (/&
518  & 0.18833348e+00,-0.26450280e-01, 0.00000000e+00/)
519 gb( 4, 6, :) = (/&
520  & 0.18833348e+00, 0.43444062e+00, 0.10000000e+01/)
521 !
522 !----- INTERVAL = 5 ----- T = 237.5
523 !
524 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
525 ga( 5, 5, :) = (/&
526  & 0.17177839e+00,-0.24145123e-01, 0.00000000e+00/)
527 gb( 5, 5, :) = (/&
528  & 0.17177839e+00, 0.40756010e+00, 0.10000000e+01/)
529 ga( 5, 6, :) = (/&
530  & 0.19100108e+00,-0.26821236e-01, 0.00000000e+00/)
531 gb( 5, 6, :) = (/&
532  & 0.19100108e+00, 0.43880316e+00, 0.10000000e+01/)
533 !
534 !----- INTERVAL = 5 ----- T = 250.0
535 !
536 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
537 ga( 6, 5, :) = (/&
538  & 0.17443933e+00,-0.24515269e-01, 0.00000000e+00/)
539 gb( 6, 5, :) = (/&
540  & 0.17443933e+00, 0.41226954e+00, 0.10000000e+01/)
541 ga( 6, 6, :) = (/&
542  & 0.19334122e+00,-0.27146657e-01, 0.00000000e+00/)
543 gb( 6, 6, :) = (/&
544  & 0.19334122e+00, 0.44258354e+00, 0.10000000e+01/)
545 !
546 !----- INTERVAL = 5 ----- T = 262.5
547 !
548 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
549 ga( 7, 5, :) = (/&
550  & 0.17683622e+00,-0.24848690e-01, 0.00000000e+00/)
551 gb( 7, 5, :) = (/&
552  & 0.17683622e+00, 0.41645142e+00, 0.10000000e+01/)
553 ga( 7, 6, :) = (/&
554  & 0.19540288e+00,-0.27433354e-01, 0.00000000e+00/)
555 gb( 7, 6, :) = (/&
556  & 0.19540288e+00, 0.44587882e+00, 0.10000000e+01/)
557 !
558 !----- INTERVAL = 5 ----- T = 275.0
559 !
560 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
561 ga( 8, 5, :) = (/&
562  & 0.17900375e+00,-0.25150210e-01, 0.00000000e+00/)
563 gb( 8, 5, :) = (/&
564  & 0.17900375e+00, 0.42018474e+00, 0.10000000e+01/)
565 ga( 8, 6, :) = (/&
566  & 0.19722732e+00,-0.27687065e-01, 0.00000000e+00/)
567 gb( 8, 6, :) = (/&
568  & 0.19722732e+00, 0.44876776e+00, 0.10000000e+01/)
569 !
570 !----- INTERVAL = 5 ----- T = 287.5
571 !
572 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
573 ga( 9, 5, :) = (/&
574  & 0.18097099e+00,-0.25423873e-01, 0.00000000e+00/)
575 gb( 9, 5, :) = (/&
576  & 0.18097099e+00, 0.42353379e+00, 0.10000000e+01/)
577 ga( 9, 6, :) = (/&
578  & 0.19884918e+00,-0.27912608e-01, 0.00000000e+00/)
579 gb( 9, 6, :) = (/&
580  & 0.19884918e+00, 0.45131451e+00, 0.10000000e+01/)
581 !
582 !----- INTERVAL = 5 ----- T = 300.0
583 !
584 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
585 ga(10, 5, :) = (/&
586  & 0.18276283e+00,-0.25673139e-01, 0.00000000e+00/)
587 gb(10, 5, :) = (/&
588  & 0.18276283e+00, 0.42655211e+00, 0.10000000e+01/)
589 ga(10, 6, :) = (/&
590  & 0.20029696e+00,-0.28113944e-01, 0.00000000e+00/)
591 gb(10, 6, :) = (/&
592  & 0.20029696e+00, 0.45357095e+00, 0.10000000e+01/)
593 !
594 !----- INTERVAL = 5 ----- T = 312.5
595 !
596 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
597 ga(11, 5, :) = (/&
598  & 0.18440117e+00,-0.25901055e-01, 0.00000000e+00/)
599 gb(11, 5, :) = (/&
600  & 0.18440117e+00, 0.42928533e+00, 0.10000000e+01/)
601 ga(11, 6, :) = (/&
602  & 0.20159300e+00,-0.28294180e-01, 0.00000000e+00/)
603 gb(11, 6, :) = (/&
604  & 0.20159300e+00, 0.45557797e+00, 0.10000000e+01/)
605 
606 
607 
608 !
609 !
610 !- H2O 800-970 + 1110-1250 cm-1
611 !
612 !--- G = 3.875E-03 ---------------
613 !
614 !----- INTERVAL = 3 ----- T = 187.5
615 !
616 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
617 ga( 1, 7, :) = (/&
618  & 0.10192131e+02, 0.80737799e+01, 0.00000000e+00/)
619 gb( 1, 7, :) = (/&
620  & 0.10192131e+02, 0.82623280e+01, 0.10000000e+01/)
621 ga( 1, 8, :) = (/&
622  & 0.92439050e+01, 0.77425778e+01, 0.00000000e+00/)
623 gb( 1, 8, :) = (/&
624  & 0.92439050e+01, 0.79342219e+01, 0.10000000e+01/)
625 !
626 !----- INTERVAL = 3 ----- T = 200.0
627 !
628 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
629 ga( 2, 7, :) = (/&
630  & 0.97258602e+01, 0.79171158e+01, 0.00000000e+00/)
631 gb( 2, 7, :) = (/&
632  & 0.97258602e+01, 0.81072291e+01, 0.10000000e+01/)
633 ga( 2, 8, :) = (/&
634  & 0.87567422e+01, 0.75443460e+01, 0.00000000e+00/)
635 gb( 2, 8, :) = (/&
636  & 0.87567422e+01, 0.77373458e+01, 0.10000000e+01/)
637 !
638 !----- INTERVAL = 3 ----- T = 212.5
639 !
640 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
641 ga( 3, 7, :) = (/&
642  & 0.92992890e+01, 0.77609605e+01, 0.00000000e+00/)
643 gb( 3, 7, :) = (/&
644  & 0.92992890e+01, 0.79523834e+01, 0.10000000e+01/)
645 ga( 3, 8, :) = (/&
646  & 0.83270144e+01, 0.73526151e+01, 0.00000000e+00/)
647 gb( 3, 8, :) = (/&
648  & 0.83270144e+01, 0.75467334e+01, 0.10000000e+01/)
649 !
650 !----- INTERVAL = 3 ----- T = 225.0
651 !
652 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
653 ga( 4, 7, :) = (/&
654  & 0.89154021e+01, 0.76087371e+01, 0.00000000e+00/)
655 gb( 4, 7, :) = (/&
656  & 0.89154021e+01, 0.78012527e+01, 0.10000000e+01/)
657 ga( 4, 8, :) = (/&
658  & 0.79528337e+01, 0.71711188e+01, 0.00000000e+00/)
659 gb( 4, 8, :) = (/&
660  & 0.79528337e+01, 0.73661786e+01, 0.10000000e+01/)
661 !
662 !----- INTERVAL = 3 ----- T = 237.5
663 !
664 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
665 ga( 5, 7, :) = (/&
666  & 0.85730084e+01, 0.74627112e+01, 0.00000000e+00/)
667 gb( 5, 7, :) = (/&
668  & 0.85730084e+01, 0.76561458e+01, 0.10000000e+01/)
669 ga( 5, 8, :) = (/&
670  & 0.76286839e+01, 0.70015571e+01, 0.00000000e+00/)
671 gb( 5, 8, :) = (/&
672  & 0.76286839e+01, 0.71974319e+01, 0.10000000e+01/)
673 !
674 !----- INTERVAL = 3 ----- T = 250.0
675 !
676 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
677 ga( 6, 7, :) = (/&
678  & 0.82685838e+01, 0.73239981e+01, 0.00000000e+00/)
679 gb( 6, 7, :) = (/&
680  & 0.82685838e+01, 0.75182174e+01, 0.10000000e+01/)
681 ga( 6, 8, :) = (/&
682  & 0.73477879e+01, 0.68442532e+01, 0.00000000e+00/)
683 gb( 6, 8, :) = (/&
684  & 0.73477879e+01, 0.70408543e+01, 0.10000000e+01/)
685 !
686 !----- INTERVAL = 3 ----- T = 262.5
687 !
688 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
689 ga( 7, 7, :) = (/&
690  & 0.79978921e+01, 0.71929934e+01, 0.00000000e+00/)
691 gb( 7, 7, :) = (/&
692  & 0.79978921e+01, 0.73878952e+01, 0.10000000e+01/)
693 ga( 7, 8, :) = (/&
694  & 0.71035818e+01, 0.66987996e+01, 0.00000000e+00/)
695 gb( 7, 8, :) = (/&
696  & 0.71035818e+01, 0.68960649e+01, 0.10000000e+01/)
697 !
698 !----- INTERVAL = 3 ----- T = 275.0
699 !
700 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
701 ga( 8, 7, :) = (/&
702  & 0.77568055e+01, 0.70697065e+01, 0.00000000e+00/)
703 gb( 8, 7, :) = (/&
704  & 0.77568055e+01, 0.72652133e+01, 0.10000000e+01/)
705 ga( 8, 8, :) = (/&
706  & 0.68903312e+01, 0.65644820e+01, 0.00000000e+00/)
707 gb( 8, 8, :) = (/&
708  & 0.68903312e+01, 0.67623672e+01, 0.10000000e+01/)
709 !
710 !----- INTERVAL = 3 ----- T = 287.5
711 !
712 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
713 ga( 9, 7, :) = (/&
714  & 0.75416266e+01, 0.69539626e+01, 0.00000000e+00/)
715 gb( 9, 7, :) = (/&
716  & 0.75416266e+01, 0.71500151e+01, 0.10000000e+01/)
717 ga( 9, 8, :) = (/&
718  & 0.67032875e+01, 0.64405267e+01, 0.00000000e+00/)
719 gb( 9, 8, :) = (/&
720  & 0.67032875e+01, 0.66389989e+01, 0.10000000e+01/)
721 !
722 !----- INTERVAL = 3 ----- T = 300.0
723 !
724 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
725 ga(10, 7, :) = (/&
726  & 0.73491694e+01, 0.68455144e+01, 0.00000000e+00/)
727 gb(10, 7, :) = (/&
728  & 0.73491694e+01, 0.70420667e+01, 0.10000000e+01/)
729 ga(10, 8, :) = (/&
730  & 0.65386461e+01, 0.63262376e+01, 0.00000000e+00/)
731 gb(10, 8, :) = (/&
732  & 0.65386461e+01, 0.65252707e+01, 0.10000000e+01/)
733 !
734 !----- INTERVAL = 3 ----- T = 312.5
735 !
736 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
737 ga(11, 7, :) = (/&
738  & 0.71767400e+01, 0.67441020e+01, 0.00000000e+00/)
739 gb(11, 7, :) = (/&
740  & 0.71767400e+01, 0.69411177e+01, 0.10000000e+01/)
741 ga(11, 8, :) = (/&
742  & 0.63934377e+01, 0.62210701e+01, 0.00000000e+00/)
743 gb(11, 8, :) = (/&
744  & 0.63934377e+01, 0.64206412e+01, 0.10000000e+01/)
745 !
746 !
747 !
748 !- H2O 970 - 1110 cm-1
749 !
750 !!-- G = 3.6E-03
751 !
752 !----- INTERVAL = 4 ----- T = 187.5
753 !
754 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
755 ga( 1, 9, :) = (/&
756  & 0.24870635e+02, 0.10542131e+02, 0.00000000e+00/)
757 gb( 1, 9, :) = (/&
758  & 0.24870635e+02, 0.10656640e+02, 0.10000000e+01/)
759 ga( 1,10, :) = (/&
760  & 0.24586283e+02, 0.10490353e+02, 0.00000000e+00/)
761 gb( 1,10, :) = (/&
762  & 0.24586283e+02, 0.10605856e+02, 0.10000000e+01/)
763 !
764 !----- INTERVAL = 4 ----- T = 200.0
765 !
766 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
767 ga( 2, 9, :) = (/&
768  & 0.24725591e+02, 0.10515895e+02, 0.00000000e+00/)
769 gb( 2, 9, :) = (/&
770  & 0.24725591e+02, 0.10630910e+02, 0.10000000e+01/)
771 ga( 2,10, :) = (/&
772  & 0.24441465e+02, 0.10463512e+02, 0.00000000e+00/)
773 gb( 2,10, :) = (/&
774  & 0.24441465e+02, 0.10579514e+02, 0.10000000e+01/)
775 !
776 !----- INTERVAL = 4 ----- T = 212.5
777 !
778 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
779 ga( 3, 9, :) = (/&
780  & 0.24600320e+02, 0.10492949e+02, 0.00000000e+00/)
781 gb( 3, 9, :) = (/&
782  & 0.24600320e+02, 0.10608399e+02, 0.10000000e+01/)
783 ga( 3,10, :) = (/&
784  & 0.24311657e+02, 0.10439183e+02, 0.00000000e+00/)
785 gb( 3,10, :) = (/&
786  & 0.24311657e+02, 0.10555632e+02, 0.10000000e+01/)
787 !
788 !----- INTERVAL = 4 ----- T = 225.0
789 !
790 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
791 ga( 4, 9, :) = (/&
792  & 0.24487300e+02, 0.10472049e+02, 0.00000000e+00/)
793 gb( 4, 9, :) = (/&
794  & 0.24487300e+02, 0.10587891e+02, 0.10000000e+01/)
795 ga( 4,10, :) = (/&
796  & 0.24196167e+02, 0.10417324e+02, 0.00000000e+00/)
797 gb( 4,10, :) = (/&
798  & 0.24196167e+02, 0.10534169e+02, 0.10000000e+01/)
799 !
800 !----- INTERVAL = 4 ----- T = 237.5
801 !
802 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
803 ga( 5, 9, :) = (/&
804  & 0.24384935e+02, 0.10452961e+02, 0.00000000e+00/)
805 gb( 5, 9, :) = (/&
806  & 0.24384935e+02, 0.10569156e+02, 0.10000000e+01/)
807 ga( 5,10, :) = (/&
808  & 0.24093406e+02, 0.10397704e+02, 0.00000000e+00/)
809 gb( 5,10, :) = (/&
810  & 0.24093406e+02, 0.10514900e+02, 0.10000000e+01/)
811 !
812 !----- INTERVAL = 4 ----- T = 250.0
813 !
814 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
815 ga( 6, 9, :) = (/&
816  & 0.24292341e+02, 0.10435562e+02, 0.00000000e+00/)
817 gb( 6, 9, :) = (/&
818  & 0.24292341e+02, 0.10552075e+02, 0.10000000e+01/)
819 ga( 6,10, :) = (/&
820  & 0.24001597e+02, 0.10380038e+02, 0.00000000e+00/)
821 gb( 6,10, :) = (/&
822  & 0.24001597e+02, 0.10497547e+02, 0.10000000e+01/)
823 !
824 !----- INTERVAL = 4 ----- T = 262.5
825 !
826 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
827 ga( 7, 9, :) = (/&
828  & 0.24208572e+02, 0.10419710e+02, 0.00000000e+00/)
829 gb( 7, 9, :) = (/&
830  & 0.24208572e+02, 0.10536510e+02, 0.10000000e+01/)
831 ga( 7,10, :) = (/&
832  & 0.23919098e+02, 0.10364052e+02, 0.00000000e+00/)
833 gb( 7,10, :) = (/&
834  & 0.23919098e+02, 0.10481842e+02, 0.10000000e+01/)
835 !
836 !----- INTERVAL = 4 ----- T = 275.0
837 !
838 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
839 ga( 8, 9, :) = (/&
840  & 0.24132642e+02, 0.10405247e+02, 0.00000000e+00/)
841 gb( 8, 9, :) = (/&
842  & 0.24132642e+02, 0.10522307e+02, 0.10000000e+01/)
843 ga( 8,10, :) = (/&
844  & 0.23844511e+02, 0.10349509e+02, 0.00000000e+00/)
845 gb( 8,10, :) = (/&
846  & 0.23844511e+02, 0.10467553e+02, 0.10000000e+01/)
847 !
848 !----- INTERVAL = 4 ----- T = 287.5
849 !
850 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
851 ga( 9, 9, :) = (/&
852  & 0.24063614e+02, 0.10392022e+02, 0.00000000e+00/)
853 gb( 9, 9, :) = (/&
854  & 0.24063614e+02, 0.10509317e+02, 0.10000000e+01/)
855 ga( 9,10, :) = (/&
856  & 0.23776708e+02, 0.10336215e+02, 0.00000000e+00/)
857 gb( 9,10, :) = (/&
858  & 0.23776708e+02, 0.10454488e+02, 0.10000000e+01/)
859 !
860 !----- INTERVAL = 4 ----- T = 300.0
861 !
862 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
863 ga(10, 9, :) = (/&
864  & 0.24000649e+02, 0.10379892e+02, 0.00000000e+00/)
865 gb(10, 9, :) = (/&
866  & 0.24000649e+02, 0.10497402e+02, 0.10000000e+01/)
867 ga(10,10, :) = (/&
868  & 0.23714816e+02, 0.10324018e+02, 0.00000000e+00/)
869 gb(10,10, :) = (/&
870  & 0.23714816e+02, 0.10442501e+02, 0.10000000e+01/)
871 !
872 !----- INTERVAL = 4 ----- T = 312.5
873 !
874 !-- INDICES FOR PADE APPROXIMATION 1 28 37 45
875 ga(11, 9, :) = (/&
876  & 0.23943021e+02, 0.10368736e+02, 0.00000000e+00/)
877 gb(11, 9, :) = (/&
878  & 0.23943021e+02, 0.10486443e+02, 0.10000000e+01/)
879 ga(11,10, :) = (/&
880  & 0.23658197e+02, 0.10312808e+02, 0.00000000e+00/)
881 gb(11,10, :) = (/&
882  & 0.23658197e+02, 0.10431483e+02, 0.10000000e+01/)
883 !
884 !
885 !
886 !- H2O 1250-1450 + 1880-2820 cm-1
887 !
888 !--- G = 0.0
889 !
890 !
891 !----- INTERVAL = 6 ----- T = 187.5
892 !
893 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
894 ga( 1,11, :) = (/&
895  & 0.11990218e+02,-0.12823142e+01, 0.00000000e+00/)
896 gb( 1,11, :) = (/&
897  & 0.11990218e+02, 0.26681588e+02, 0.10000000e+01/)
898 ga( 1,12, :) = (/&
899  & 0.79709806e+01,-0.74805226e+00, 0.00000000e+00/)
900 gb( 1,12, :) = (/&
901  & 0.79709806e+01, 0.18377807e+02, 0.10000000e+01/)
902 !
903 !----- INTERVAL = 6 ----- T = 200.0
904 !
905 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
906 ga( 2,11, :) = (/&
907  & 0.10904073e+02,-0.10571588e+01, 0.00000000e+00/)
908 gb( 2,11, :) = (/&
909  & 0.10904073e+02, 0.24728346e+02, 0.10000000e+01/)
910 ga( 2,12, :) = (/&
911  & 0.75400737e+01,-0.56252739e+00, 0.00000000e+00/)
912 gb( 2,12, :) = (/&
913  & 0.75400737e+01, 0.17643148e+02, 0.10000000e+01/)
914 !
915 !----- INTERVAL = 6 ----- T = 212.5
916 !
917 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
918 ga( 3,11, :) = (/&
919  & 0.89126838e+01,-0.74864953e+00, 0.00000000e+00/)
920 gb( 3,11, :) = (/&
921  & 0.89126838e+01, 0.20551342e+02, 0.10000000e+01/)
922 ga( 3,12, :) = (/&
923  & 0.81804377e+01,-0.46188072e+00, 0.00000000e+00/)
924 gb( 3,12, :) = (/&
925  & 0.81804377e+01, 0.19296161e+02, 0.10000000e+01/)
926 !
927 !----- INTERVAL = 6 ----- T = 225.0
928 !
929 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
930 ga( 4,11, :) = (/&
931  & 0.85622405e+01,-0.58705980e+00, 0.00000000e+00/)
932 gb( 4,11, :) = (/&
933  & 0.85622405e+01, 0.19955244e+02, 0.10000000e+01/)
934 ga( 4,12, :) = (/&
935  & 0.10564339e+02,-0.40712065e+00, 0.00000000e+00/)
936 gb( 4,12, :) = (/&
937  & 0.10564339e+02, 0.24951120e+02, 0.10000000e+01/)
938 !
939 !----- INTERVAL = 6 ----- T = 237.5
940 !
941 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
942 ga( 5,11, :) = (/&
943  & 0.94892164e+01,-0.49305772e+00, 0.00000000e+00/)
944 gb( 5,11, :) = (/&
945  & 0.94892164e+01, 0.22227100e+02, 0.10000000e+01/)
946 ga( 5,12, :) = (/&
947  & 0.46896789e+02,-0.15295996e+01, 0.00000000e+00/)
948 gb( 5,12, :) = (/&
949  & 0.46896789e+02, 0.10957372e+03, 0.10000000e+01/)
950 !
951 !----- INTERVAL = 6 ----- T = 250.0
952 !
953 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
954 ga( 6,11, :) = (/&
955  & 0.13580937e+02,-0.51461431e+00, 0.00000000e+00/)
956 gb( 6,11, :) = (/&
957  & 0.13580937e+02, 0.31770288e+02, 0.10000000e+01/)
958 ga( 6,12, :) = (/&
959  &-0.30926524e+01, 0.43555255e+00, 0.00000000e+00/)
960 gb( 6,12, :) = (/&
961  &-0.30926524e+01,-0.67432659e+01, 0.10000000e+01/)
962 !
963 !----- INTERVAL = 6 ----- T = 262.5
964 !
965 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
966 ga( 7,11, :) = (/&
967  &-0.32050918e+03, 0.12373350e+02, 0.00000000e+00/)
968 gb( 7,11, :) = (/&
969  &-0.32050918e+03,-0.74061287e+03, 0.10000000e+01/)
970 ga( 7,12, :) = (/&
971  & 0.85742941e+00, 0.50380874e+00, 0.00000000e+00/)
972 gb( 7,12, :) = (/&
973  & 0.85742941e+00, 0.24550746e+01, 0.10000000e+01/)
974 !
975 !----- INTERVAL = 6 ----- T = 275.0
976 !
977 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
978 ga( 8,11, :) = (/&
979  &-0.37133165e+01, 0.44809588e+00, 0.00000000e+00/)
980 gb( 8,11, :) = (/&
981  &-0.37133165e+01,-0.81329826e+01, 0.10000000e+01/)
982 ga( 8,12, :) = (/&
983  & 0.19164038e+01, 0.68537352e+00, 0.00000000e+00/)
984 gb( 8,12, :) = (/&
985  & 0.19164038e+01, 0.49089917e+01, 0.10000000e+01/)
986 !
987 !----- INTERVAL = 6 ----- T = 287.5
988 !
989 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
990 ga( 9,11, :) = (/&
991  & 0.18890836e+00, 0.46548918e+00, 0.00000000e+00/)
992 gb( 9,11, :) = (/&
993  & 0.18890836e+00, 0.90279822e+00, 0.10000000e+01/)
994 ga( 9,12, :) = (/&
995  & 0.23513199e+01, 0.89437630e+00, 0.00000000e+00/)
996 gb( 9,12, :) = (/&
997  & 0.23513199e+01, 0.59008712e+01, 0.10000000e+01/)
998 !
999 !----- INTERVAL = 6 ----- T = 300.0
1000 !
1001 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
1002 ga(10,11, :) = (/&
1003  & 0.14209226e+01, 0.59121475e+00, 0.00000000e+00/)
1004 gb(10,11, :) = (/&
1005  & 0.14209226e+01, 0.37532746e+01, 0.10000000e+01/)
1006 ga(10,12, :) = (/&
1007  & 0.25566644e+01, 0.11127003e+01, 0.00000000e+00/)
1008 gb(10,12, :) = (/&
1009  & 0.25566644e+01, 0.63532616e+01, 0.10000000e+01/)
1010 !
1011 !----- INTERVAL = 6 ----- T = 312.5
1012 !
1013 !-- INDICES FOR PADE APPROXIMATION 1 35 40 45
1014 ga(11,11, :) = (/&
1015  & 0.19817679e+01, 0.74676119e+00, 0.00000000e+00/)
1016 gb(11,11, :) = (/&
1017  & 0.19817679e+01, 0.50437916e+01, 0.10000000e+01/)
1018 ga(11,12, :) = (/&
1019  & 0.26555181e+01, 0.13329782e+01, 0.00000000e+00/)
1020 gb(11,12, :) = (/&
1021  & 0.26555181e+01, 0.65558627e+01, 0.10000000e+01/)
1022 !
1023 !
1024 
1025 !
1026 !
1027 !
1028 !-- CO2 -- INT.2 -- 500-800 CM-1 --- FROM ABS225 ----------------------
1029 !
1030 !
1031 !
1032 !-- FIU = 0.8 + MAX(0.35,(7-IU)*0.9) , X/T, 9
1033 !
1034 !----- INTERVAL = 2 ----- T = 187.5
1035 !
1036 !-- INDICES FOR PADE APPROXIMATION 1 30 38 45
1037 ga( 1,13, :) = (/&
1038  & 0.87668459e-01, 0.13845511e+01, 0.00000000e+00/)
1039 gb( 1,13, :) = (/&
1040  & 0.87668459e-01, 0.23203798e+01, 0.10000000e+01/)
1041 ga( 1,14, :) = (/&
1042  & 0.74878820e-01, 0.11718758e+01, 0.00000000e+00/)
1043 gb( 1,14, :) = (/&
1044  & 0.74878820e-01, 0.20206726e+01, 0.10000000e+01/)
1045 !
1046 !----- INTERVAL = 2 ----- T = 200.0
1047 !
1048 !-- INDICES FOR PADE APPROXIMATION 1 30 38 45
1049 ga( 2,13, :) = (/&
1050  & 0.83754276e-01, 0.13187042e+01, 0.00000000e+00/)
1051 gb( 2,13, :) = (/&
1052  & 0.83754276e-01, 0.22288925e+01, 0.10000000e+01/)
1053 ga( 2,14, :) = (/&
1054  & 0.71650966e-01, 0.11216131e+01, 0.00000000e+00/)
1055 gb( 2,14, :) = (/&
1056  & 0.71650966e-01, 0.19441824e+01, 0.10000000e+01/)
1057 !
1058 !----- INTERVAL = 2 ----- T = 212.5
1059 !
1060 !-- INDICES FOR PADE APPROXIMATION 1 30 38 45
1061 ga( 3,13, :) = (/&
1062  & 0.80460283e-01, 0.12644396e+01, 0.00000000e+00/)
1063 gb( 3,13, :) = (/&
1064  & 0.80460283e-01, 0.21515593e+01, 0.10000000e+01/)
1065 ga( 3,14, :) = (/&
1066  & 0.68979615e-01, 0.10809473e+01, 0.00000000e+00/)
1067 gb( 3,14, :) = (/&
1068  & 0.68979615e-01, 0.18807257e+01, 0.10000000e+01/)
1069 !
1070 !----- INTERVAL = 2 ----- T = 225.0
1071 !
1072 !-- INDICES FOR PADE APPROXIMATION 1 30 38 45
1073 ga( 4,13, :) = (/&
1074  & 0.77659686e-01, 0.12191543e+01, 0.00000000e+00/)
1075 gb( 4,13, :) = (/&
1076  & 0.77659686e-01, 0.20855896e+01, 0.10000000e+01/)
1077 ga( 4,14, :) = (/&
1078  & 0.66745345e-01, 0.10476396e+01, 0.00000000e+00/)
1079 gb( 4,14, :) = (/&
1080  & 0.66745345e-01, 0.18275618e+01, 0.10000000e+01/)
1081 !
1082 !----- INTERVAL = 2 ----- T = 237.5
1083 !
1084 !-- INDICES FOR PADE APPROXIMATION 1 30 38 45
1085 ga( 5,13, :) = (/&
1086  & 0.75257056e-01, 0.11809511e+01, 0.00000000e+00/)
1087 gb( 5,13, :) = (/&
1088  & 0.75257056e-01, 0.20288489e+01, 0.10000000e+01/)
1089 ga( 5,14, :) = (/&
1090  & 0.64857571e-01, 0.10200373e+01, 0.00000000e+00/)
1091 gb( 5,14, :) = (/&
1092  & 0.64857571e-01, 0.17825910e+01, 0.10000000e+01/)
1093 !
1094 !----- INTERVAL = 2 ----- T = 250.0
1095 !
1096 !-- INDICES FOR PADE APPROXIMATION 1 30 38 45
1097 ga( 6,13, :) = (/&
1098  & 0.73179175e-01, 0.11484154e+01, 0.00000000e+00/)
1099 gb( 6,13, :) = (/&
1100  & 0.73179175e-01, 0.19796791e+01, 0.10000000e+01/)
1101 ga( 6,14, :) = (/&
1102  & 0.63248495e-01, 0.99692726e+00, 0.00000000e+00/)
1103 gb( 6,14, :) = (/&
1104  & 0.63248495e-01, 0.17442308e+01, 0.10000000e+01/)
1105 !
1106 !----- INTERVAL = 2 ----- T = 262.5
1107 !
1108 !-- INDICES FOR PADE APPROXIMATION 1 30 38 45
1109 ga( 7,13, :) = (/&
1110  & 0.71369063e-01, 0.11204723e+01, 0.00000000e+00/)
1111 gb( 7,13, :) = (/&
1112  & 0.71369063e-01, 0.19367778e+01, 0.10000000e+01/)
1113 ga( 7,14, :) = (/&
1114  & 0.61866970e-01, 0.97740923e+00, 0.00000000e+00/)
1115 gb( 7,14, :) = (/&
1116  & 0.61866970e-01, 0.17112809e+01, 0.10000000e+01/)
1117 !
1118 !----- INTERVAL = 2 ----- T = 275.0
1119 !
1120 !-- INDICES FOR PADE APPROXIMATION 1 30 38 45
1121 ga( 8,13, :) = (/&
1122  & 0.69781812e-01, 0.10962918e+01, 0.00000000e+00/)
1123 gb( 8,13, :) = (/&
1124  & 0.69781812e-01, 0.18991112e+01, 0.10000000e+01/)
1125 ga( 8,14, :) = (/&
1126  & 0.60673632e-01, 0.96080188e+00, 0.00000000e+00/)
1127 gb( 8,14, :) = (/&
1128  & 0.60673632e-01, 0.16828137e+01, 0.10000000e+01/)
1129 !
1130 !----- INTERVAL = 2 ----- T = 287.5
1131 !
1132 !-- INDICES FOR PADE APPROXIMATION 1 30 38 45
1133 ga( 9,13, :) = (/&
1134  & 0.68381606e-01, 0.10752229e+01, 0.00000000e+00/)
1135 gb( 9,13, :) = (/&
1136  & 0.68381606e-01, 0.18658501e+01, 0.10000000e+01/)
1137 ga( 9,14, :) = (/&
1138  & 0.59637277e-01, 0.94657562e+00, 0.00000000e+00/)
1139 gb( 9,14, :) = (/&
1140  & 0.59637277e-01, 0.16580908e+01, 0.10000000e+01/)
1141 !
1142 !----- INTERVAL = 2 ----- T = 300.0
1143 !
1144 !-- INDICES FOR PADE APPROXIMATION 1 30 38 45
1145 ga(10,13, :) = (/&
1146  & 0.67139539e-01, 0.10567474e+01, 0.00000000e+00/)
1147 gb(10,13, :) = (/&
1148  & 0.67139539e-01, 0.18363226e+01, 0.10000000e+01/)
1149 ga(10,14, :) = (/&
1150  & 0.58732178e-01, 0.93430511e+00, 0.00000000e+00/)
1151 gb(10,14, :) = (/&
1152  & 0.58732178e-01, 0.16365014e+01, 0.10000000e+01/)
1153 !
1154 !----- INTERVAL = 2 ----- T = 312.5
1155 !
1156 !-- INDICES FOR PADE APPROXIMATION 1 30 38 45
1157 ga(11,13, :) = (/&
1158  & 0.66032012e-01, 0.10404465e+01, 0.00000000e+00/)
1159 gb(11,13, :) = (/&
1160  & 0.66032012e-01, 0.18099779e+01, 0.10000000e+01/)
1161 ga(11,14, :) = (/&
1162  & 0.57936092e-01, 0.92363528e+00, 0.00000000e+00/)
1163 gb(11,14, :) = (/&
1164  & 0.57936092e-01, 0.16175164e+01, 0.10000000e+01/)
1165 !
1166 !
1167 !
1168 !
1169 !
1170 !
1171 !
1172 !
1173 !
1174 !
1175 !-- CARBON DIOXIDE LINES IN THE WINDOW REGION (800-1250 CM-1)
1176 !
1177 !
1178 !-- G = 0.0
1179 !
1180 !
1181 !----- INTERVAL = 4 ----- T = 187.5
1182 !
1183 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
1184 ga( 1,15, :) = (/&
1185  & 0.13230067e+02, 0.22042132e+02, 0.00000000e+00/)
1186 gb( 1,15, :) = (/&
1187  & 0.13230067e+02, 0.22051750e+02, 0.10000000e+01/)
1188 ga( 1,16, :) = (/&
1189  & 0.13183816e+02, 0.22169501e+02, 0.00000000e+00/)
1190 gb( 1,16, :) = (/&
1191  & 0.13183816e+02, 0.22178972e+02, 0.10000000e+01/)
1192 !
1193 !----- INTERVAL = 4 ----- T = 200.0
1194 !
1195 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
1196 ga( 2,15, :) = (/&
1197  & 0.13213564e+02, 0.22107298e+02, 0.00000000e+00/)
1198 gb( 2,15, :) = (/&
1199  & 0.13213564e+02, 0.22116850e+02, 0.10000000e+01/)
1200 ga( 2,16, :) = (/&
1201  & 0.13189991e+02, 0.22270075e+02, 0.00000000e+00/)
1202 gb( 2,16, :) = (/&
1203  & 0.13189991e+02, 0.22279484e+02, 0.10000000e+01/)
1204 !
1205 !----- INTERVAL = 4 ----- T = 212.5
1206 !
1207 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
1208 ga( 3,15, :) = (/&
1209  & 0.13209140e+02, 0.22180915e+02, 0.00000000e+00/)
1210 gb( 3,15, :) = (/&
1211  & 0.13209140e+02, 0.22190410e+02, 0.10000000e+01/)
1212 ga( 3,16, :) = (/&
1213  & 0.13209485e+02, 0.22379193e+02, 0.00000000e+00/)
1214 gb( 3,16, :) = (/&
1215  & 0.13209485e+02, 0.22388551e+02, 0.10000000e+01/)
1216 !
1217 !----- INTERVAL = 4 ----- T = 225.0
1218 !
1219 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
1220 ga( 4,15, :) = (/&
1221  & 0.13213894e+02, 0.22259478e+02, 0.00000000e+00/)
1222 gb( 4,15, :) = (/&
1223  & 0.13213894e+02, 0.22268925e+02, 0.10000000e+01/)
1224 ga( 4,16, :) = (/&
1225  & 0.13238789e+02, 0.22492992e+02, 0.00000000e+00/)
1226 gb( 4,16, :) = (/&
1227  & 0.13238789e+02, 0.22502309e+02, 0.10000000e+01/)
1228 !
1229 !----- INTERVAL = 4 ----- T = 237.5
1230 !
1231 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
1232 ga( 5,15, :) = (/&
1233  & 0.13225963e+02, 0.22341039e+02, 0.00000000e+00/)
1234 gb( 5,15, :) = (/&
1235  & 0.13225963e+02, 0.22350445e+02, 0.10000000e+01/)
1236 ga( 5,16, :) = (/&
1237  & 0.13275017e+02, 0.22608508e+02, 0.00000000e+00/)
1238 gb( 5,16, :) = (/&
1239  & 0.13275017e+02, 0.22617792e+02, 0.10000000e+01/)
1240 !
1241 !----- INTERVAL = 4 ----- T = 250.0
1242 !
1243 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
1244 ga( 6,15, :) = (/&
1245  & 0.13243806e+02, 0.22424247e+02, 0.00000000e+00/)
1246 gb( 6,15, :) = (/&
1247  & 0.13243806e+02, 0.22433617e+02, 0.10000000e+01/)
1248 ga( 6,16, :) = (/&
1249  & 0.13316096e+02, 0.22723843e+02, 0.00000000e+00/)
1250 gb( 6,16, :) = (/&
1251  & 0.13316096e+02, 0.22733099e+02, 0.10000000e+01/)
1252 !
1253 !----- INTERVAL = 4 ----- T = 262.5
1254 !
1255 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
1256 ga( 7,15, :) = (/&
1257  & 0.13266104e+02, 0.22508089e+02, 0.00000000e+00/)
1258 gb( 7,15, :) = (/&
1259  & 0.13266104e+02, 0.22517429e+02, 0.10000000e+01/)
1260 ga( 7,16, :) = (/&
1261  & 0.13360555e+02, 0.22837837e+02, 0.00000000e+00/)
1262 gb( 7,16, :) = (/&
1263  & 0.13360555e+02, 0.22847071e+02, 0.10000000e+01/)
1264 !
1265 !----- INTERVAL = 4 ----- T = 275.0
1266 !
1267 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
1268 ga( 8,15, :) = (/&
1269  & 0.13291782e+02, 0.22591771e+02, 0.00000000e+00/)
1270 gb( 8,15, :) = (/&
1271  & 0.13291782e+02, 0.22601086e+02, 0.10000000e+01/)
1272 ga( 8,16, :) = (/&
1273  & 0.13407324e+02, 0.22949751e+02, 0.00000000e+00/)
1274 gb( 8,16, :) = (/&
1275  & 0.13407324e+02, 0.22958967e+02, 0.10000000e+01/)
1276 !
1277 !----- INTERVAL = 4 ----- T = 287.5
1278 !
1279 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
1280 ga( 9,15, :) = (/&
1281  & 0.13319961e+02, 0.22674661e+02, 0.00000000e+00/)
1282 gb( 9,15, :) = (/&
1283  & 0.13319961e+02, 0.22683956e+02, 0.10000000e+01/)
1284 ga( 9,16, :) = (/&
1285  & 0.13455544e+02, 0.23059032e+02, 0.00000000e+00/)
1286 gb( 9,16, :) = (/&
1287  & 0.13455544e+02, 0.23068234e+02, 0.10000000e+01/)
1288 !
1289 !----- INTERVAL = 4 ----- T = 300.0
1290 !
1291 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
1292 ga(10,15, :) = (/&
1293  & 0.13349927e+02, 0.22756246e+02, 0.00000000e+00/)
1294 gb(10,15, :) = (/&
1295  & 0.13349927e+02, 0.22765522e+02, 0.10000000e+01/)
1296 ga(10,16, :) = (/&
1297  & 0.13504450e+02, 0.23165146e+02, 0.00000000e+00/)
1298 gb(10,16, :) = (/&
1299  & 0.13504450e+02, 0.23174336e+02, 0.10000000e+01/)
1300 !
1301 !----- INTERVAL = 4 ----- T = 312.5
1302 !
1303 !-- INDICES FOR PADE APPROXIMATION 1 15 29 45
1304 ga(11,15, :) = (/&
1305  & 0.13381108e+02, 0.22836093e+02, 0.00000000e+00/)
1306 gb(11,15, :) = (/&
1307  & 0.13381108e+02, 0.22845354e+02, 0.10000000e+01/)
1308 ga(11,16, :) = (/&
1309  & 0.13553282e+02, 0.23267456e+02, 0.00000000e+00/)
1310 gb(11,16, :) = (/&
1311  & 0.13553282e+02, 0.23276638e+02, 0.10000000e+01/)
1312 !
1313 RETURN
1314 END SUBROUTINE suolw
suolw
subroutine suolw
Definition: suolwn.F90:2
yoeolw
Definition: yoeolw.F90:1