/[lmdze]/trunk/phylmd/Radlwsw/lwb.f90
ViewVC logotype

Annotation of /trunk/phylmd/Radlwsw/lwb.f90

Parent Directory Parent Directory | Revision Log Revision Log

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