/[lmdze]/trunk/libf/phylmd/Radlwsw/lwb.f
ViewVC logotype

Annotation of /trunk/libf/phylmd/Radlwsw/lwb.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 71 - (hide annotations)
Mon Jul 8 18:12:18 2013 UTC (10 years, 10 months ago) by guez
File size: 49691 byte(s) 
No reason to call inidissip in ce0l.

In inidissip, set random seed to 1 beacuse PGI compiler does not
accept all zeros.

dq was computed needlessly in caladvtrac. Arguments masse and dq of
calfis not used.

Replaced real*8 by double precision.

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