source: LMDZ5/trunk/libf/phymar/suolwn.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: 39.3 KB
Line 
1SUBROUTINE 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
50USE 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
58IMPLICIT NONE
59
60
61!      ----------------------------------------------------------------
62!
63!*       1.    SET VALUES.
64!              -----------
65!
66NG1=2
67NG1P1=NG1+1
68NISP=6
69NIPD=8
70NIPD2=2*NIPD
71NTR=11
72NTRA=15
73NUA=24
74!
75!     
76!--  COEFFICIENTS FOR THE PSEUDO-VOIGT LINE PROFILE (PA)
77!
78RVGCO2= 60.
79RVGH2O= 30.
80RVGO3 =400.
81!
82!
83!-- REFERENCE TEMPERATURES FOR THE PLANCK FUNCTIONS ----- TSTAND
84!          TEMPERATURE DEPENDENCE OF THE ABSORPTION ----- TREF
85TSTAND= 250.00
86TREF  = 250.00
87!
88
89!-- ROOTS AND WEIGHTS FOR THE 2-POINT GAUSSIAN QUADRATURE
90RT1 = (/ -0.577350269, +0.577350269 /)
91WG1 = (/  1.         ,  1.          /)
92
93!-- TEMPERATURE STEP AND MAXIMUM INDEX FOR THE TRANSMISSIVITIES
94TSTP  = 12.5
95MXIXT = 10
96
97!
98!-- REFERENCE TEMPERATURE FOR THE TRANSMISSIVITIES
99
100TINTP = (/ 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!
106XP( :, 1) = (/&
107 &0.46430621E+02, 0.12928299E+03, 0.20732648E+03,&
108 &0.31398411E+03, 0.18373177E+03,-0.11412303E+03/)
109XP( :, 2) = (/&
110 &0.73604774E+02, 0.27887914E+03, 0.27076947E+03,&
111 &-.57322111E+02,-0.64742459E+02, 0.87238280E+02/)
112XP( :, 3) = (/&
113 &0.37050866E+02, 0.20498759E+03, 0.37558029E+03,&
114 &0.17401171E+03,-0.13350302E+03,-0.37651795E+02/)
115XP( :, 4) = (/&
116 &0.14930141E+02, 0.89161160E+02, 0.17793062E+03,&
117 &0.93433860E+02,-0.70646020E+02,-0.26373150E+02/)
118XP( :, 5) = (/&
119 &0.40386780E+02, 0.10855270E+03, 0.50755010E+02,&
120&-0.31496190E+02, 0.12791300E+00, 0.18017770E+01/)
121XP( :, 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
132AT(1,:)= (/ 0.298199E-02,-.394023E-03,0.319566E-04 /)
133BT(1,:)= (/-0.106432E-04,0.660324E-06,0.174356E-06 /)
134!
135!-- H2O --  500 -  800
136AT(2,:)= (/ 0.143676E-01,0.366501E-02,-.160822E-02 /)
137BT(2,:)= (/-0.553979E-04,-.101701E-04,0.920868E-05 /)
138!
139!-- H2O --  800 - 970 + 1110 - 1250
140AT(3,:)= (/ 0.197861E-01,0.315541E-02,-.174547E-02 /)
141BT(3,:)= (/-0.877012E-04,0.513302E-04,0.523138E-06 /)
142!
143!-- H2O --  970 - 1110
144AT(4,:)= (/ 0.289560E-01,-.208807E-02,-.121943E-02 /)
145BT(4,:)= (/-0.165960E-03,0.157704E-03,-.146427E-04 /)
146!
147!-- H2O --  350 -  500
148AT(5,:)= (/ 0.103800E-01,0.436296E-02,-.161431E-02 /)
149BT(5,:)= (/ -.276744E-04,-.327381E-04,0.127646E-04 /)
150!
151!-- H2O -- 1250 - 1450 + 1880 - 2820
152AT(6,:)= (/ 0.868859E-02,-.972752E-03,0.000000E-00 /)
153BT(6,:)= (/ -.278412E-04,-.713940E-06,0.117469E-05 /)
154!
155
156!-------CO2
157
158AT(7,:)= (/&
159     & 0.250073E-03,0.455875E-03,0.109242E-03 /)
160BT(7,:)= (/&
161     & 0.199846E-05,-.216313E-05,0.175991E-06 /)
162AT(8,:)= (/&
163     & 0.307423E-01,0.110879E-02,-.322172E-03 /)
164BT(8,:)= (/&
165     &-0.108482E-03,0.258096E-05,-.814575E-06 /)
166!
167!--  COEFFICIENTS FOR OZONE ABSORPTION
168!
169O1H = 2230.
170O2H =  100.
171RPIALF0= 2.
172OCT(1)= -.326E-03
173OCT(2)= -.102E-05
174OCT(3)=  .137E-02
175OCT(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
200GA( 1, 1, :) = (/&
201     & 0.63499072E-02,-0.99506586E-03, 0.00000000E+00/)
202GB( 1, 1, :) = (/&
203     & 0.63499072E-02, 0.97222852E-01, 0.10000000E+01/)
204GA( 1, 2, :) = (/&
205     & 0.77266491E-02,-0.11661515E-02, 0.00000000E+00/)
206GB( 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
212GA( 2, 1, :) = (/&
213     & 0.65566348E-02,-0.10184169E-02, 0.00000000E+00/)
214GB( 2, 1, :) = (/&
215     & 0.65566348E-02, 0.98862238E-01, 0.10000000E+01/)
216GA( 2, 2, :) = (/&
217     & 0.81323287E-02,-0.11886130E-02, 0.00000000E+00/)
218GB( 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
224GA( 3, 1, :) = (/&
225     & 0.67849730E-02,-0.10404730E-02, 0.00000000E+00/)
226GB( 3, 1, :) = (/&
227     & 0.67849730E-02, 0.10061504E+00, 0.10000000E+01/)
228GA( 3, 2, :) = (/&
229     & 0.86507620E-02,-0.12139929E-02, 0.00000000E+00/)
230GB( 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
236GA( 4, 1, :) = (/&
237     & 0.70481947E-02,-0.10621792E-02, 0.00000000E+00/)
238GB( 4, 1, :) = (/&
239     & 0.70481947E-02, 0.10256222E+00, 0.10000000E+01/)
240GA( 4, 2, :) = (/&
241     & 0.92776391E-02,-0.12445811E-02, 0.00000000E+00/)
242GB( 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
248GA( 5, 1, :) = (/&
249     & 0.73585943E-02,-0.10847662E-02, 0.00000000E+00/)
250GB( 5, 1, :) = (/&
251     & 0.73585943E-02, 0.10475952E+00, 0.10000000E+01/)
252GA( 5, 2, :) = (/&
253     & 0.99806312E-02,-0.12807672E-02, 0.00000000E+00/)
254GB( 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
260GA( 6, 1, :) = (/&
261     & 0.77242818E-02,-0.11094726E-02, 0.00000000E+00/)
262GB( 6, 1, :) = (/&
263     & 0.77242818E-02, 0.10720986E+00, 0.10000000E+01/)
264GA( 6, 2, :) = (/&
265     & 0.10709803E-01,-0.13208251E-02, 0.00000000E+00/)
266GB( 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
272GA( 7, 1, :) = (/&
273     & 0.81472693E-02,-0.11372949E-02, 0.00000000E+00/)
274GB( 7, 1, :) = (/&
275     & 0.81472693E-02, 0.10985370E+00, 0.10000000E+01/)
276GA( 7, 2, :) = (/&
277     & 0.11414739E-01,-0.13619034E-02, 0.00000000E+00/)
278GB( 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
284GA( 8, 1, :) = (/&
285     & 0.86227527E-02,-0.11687683E-02, 0.00000000E+00/)
286GB( 8, 1, :) = (/&
287     & 0.86227527E-02, 0.11257633E+00, 0.10000000E+01/)
288GA( 8, 2, :) = (/&
289     & 0.12058772E-01,-0.14014165E-02, 0.00000000E+00/)
290GB( 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
296GA( 9, 1, :) = (/&
297     & 0.91396814E-02,-0.12038314E-02, 0.00000000E+00/)
298GB( 9, 1, :) = (/&
299     & 0.91396814E-02, 0.11522980E+00, 0.10000000E+01/)
300GA( 9, 2, :) = (/&
301     & 0.12623992E-01,-0.14378639E-02, 0.00000000E+00/)
302GB( 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
308GA(10, 1, :) = (/&
309     & 0.96825438E-02,-0.12418367E-02, 0.00000000E+00/)
310GB(10, 1, :) = (/&
311     & 0.96825438E-02, 0.11766343E+00, 0.10000000E+01/)
312GA(10, 2, :) = (/&
313     & 0.13108146E-01,-0.14708488E-02, 0.00000000E+00/)
314GB(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
320GA(11, 1, :) = (/&
321     & 0.10233955E-01,-0.12817135E-02, 0.00000000E+00/)
322GB(11, 1, :) = (/&
323     & 0.10233955E-01, 0.11975320E+00, 0.10000000E+01/)
324GA(11, 2, :) = (/&
325     & 0.13518390E-01,-0.15006791E-02, 0.00000000E+00/)
326GB(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
338GA( 1, 3, :) = (/&
339     & 0.11644593E+01, 0.41243390E+00, 0.00000000E+00/)
340GB( 1, 3, :) = (/&
341     & 0.11644593E+01, 0.10346097E+01, 0.10000000E+01/)
342GA( 1, 4, :) = (/&
343     & 0.12006968E+01, 0.48318936E+00, 0.00000000E+00/)
344GB( 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
350GA( 2, 3, :) = (/&
351     & 0.11747203E+01, 0.43407282E+00, 0.00000000E+00/)
352GB( 2, 3, :) = (/&
353     & 0.11747203E+01, 0.10433655E+01, 0.10000000E+01/)
354GA( 2, 4, :) = (/&
355     & 0.12108196E+01, 0.50501827E+00, 0.00000000E+00/)
356GB( 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
362GA( 3, 3, :) = (/&
363     & 0.11837872E+01, 0.45331413E+00, 0.00000000E+00/)
364GB( 3, 3, :) = (/&
365     & 0.11837872E+01, 0.10511933E+01, 0.10000000E+01/)
366GA( 3, 4, :) = (/&
367     & 0.12196717E+01, 0.52409502E+00, 0.00000000E+00/)
368GB( 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
374GA( 4, 3, :) = (/&
375     & 0.11918561E+01, 0.47048604E+00, 0.00000000E+00/)
376GB( 4, 3, :) = (/&
377     & 0.11918561E+01, 0.10582150E+01, 0.10000000E+01/)
378GA( 4, 4, :) = (/&
379     & 0.12274493E+01, 0.54085277E+00, 0.00000000E+00/)
380GB( 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
386GA( 5, 3, :) = (/&
387     & 0.11990757E+01, 0.48586286E+00, 0.00000000E+00/)
388GB( 5, 3, :) = (/&
389     & 0.11990757E+01, 0.10645317E+01, 0.10000000E+01/)
390GA( 5, 4, :) = (/&
391     & 0.12343189E+01, 0.55565422E+00, 0.00000000E+00/)
392GB( 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
398GA( 6, 3, :) = (/&
399     & 0.12055643E+01, 0.49968044E+00, 0.00000000E+00/)
400GB( 6, 3, :) = (/&
401     & 0.12055643E+01, 0.10702313E+01, 0.10000000E+01/)
402GA( 6, 4, :) = (/&
403     & 0.12404147E+01, 0.56878618E+00, 0.00000000E+00/)
404GB( 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
410GA( 7, 3, :) = (/&
411     & 0.12114186E+01, 0.51214132E+00, 0.00000000E+00/)
412GB( 7, 3, :) = (/&
413     & 0.12114186E+01, 0.10753907E+01, 0.10000000E+01/)
414GA( 7, 4, :) = (/&
415     & 0.12458431E+01, 0.58047395E+00, 0.00000000E+00/)
416GB( 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
422GA( 8, 3, :) = (/&
423     & 0.12167192E+01, 0.52341830E+00, 0.00000000E+00/)
424GB( 8, 3, :) = (/&
425     & 0.12167192E+01, 0.10800762E+01, 0.10000000E+01/)
426GA( 8, 4, :) = (/&
427     & 0.12506907E+01, 0.59089894E+00, 0.00000000E+00/)
428GB( 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
434GA( 9, 3, :) = (/&
435     & 0.12215344E+01, 0.53365803E+00, 0.00000000E+00/)
436GB( 9, 3, :) = (/&
437     & 0.12215344E+01, 0.10843446E+01, 0.10000000E+01/)
438GA( 9, 4, :) = (/&
439     & 0.12550299E+01, 0.60021475E+00, 0.00000000E+00/)
440GB( 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
446GA(10, 3, :) = (/&
447     & 0.12259226E+01, 0.54298448E+00, 0.00000000E+00/)
448GB(10, 3, :) = (/&
449     & 0.12259226E+01, 0.10882439E+01, 0.10000000E+01/)
450GA(10, 4, :) = (/&
451     & 0.12589256E+01, 0.60856112E+00, 0.00000000E+00/)
452GB(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
458GA(11, 3, :) = (/&
459     & 0.12299344E+01, 0.55150227E+00, 0.00000000E+00/)
460GB(11, 3, :) = (/&
461     & 0.12299344E+01, 0.10918144E+01, 0.10000000E+01/)
462GA(11, 4, :) = (/&
463     & 0.12624402E+01, 0.61607594E+00, 0.00000000E+00/)
464GB(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
477GA( 1, 5, :) = (/&
478     & 0.15750172E+00,-0.22159303E-01, 0.00000000E+00/)
479GB( 1, 5, :) = (/&
480     & 0.15750172E+00, 0.38103212E+00, 0.10000000E+01/)
481GA( 1, 6, :) = (/&
482     & 0.17770551E+00,-0.24972399E-01, 0.00000000E+00/)
483GB( 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
489GA( 2, 5, :) = (/&
490     & 0.16174076E+00,-0.22748917E-01, 0.00000000E+00/)
491GB( 2, 5, :) = (/&
492     & 0.16174076E+00, 0.38913800E+00, 0.10000000E+01/)
493GA( 2, 6, :) = (/&
494     & 0.18176757E+00,-0.25537247E-01, 0.00000000E+00/)
495GB( 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
501GA( 3, 5, :) = (/&
502     & 0.16548628E+00,-0.23269898E-01, 0.00000000E+00/)
503GB( 3, 5, :) = (/&
504     & 0.16548628E+00, 0.39613651E+00, 0.10000000E+01/)
505GA( 3, 6, :) = (/&
506     & 0.18527967E+00,-0.26025624E-01, 0.00000000E+00/)
507GB( 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
513GA( 4, 5, :) = (/&
514     & 0.16881124E+00,-0.23732392E-01, 0.00000000E+00/)
515GB( 4, 5, :) = (/&
516     & 0.16881124E+00, 0.40222421E+00, 0.10000000E+01/)
517GA( 4, 6, :) = (/&
518     & 0.18833348E+00,-0.26450280E-01, 0.00000000E+00/)
519GB( 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
525GA( 5, 5, :) = (/&
526     & 0.17177839E+00,-0.24145123E-01, 0.00000000E+00/)
527GB( 5, 5, :) = (/&
528     & 0.17177839E+00, 0.40756010E+00, 0.10000000E+01/)
529GA( 5, 6, :) = (/&
530     & 0.19100108E+00,-0.26821236E-01, 0.00000000E+00/)
531GB( 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
537GA( 6, 5, :) = (/&
538     & 0.17443933E+00,-0.24515269E-01, 0.00000000E+00/)
539GB( 6, 5, :) = (/&
540     & 0.17443933E+00, 0.41226954E+00, 0.10000000E+01/)
541GA( 6, 6, :) = (/&
542     & 0.19334122E+00,-0.27146657E-01, 0.00000000E+00/)
543GB( 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
549GA( 7, 5, :) = (/&
550     & 0.17683622E+00,-0.24848690E-01, 0.00000000E+00/)
551GB( 7, 5, :) = (/&
552     & 0.17683622E+00, 0.41645142E+00, 0.10000000E+01/)
553GA( 7, 6, :) = (/&
554     & 0.19540288E+00,-0.27433354E-01, 0.00000000E+00/)
555GB( 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
561GA( 8, 5, :) = (/&
562     & 0.17900375E+00,-0.25150210E-01, 0.00000000E+00/)
563GB( 8, 5, :) = (/&
564     & 0.17900375E+00, 0.42018474E+00, 0.10000000E+01/)
565GA( 8, 6, :) = (/&
566     & 0.19722732E+00,-0.27687065E-01, 0.00000000E+00/)
567GB( 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
573GA( 9, 5, :) = (/&
574     & 0.18097099E+00,-0.25423873E-01, 0.00000000E+00/)
575GB( 9, 5, :) = (/&
576     & 0.18097099E+00, 0.42353379E+00, 0.10000000E+01/)
577GA( 9, 6, :) = (/&
578     & 0.19884918E+00,-0.27912608E-01, 0.00000000E+00/)
579GB( 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
585GA(10, 5, :) = (/&
586     & 0.18276283E+00,-0.25673139E-01, 0.00000000E+00/)
587GB(10, 5, :) = (/&
588     & 0.18276283E+00, 0.42655211E+00, 0.10000000E+01/)
589GA(10, 6, :) = (/&
590     & 0.20029696E+00,-0.28113944E-01, 0.00000000E+00/)
591GB(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
597GA(11, 5, :) = (/&
598     & 0.18440117E+00,-0.25901055E-01, 0.00000000E+00/)
599GB(11, 5, :) = (/&
600     & 0.18440117E+00, 0.42928533E+00, 0.10000000E+01/)
601GA(11, 6, :) = (/&
602     & 0.20159300E+00,-0.28294180E-01, 0.00000000E+00/)
603GB(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
617GA( 1, 7, :) = (/&
618     & 0.10192131E+02, 0.80737799E+01, 0.00000000E+00/)
619GB( 1, 7, :) = (/&
620     & 0.10192131E+02, 0.82623280E+01, 0.10000000E+01/)
621GA( 1, 8, :) = (/&
622     & 0.92439050E+01, 0.77425778E+01, 0.00000000E+00/)
623GB( 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
629GA( 2, 7, :) = (/&
630     & 0.97258602E+01, 0.79171158E+01, 0.00000000E+00/)
631GB( 2, 7, :) = (/&
632     & 0.97258602E+01, 0.81072291E+01, 0.10000000E+01/)
633GA( 2, 8, :) = (/&
634     & 0.87567422E+01, 0.75443460E+01, 0.00000000E+00/)
635GB( 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
641GA( 3, 7, :) = (/&
642     & 0.92992890E+01, 0.77609605E+01, 0.00000000E+00/)
643GB( 3, 7, :) = (/&
644     & 0.92992890E+01, 0.79523834E+01, 0.10000000E+01/)
645GA( 3, 8, :) = (/&
646     & 0.83270144E+01, 0.73526151E+01, 0.00000000E+00/)
647GB( 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
653GA( 4, 7, :) = (/&
654     & 0.89154021E+01, 0.76087371E+01, 0.00000000E+00/)
655GB( 4, 7, :) = (/&
656     & 0.89154021E+01, 0.78012527E+01, 0.10000000E+01/)
657GA( 4, 8, :) = (/&
658     & 0.79528337E+01, 0.71711188E+01, 0.00000000E+00/)
659GB( 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
665GA( 5, 7, :) = (/&
666     & 0.85730084E+01, 0.74627112E+01, 0.00000000E+00/)
667GB( 5, 7, :) = (/&
668     & 0.85730084E+01, 0.76561458E+01, 0.10000000E+01/)
669GA( 5, 8, :) = (/&
670     & 0.76286839E+01, 0.70015571E+01, 0.00000000E+00/)
671GB( 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
677GA( 6, 7, :) = (/&
678     & 0.82685838E+01, 0.73239981E+01, 0.00000000E+00/)
679GB( 6, 7, :) = (/&
680     & 0.82685838E+01, 0.75182174E+01, 0.10000000E+01/)
681GA( 6, 8, :) = (/&
682     & 0.73477879E+01, 0.68442532E+01, 0.00000000E+00/)
683GB( 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
689GA( 7, 7, :) = (/&
690     & 0.79978921E+01, 0.71929934E+01, 0.00000000E+00/)
691GB( 7, 7, :) = (/&
692     & 0.79978921E+01, 0.73878952E+01, 0.10000000E+01/)
693GA( 7, 8, :) = (/&
694     & 0.71035818E+01, 0.66987996E+01, 0.00000000E+00/)
695GB( 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
701GA( 8, 7, :) = (/&
702     & 0.77568055E+01, 0.70697065E+01, 0.00000000E+00/)
703GB( 8, 7, :) = (/&
704     & 0.77568055E+01, 0.72652133E+01, 0.10000000E+01/)
705GA( 8, 8, :) = (/&
706     & 0.68903312E+01, 0.65644820E+01, 0.00000000E+00/)
707GB( 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
713GA( 9, 7, :) = (/&
714     & 0.75416266E+01, 0.69539626E+01, 0.00000000E+00/)
715GB( 9, 7, :) = (/&
716     & 0.75416266E+01, 0.71500151E+01, 0.10000000E+01/)
717GA( 9, 8, :) = (/&
718     & 0.67032875E+01, 0.64405267E+01, 0.00000000E+00/)
719GB( 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
725GA(10, 7, :) = (/&
726     & 0.73491694E+01, 0.68455144E+01, 0.00000000E+00/)
727GB(10, 7, :) = (/&
728     & 0.73491694E+01, 0.70420667E+01, 0.10000000E+01/)
729GA(10, 8, :) = (/&
730     & 0.65386461E+01, 0.63262376E+01, 0.00000000E+00/)
731GB(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
737GA(11, 7, :) = (/&
738     & 0.71767400E+01, 0.67441020E+01, 0.00000000E+00/)
739GB(11, 7, :) = (/&
740     & 0.71767400E+01, 0.69411177E+01, 0.10000000E+01/)
741GA(11, 8, :) = (/&
742     & 0.63934377E+01, 0.62210701E+01, 0.00000000E+00/)
743GB(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
755GA( 1, 9, :) = (/&
756     & 0.24870635E+02, 0.10542131E+02, 0.00000000E+00/)
757GB( 1, 9, :) = (/&
758     & 0.24870635E+02, 0.10656640E+02, 0.10000000E+01/)
759GA( 1,10, :) = (/&
760     & 0.24586283E+02, 0.10490353E+02, 0.00000000E+00/)
761GB( 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
767GA( 2, 9, :) = (/&
768     & 0.24725591E+02, 0.10515895E+02, 0.00000000E+00/)
769GB( 2, 9, :) = (/&
770     & 0.24725591E+02, 0.10630910E+02, 0.10000000E+01/)
771GA( 2,10, :) = (/&
772     & 0.24441465E+02, 0.10463512E+02, 0.00000000E+00/)
773GB( 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
779GA( 3, 9, :) = (/&
780     & 0.24600320E+02, 0.10492949E+02, 0.00000000E+00/)
781GB( 3, 9, :) = (/&
782     & 0.24600320E+02, 0.10608399E+02, 0.10000000E+01/)
783GA( 3,10, :) = (/&
784     & 0.24311657E+02, 0.10439183E+02, 0.00000000E+00/)
785GB( 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
791GA( 4, 9, :) = (/&
792     & 0.24487300E+02, 0.10472049E+02, 0.00000000E+00/)
793GB( 4, 9, :) = (/&
794     & 0.24487300E+02, 0.10587891E+02, 0.10000000E+01/)
795GA( 4,10, :) = (/&
796     & 0.24196167E+02, 0.10417324E+02, 0.00000000E+00/)
797GB( 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
803GA( 5, 9, :) = (/&
804     & 0.24384935E+02, 0.10452961E+02, 0.00000000E+00/)
805GB( 5, 9, :) = (/&
806     & 0.24384935E+02, 0.10569156E+02, 0.10000000E+01/)
807GA( 5,10, :) = (/&
808     & 0.24093406E+02, 0.10397704E+02, 0.00000000E+00/)
809GB( 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
815GA( 6, 9, :) = (/&
816     & 0.24292341E+02, 0.10435562E+02, 0.00000000E+00/)
817GB( 6, 9, :) = (/&
818     & 0.24292341E+02, 0.10552075E+02, 0.10000000E+01/)
819GA( 6,10, :) = (/&
820     & 0.24001597E+02, 0.10380038E+02, 0.00000000E+00/)
821GB( 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
827GA( 7, 9, :) = (/&
828     & 0.24208572E+02, 0.10419710E+02, 0.00000000E+00/)
829GB( 7, 9, :) = (/&
830     & 0.24208572E+02, 0.10536510E+02, 0.10000000E+01/)
831GA( 7,10, :) = (/&
832     & 0.23919098E+02, 0.10364052E+02, 0.00000000E+00/)
833GB( 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
839GA( 8, 9, :) = (/&
840     & 0.24132642E+02, 0.10405247E+02, 0.00000000E+00/)
841GB( 8, 9, :) = (/&
842     & 0.24132642E+02, 0.10522307E+02, 0.10000000E+01/)
843GA( 8,10, :) = (/&
844     & 0.23844511E+02, 0.10349509E+02, 0.00000000E+00/)
845GB( 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
851GA( 9, 9, :) = (/&
852     & 0.24063614E+02, 0.10392022E+02, 0.00000000E+00/)
853GB( 9, 9, :) = (/&
854     & 0.24063614E+02, 0.10509317E+02, 0.10000000E+01/)
855GA( 9,10, :) = (/&
856     & 0.23776708E+02, 0.10336215E+02, 0.00000000E+00/)
857GB( 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
863GA(10, 9, :) = (/&
864     & 0.24000649E+02, 0.10379892E+02, 0.00000000E+00/)
865GB(10, 9, :) = (/&
866     & 0.24000649E+02, 0.10497402E+02, 0.10000000E+01/)
867GA(10,10, :) = (/&
868     & 0.23714816E+02, 0.10324018E+02, 0.00000000E+00/)
869GB(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
875GA(11, 9, :) = (/&
876     & 0.23943021E+02, 0.10368736E+02, 0.00000000E+00/)
877GB(11, 9, :) = (/&
878     & 0.23943021E+02, 0.10486443E+02, 0.10000000E+01/)
879GA(11,10, :) = (/&
880     & 0.23658197E+02, 0.10312808E+02, 0.00000000E+00/)
881GB(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
894GA( 1,11, :) = (/&
895     & 0.11990218E+02,-0.12823142E+01, 0.00000000E+00/)
896GB( 1,11, :) = (/&
897     & 0.11990218E+02, 0.26681588E+02, 0.10000000E+01/)
898GA( 1,12, :) = (/&
899     & 0.79709806E+01,-0.74805226E+00, 0.00000000E+00/)
900GB( 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
906GA( 2,11, :) = (/&
907     & 0.10904073E+02,-0.10571588E+01, 0.00000000E+00/)
908GB( 2,11, :) = (/&
909     & 0.10904073E+02, 0.24728346E+02, 0.10000000E+01/)
910GA( 2,12, :) = (/&
911     & 0.75400737E+01,-0.56252739E+00, 0.00000000E+00/)
912GB( 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
918GA( 3,11, :) = (/&
919     & 0.89126838E+01,-0.74864953E+00, 0.00000000E+00/)
920GB( 3,11, :) = (/&
921     & 0.89126838E+01, 0.20551342E+02, 0.10000000E+01/)
922GA( 3,12, :) = (/&
923     & 0.81804377E+01,-0.46188072E+00, 0.00000000E+00/)
924GB( 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
930GA( 4,11, :) = (/&
931     & 0.85622405E+01,-0.58705980E+00, 0.00000000E+00/)
932GB( 4,11, :) = (/&
933     & 0.85622405E+01, 0.19955244E+02, 0.10000000E+01/)
934GA( 4,12, :) = (/&
935     & 0.10564339E+02,-0.40712065E+00, 0.00000000E+00/)
936GB( 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
942GA( 5,11, :) = (/&
943     & 0.94892164E+01,-0.49305772E+00, 0.00000000E+00/)
944GB( 5,11, :) = (/&
945     & 0.94892164E+01, 0.22227100E+02, 0.10000000E+01/)
946GA( 5,12, :) = (/&
947     & 0.46896789E+02,-0.15295996E+01, 0.00000000E+00/)
948GB( 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
954GA( 6,11, :) = (/&
955     & 0.13580937E+02,-0.51461431E+00, 0.00000000E+00/)
956GB( 6,11, :) = (/&
957     & 0.13580937E+02, 0.31770288E+02, 0.10000000E+01/)
958GA( 6,12, :) = (/&
959     &-0.30926524E+01, 0.43555255E+00, 0.00000000E+00/)
960GB( 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
966GA( 7,11, :) = (/&
967     &-0.32050918E+03, 0.12373350E+02, 0.00000000E+00/)
968GB( 7,11, :) = (/&
969     &-0.32050918E+03,-0.74061287E+03, 0.10000000E+01/)
970GA( 7,12, :) = (/&
971     & 0.85742941E+00, 0.50380874E+00, 0.00000000E+00/)
972GB( 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
978GA( 8,11, :) = (/&
979     &-0.37133165E+01, 0.44809588E+00, 0.00000000E+00/)
980GB( 8,11, :) = (/&
981     &-0.37133165E+01,-0.81329826E+01, 0.10000000E+01/)
982GA( 8,12, :) = (/&
983     & 0.19164038E+01, 0.68537352E+00, 0.00000000E+00/)
984GB( 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
990GA( 9,11, :) = (/&
991     & 0.18890836E+00, 0.46548918E+00, 0.00000000E+00/)
992GB( 9,11, :) = (/&
993     & 0.18890836E+00, 0.90279822E+00, 0.10000000E+01/)
994GA( 9,12, :) = (/&
995     & 0.23513199E+01, 0.89437630E+00, 0.00000000E+00/)
996GB( 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
1002GA(10,11, :) = (/&
1003     & 0.14209226E+01, 0.59121475E+00, 0.00000000E+00/)
1004GB(10,11, :) = (/&
1005     & 0.14209226E+01, 0.37532746E+01, 0.10000000E+01/)
1006GA(10,12, :) = (/&
1007     & 0.25566644E+01, 0.11127003E+01, 0.00000000E+00/)
1008GB(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
1014GA(11,11, :) = (/&
1015     & 0.19817679E+01, 0.74676119E+00, 0.00000000E+00/)
1016GB(11,11, :) = (/&
1017     & 0.19817679E+01, 0.50437916E+01, 0.10000000E+01/)
1018GA(11,12, :) = (/&
1019     & 0.26555181E+01, 0.13329782E+01, 0.00000000E+00/)
1020GB(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
1037GA( 1,13, :) = (/&
1038     & 0.87668459E-01, 0.13845511E+01, 0.00000000E+00/)
1039GB( 1,13, :) = (/&
1040     & 0.87668459E-01, 0.23203798E+01, 0.10000000E+01/)
1041GA( 1,14, :) = (/&
1042     & 0.74878820E-01, 0.11718758E+01, 0.00000000E+00/)
1043GB( 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
1049GA( 2,13, :) = (/&
1050     & 0.83754276E-01, 0.13187042E+01, 0.00000000E+00/)
1051GB( 2,13, :) = (/&
1052     & 0.83754276E-01, 0.22288925E+01, 0.10000000E+01/)
1053GA( 2,14, :) = (/&
1054     & 0.71650966E-01, 0.11216131E+01, 0.00000000E+00/)
1055GB( 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
1061GA( 3,13, :) = (/&
1062     & 0.80460283E-01, 0.12644396E+01, 0.00000000E+00/)
1063GB( 3,13, :) = (/&
1064     & 0.80460283E-01, 0.21515593E+01, 0.10000000E+01/)
1065GA( 3,14, :) = (/&
1066     & 0.68979615E-01, 0.10809473E+01, 0.00000000E+00/)
1067GB( 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
1073GA( 4,13, :) = (/&
1074     & 0.77659686E-01, 0.12191543E+01, 0.00000000E+00/)
1075GB( 4,13, :) = (/&
1076     & 0.77659686E-01, 0.20855896E+01, 0.10000000E+01/)
1077GA( 4,14, :) = (/&
1078     & 0.66745345E-01, 0.10476396E+01, 0.00000000E+00/)
1079GB( 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
1085GA( 5,13, :) = (/&
1086     & 0.75257056E-01, 0.11809511E+01, 0.00000000E+00/)
1087GB( 5,13, :) = (/&
1088     & 0.75257056E-01, 0.20288489E+01, 0.10000000E+01/)
1089GA( 5,14, :) = (/&
1090     & 0.64857571E-01, 0.10200373E+01, 0.00000000E+00/)
1091GB( 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
1097GA( 6,13, :) = (/&
1098     & 0.73179175E-01, 0.11484154E+01, 0.00000000E+00/)
1099GB( 6,13, :) = (/&
1100     & 0.73179175E-01, 0.19796791E+01, 0.10000000E+01/)
1101GA( 6,14, :) = (/&
1102     & 0.63248495E-01, 0.99692726E+00, 0.00000000E+00/)
1103GB( 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
1109GA( 7,13, :) = (/&
1110     & 0.71369063E-01, 0.11204723E+01, 0.00000000E+00/)
1111GB( 7,13, :) = (/&
1112     & 0.71369063E-01, 0.19367778E+01, 0.10000000E+01/)
1113GA( 7,14, :) = (/&
1114     & 0.61866970E-01, 0.97740923E+00, 0.00000000E+00/)
1115GB( 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
1121GA( 8,13, :) = (/&
1122     & 0.69781812E-01, 0.10962918E+01, 0.00000000E+00/)
1123GB( 8,13, :) = (/&
1124     & 0.69781812E-01, 0.18991112E+01, 0.10000000E+01/)
1125GA( 8,14, :) = (/&
1126     & 0.60673632E-01, 0.96080188E+00, 0.00000000E+00/)
1127GB( 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
1133GA( 9,13, :) = (/&
1134     & 0.68381606E-01, 0.10752229E+01, 0.00000000E+00/)
1135GB( 9,13, :) = (/&
1136     & 0.68381606E-01, 0.18658501E+01, 0.10000000E+01/)
1137GA( 9,14, :) = (/&
1138     & 0.59637277E-01, 0.94657562E+00, 0.00000000E+00/)
1139GB( 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
1145GA(10,13, :) = (/&
1146     & 0.67139539E-01, 0.10567474E+01, 0.00000000E+00/)
1147GB(10,13, :) = (/&
1148     & 0.67139539E-01, 0.18363226E+01, 0.10000000E+01/)
1149GA(10,14, :) = (/&
1150     & 0.58732178E-01, 0.93430511E+00, 0.00000000E+00/)
1151GB(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
1157GA(11,13, :) = (/&
1158     & 0.66032012E-01, 0.10404465E+01, 0.00000000E+00/)
1159GB(11,13, :) = (/&
1160     & 0.66032012E-01, 0.18099779E+01, 0.10000000E+01/)
1161GA(11,14, :) = (/&
1162     & 0.57936092E-01, 0.92363528E+00, 0.00000000E+00/)
1163GB(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
1184GA( 1,15, :) = (/&
1185     & 0.13230067E+02, 0.22042132E+02, 0.00000000E+00/)
1186GB( 1,15, :) = (/&
1187     & 0.13230067E+02, 0.22051750E+02, 0.10000000E+01/)
1188GA( 1,16, :) = (/&
1189     & 0.13183816E+02, 0.22169501E+02, 0.00000000E+00/)
1190GB( 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
1196GA( 2,15, :) = (/&
1197     & 0.13213564E+02, 0.22107298E+02, 0.00000000E+00/)
1198GB( 2,15, :) = (/&
1199     & 0.13213564E+02, 0.22116850E+02, 0.10000000E+01/)
1200GA( 2,16, :) = (/&
1201     & 0.13189991E+02, 0.22270075E+02, 0.00000000E+00/)
1202GB( 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
1208GA( 3,15, :) = (/&
1209     & 0.13209140E+02, 0.22180915E+02, 0.00000000E+00/)
1210GB( 3,15, :) = (/&
1211     & 0.13209140E+02, 0.22190410E+02, 0.10000000E+01/)
1212GA( 3,16, :) = (/&
1213     & 0.13209485E+02, 0.22379193E+02, 0.00000000E+00/)
1214GB( 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
1220GA( 4,15, :) = (/&
1221     & 0.13213894E+02, 0.22259478E+02, 0.00000000E+00/)
1222GB( 4,15, :) = (/&
1223     & 0.13213894E+02, 0.22268925E+02, 0.10000000E+01/)
1224GA( 4,16, :) = (/&
1225     & 0.13238789E+02, 0.22492992E+02, 0.00000000E+00/)
1226GB( 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
1232GA( 5,15, :) = (/&
1233     & 0.13225963E+02, 0.22341039E+02, 0.00000000E+00/)
1234GB( 5,15, :) = (/&
1235     & 0.13225963E+02, 0.22350445E+02, 0.10000000E+01/)
1236GA( 5,16, :) = (/&
1237     & 0.13275017E+02, 0.22608508E+02, 0.00000000E+00/)
1238GB( 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
1244GA( 6,15, :) = (/&
1245     & 0.13243806E+02, 0.22424247E+02, 0.00000000E+00/)
1246GB( 6,15, :) = (/&
1247     & 0.13243806E+02, 0.22433617E+02, 0.10000000E+01/)
1248GA( 6,16, :) = (/&
1249     & 0.13316096E+02, 0.22723843E+02, 0.00000000E+00/)
1250GB( 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
1256GA( 7,15, :) = (/&
1257     & 0.13266104E+02, 0.22508089E+02, 0.00000000E+00/)
1258GB( 7,15, :) = (/&
1259     & 0.13266104E+02, 0.22517429E+02, 0.10000000E+01/)
1260GA( 7,16, :) = (/&
1261     & 0.13360555E+02, 0.22837837E+02, 0.00000000E+00/)
1262GB( 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
1268GA( 8,15, :) = (/&
1269     & 0.13291782E+02, 0.22591771E+02, 0.00000000E+00/)
1270GB( 8,15, :) = (/&
1271     & 0.13291782E+02, 0.22601086E+02, 0.10000000E+01/)
1272GA( 8,16, :) = (/&
1273     & 0.13407324E+02, 0.22949751E+02, 0.00000000E+00/)
1274GB( 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
1280GA( 9,15, :) = (/&
1281     & 0.13319961E+02, 0.22674661E+02, 0.00000000E+00/)
1282GB( 9,15, :) = (/&
1283     & 0.13319961E+02, 0.22683956E+02, 0.10000000E+01/)
1284GA( 9,16, :) = (/&
1285     & 0.13455544E+02, 0.23059032E+02, 0.00000000E+00/)
1286GB( 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
1292GA(10,15, :) = (/&
1293     & 0.13349927E+02, 0.22756246E+02, 0.00000000E+00/)
1294GB(10,15, :) = (/&
1295     & 0.13349927E+02, 0.22765522E+02, 0.10000000E+01/)
1296GA(10,16, :) = (/&
1297     & 0.13504450E+02, 0.23165146E+02, 0.00000000E+00/)
1298GB(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
1304GA(11,15, :) = (/&
1305     & 0.13381108E+02, 0.22836093E+02, 0.00000000E+00/)
1306GB(11,15, :) = (/&
1307     & 0.13381108E+02, 0.22845354E+02, 0.10000000E+01/)
1308GA(11,16, :) = (/&
1309     & 0.13553282E+02, 0.23267456E+02, 0.00000000E+00/)
1310GB(11,16, :) = (/&
1311     & 0.13553282E+02, 0.23276638E+02, 0.10000000E+01/)
1312!     
1313RETURN
1314END SUBROUTINE SUOLW
Note: See TracBrowser for help on using the repository browser.