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

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

Parent Directory Parent Directory | Revision Log Revision Log

Revision 155 - (hide annotations)
Wed Jul 8 17:03:45 2015 UTC (8 years, 11 months ago) by guez
File size: 50320 byte(s) 
Do not write any longer to startphy.nc nor read from restartphy.nc the
NetCDF variable ALBLW: it was the same than ALBE. ALBE was for the
visible and ALBLW for the near infrared. In physiq, use only variables
falbe and albsol, removed falblw and albsollw. See revision 888 of
LMDZ.

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