/[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 24 - (hide annotations)
Wed Mar 3 13:23:49 2010 UTC (14 years, 2 months ago) by guez
File size: 49461 byte(s) 
Created directory "phylmd/Radlwsw". Split "radlwsw.f" in files
containing a single procedure.

Removed variable "itaufinp1" in "leapfrog".
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     REAL*8 PDT0(KDLON)
61     REAL*8 PTAVE(KDLON,KFLEV)
62     REAL*8 PTL(KDLON,KFLEV+1)
63     C
64     REAL*8 PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF LEVEL PLANCK FUNCTION
65     REAL*8 PBINT(KDLON,KFLEV+1) ! HALF LEVEL PLANCK FUNCTION
66     REAL*8 PBSUIN(KDLON) ! SURFACE PLANCK FUNCTION
67     REAL*8 PBSUR(KDLON,Ninter) ! SURFACE SPECTRAL PLANCK FUNCTION
68     REAL*8 PBTOP(KDLON,Ninter) ! TOP SPECTRAL PLANCK FUNCTION
69     REAL*8 PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
70     REAL*8 PGA(KDLON,8,2,KFLEV) ! dB/dT-weighted LAYER PADE APPROXIMANTS
71     REAL*8 PGB(KDLON,8,2,KFLEV) ! dB/dT-weighted LAYER PADE APPROXIMANTS
72     REAL*8 PGASUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
73     REAL*8 PGBSUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
74     REAL*8 PGATOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
75     REAL*8 PGBTOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
76     C
77     C-------------------------------------------------------------------------
78     C* LOCAL VARIABLES:
79     INTEGER INDB(KDLON),INDS(KDLON)
80     REAL*8 ZBLAY(KDLON,KFLEV),ZBLEV(KDLON,KFLEV+1)
81     REAL*8 ZRES(KDLON),ZRES2(KDLON),ZTI(KDLON),ZTI2(KDLON)
82     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     REAL*8 zdsto1, zdstox, zdst1, zdstx
88     c
89     C* Quelques parametres:
90     REAL*8 TSTAND
91     PARAMETER (TSTAND=250.0)
92     REAL*8 TSTP
93     PARAMETER (TSTP=12.5)
94     INTEGER MXIXT
95     PARAMETER (MXIXT=10)
96     C
97     C* Used Data Block:
98     REAL*8 TINTP(11)
99     SAVE TINTP
100     REAL*8 GA(11,16,3), GB(11,16,3)
101     SAVE GA, GB
102     REAL*8 XP(6,6)
103     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