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

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

Parent Directory Parent Directory | Revision Log Revision Log

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

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

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

Replaced real*8 by double precision.

Pass arrays with inverted order of vertical levels to conflx instead
of creating local variables for this inside conflx.
1 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 DOUBLE PRECISION PDT0(KDLON)
61 DOUBLE PRECISION PTAVE(KDLON,KFLEV)
62 DOUBLE PRECISION PTL(KDLON,KFLEV+1)
63 C
64 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 C
77 C-------------------------------------------------------------------------
78 C* LOCAL VARIABLES:
79 INTEGER INDB(KDLON),INDS(KDLON)
80 DOUBLE PRECISION ZBLAY(KDLON,KFLEV),ZBLEV(KDLON,KFLEV+1)
81 DOUBLE PRECISION 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 DOUBLE PRECISION zdsto1, zdstox, zdst1, zdstx
88 c
89 C* Quelques parametres:
90 DOUBLE PRECISION TSTAND
91 PARAMETER (TSTAND=250.0)
92 DOUBLE PRECISION TSTP
93 PARAMETER (TSTP=12.5)
94 INTEGER MXIXT
95 PARAMETER (MXIXT=10)
96 C
97 C* Used Data Block:
98 DOUBLE PRECISION TINTP(11)
99 SAVE TINTP
100 DOUBLE PRECISION GA(11,16,3), GB(11,16,3)
101 SAVE GA, GB
102 DOUBLE PRECISION 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