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

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