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