This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / GEANT321 / gphys / gphxin.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1995/10/24 10:21:30  cernlib
6 * Geant
7 *
8 *
9 #include "geant321/pilot.h"
10 *CMZ :  3.21/02 29/03/94  15.41.22  by  S.Giani
11 *-- Author :
12       SUBROUTINE GPHXIN
13 C.
14 C.    ******************************************************************
15 C.    *                                                                *
16 C     *       Initialize x-section constants for                       *
17 C.    *       the photoelectric effect                                 *
18 C.    *                                                                *
19 C.    *       The constants are from SANDIA REPORT SAND87-0070 UC-34   *
20 C.    *       by F. Biggs and R. Lighthill                             *
21 C.    *                                                                *
22 C.    *    ==>Called by : GPHINI                                       *
23 C.    *       Author    J. Chwastowski                                 *
24 C.    *                                                                *
25 C.    ******************************************************************
26 C.
27       IMPLICIT NONE
28 #include "geant321/gcpmxz.inc"
29 #include "geant321/gcphxs.inc"
30 #include "geant321/gc10ev.inc"
31 #include "geant321/gconsp.inc"
32       INTEGER I,J,K
33       REAL TMP,TMP1
34       CHARACTER*6 CTMP(MAXINT,MAXELZ)
35       DIMENSION TMP(MAXPOW,MAXINT,MAXELZ),TMP1(MAXELZ)
36       DATA (CTMP(I,  1),I=1,MAXINT)
37      +/ 'K     ','.1    ','.8    ','4.    ','20.   ','100.  '
38      +,'500.  ','I     ',5*'      ' /
39       DATA (CTMP(I,  2),I=1,MAXINT)
40      +/ 'K     ','.15   ','.8    ','4.    ','20.   ','100.  '
41      +,'500.  ','I     ',5*'      ' /
42       DATA (CTMP(I,  3),I=1,MAXINT)
43      +/ 'K     ','.8    ','4.    ','20.   ','100.  ','500.  '
44      +,'I     ',6*'      ' /
45       DATA (CTMP(I,  4),I=1,MAXINT)
46      +/ 'K     ','.8    ','4.    ','20.   ','100.  ','500.  '
47      +,'I     ',6*'      ' /
48       DATA (CTMP(I,  5),I=1,MAXINT)
49      +/ 'K     ','.8    ','4.    ','20.   ','100.  ','500.  '
50      +,'I     ',6*'      ' /
51       DATA (CTMP(I,  6),I=1,MAXINT)
52      +/ '.0457 ','K     ','.8    ','4.    ','20.   ','100.  '
53      +,'500.  ','I     ',5*'      ' /
54       DATA (CTMP(I,  7),I=1,MAXINT)
55      +/ '.0404 ','K     ','.8    ','4.    ','20.   ','100.  '
56      +,'500.  ','I     ',5*'      ' /
57       DATA (CTMP(I,  8),I=1,MAXINT)
58      +/ '.0483 ','K     ','4.    ','20.   ','100.  ','500.  '
59      +,'I     ',6*'      ' /
60       DATA (CTMP(I,  9),I=1,MAXINT)
61      +/ '.0724 ','K     ','4.    ','20.   ','100.  ','500.  '
62      +,'I     ',6*'      ' /
63       DATA (CTMP(I, 10),I=1,MAXINT)
64      +/ 'L3    ','L1    ','.124  ','K     ','4.    ','20.   '
65      +,'100.  ','500.  ','I     ',4*'      ' /
66       DATA (CTMP(I, 11),I=1,MAXINT)
67      +/ 'L3    ','L1    ','K     ','20.   ','100.  ','500.  '
68      +,'I     ',6*'      ' /
69       DATA (CTMP(I, 12),I=1,MAXINT)
70      +/ 'L3    ','.123  ','K     ','20.   ','100.  ','500.  '
71      +,'I     ',6*'      ' /
72       DATA (CTMP(I, 13),I=1,MAXINT)
73      +/ '.0159 ','L3    ','L1    ','K     ','20.   ','100.  '
74      +,'500.  ','I     ',5*'      ' /
75       DATA (CTMP(I, 14),I=1,MAXINT)
76      +/ '.0305 ','L3    ','K     ','20.   ','100.  ','500.  '
77      +,'I     ',6*'      ' /
78       DATA (CTMP(I, 15),I=1,MAXINT)
79      +/ '.0306 ','L3    ','K     ','20.   ','100.  ','500.  '
80      +,'I     ',6*'      ' /
81       DATA (CTMP(I, 16),I=1,MAXINT)
82      +/ '.0341 ','L3    ','K     ','20.   ','100.  ','500.  '
83      +,'I     ',6*'      ' /
84       DATA (CTMP(I, 17),I=1,MAXINT)
85      +/ '.0436 ','L3    ','K     ','20.   ','100.  ','500.  '
86      +,'I     ',6*'      ' /
87       DATA (CTMP(I, 18),I=1,MAXINT)
88      +/ '.0362 ','L3    ','K     ','20.   ','100.  ','500.  '
89      +,'I     ',6*'      ' /
90       DATA (CTMP(I, 19),I=1,MAXINT)
91      +/ '.0305 ','L3    ','K     ','20.   ','100.  ','500.  '
92      +,'I     ',6*'      ' /
93       DATA (CTMP(I, 20),I=1,MAXINT)
94      +/ '.0305 ','L3    ','K     ','20.   ','100.  ','500.  '
95      +,'I     ',6*'      ' /
96       DATA (CTMP(I, 21),I=1,MAXINT)
97      +/ 'M3    ','L3    ','K     ','100.  ','500.  ','I     '
98      +,7*'      ' /
99       DATA (CTMP(I, 22),I=1,MAXINT)
100      +/ '.1    ','L3    ','K     ','100.  ','500.  ','I     '
101      +,7*'      ' /
102       DATA (CTMP(I, 23),I=1,MAXINT)
103      +/ '.12   ','L3    ','L1    ','K     ','100.  ','500.  '
104      +,'I     ',6*'      ' /
105       DATA (CTMP(I, 24),I=1,MAXINT)
106      +/ '.106  ','L3    ','K     ','100.  ','500.  ','I     '
107      +,7*'      ' /
108       DATA (CTMP(I, 25),I=1,MAXINT)
109      +/ '.1    ','L3    ','K     ','100.  ','500.  ','I     '
110      +,7*'      ' /
111       DATA (CTMP(I, 26),I=1,MAXINT)
112      +/ '.024  ','.12   ','L3    ','K     ','100.  ','500.  '
113      +,'I     ',6*'      ' /
114       DATA (CTMP(I, 27),I=1,MAXINT)
115      +/ '.0308 ','.093  ','L3    ','K     ','100.  ','500.  '
116      +,'I     ',6*'      ' /
117       DATA (CTMP(I, 28),I=1,MAXINT)
118      +/ '.1    ','L3    ','L2    ','L1    ','K     ','100.  '
119      +,'500.  ','I     ',5*'      ' /
120       DATA (CTMP(I, 29),I=1,MAXINT)
121      +/ '.1    ','L3    ','L1    ','K     ','100.  ','500.  '
122      +,'I     ',6*'      ' /
123       DATA (CTMP(I, 30),I=1,MAXINT)
124      +/ '.0307 ','.109  ','L3    ','L2    ','L1    ','K     '
125      +,'100.  ','500.  ','I     ',4*'      ' /
126       DATA (CTMP(I, 31),I=1,MAXINT)
127      +/ '.0278 ','.086  ','L3    ','L2    ','L1    ','K     '
128      +,'100.  ','500.  ','I     ',4*'      ' /
129       DATA (CTMP(I, 32),I=1,MAXINT)
130      +/ '.0305 ','.1    ','L3    ','L2    ','L1    ','K     '
131      +,'100.  ','500.  ','I     ',4*'      ' /
132       DATA (CTMP(I, 33),I=1,MAXINT)
133      +/ '.0371 ','.1    ','L3    ','L2    ','L1    ','K     '
134      +,'100.  ','500.  ','I     ',4*'      ' /
135       DATA (CTMP(I, 34),I=1,MAXINT)
136      +/ '.0371 ','.114  ','L3    ','L2    ','L1    ','K     '
137      +,'100.  ','500.  ','I     ',4*'      ' /
138       DATA (CTMP(I, 35),I=1,MAXINT)
139      +/ '.035  ','.1    ','.5    ','L3    ','L2    ','L1    '
140      +,'K     ','100.  ','500.  ','I     ',3*'      ' /
141       DATA (CTMP(I, 36),I=1,MAXINT)
142      +/ '.0305 ','.097  ','.3    ','L3    ','L2    ','L1    '
143      +,'K     ','100.  ','500.  ','I     ',3*'      ' /
144       DATA (CTMP(I, 37),I=1,MAXINT)
145      +/ '.0305 ','M5    ','.3    ','L3    ','L2    ','L1    '
146      +,'K     ','100.  ','500.  ','I     ',3*'      ' /
147       DATA (CTMP(I, 38),I=1,MAXINT)
148      +/ '.0307 ','M5    ','.3    ','L3    ','L2    ','L1    '
149      +,'K     ','100.  ','500.  ','I     ',3*'      ' /
150       DATA (CTMP(I, 39),I=1,MAXINT)
151      +/ '.0232 ','.1677 ','.8    ','L3    ','L2    ','L1    '
152      +,'K     ','100.  ','500.  ','I     ',3*'      ' /
153       DATA (CTMP(I, 40),I=1,MAXINT)
154      +/ 'M5    ','.5    ','L3    ','L2    ','L1    ','K     '
155      +,'100.  ','500.  ','I     ',4*'      ' /
156       DATA (CTMP(I, 41),I=1,MAXINT)
157      +/ '.0297 ','.2083 ','.9    ','L3    ','L2    ','L1    '
158      +,'K     ','100.  ','500.  ','I     ',3*'      ' /
159       DATA (CTMP(I, 42),I=1,MAXINT)
160      +/ '.03   ','M5    ','.8    ','L3    ','L2    ','L1    '
161      +,'K     ','100.  ','500.  ','I     ',3*'      ' /
162       DATA (CTMP(I, 43),I=1,MAXINT)
163      +/ '.068  ','M5    ','M3    ','L3    ','L2    ','L1    '
164      +,'K     ','100.  ','500.  ','I     ',3*'      ' /
165       DATA (CTMP(I, 44),I=1,MAXINT)
166      +/ '.0716 ','M5    ','M3    ','L3    ','L2    ','L1    '
167      +,'K     ','100.  ','500.  ','I     ',3*'      ' /
168       DATA (CTMP(I, 45),I=1,MAXINT)
169      +/ '.072  ','M5    ','M3    ','L3    ','L2    ','L1    '
170      +,'K     ','100.  ','500.  ','I     ',3*'      ' /
171       DATA (CTMP(I, 46),I=1,MAXINT)
172      +/ '.07   ','M5    ','M3    ','L3    ','L2    ','L1    '
173      +,'K     ','100.  ','500.  ','I     ',3*'      ' /
174       DATA (CTMP(I, 47),I=1,MAXINT)
175      +/ '.086  ','M5    ','M3    ','L3    ','L2    ','L1    '
176      +,'K     ','100.  ','500.  ','I     ',3*'      ' /
177       DATA (CTMP(I, 48),I=1,MAXINT)
178      +/ '.03   ','.098  ','M5    ','M3    ','L3    ','L2    '
179      +,'L1    ','K     ','100.  ','500.  ','I     ',2*'      ' /
180       DATA (CTMP(I, 49),I=1,MAXINT)
181      +/ '.031  ','.103  ','M5    ','M2    ','L3    ','L2    '
182      +,'L1    ','K     ','100.  ','500.  ','I     ',2*'      ' /
183       DATA (CTMP(I, 50),I=1,MAXINT)
184      +/ '.0318 ','.093  ','M5    ','M3    ','L3    ','L2    '
185      +,'L1    ','K     ','100.  ','500.  ','I     ',2*'      ' /
186       DATA (CTMP(I, 51),I=1,MAXINT)
187      +/ '.0304 ','.096  ','M5    ','M3    ','L3    ','L2    '
188      +,'L1    ','K     ','100.  ','500.  ','I     ',2*'      ' /
189       DATA (CTMP(I, 52),I=1,MAXINT)
190      +/ '.0528 ','.106  ','M5    ','M1    ','L3    ','L2    '
191      +,'L1    ','K     ','100.  ','500.  ','I     ',2*'      ' /
192       DATA (CTMP(I, 53),I=1,MAXINT)
193      +/ '.034  ','.103  ','M5    ','M2    ','M1    ','L3    '
194      +,'L2    ','L1    ','K     ','100.  ','500.  ','I     ','      ' /
195       DATA (CTMP(I, 54),I=1,MAXINT)
196      +/ '.0515 ','.1073 ','M5    ','M3    ','M1    ','L3    '
197      +,'L2    ','L1    ','K     ','100.  ','500.  ','I     ','      ' /
198       DATA (CTMP(I, 55),I=1,MAXINT)
199      +/ '.0704 ','.111  ','M5    ','M2    ','L3    ','L2    '
200      +,'L1    ','K     ','100.  ','500.  ','I     ',2*'      ' /
201       DATA (CTMP(I, 56),I=1,MAXINT)
202      +/ '.0305 ','N5    ','M5    ','M2    ','L3    ','L2    '
203      +,'L1    ','K     ','100.  ','500.  ','I     ',2*'      ' /
204       DATA (CTMP(I, 57),I=1,MAXINT)
205      +/ '.04   ','.13   ','M5    ','M3    ','M2    ','L3    '
206      +,'L2    ','L1    ','K     ','100.  ','500.  ','I     ','      ' /
207       DATA (CTMP(I, 58),I=1,MAXINT)
208      +/ '.053  ','N5    ','M5    ','M3    ','M2    ','M1    '
209      +,'L3    ','L2    ','L1    ','K     ','100.  ','500.  ','I     ' /
210       DATA (CTMP(I, 59),I=1,MAXINT)
211      +/ '.04   ','N5    ','M5    ','M3    ','M2    ','M1    '
212      +,'L3    ','L2    ','L1    ','K     ','100.  ','500.  ','I     ' /
213       DATA (CTMP(I, 60),I=1,MAXINT)
214      +/ '.0305 ','N5    ','M5    ','M3    ','M1    ','L3    '
215      +,'L2    ','L1    ','K     ','100.  ','500.  ','I     ','      ' /
216       DATA (CTMP(I, 61),I=1,MAXINT)
217      +/ '.0305 ','N5    ','M5    ','M3    ','M2    ','M1    '
218      +,'L3    ','L2    ','L1    ','K     ','100.  ','500.  ','I     ' /
219       DATA (CTMP(I, 62),I=1,MAXINT)
220      +/ '.0305 ','N5    ','M5    ','M3    ','L3    ','L2    '
221      +,'L1    ','K     ','100.  ','500.  ','I     ',2*'      ' /
222       DATA (CTMP(I, 63),I=1,MAXINT)
223      +/ '.0305 ','N5    ','M5    ','M3    ','M2    ','M1    '
224      +,'L3    ','L2    ','L1    ','K     ','100.  ','500.  ','I     ' /
225       DATA (CTMP(I, 64),I=1,MAXINT)
226      +/ '.0305 ','N5    ','M5    ','M3    ','M2    ','M1    '
227      +,'L3    ','L2    ','L1    ','K     ','500.  ','I     ','      ' /
228       DATA (CTMP(I, 65),I=1,MAXINT)
229      +/ '.0305 ','N5    ','M5    ','M3    ','M2    ','M1    '
230      +,'L3    ','L2    ','L1    ','K     ','500.  ','I     ','      ' /
231       DATA (CTMP(I, 66),I=1,MAXINT)
232      +/ '.0305 ','N5    ','M5    ','M3    ','M2    ','M1    '
233      +,'L3    ','L2    ','L1    ','K     ','500.  ','I     ','      ' /
234       DATA (CTMP(I, 67),I=1,MAXINT)
235      +/ '.0426 ','N5    ','M5    ','M3    ','M2    ','M1    '
236      +,'L3    ','L2    ','L1    ','K     ','500.  ','I     ','      ' /
237       DATA (CTMP(I, 68),I=1,MAXINT)
238      +/ '.0322 ','N5    ','M5    ','M3    ','M2    ','M1    '
239      +,'L3    ','L2    ','L1    ','K     ','500.  ','I     ','      ' /
240       DATA (CTMP(I, 69),I=1,MAXINT)
241      +/ '.031  ','N5    ','M5    ','M3    ','L3    ','L2    '
242      +,'L1    ','K     ','500.  ','I     ',3*'      ' /
243       DATA (CTMP(I, 70),I=1,MAXINT)
244      +/ '.075  ','N5    ','M5    ','M3    ','L3    ','L2    '
245      +,'L1    ','K     ','500.  ','I     ',3*'      ' /
246       DATA (CTMP(I, 71),I=1,MAXINT)
247      +/ '.09   ','N5    ','M5    ','M3    ','L3    ','L2    '
248      +,'L1    ','K     ','500.  ','I     ',3*'      ' /
249       DATA (CTMP(I, 72),I=1,MAXINT)
250      +/ '.0512 ','N5    ','M5    ','M3    ','L3    ','L2    '
251      +,'L1    ','K     ','500.  ','I     ',3*'      ' /
252       DATA (CTMP(I, 73),I=1,MAXINT)
253      +/ '.0233 ','.0428 ','N5    ','M5    ','M3    ','L3    '
254      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
255       DATA (CTMP(I, 74),I=1,MAXINT)
256      +/ '.031  ','.0493 ','N5    ','M5    ','M3    ','L3    '
257      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
258       DATA (CTMP(I, 75),I=1,MAXINT)
259      +/ '.0384 ','N5    ','M5    ','M3    ','L3    ','L2    '
260      +,'L1    ','K     ','500.  ','I     ',3*'      ' /
261       DATA (CTMP(I, 76),I=1,MAXINT)
262      +/ '.0464 ','N5    ','M5    ','M3    ','L3    ','L2    '
263      +,'L1    ','K     ','500.  ','I     ',3*'      ' /
264       DATA (CTMP(I, 77),I=1,MAXINT)
265      +/ '.048  ','N5    ','M5    ','M3    ','L3    ','L2    '
266      +,'L1    ','K     ','500.  ','I     ',3*'      ' /
267       DATA (CTMP(I, 78),I=1,MAXINT)
268      +/ '.0266 ','.0735 ','N5    ','M5    ','M3    ','L3    '
269      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
270       DATA (CTMP(I, 79),I=1,MAXINT)
271      +/ '.07   ','N5    ','M5    ','M3    ','L3    ','L2    '
272      +,'L1    ','K     ','500.  ','I     ',3*'      ' /
273       DATA (CTMP(I, 80),I=1,MAXINT)
274      +/ '.09   ','N5    ','M5    ','M3    ','L3    ','L2    '
275      +,'L1    ','K     ','500.  ','I     ',3*'      ' /
276       DATA (CTMP(I, 81),I=1,MAXINT)
277      +/ '.0313 ','.0985 ','.386  ','M5    ','M3    ','L3    '
278      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
279       DATA (CTMP(I, 82),I=1,MAXINT)
280      +/ '.0311 ','.095  ','.413  ','M5    ','M3    ','L3    '
281      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
282       DATA (CTMP(I, 83),I=1,MAXINT)
283      +/ 'O5    ','O4    ','.0379 ','.1136 ','.362  ','M5    '
284      +,'M3    ','L3    ','L2    ','L1    ','K     ','500.  ','I     ' /
285       DATA (CTMP(I, 84),I=1,MAXINT)
286      +/ '.0305 ','.07   ','.2415 ','M5    ','M3    ','L3    '
287      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
288       DATA (CTMP(I, 85),I=1,MAXINT)
289      +/ '.0305 ','.072  ','.244  ','M5    ','M3    ','L3    '
290      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
291       DATA (CTMP(I, 86),I=1,MAXINT)
292      +/ '.0466 ','.093  ','.305  ','M5    ','M3    ','L3    '
293      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
294       DATA (CTMP(I, 87),I=1,MAXINT)
295      +/ '.0494 ','.196  ','.525  ','M5    ','M3    ','L3    '
296      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
297       DATA (CTMP(I, 88),I=1,MAXINT)
298      +/ '.0579 ','.15   ','.396  ','M5    ','M3    ','L3    '
299      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
300       DATA (CTMP(I, 89),I=1,MAXINT)
301      +/ '.0712 ','.216  ','.375  ','M5    ','M3    ','L3    '
302      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
303       DATA (CTMP(I, 90),I=1,MAXINT)
304      +/ 'O5    ','N5    ','M5    ','M3    ','L3    ','L2    '
305      +,'L1    ','K     ','500.  ','I     ',3*'      ' /
306       DATA (CTMP(I, 91),I=1,MAXINT)
307      +/ 'O5    ','.371  ','M5    ','M3    ','L3    ','L2    '
308      +,'L1    ','K     ','500.  ','I     ',3*'      ' /
309       DATA (CTMP(I, 92),I=1,MAXINT)
310      +/ '.0171 ','O5    ','.617  ','M5    ','M3    ','L3    '
311      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
312       DATA (CTMP(I, 93),I=1,MAXINT)
313      +/ 'O5    ','.431  ','M5    ','M3    ','L3    ','L2    '
314      +,'L1    ','K     ','500.  ','I     ',3*'      ' /
315       DATA (CTMP(I, 94),I=1,MAXINT)
316      +/ '.0306 ','O5    ','.405  ','M5    ','M3    ','L3    '
317      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
318       DATA (CTMP(I, 95),I=1,MAXINT)
319      +/ 'O5    ','.498  ','N5    ','M5    ','M3    ','L3    '
320      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
321       DATA (CTMP(I, 96),I=1,MAXINT)
322      +/ 'O5    ','.498  ','N5    ','M5    ','M3    ','L3    '
323      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
324       DATA (CTMP(I, 97),I=1,MAXINT)
325      +/ 'O5    ','.531  ','N5    ','M5    ','M3    ','L3    '
326      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
327       DATA (CTMP(I, 98),I=1,MAXINT)
328      +/ 'O5    ','.548  ','N5    ','M5    ','M3    ','L3    '
329      +,'L2    ','L1    ','K     ','500.  ','I     ',2*'      ' /
330       DATA (CTMP(I, 99),I=1,MAXINT)
331      +/ 'N5    ','M5    ','M4    ','M3    ','M2    ','M1    '
332      +,'L3    ','L2    ','L1    ','K     ','500.  ','I     ','      ' /
333       DATA (CTMP(I,100),I=1,MAXINT)
334      +/ 'N5    ','M5    ','M4    ','M3    ','M2    ','M1    '
335      +,'L3    ','L2    ','L1    ','K     ','500.  ','I     ','      ' /
336 C
337 C For HYDROGEN the parametrization is changed below K shell
338 C  Orig. SANDIA     +   0.1000E-07,  0.0000E+00,  0.0000E+00,  0.0000E+00
339 C  New param        +   0.0000E+00,  0.0000E+00,  0.0000E+00,  0.0000E+00
340 C
341       DATA ((TMP(I,J,  1),I=1,MAXPOW),J=1,MAXINT) /
342      +   0.0000E+00,  0.0000E+00,  0.0000E+00,  0.0000E+00
343      + , -0.6383E+02, -0.6446E+01,  0.1317E+02, -0.5045E-01
344      + ,  0.3051E+01, -0.7818E+01,  0.1144E+02,  0.6959E-01
345      + ,  0.7636E-01, -0.9406E+00,  0.6144E+01,  0.1425E+01
346      + ,  0.1180E-02, -0.8230E-01,  0.2880E+01,  0.5534E+01
347      + ,  0.1620E-04, -0.5610E-02,  0.1214E+01,  0.1761E+02
348      + ,  0.1034E-05, -0.4114E-03,  0.6287E+00,  0.3927E+02
349      + ,  0.4599E-06,  0.5006E-03, -0.1425E-01,  0.1960E+03,  20*0.0 /
350       DATA ((TMP(I,J,  2),I=1,MAXPOW),J=1,MAXINT) /
351      +   0.1000E-07,  0.0000E+00,  0.0000E+00,  0.0000E+00
352      + , -0.1196E+04,  0.4759E+03,  0.1979E+02, -0.3540E+00
353      + ,  0.1857E+02, -0.4079E+02,  0.8609E+02, -0.2618E+01
354      + ,  0.1606E+01, -0.1726E+02,  0.8101E+02, -0.4928E+01
355      + ,  0.3663E-02, -0.3782E+00,  0.2288E+02,  0.6006E+02
356      + ,  0.5934E-03, -0.1109E+00,  0.1388E+02,  0.1680E+03
357      + ,  0.1120E-04, -0.3894E-02,  0.6136E+01,  0.4354E+03
358      + ,  0.4921E-05,  0.1493E-02,  0.4725E+01,  0.7128E+03,  20*0.0 /
359       DATA ((TMP(I,J,  3),I=1,MAXPOW),J=1,MAXINT) /
360      +  -0.2418E+03,  0.3530E+02, -0.2699E+00,  0.5580E-03
361      + , -0.1474E+03,  0.1578E+03,  0.2036E+03, -0.8521E+01
362      + ,  0.4375E+01, -0.5094E+02,  0.3053E+03, -0.3064E+02
363      + ,  0.4040E-01, -0.2930E+01,  0.1293E+03,  0.1828E+03
364      + ,  0.3139E-02, -0.6107E+00,  0.7342E+02,  0.6852E+03
365      + ,  0.4823E-04, -0.1161E-01,  0.2898E+02,  0.2129E+04
366      + ,  0.2128E-04,  0.1951E-01,  0.1009E+02,  0.7144E+04,  24*0.0 /
367       DATA ((TMP(I,J,  4),I=1,MAXPOW),J=1,MAXINT) /
368      +  -0.2785E+03,  0.1252E+03, -0.1912E+01,  0.8371E-02
369      + , -0.4798E+03,  0.6302E+03,  0.4187E+03, -0.2783E+02
370      + ,  0.1037E+02, -0.1280E+03,  0.8807E+03, -0.1646E+03
371      + ,  0.1018E+00, -0.8264E+01,  0.4158E+03,  0.4366E+03
372      + ,  0.5893E-02, -0.1595E+01,  0.2644E+03,  0.1566E+04
373      + ,  0.2129E-03, -0.7972E-01,  0.1155E+03,  0.6065E+04
374      + ,  0.7797E-04,  0.7695E-01,  0.3396E+02,  0.2439E+05,  24*0.0 /
375       DATA ((TMP(I,J,  5),I=1,MAXPOW),J=1,MAXINT) /
376      +  -0.3387E+03,  0.1923E+03, -0.2742E+01,  0.1204E-01
377      + , -0.9943E+03,  0.1816E+04,  0.3939E+03, -0.1755E+02
378      + ,  0.3689E+01, -0.8834E+02,  0.1525E+04, -0.2145E+03
379      + ,  0.6447E+00, -0.3953E+02,  0.1281E+04,  0.1768E+03
380      + ,  0.1852E-02, -0.1358E+01,  0.5252E+03,  0.5160E+04
381      + ,  0.4282E-03, -0.1297E+00,  0.2892E+03,  0.1601E+05
382      + ,  0.1965E-03,  0.1947E+00,  0.1001E+03,  0.5902E+05,  24*0.0 /
383       DATA ((TMP(I,J,  6),I=1,MAXPOW),J=1,MAXINT) /
384      +   0.5704E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
385      + , -0.3935E+03,  0.3219E+03, -0.6549E+01,  0.2086E+00
386      + , -0.9022E+03,  0.1760E+04,  0.1549E+04, -0.2280E+03
387      + , -0.7363E+01, -0.1537E+01,  0.2672E+04, -0.4482E+03
388      + ,  0.1640E+01, -0.9428E+02,  0.2872E+04, -0.5583E+03
389      + , -0.3742E-02, -0.1232E+01,  0.9489E+03,  0.1406E+05
390      + ,  0.6760E-03, -0.1164E+00,  0.6566E+03,  0.2739E+05
391      + ,  0.4503E-03,  0.4692E+00,  0.2183E+03,  0.1308E+06,  20*0.0 /
392       DATA ((TMP(I,J,  7),I=1,MAXPOW),J=1,MAXINT) /
393      +   0.1010E+05,  0.0000E+00,  0.0000E+00,  0.0000E+00
394      + , -0.3622E+03,  0.3873E+03,  0.1244E+02, -0.4452E+00
395      + , -0.2338E+04,  0.5732E+04, -0.2082E+03,  0.1482E+03
396      + , -0.4940E+01, -0.8442E+02,  0.4620E+04, -0.1186E+04
397      + ,  0.2019E+01, -0.1249E+03,  0.4609E+04, -0.9421E+03
398      + ,  0.1709E-01, -0.8196E+01,  0.2345E+04,  0.1369E+05
399      + ,  0.1872E-02, -0.6732E+00,  0.1282E+04,  0.5700E+05
400      + ,  0.8122E-03,  0.8364E+00,  0.4410E+03,  0.2358E+06,  20*0.0 /
401       DATA ((TMP(I,J,  8),I=1,MAXPOW),J=1,MAXINT) /
402      +   0.1144E+05,  0.0000E+00,  0.0000E+00,  0.0000E+00
403      + , -0.2863E+03,  0.4086E+03,  0.4436E+02, -0.1782E+01
404      + , -0.7181E+02,  0.4748E+03,  0.5542E+04, -0.1363E+04
405      + ,  0.2745E+01, -0.1747E+03,  0.7159E+04, -0.2213E+04
406      + ,  0.3774E-01, -0.1559E+02,  0.4045E+04,  0.1810E+05
407      + ,  0.3169E-02, -0.1146E+01,  0.2194E+04,  0.9131E+05
408      + ,  0.1367E-02,  0.1473E+01,  0.7214E+03,  0.4048E+06,  24*0.0 /
409       DATA ((TMP(I,J,  9),I=1,MAXPOW),J=1,MAXINT) /
410      +   0.1129E+05,  0.0000E+00,  0.0000E+00,  0.0000E+00
411      + , -0.2676E+03,  0.4767E+03,  0.8767E+02, -0.4504E+01
412      + , -0.8907E+02,  0.6850E+03,  0.7217E+04, -0.2022E+04
413      + ,  0.2813E+01, -0.2020E+03,  0.9903E+04, -0.4454E+04
414      + ,  0.6097E-01, -0.2435E+02,  0.6099E+04,  0.2258E+05
415      + ,  0.4857E-02, -0.1806E+01,  0.3335E+04,  0.1285E+06
416      + ,  0.2041E-02,  0.2257E+01,  0.1084E+04,  0.5997E+06,  24*0.0 /
417       DATA ((TMP(I,J, 10),I=1,MAXPOW),J=1,MAXINT) /
418      +   0.1000E-05,  0.0000E+00,  0.0000E+00,  0.0000E+00
419      + ,  0.8235E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
420      + , -0.6059E+04,  0.3242E+04, -0.1796E+03,  0.3141E+01
421      + , -0.1940E+03,  0.4340E+03,  0.2357E+03, -0.1666E+02
422      + , -0.8680E+02,  0.7434E+03,  0.1077E+05, -0.3822E+04
423      + ,  0.3167E+01, -0.2491E+03,  0.1439E+05, -0.8187E+04
424      + ,  0.1288E+00, -0.4471E+02,  0.9862E+04,  0.2497E+05
425      + ,  0.5202E-02,  0.7337E-01,  0.4614E+04,  0.2256E+06
426      + ,  0.3087E-02,  0.3402E+01,  0.1569E+04,  0.1147E+07,  16*0.0 /
427       DATA ((TMP(I,J, 11),I=1,MAXPOW),J=1,MAXINT) /
428      +   0.5429E+02,  0.0000E+00,  0.0000E+00,  0.0000E+00
429      + ,  0.8874E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
430      + , -0.5565E+03,  0.1000E+04,  0.1062E+03, -0.6846E+01
431      + ,  0.1732E+01, -0.1853E+03,  0.1712E+05, -0.8415E+04
432      + ,  0.2726E+00, -0.8218E+02,  0.1475E+05,  0.9467E+04
433      + ,  0.2083E-02,  0.4548E+01,  0.5521E+04,  0.3358E+06
434      + ,  0.4199E-02,  0.4808E+01,  0.1793E+04,  0.1898E+07,  24*0.0 /
435       DATA ((TMP(I,J, 12),I=1,MAXPOW),J=1,MAXINT) /
436      +  -0.5574E+03,  0.8687E+02, -0.1493E+01,  0.8660E-02
437      + , -0.2726E+05,  0.9700E+04, -0.6921E+03,  0.1446E+02
438      + , -0.6179E+03,  0.1379E+04,  0.2321E+03, -0.2291E+02
439      + , -0.6648E+00, -0.1084E+03,  0.2219E+05, -0.1241E+05
440      + ,  0.4080E+00, -0.1205E+03,  0.2139E+05, -0.1143E+03
441      + ,  0.8338E-02,  0.2236E+01,  0.8819E+04,  0.4284E+06
442      + ,  0.6233E-02,  0.7344E+01,  0.2662E+04,  0.2460E+07,  24*0.0 /
443       DATA ((TMP(I,J, 13),I=1,MAXPOW),J=1,MAXINT) /
444      +  -0.1654E+05,  0.1585E+03,  0.3907E+01, -0.3383E-01
445      + ,  0.1122E+04, -0.4015E+02,  0.6623E+00, -0.2813E-02
446      + ,  0.2390E+05, -0.6953E+03, -0.7978E+02,  0.1974E+01
447      + , -0.5284E+03,  0.1399E+04,  0.4360E+03, -0.4747E+02
448      + , -0.3674E+01, -0.1622E+02,  0.2732E+05, -0.1752E+05
449      + ,  0.4158E+00, -0.1351E+03,  0.2716E+05,  0.3723E+03
450      + ,  0.1125E-01,  0.2747E+01,  0.1174E+05,  0.5695E+06
451      + ,  0.8505E-02,  0.1041E+02,  0.4429E+04,  0.2413E+07,  20*0.0 /
452       DATA ((TMP(I,J, 14),I=1,MAXPOW),J=1,MAXINT) /
453      +   0.4575E+03,  0.0000E+00,  0.0000E+00,  0.0000E+00
454      + , -0.1168E+03,  0.1269E+03, -0.5506E+01,  0.6616E-01
455      + , -0.6451E+03,  0.1836E+04,  0.5350E+03, -0.6436E+02
456      + , -0.4056E+01, -0.1089E+02,  0.3544E+05, -0.2595E+05
457      + ,  0.4945E+00, -0.1682E+03,  0.3627E+05, -0.1605E+05
458      + ,  0.2716E-01, -0.5428E+01,  0.1775E+05,  0.6765E+06
459      + ,  0.1127E-01,  0.1376E+02,  0.4871E+04,  0.4477E+07,  24*0.0 /
460       DATA ((TMP(I,J, 15),I=1,MAXPOW),J=1,MAXINT) /
461      +   0.4728E+03,  0.0000E+00,  0.0000E+00,  0.0000E+00
462      + , -0.2227E+03,  0.2204E+03, -0.1122E+02,  0.1569E+00
463      + , -0.5700E+03,  0.1909E+04,  0.7651E+03, -0.1123E+03
464      + , -0.6098E+01,  0.7712E+02,  0.4196E+05, -0.3437E+05
465      + ,  0.5710E+00, -0.2001E+03,  0.4505E+05, -0.3853E+05
466      + ,  0.3613E-01, -0.8426E+01,  0.2277E+05,  0.8080E+06
467      + ,  0.1421E-01,  0.1770E+02,  0.6221E+04,  0.5554E+07,  24*0.0 /
468       DATA ((TMP(I,J, 16),I=1,MAXPOW),J=1,MAXINT) /
469      +   0.5368E+03,  0.0000E+00,  0.0000E+00,  0.0000E+00
470      + , -0.2097E+03,  0.3174E+03, -0.1815E+02,  0.2796E+00
471      + , -0.4885E+03,  0.1976E+04,  0.1206E+04, -0.1934E+03
472      + , -0.1021E+02,  0.2346E+03,  0.5164E+05, -0.4707E+05
473      + ,  0.7467E+00, -0.2534E+03,  0.5800E+05, -0.6691E+05
474      + ,  0.6338E-01, -0.3443E+02,  0.4058E+05,  0.2233E+06
475      + ,  0.1867E-01,  0.2398E+02,  0.6706E+04,  0.8358E+07,  24*0.0 /
476       DATA ((TMP(I,J, 17),I=1,MAXPOW),J=1,MAXINT) /
477      +   0.7808E+03,  0.0000E+00,  0.0000E+00,  0.0000E+00
478      + , -0.6163E+02,  0.3934E+03, -0.2563E+02,  0.4396E+00
479      + , -0.2147E+03,  0.1269E+04,  0.2099E+04, -0.3705E+03
480      + , -0.1192E+02,  0.3253E+03,  0.5927E+05, -0.5980E+05
481      + ,  0.8789E+00, -0.2920E+03,  0.6860E+05, -0.1018E+06
482      + ,  0.6556E-01, -0.2175E+02,  0.3895E+05,  0.9737E+06
483      + ,  0.2236E-01,  0.2884E+02,  0.9717E+04,  0.9001E+07,  24*0.0 /
484       DATA ((TMP(I,J, 18),I=1,MAXPOW),J=1,MAXINT) /
485      +   0.1275E+05,  0.0000E+00,  0.0000E+00,  0.0000E+00
486      + , -0.1674E+04,  0.1047E+04, -0.9010E+02,  0.2166E+01
487      + , -0.2244E+03,  0.1388E+04,  0.2571E+04, -0.5113E+03
488      + , -0.1446E+02,  0.4359E+03,  0.6578E+05, -0.7284E+05
489      + ,  0.8786E+00, -0.3001E+03,  0.7682E+05, -0.1219E+06
490      + ,  0.6779E-01, -0.2695E+02,  0.5069E+05,  0.5281E+06
491      + ,  0.2598E-01,  0.3454E+02,  0.1059E+05,  0.1068E+08,  24*0.0 /
492       DATA ((TMP(I,J, 19),I=1,MAXPOW),J=1,MAXINT) /
493      +   0.3111E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
494      + ,  0.6162E+02,  0.5770E+03, -0.4497E+02,  0.9218E+00
495      + , -0.3601E+03,  0.2287E+04,  0.2513E+04, -0.4591E+03
496      + , -0.2113E+02,  0.7161E+03,  0.8183E+05, -0.9859E+05
497      + ,  0.1048E+01, -0.3644E+03,  0.9844E+05, -0.1760E+06
498      + ,  0.1040E+00, -0.3553E+02,  0.6096E+05,  0.1226E+07
499      + ,  0.3453E-01,  0.4621E+02,  0.1669E+05,  0.1238E+08,  24*0.0 /
500       DATA ((TMP(I,J, 20),I=1,MAXPOW),J=1,MAXINT) /
501      +   0.3029E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
502      + ,  0.3816E+03,  0.5887E+03, -0.4888E+02,  0.1019E+01
503      + , -0.1054E+03,  0.1120E+04,  0.5115E+04, -0.1325E+04
504      + , -0.3015E+02,  0.1141E+04,  0.9459E+05, -0.1195E+06
505      + ,  0.1408E+01, -0.4553E+03,  0.1206E+06, -0.2531E+06
506      + ,  0.1204E+00, -0.5356E+02,  0.8779E+05,  0.6252E+06
507      + ,  0.4178E-01,  0.5761E+02,  0.1377E+05,  0.2019E+08,  24*0.0 /
508       DATA ((TMP(I,J, 21),I=1,MAXPOW),J=1,MAXINT) /
509      +   0.1121E+04,  0.7413E+01, -0.3020E+00,  0.1510E-02
510      + ,  0.2970E+03,  0.6359E+03, -0.4935E+02,  0.1044E+01
511      + , -0.1507E+03,  0.1501E+04,  0.5139E+04, -0.1442E+04
512      + ,  0.4615E+00, -0.2614E+03,  0.1217E+06, -0.2057E+06
513      + ,  0.1509E+00, -0.6913E+02,  0.9253E+05,  0.1096E+07
514      + ,  0.4833E-01,  0.6880E+02,  0.2762E+05,  0.1290E+08,  28*0.0 /
515       DATA ((TMP(I,J, 22),I=1,MAXPOW),J=1,MAXINT) /
516      +   0.3596E+04, -0.6877E+02,  0.9609E+00, -0.4754E-02
517      + , -0.3132E+03,  0.1024E+04, -0.1008E+03,  0.3096E+01
518      + , -0.4626E+02,  0.8678E+03,  0.7308E+04, -0.2383E+04
519      + ,  0.2074E+00, -0.2292E+03,  0.1371E+06, -0.2489E+06
520      + ,  0.1972E+00, -0.9550E+02,  0.1106E+06,  0.1067E+07
521      + ,  0.5606E-01,  0.8140E+02,  0.3106E+05,  0.1567E+08,  28*0.0 /
522       DATA ((TMP(I,J, 23),I=1,MAXPOW),J=1,MAXINT) /
523      +   0.4117E+04, -0.1054E+03,  0.1531E+01, -0.7725E-02
524      + , -0.4052E+03,  0.1169E+04, -0.1212E+03,  0.4148E+01
525      + ,  0.9968E+04, -0.1796E+04, -0.1204E+04,  0.1481E+04
526      + , -0.1429E+03,  0.1632E+04,  0.7203E+04, -0.2283E+04
527      + ,  0.3940E+00, -0.2733E+03,  0.1564E+06, -0.3184E+06
528      + ,  0.2007E+00, -0.8911E+02,  0.1254E+06,  0.1137E+07
529      + ,  0.6449E-01,  0.9483E+02,  0.3584E+05,  0.1821E+08,  24*0.0 /
530       DATA ((TMP(I,J, 24),I=1,MAXPOW),J=1,MAXINT) /
531      +   0.5069E+04, -0.7968E+02,  0.8395E+00, -0.4494E-02
532      + , -0.4404E+03,  0.1301E+04, -0.1166E+03,  0.3481E+01
533      + ,  0.3932E+02,  0.1537E+03,  0.1228E+05, -0.4984E+04
534      + ,  0.2740E+00, -0.2665E+03,  0.1818E+06, -0.3997E+06
535      + ,  0.2727E+00, -0.1288E+03,  0.1543E+06,  0.9755E+06
536      + ,  0.7660E-01,  0.1142E+03,  0.4290E+05,  0.2161E+08,  28*0.0 /
537       DATA ((TMP(I,J, 25),I=1,MAXPOW),J=1,MAXINT) /
538      +   0.5692E+04, -0.1278E+03,  0.1150E+01, -0.3774E-02
539      + , -0.4901E+03,  0.1450E+04, -0.1365E+03,  0.4134E+01
540      + , -0.3065E+02,  0.7209E+03,  0.1275E+05, -0.5462E+04
541      + ,  0.2930E+00, -0.2703E+03,  0.2034E+06, -0.4886E+06
542      + ,  0.2940E+00, -0.1342E+03,  0.1761E+06,  0.8747E+06
543      + ,  0.8736E-01,  0.1322E+03,  0.4952E+05,  0.2436E+08,  28*0.0 /
544       DATA ((TMP(I,J, 26),I=1,MAXPOW),J=1,MAXINT) /
545      +  -0.3328E+04,  0.2590E+03, -0.3812E+01,  0.1679E-01
546      + ,  0.4132E+04,  0.3963E+03, -0.2505E+02,  0.3437E+00
547      + , -0.2683E+02,  0.1250E+04, -0.7321E+02,  0.1137E+01
548      + , -0.8128E+01,  0.6550E+03,  0.1511E+05, -0.7079E+04
549      + ,  0.2151E+00, -0.2680E+03,  0.2341E+06, -0.6063E+06
550      + ,  0.3511E+00, -0.1633E+03,  0.2090E+06,  0.7131E+06
551      + ,  0.1026E+00,  0.1572E+03,  0.5901E+05,  0.2806E+08,  24*0.0 /
552       DATA ((TMP(I,J, 27),I=1,MAXPOW),J=1,MAXINT) /
553      +   0.2074E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
554      + ,  0.6932E+04, -0.6765E-01, -0.9182E+01,  0.1423E+00
555      + , -0.8228E+03,  0.1982E+04, -0.1782E+03,  0.4937E+01
556      + ,  0.2092E+02,  0.3311E+03,  0.1820E+05, -0.9303E+04
557      + ,  0.4947E-01, -0.2478E+03,  0.2586E+06, -0.7185E+06
558      + ,  0.3588E+00, -0.1492E+03,  0.2296E+06,  0.8856E+06
559      + ,  0.1154E+00,  0.1796E+03,  0.6740E+05,  0.3100E+08,  24*0.0 /
560       DATA ((TMP(I,J, 28),I=1,MAXPOW),J=1,MAXINT) /
561      +   0.8515E+04, -0.2121E+03,  0.2132E+01, -0.7667E-02
562      + , -0.1021E+04,  0.2388E+04, -0.2361E+03,  0.7261E+01
563      + ,  0.0000E+00,  0.6600E+04,  0.0000E+00,  0.0000E+00
564      + ,  0.0000E+00,  0.9810E+04,  0.0000E+00,  0.0000E+00
565      + , -0.1175E+03,  0.2171E+04,  0.1514E+05, -0.6042E+04
566      + ,  0.1685E-01, -0.2522E+03,  0.3005E+06, -0.9000E+06
567      + ,  0.4681E+00, -0.2085E+03,  0.2782E+06,  0.4399E+06
568      + ,  0.1364E+00,  0.2148E+03,  0.8039E+05,  0.3626E+08,  20*0.0 /
569       DATA ((TMP(I,J, 29),I=1,MAXPOW),J=1,MAXINT) /
570      +   0.9992E+04, -0.3780E+03,  0.5657E+01, -0.2788E-01
571      + , -0.8896E+03,  0.2391E+04, -0.1986E+03,  0.3649E+01
572      + , -0.5089E+03,  0.4139E+04,  0.7468E+04, -0.6981E+03
573      + , -0.9173E+02,  0.1862E+04,  0.1777E+05, -0.8202E+04
574      + ,  0.7440E-02, -0.2371E+03,  0.3202E+06, -0.1028E+07
575      + ,  0.4729E+00, -0.2055E+03,  0.3000E+06,  0.2182E+06
576      + ,  0.1472E+00,  0.2345E+03,  0.8644E+05,  0.3990E+08,  24*0.0 /
577       DATA ((TMP(I,J, 30),I=1,MAXPOW),J=1,MAXINT) /
578      +   0.1479E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
579      + ,  0.1325E+05, -0.6762E+03,  0.1188E+02, -0.6785E-01
580      + , -0.1330E+04,  0.3104E+04, -0.3092E+03,  0.8946E+01
581      + , -0.2807E+05,  0.1091E+06, -0.1511E+06,  0.7897E+05
582      + , -0.3057E+04,  0.1385E+05,  0.3780E+03, -0.2163E+03
583      + , -0.8018E+02,  0.1794E+04,  0.2095E+05, -0.1070E+05
584      + , -0.1158E+01, -0.3871E+02,  0.3499E+06, -0.1154E+07
585      + ,  0.4820E+00, -0.1832E+03,  0.3296E+06,  0.6802E+06
586      + ,  0.1669E+00,  0.2696E+03,  0.1000E+06,  0.4382E+08,  16*0.0 /
587       DATA ((TMP(I,J, 31),I=1,MAXPOW),J=1,MAXINT) /
588      +   0.1479E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
589      + ,  0.1211E+05, -0.7144E+03,  0.1535E+02, -0.1030E+00
590      + , -0.1199E+04,  0.3215E+04, -0.3594E+03,  0.1152E+02
591      + , -0.1457E+05,  0.6363E+05, -0.1011E+06,  0.6416E+05
592      + ,  0.3470E+04, -0.8674E+04,  0.2548E+05, -0.8411E+04
593      + , -0.9431E+02,  0.2078E+04,  0.2172E+05, -0.1101E+05
594      + , -0.1253E+01, -0.2897E+02,  0.3783E+06, -0.1357E+07
595      + ,  0.5591E+00, -0.2217E+03,  0.3625E+06,  0.3392E+06
596      + ,  0.1949E+00,  0.2371E+03,  0.1767E+06,  0.2400E+08,  16*0.0 /
597       DATA ((TMP(I,J, 32),I=1,MAXPOW),J=1,MAXINT) /
598      +   0.4941E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
599      + ,  0.1269E+05, -0.9339E+03,  0.2767E+02, -0.1952E+00
600      + , -0.1358E+04,  0.3695E+04, -0.4569E+03,  0.1599E+02
601      + , -0.7402E+04,  0.3636E+05, -0.6670E+05,  0.5346E+05
602      + , -0.3400E+03,  0.3232E+04,  0.1163E+05, -0.1089E+04
603      + , -0.5940E+02,  0.1646E+04,  0.2621E+05, -0.1507E+05
604      + , -0.2025E+01,  0.1283E+03,  0.4069E+06, -0.1510E+07
605      + ,  0.6478E+00, -0.2603E+03,  0.4045E+06, -0.5153E+05
606      + ,  0.2155E+00,  0.2609E+03,  0.2074E+06,  0.2220E+08,  16*0.0 /
607       DATA ((TMP(I,J, 33),I=1,MAXPOW),J=1,MAXINT) /
608      +   0.2547E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
609      + ,  0.2402E+05, -0.3110E+04,  0.1364E+03, -0.1877E+01
610      + , -0.1536E+04,  0.4198E+04, -0.5497E+03,  0.1920E+02
611      + , -0.4525E+04,  0.2468E+05, -0.5130E+05,  0.5102E+05
612      + , -0.3858E+03,  0.1013E+05,  0.3603E+04, -0.7657E+03
613      + , -0.5339E+02,  0.1649E+04,  0.2950E+05, -0.1800E+05
614      + , -0.1655E+01,  0.8588E+02,  0.4504E+06, -0.1795E+07
615      + ,  0.7880E+00, -0.3490E+03,  0.4641E+06, -0.1258E+07
616      + ,  0.2393E+00,  0.2906E+03,  0.2360E+06,  0.2145E+08,  16*0.0 /
617       DATA ((TMP(I,J, 34),I=1,MAXPOW),J=1,MAXINT) /
618      +   0.1577E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
619      + ,  0.1806E+05, -0.2230E+04,  0.9408E+02, -0.1263E+01
620      + , -0.1610E+04,  0.4614E+04, -0.6392E+03,  0.2252E+02
621      + , -0.6637E+03,  0.4194E+04, -0.1104E+05,  0.2806E+05
622      + , -0.4265E+03,  0.9660E+04,  0.5343E+04, -0.1174E+04
623      + , -0.3968E+02,  0.1459E+04,  0.3327E+05, -0.2229E+05
624      + , -0.3461E+01,  0.4401E+03,  0.4690E+06, -0.1943E+07
625      + ,  0.8564E+00, -0.3723E+03,  0.5020E+06, -0.1430E+07
626      + ,  0.2590E+00,  0.3237E+03,  0.2541E+06,  0.2321E+08,  16*0.0 /
627       DATA ((TMP(I,J, 35),I=1,MAXPOW),J=1,MAXINT) /
628      +   0.2129E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
629      + ,  0.1451E+05, -0.1941E+04,  0.8744E+02, -0.1214E+01
630      + , -0.3534E+04,  0.6788E+04, -0.1134E+04,  0.5165E+02
631      + , -0.3347E+03,  0.2195E+04,  0.1033E+04, -0.2839E+03
632      + ,  0.1461E+03, -0.1050E+04,  0.3457E+04,  0.2000E+05
633      + , -0.4318E+03,  0.1070E+05,  0.6548E+04, -0.1535E+04
634      + , -0.4033E+02,  0.1607E+04,  0.3681E+05, -0.2478E+05
635      + , -0.3955E+01,  0.5482E+03,  0.5202E+06, -0.2262E+07
636      + ,  0.9718E+00, -0.4209E+03,  0.5662E+06, -0.2099E+07
637      + ,  0.2907E+00,  0.3653E+03,  0.2908E+06,  0.2417E+08,  12*0.0 /
638       DATA ((TMP(I,J, 36),I=1,MAXPOW),J=1,MAXINT) /
639      +   0.3691E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
640      + ,  0.2024E+04, -0.2083E+03,  0.4439E+01,  0.1056E+00
641      + ,  0.7648E+03,  0.5257E+04, -0.9960E+03,  0.4686E+02
642      + , -0.2451E+03,  0.2199E+04,  0.1217E+04, -0.3120E+03
643      + ,  0.3611E+03, -0.2974E+04,  0.1308E+05,  0.1268E+05
644      + , -0.4163E+03,  0.1017E+05,  0.7535E+04, -0.1891E+04
645      + , -0.9924E+00,  0.8455E-01,  0.4446E+05, -0.3490E+05
646      + , -0.3076E+01,  0.4257E+03,  0.5613E+06, -0.2587E+07
647      + ,  0.1050E+01, -0.4435E+03,  0.6114E+06, -0.3028E+07
648      + ,  0.3136E+00,  0.3952E+03,  0.3252E+06,  0.2248E+08,  12*0.0 /
649       DATA ((TMP(I,J, 37),I=1,MAXPOW),J=1,MAXINT) /
650      +   0.4606E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
651      + ,  0.1441E+04, -0.1761E+03,  0.6821E+01,  0.4470E-01
652      + , -0.6763E+03,  0.6772E+04, -0.1372E+04,  0.7035E+02
653      + , -0.3457E+03,  0.2591E+04,  0.1242E+04, -0.3481E+03
654      + ,  0.3524E+03, -0.3194E+04,  0.1657E+05,  0.1224E+05
655      + , -0.4025E+03,  0.1067E+05,  0.8396E+04, -0.2275E+04
656      + ,  0.5234E+01,  0.7494E+03,  0.5033E+05, -0.4276E+05
657      + , -0.2207E+01,  0.2653E+03,  0.6269E+06, -0.3163E+07
658      + ,  0.1101E+01, -0.4280E+03,  0.6664E+06, -0.3481E+07
659      + ,  0.3467E+00,  0.4515E+03,  0.3525E+06,  0.2783E+08,  12*0.0 /
660       DATA ((TMP(I,J, 38),I=1,MAXPOW),J=1,MAXINT) /
661      +   0.4850E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
662      + ,  0.7854E+03, -0.6513E+02,  0.2134E+01,  0.1138E+00
663      + , -0.2484E+04,  0.8950E+04, -0.2019E+04,  0.1171E+03
664      + , -0.3209E+03,  0.2818E+04,  0.1344E+04, -0.3984E+03
665      + ,  0.2746E+03, -0.2830E+04,  0.2149E+05,  0.7655E+04
666      + , -0.4181E+03,  0.1025E+05,  0.1159E+05, -0.3199E+04
667      + , -0.4367E+02,  0.1966E+04,  0.4835E+05, -0.3788E+05
668      + , -0.5177E+00, -0.2546E+02,  0.6941E+06, -0.3733E+07
669      + ,  0.1296E+01, -0.5192E+03,  0.7384E+06, -0.5040E+07
670      + ,  0.3800E+00,  0.4950E+03,  0.4110E+06,  0.1955E+08,  12*0.0 /
671       DATA ((TMP(I,J, 39),I=1,MAXPOW),J=1,MAXINT) /
672      +   0.3854E+05, -0.1366E+04,  0.1672E+02, -0.6795E-01
673      + ,  0.2849E+04, -0.2828E+03,  0.1529E+02, -0.1723E+00
674      + , -0.6721E+04,  0.1290E+05, -0.3111E+04,  0.1985E+03
675      + ,  0.3585E+03, -0.4223E+03,  0.6620E+04, -0.2706E+04
676      + , -0.8813E+02,  0.1840E+04,  0.2198E+05, -0.1440E+04
677      + , -0.4141E+03,  0.1043E+05,  0.1385E+05, -0.4020E+04
678      + , -0.1282E+02,  0.1204E+04,  0.6028E+05, -0.5752E+05
679      + ,  0.4265E+00, -0.1646E+03,  0.7684E+06, -0.4421E+07
680      + ,  0.1334E+01, -0.4917E+03,  0.8066E+06, -0.5879E+07
681      + ,  0.4194E+00,  0.5651E+03,  0.4353E+06,  0.2981E+08,  12*0.0 /
682       DATA ((TMP(I,J, 40),I=1,MAXPOW),J=1,MAXINT) /
683      +   0.2993E+03,  0.1449E+03, -0.3012E+01,  0.1713E-01
684      + , -0.7040E+04,  0.1434E+05, -0.3926E+04,  0.2887E+03
685      + , -0.1647E+03,  0.2296E+04,  0.3166E+04, -0.1127E+04
686      + , -0.2164E+03,  0.3512E+04,  0.2053E+05, -0.3382E+04
687      + , -0.3996E+03,  0.1064E+05,  0.1548E+05, -0.4776E+04
688      + , -0.3011E+02,  0.1794E+04,  0.6153E+05, -0.5699E+05
689      + , -0.5900E+00,  0.1945E+03,  0.8031E+06, -0.4579E+07
690      + ,  0.1474E+01, -0.5530E+03,  0.8907E+06, -0.7927E+07
691      + ,  0.4562E+00,  0.6248E+03,  0.4760E+06,  0.3216E+08,  16*0.0 /
692       DATA ((TMP(I,J, 41),I=1,MAXPOW),J=1,MAXINT) /
693      +   0.3251E+04, -0.6844E+02,  0.1575E+01, -0.1094E-01
694      + ,  0.2764E+04, -0.4897E+03,  0.4237E+02, -0.8372E+00
695      + , -0.2018E+04,  0.8543E+04, -0.1769E+04,  0.2770E+02
696      + ,  0.2150E+03,  0.4490E+03,  0.6409E+04, -0.2429E+04
697      + , -0.2981E+03,  0.6188E+04,  0.1735E+05, -0.4420E+04
698      + , -0.3923E+03,  0.1102E+05,  0.1751E+05, -0.5733E+04
699      + , -0.1658E+02,  0.1435E+04,  0.7050E+05, -0.7037E+05
700      + , -0.2292E+01,  0.5516E+03,  0.8555E+06, -0.5018E+07
701      + ,  0.1599E+01, -0.5764E+03,  0.9643E+06, -0.8517E+07
702      + ,  0.4986E+00,  0.7224E+03,  0.4910E+06,  0.4104E+08,  12*0.0 /
703       DATA ((TMP(I,J, 42),I=1,MAXPOW),J=1,MAXINT) /
704      +  -0.3508E+04,  0.3050E+03, -0.4620E+01,  0.1998E-01
705      + ,  0.3053E+04, -0.6192E+03,  0.5104E+02, -0.9955E+00
706      + , -0.7786E+04,  0.1722E+05, -0.5215E+04,  0.3992E+03
707      + , -0.1114E+03,  0.1895E+04,  0.5405E+04, -0.2261E+04
708      + , -0.3027E+03,  0.7576E+04,  0.1570E+05, -0.4795E+04
709      + , -0.3784E+03,  0.1099E+05,  0.1983E+05, -0.6801E+04
710      + , -0.1840E+02,  0.1552E+04,  0.7635E+05, -0.8181E+05
711      + , -0.2256E+01,  0.7317E+03,  0.9005E+06, -0.5349E+07
712      + ,  0.1540E+01, -0.4834E+03,  0.1021E+07, -0.9300E+07
713      + ,  0.5359E+00,  0.7662E+03,  0.5243E+06,  0.5200E+08,  12*0.0 /
714       DATA ((TMP(I,J, 43),I=1,MAXPOW),J=1,MAXINT) /
715      +   0.2680E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
716      + , -0.7592E+03,  0.8677E+03, -0.1346E+03,  0.6198E+01
717      + , -0.3269E+05,  0.4110E+05, -0.1328E+05,  0.1384E+04
718      + , -0.3863E+03,  0.3458E+04,  0.3500E+04, -0.1341E+04
719      + , -0.2871E+03,  0.7198E+04,  0.1867E+05, -0.5767E+04
720      + , -0.1101E+03,  0.1625E+05,  0.3183E+04, -0.1491E+04
721      + , -0.2936E+02,  0.1924E+04,  0.7955E+05, -0.8669E+05
722      + ,  0.5721E+00,  0.7164E+02,  0.1004E+07, -0.6678E+07
723      + ,  0.1769E+01, -0.5691E+03,  0.1097E+07, -0.1069E+08
724      + ,  0.5748E+00,  0.8119E+03,  0.6113E+06,  0.3605E+08,  12*0.0 /
725       DATA ((TMP(I,J, 44),I=1,MAXPOW),J=1,MAXINT) /
726      +   0.2116E+04,  0.1001E+03, -0.2614E+01,  0.1466E-01
727      + ,  0.1947E+03,  0.7270E+03, -0.1490E+03,  0.8092E+01
728      + , -0.2141E+05,  0.2872E+05, -0.8592E+04,  0.8312E+03
729      + , -0.2430E+01,  0.2793E+04,  0.5237E+04, -0.2042E+04
730      + , -0.2831E+03,  0.7880E+04,  0.1913E+05, -0.6518E+04
731      + , -0.3547E+03,  0.1093E+05,  0.2575E+05, -0.9586E+04
732      + , -0.1311E+02,  0.1428E+04,  0.9237E+05, -0.1148E+06
733      + ,  0.1007E-01,  0.1277E+03,  0.1086E+07, -0.7755E+07
734      + ,  0.1633E+01, -0.4030E+03,  0.1144E+07, -0.9817E+07
735      + ,  0.6218E+00,  0.9117E+03,  0.5873E+06,  0.6614E+08,  12*0.0 /
736       DATA ((TMP(I,J, 45),I=1,MAXPOW),J=1,MAXINT) /
737      +   0.2535E+04,  0.1329E+03, -0.3193E+01,  0.1648E-01
738      + , -0.1037E+04,  0.1357E+04, -0.2412E+03,  0.1212E+02
739      + ,  0.1155E+04,  0.3456E+04,  0.8529E+03, -0.3075E+03
740      + , -0.4773E+03,  0.3678E+04,  0.5261E+04, -0.2225E+04
741      + , -0.2665E+03,  0.7197E+04,  0.2390E+05, -0.7946E+04
742      + , -0.3285E+03,  0.9342E+04,  0.3454E+05, -0.1203E+05
743      + , -0.2532E+02,  0.2009E+04,  0.9447E+05, -0.1115E+06
744      + , -0.3043E+01,  0.9079E+03,  0.1126E+07, -0.8086E+07
745      + ,  0.1870E+01, -0.4570E+03,  0.1251E+07, -0.1191E+08
746      + ,  0.6726E+00,  0.9644E+03,  0.7281E+06,  0.4409E+08,  12*0.0 /
747       DATA ((TMP(I,J, 46),I=1,MAXPOW),J=1,MAXINT) /
748      +   0.1095E+05, -0.2552E+03,  0.1914E+01, -0.3916E-02
749      + ,  0.2418E+04,  0.9664E+02, -0.1050E+03,  0.9047E+01
750      + , -0.1605E+05,  0.3107E+05, -0.1240E+05,  0.1658E+04
751      + , -0.4705E+03,  0.3846E+04,  0.5959E+04, -0.2639E+04
752      + , -0.2530E+03,  0.7129E+04,  0.2620E+05, -0.9068E+04
753      + , -0.3225E+03,  0.1001E+05,  0.3550E+05, -0.1354E+05
754      + ,  0.1642E+02,  0.4603E+03,  0.1181E+06, -0.1737E+06
755      + , -0.2888E+01,  0.8714E+03,  0.1203E+07, -0.9306E+07
756      + ,  0.1857E+01, -0.3662E+03,  0.1308E+07, -0.1219E+08
757      + ,  0.7151E+00,  0.1066E+04,  0.7318E+06,  0.6061E+08,  12*0.0 /
758       DATA ((TMP(I,J, 47),I=1,MAXPOW),J=1,MAXINT) /
759      +   0.1054E+05, -0.1973E+03,  0.1033E+01, -0.6574E-03
760      + , -0.4089E+04,  0.3916E+04, -0.8005E+03,  0.4841E+02
761      + , -0.2951E+04,  0.6256E+04,  0.1909E+04, -0.8620E+03
762      + , -0.2321E+02,  0.1892E+04,  0.9451E+04, -0.4150E+04
763      + , -0.2547E+03,  0.7851E+04,  0.2733E+05, -0.1038E+05
764      + , -0.3170E+03,  0.1045E+05,  0.3850E+05, -0.1560E+05
765      + , -0.1817E+01,  0.1265E+04,  0.1192E+06, -0.1684E+06
766      + , -0.4061E+01,  0.1339E+04,  0.1260E+07, -0.9752E+07
767      + ,  0.2065E+01, -0.4444E+03,  0.1433E+07, -0.1533E+08
768      + ,  0.7734E+00,  0.1157E+04,  0.8006E+06,  0.6214E+08,  12*0.0 /
769       DATA ((TMP(I,J, 48),I=1,MAXPOW),J=1,MAXINT) /
770      +   0.1876E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
771      + ,  0.8550E+04, -0.2785E+03,  0.2021E+01,  0.1032E-01
772      + ,  0.4933E+03,  0.1167E+04, -0.3452E+03,  0.2759E+02
773      + ,  0.2088E+04,  0.3036E+04,  0.1247E+04, -0.2691E+03
774      + ,  0.7508E+01,  0.2031E+04,  0.1014E+05, -0.4621E+04
775      + , -0.2454E+03,  0.8162E+04,  0.2822E+05, -0.1156E+05
776      + , -0.2677E+03,  0.8244E+04,  0.4935E+05, -0.1791E+05
777      + , -0.2306E+02,  0.2129E+04,  0.1161E+06, -0.1547E+06
778      + , -0.8206E+01,  0.2269E+04,  0.1268E+07, -0.9746E+07
779      + ,  0.2254E+01, -0.4904E+03,  0.1506E+07, -0.1642E+08
780      + ,  0.8123E+00,  0.1229E+04,  0.8686E+06,  0.5293E+08,  8*0.0 /
781       DATA ((TMP(I,J, 49),I=1,MAXPOW),J=1,MAXINT) /
782      +   0.1702E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
783      + ,  0.2333E+05, -0.1510E+04,  0.3498E+02, -0.2775E+00
784      + , -0.4354E+04,  0.4933E+04, -0.1221E+04,  0.9073E+02
785      + , -0.2048E+04,  0.9964E+04, -0.2857E+04,  0.7223E+03
786      + , -0.1763E+03,  0.3053E+04,  0.9756E+04, -0.4718E+04
787      + , -0.2324E+03,  0.7725E+04,  0.3333E+05, -0.1361E+05
788      + , -0.2675E+03,  0.8828E+04,  0.5249E+05, -0.2065E+05
789      + , -0.2149E+02,  0.2143E+04,  0.1256E+06, -0.1775E+06
790      + , -0.2097E+01,  0.1045E+04,  0.1428E+07, -0.1216E+08
791      + ,  0.2233E+01, -0.3526E+03,  0.1578E+07, -0.1747E+08
792      + ,  0.8685E+00,  0.1331E+04,  0.9171E+06,  0.6234E+08,  8*0.0 /
793       DATA ((TMP(I,J, 50),I=1,MAXPOW),J=1,MAXINT) /
794      +   0.1857E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
795      + ,  0.3071E+05, -0.2069E+04,  0.4886E+02, -0.3894E+00
796      + ,  0.6119E+02,  0.1769E+04, -0.5287E+03,  0.4438E+02
797      + ,  0.2325E+04,  0.2905E+04,  0.1395E+04, -0.1160E+03
798      + , -0.8411E+03,  0.7174E+04,  0.3429E+04, -0.1945E+04
799      + , -0.2078E+03,  0.6795E+04,  0.3953E+05, -0.1546E+05
800      + , -0.2735E+03,  0.1006E+05,  0.5154E+05, -0.2316E+05
801      + ,  0.1409E+02,  0.5245E+03,  0.1554E+06, -0.2738E+06
802      + , -0.6053E+01,  0.2153E+04,  0.1430E+07, -0.1218E+08
803      + ,  0.2434E+01, -0.4441E+03,  0.1694E+07, -0.2116E+08
804      + ,  0.9156E+00,  0.1438E+04,  0.9506E+06,  0.6990E+08,  8*0.0 /
805       DATA ((TMP(I,J, 51),I=1,MAXPOW),J=1,MAXINT) /
806      +   0.7076E+03,  0.0000E+00,  0.0000E+00,  0.0000E+00
807      + ,  0.2538E+05, -0.2014E+04,  0.5518E+02, -0.5093E+00
808      + , -0.2753E+04,  0.3809E+04, -0.9187E+03,  0.6404E+02
809      + ,  0.4191E+04,  0.1201E+04,  0.6138E+03,  0.8057E+03
810      + ,  0.9566E+02,  0.1292E+04,  0.1534E+05, -0.8151E+04
811      + , -0.2154E+03,  0.7884E+04,  0.3884E+05, -0.1742E+05
812      + , -0.2413E+03,  0.8640E+04,  0.6267E+05, -0.2629E+05
813      + , -0.2300E+02,  0.2351E+04,  0.1410E+06, -0.2130E+06
814      + , -0.9816E+01,  0.3157E+04,  0.1443E+07, -0.1193E+08
815      + ,  0.2470E+01, -0.3671E+03,  0.1779E+07, -0.2260E+08
816      + ,  0.9713E+00,  0.1520E+04,  0.1016E+07,  0.7466E+08,  8*0.0 /
817       DATA ((TMP(I,J, 52),I=1,MAXPOW),J=1,MAXINT) /
818      +   0.9416E+03,  0.0000E+00,  0.0000E+00,  0.0000E+00
819      + , -0.5818E+05,  0.1470E+05, -0.1039E+04,  0.2257E+02
820      + , -0.2734E+04,  0.3947E+04, -0.9985E+03,  0.7313E+02
821      + , -0.1058E+06,  0.2917E+06, -0.2491E+06,  0.7124E+05
822      + , -0.4397E+02,  0.2059E+04,  0.1571E+05, -0.9021E+04
823      + , -0.1946E+03,  0.7129E+04,  0.4399E+05, -0.1939E+05
824      + , -0.2372E+03,  0.9166E+04,  0.6308E+05, -0.2911E+05
825      + ,  0.1086E+01,  0.1149E+04,  0.1640E+06, -0.2950E+06
826      + , -0.3401E+01,  0.1832E+04,  0.1571E+07, -0.1426E+08
827      + ,  0.2389E+01, -0.1733E+03,  0.1798E+07, -0.2273E+08
828      + ,  0.1007E+01,  0.1586E+04,  0.1076E+07,  0.7154E+08,  8*0.0 /
829       DATA ((TMP(I,J, 53),I=1,MAXPOW),J=1,MAXINT) /
830      +   0.7020E+03,  0.0000E+00,  0.0000E+00,  0.0000E+00
831      + ,  0.2410E+05, -0.2114E+04,  0.6324E+02, -0.6266E+00
832      + , -0.2554E+04,  0.4115E+04, -0.1059E+04,  0.7800E+02
833      + , -0.4564E+05,  0.1181E+06, -0.9050E+05,  0.2490E+05
834      + ,  0.1940E+05, -0.1268E+05, -0.1443E+05,  0.1681E+05
835      + , -0.1016E+04,  0.1150E+05, -0.8322E+04,  0.1024E+05
836      + , -0.1833E+03,  0.6763E+04,  0.5201E+05, -0.2240E+05
837      + , -0.2150E+03,  0.8275E+04,  0.7525E+05, -0.3297E+05
838      + , -0.2566E+02,  0.2615E+04,  0.1595E+06, -0.2556E+06
839      + , -0.3670E+01,  0.1872E+04,  0.1719E+07, -0.1646E+08
840      + ,  0.2225E+01,  0.1829E+03,  0.1881E+07, -0.2157E+08
841      + ,  0.1098E+01,  0.1755E+04,  0.1154E+07,  0.8971E+08,  4*0.0 /
842       DATA ((TMP(I,J, 54),I=1,MAXPOW),J=1,MAXINT) /
843      +   0.6262E+03,  0.0000E+00,  0.0000E+00,  0.0000E+00
844      + ,  0.7578E+05, -0.1206E+05,  0.6503E+03, -0.1177E+02
845      + , -0.1456E+04,  0.3607E+04, -0.1050E+04,  0.8562E+02
846      + ,  0.5494E+05, -0.4316E+05, -0.3353E+05,  0.3122E+05
847      + , -0.2422E+05,  0.6454E+05, -0.3914E+05,  0.8599E+04
848      + , -0.6606E+03,  0.7269E+04,  0.6119E+04, -0.3383E+04
849      + , -0.1565E+03,  0.5749E+04,  0.5985E+05, -0.2396E+05
850      + , -0.2236E+03,  0.9444E+04,  0.7463E+05, -0.3742E+05
851      + , -0.2010E+02,  0.2400E+04,  0.1720E+06, -0.2977E+06
852      + , -0.1622E+01,  0.1830E+04,  0.1780E+07, -0.1721E+08
853      + ,  0.2082E+01,  0.3625E+03,  0.1962E+07, -0.2448E+08
854      + ,  0.1148E+01,  0.1892E+04,  0.1133E+07,  0.1244E+09,  4*0.0 /
855       DATA ((TMP(I,J, 55),I=1,MAXPOW),J=1,MAXINT) /
856      +   0.6813E+03,  0.0000E+00,  0.0000E+00,  0.0000E+00
857      + ,  0.1642E+06, -0.2922E+05,  0.1724E+04, -0.3361E+02
858      + , -0.2400E+04,  0.4936E+04, -0.1531E+04,  0.1354E+03
859      + ,  0.1295E+05, -0.7423E+04, -0.6620E+04,  0.9125E+04
860      + ,  0.6121E+03, -0.1572E+04,  0.2712E+05, -0.1563E+05
861      + , -0.8402E+02,  0.2872E+04,  0.7699E+05, -0.1853E+05
862      + , -0.2254E+03,  0.1010E+05,  0.7940E+05, -0.4284E+05
863      + , -0.1987E+01,  0.1400E+04,  0.2040E+06, -0.4166E+06
864      + , -0.4238E+01,  0.2444E+04,  0.1873E+07, -0.1899E+08
865      + ,  0.2278E+01,  0.4044E+03,  0.2085E+07, -0.2635E+08
866      + ,  0.1226E+01,  0.2030E+04,  0.1249E+07,  0.1169E+09,  8*0.0 /
867       DATA ((TMP(I,J, 56),I=1,MAXPOW),J=1,MAXINT) /
868      +   0.5765E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
869      + ,  0.1759E+04, -0.2331E+03,  0.9598E+01,  0.3960E-01
870      + , -0.4276E+03,  0.2995E+04, -0.9837E+03,  0.8771E+02
871      + ,  0.6154E+04,  0.1190E+04, -0.2858E+03,  0.1593E+04
872      + ,  0.2420E+03, -0.4033E+03,  0.3138E+05, -0.2304E+05
873      + , -0.1084E+03,  0.4043E+04,  0.7718E+05, -0.2532E+05
874      + , -0.1836E+03,  0.7952E+04,  0.9545E+05, -0.4532E+05
875      + ,  0.8708E+01,  0.7581E+03,  0.2258E+06, -0.5011E+06
876      + , -0.5846E+01,  0.3162E+04,  0.1902E+07, -0.1948E+08
877      + ,  0.2590E+01,  0.2297E+03,  0.2236E+07, -0.3193E+08
878      + ,  0.1281E+01,  0.2144E+04,  0.1303E+07,  0.1193E+09,  8*0.0 /
879       DATA ((TMP(I,J, 57),I=1,MAXPOW),J=1,MAXINT) /
880      +   0.4514E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
881      + ,  0.2741E+05, -0.4218E+04,  0.2117E+03, -0.3183E+01
882      + , -0.1231E+04,  0.4053E+04, -0.1329E+04,  0.1203E+03
883      + ,  0.2431E+04,  0.3780E+04,  0.2792E+04, -0.2167E+03
884      + , -0.4024E+03,  0.8837E+04,  0.2398E+04, -0.4154E+03
885      + , -0.2694E+03,  0.4549E+04,  0.1944E+05, -0.1264E+05
886      + , -0.9085E+02,  0.3459E+04,  0.8639E+05, -0.2558E+05
887      + , -0.1485E+03,  0.6357E+04,  0.1121E+06, -0.4626E+05
888      + , -0.6527E+01,  0.1901E+04,  0.2209E+06, -0.4602E+06
889      + , -0.7588E+01,  0.3685E+04,  0.1991E+07, -0.2054E+08
890      + ,  0.2727E+01,  0.3137E+03,  0.2346E+07, -0.3362E+08
891      + ,  0.1365E+01,  0.2298E+04,  0.1393E+07,  0.1220E+09,  4*0.0 /
892       DATA ((TMP(I,J, 58),I=1,MAXPOW),J=1,MAXINT) /
893      +   0.1542E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
894      + ,  0.2448E+05, -0.5012E+04,  0.3411E+03, -0.7413E+01
895      + , -0.1622E+04,  0.4662E+04, -0.1505E+04,  0.1385E+03
896      + ,  0.5551E+05, -0.1769E+06,  0.2024E+06, -0.7125E+05
897      + ,  0.4656E+04,  0.6580E+04, -0.3256E+04,  0.8079E+03
898      + , -0.2828E+03,  0.1019E+05,  0.2005E+04, -0.4180E+03
899      + ,  0.3264E+03, -0.7440E+03,  0.3609E+05,  0.2685E+05
900      + , -0.7263E+02,  0.2822E+04,  0.9606E+05, -0.2458E+05
901      + , -0.1695E+03,  0.8899E+04,  0.1139E+06, -0.5644E+05
902      + ,  0.1939E+02,  0.3799E+03,  0.2661E+06, -0.6583E+06
903      + , -0.1029E+02,  0.4409E+04,  0.2101E+07, -0.2280E+08
904      + ,  0.2723E+01,  0.4916E+03,  0.2494E+07, -0.3595E+08
905      + ,  0.1464E+01,  0.2444E+04,  0.1486E+07,  0.1376E+09
906      +/
907       DATA ((TMP(I,J, 59),I=1,MAXPOW),J=1,MAXINT) /
908      +   0.3935E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
909      + ,  0.5885E+02,  0.5060E+02,  0.1142E+02, -0.2858E+00
910      + , -0.2475E+04,  0.8072E+04, -0.2077E+04,  0.2099E+03
911      + , -0.6489E+04,  0.2485E+05, -0.1060E+05,  0.2589E+04
912      + , -0.5349E+03,  0.8049E+04,  0.5499E+04, -0.9313E+03
913      + , -0.3829E+03,  0.1003E+05,  0.3437E+04, -0.7202E+03
914      + ,  0.1158E+03,  0.1754E+04,  0.2978E+05, -0.2081E+05
915      + , -0.4523E+02,  0.1773E+04,  0.1075E+06, -0.1884E+05
916      + , -0.1401E+03,  0.6504E+04,  0.1315E+06, -0.5758E+05
917      + ,  0.8492E+01,  0.8664E+03,  0.2819E+06, -0.7256E+06
918      + , -0.9126E+01,  0.4223E+04,  0.2290E+07, -0.2709E+08
919      + ,  0.3207E+01,  0.3424E+03,  0.2697E+07, -0.4123E+08
920      + ,  0.1557E+01,  0.2710E+04,  0.1566E+07,  0.1383E+09
921      +/
922       DATA ((TMP(I,J, 60),I=1,MAXPOW),J=1,MAXINT) /
923      +   0.2269E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
924      + , -0.1142E+04,  0.4035E+03, -0.1726E+02,  0.2479E+00
925      + , -0.4138E+03,  0.3657E+04, -0.1265E+04,  0.1287E+03
926      + ,  0.1241E+05, -0.6608E+04, -0.8545E+04,  0.1174E+05
927      + ,  0.2239E+03,  0.1119E+05, -0.1181E+04,  0.2885E+03
928      + , -0.3530E+03,  0.6609E+04,  0.1641E+05, -0.7462E+04
929      + , -0.6644E+02,  0.2796E+04,  0.1098E+06, -0.2924E+05
930      + , -0.1114E+03,  0.5206E+04,  0.1457E+06, -0.5559E+05
931      + ,  0.8287E+01,  0.9600E+03,  0.2935E+06, -0.7637E+06
932      + , -0.9944E+01,  0.4886E+04,  0.2320E+07, -0.2702E+08
933      + ,  0.3820E+01,  0.8492E+02,  0.2868E+07, -0.4750E+08
934      + ,  0.1633E+01,  0.2833E+04,  0.1759E+07,  0.9295E+08,  4*0.0 /
935       DATA ((TMP(I,J, 61),I=1,MAXPOW),J=1,MAXINT) /
936      +   0.2285E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
937      + ,  0.3175E+04, -0.4311E+03,  0.2911E+02, -0.5124E+00
938      + , -0.7108E+03,  0.4102E+04, -0.1397E+04,  0.1482E+03
939      + , -0.5546E+03,  0.7912E+04,  0.4958E+04, -0.7606E+03
940      + , -0.5157E+03,  0.8937E+04,  0.6066E+04, -0.1148E+04
941      + , -0.4607E+03,  0.9760E+04,  0.6028E+04, -0.1298E+04
942      + , -0.1912E+03,  0.5138E+04,  0.2120E+05, -0.1345E+05
943      + , -0.7165E+02,  0.3169E+04,  0.1166E+06, -0.3499E+05
944      + , -0.2589E+03,  0.8215E+04, -0.1102E+06,  0.1824E+07
945      + , -0.4251E+02,  0.4740E+04,  0.2377E+06, -0.4279E+06
946      + ,  0.9003E-01,  0.2739E+04,  0.2641E+07, -0.3336E+08
947      + ,  0.2857E+01,  0.1090E+04,  0.2888E+07, -0.4432E+08
948      + ,  0.1742E+01,  0.3069E+04,  0.1744E+07,  0.1719E+09
949      +/
950       DATA ((TMP(I,J, 62),I=1,MAXPOW),J=1,MAXINT) /
951      +   0.1754E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
952      + ,  0.1174E+04,  0.1370E+03, -0.4815E+01,  0.3574E-01
953      + , -0.1315E+04,  0.5218E+04, -0.1900E+04,  0.2181E+03
954      + , -0.5499E+03,  0.7881E+04,  0.5664E+04, -0.8952E+03
955      + , -0.1126E+03,  0.3912E+04,  0.3230E+05, -0.2774E+05
956      + , -0.3580E+02,  0.1576E+04,  0.1290E+06, -0.2201E+05
957      + , -0.6400E+02,  0.3084E+04,  0.1751E+06, -0.4595E+05
958      + ,  0.7096E+01,  0.1071E+04,  0.3267E+06, -0.9193E+06
959      + , -0.1082E+02,  0.5919E+04,  0.2475E+07, -0.2997E+08
960      + ,  0.3844E+01,  0.4077E+03,  0.3137E+07, -0.5576E+08
961      + ,  0.1799E+01,  0.3181E+04,  0.1898E+07,  0.1263E+09,  8*0.0 /
962       DATA ((TMP(I,J, 63),I=1,MAXPOW),J=1,MAXINT) /
963      +   0.2047E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
964      + ,  0.4276E+04, -0.4382E+03,  0.2276E+02, -0.3499E+00
965      + , -0.4678E+03,  0.4034E+04, -0.1326E+04,  0.1388E+03
966      + , -0.5507E+03,  0.7936E+04,  0.6518E+04, -0.1058E+04
967      + , -0.4589E+03,  0.1019E+05,  0.5669E+04, -0.1234E+04
968      + , -0.4913E+03,  0.8740E+04,  0.1016E+05, -0.2142E+04
969      + , -0.2288E+03,  0.5934E+04,  0.2154E+05, -0.1326E+05
970      + , -0.2466E+02,  0.1115E+04,  0.1400E+06, -0.1768E+05
971      + , -0.6266E+02,  0.3143E+04,  0.1878E+06, -0.5089E+05
972      + , -0.4306E+01,  0.2026E+04,  0.3280E+06, -0.9012E+06
973      + , -0.7610E+01,  0.5488E+04,  0.2642E+07, -0.3287E+08
974      + ,  0.3783E+01,  0.7655E+03,  0.3245E+07, -0.5731E+08
975      + ,  0.1904E+01,  0.3400E+04,  0.2020E+07,  0.1315E+09
976      +/
977       DATA ((TMP(I,J, 64),I=1,MAXPOW),J=1,MAXINT) /
978      +   0.1620E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
979      + ,  0.1627E+04,  0.8320E+02, -0.5112E+01,  0.7828E-01
980      + , -0.2404E+04,  0.7220E+04, -0.2684E+04,  0.2976E+03
981      + , -0.5394E+03,  0.7849E+04,  0.7295E+04, -0.1216E+04
982      + , -0.4486E+03,  0.1029E+05,  0.6091E+04, -0.1380E+04
983      + , -0.4679E+03,  0.9424E+04,  0.9571E+04, -0.2201E+04
984      + , -0.3017E+03,  0.6938E+04,  0.1944E+05, -0.1104E+05
985      + , -0.5768E+02,  0.2841E+04,  0.1366E+06, -0.4131E+05
986      + , -0.3358E+02,  0.1701E+04,  0.2029E+06, -0.3312E+05
987      + , -0.8668E+01,  0.2379E+04,  0.3331E+06, -0.9285E+06
988      + ,  0.3209E+01,  0.1579E+04,  0.3144E+07, -0.4804E+08
989      + ,  0.1966E+01,  0.3553E+04,  0.2103E+07,  0.1345E+09,  4*0.0 /
990       DATA ((TMP(I,J, 65),I=1,MAXPOW),J=1,MAXINT) /
991      +   0.1879E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
992      + ,  0.3021E+04, -0.1618E+03,  0.9979E+01, -0.1862E+00
993      + , -0.1406E+04,  0.5621E+04, -0.1912E+04,  0.2075E+03
994      + , -0.5345E+03,  0.7632E+04,  0.8660E+04, -0.1449E+04
995      + , -0.4988E+03,  0.9401E+04,  0.9226E+04, -0.2000E+04
996      + , -0.4700E+03,  0.9391E+04,  0.1117E+05, -0.2613E+04
997      + , -0.1902E+03,  0.5724E+04,  0.2615E+05, -0.1782E+05
998      + , -0.1568E+03,  0.1116E+05,  0.8757E+05, -0.7647E+05
999      + , -0.4533E+02,  0.2421E+04,  0.2126E+06, -0.4859E+05
1000      + , -0.4721E+01,  0.2127E+04,  0.3644E+06, -0.1109E+07
1001      + ,  0.3903E+01,  0.1250E+04,  0.3405E+07, -0.5637E+08
1002      + ,  0.2077E+01,  0.3763E+04,  0.2261E+07,  0.1115E+09,  4*0.0 /
1003       DATA ((TMP(I,J, 66),I=1,MAXPOW),J=1,MAXINT) /
1004      +   0.3294E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1005      + ,  0.3167E+04, -0.4242E+02,  0.5269E+01, -0.1178E+00
1006      + , -0.1109E+04,  0.5234E+04, -0.1650E+04,  0.1694E+03
1007      + , -0.5242E+03,  0.7751E+04,  0.9341E+04, -0.1624E+04
1008      + , -0.4635E+03,  0.1050E+05,  0.8066E+04, -0.1930E+04
1009      + , -0.4608E+03,  0.9448E+04,  0.1216E+05, -0.2942E+04
1010      + , -0.1077E+03,  0.4675E+04,  0.3253E+05, -0.2519E+05
1011      + , -0.6326E+02,  0.3407E+04,  0.1499E+06, -0.5526E+05
1012      + , -0.7262E+02,  0.4176E+04,  0.2130E+06, -0.8026E+05
1013      + , -0.1800E+02,  0.3374E+04,  0.3526E+06, -0.9982E+06
1014      + ,  0.4112E+01,  0.1344E+04,  0.3520E+07, -0.5867E+08
1015      + ,  0.2166E+01,  0.3928E+04,  0.2396E+07,  0.1007E+09,  4*0.0 /
1016       DATA ((TMP(I,J, 67),I=1,MAXPOW),J=1,MAXINT) /
1017      +   0.1647E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1018      + ,  0.2949E+04, -0.5473E+02, -0.1313E+01,  0.5447E-01
1019      + , -0.1961E+04,  0.6859E+04, -0.2353E+04,  0.2538E+03
1020      + , -0.5221E+03,  0.8169E+04,  0.9852E+04, -0.1807E+04
1021      + , -0.4878E+03,  0.9589E+04,  0.1119E+05, -0.2580E+04
1022      + , -0.4549E+03,  0.9349E+04,  0.1378E+05, -0.3393E+04
1023      + , -0.1379E+03,  0.5032E+04,  0.3428E+05, -0.2787E+05
1024      + , -0.4939E+02,  0.2705E+04,  0.1628E+06, -0.5056E+05
1025      + , -0.4426E+02,  0.2555E+04,  0.2349E+06, -0.5948E+05
1026      + , -0.2988E+01,  0.2094E+04,  0.4072E+06, -0.1335E+07
1027      + ,  0.5131E+01,  0.5023E+03,  0.4000E+07, -0.7974E+08
1028      + ,  0.2274E+01,  0.4258E+04,  0.2387E+07,  0.1450E+09,  4*0.0 /
1029       DATA ((TMP(I,J, 68),I=1,MAXPOW),J=1,MAXINT) /
1030      +   0.2419E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1031      + ,  0.3178E+04, -0.8940E+02,  0.5741E+01, -0.1173E+00
1032      + , -0.1598E+04,  0.6371E+04, -0.2015E+04,  0.1951E+03
1033      + , -0.5148E+03,  0.8464E+04,  0.1038E+05, -0.1998E+04
1034      + , -0.4781E+03,  0.8546E+04,  0.1463E+05, -0.3195E+04
1035      + , -0.4470E+03,  0.9490E+04,  0.1489E+05, -0.3799E+04
1036      + , -0.2279E+03,  0.6501E+04,  0.3083E+05, -0.2363E+05
1037      + , -0.5162E+02,  0.2947E+04,  0.1710E+06, -0.5828E+05
1038      + , -0.5360E+02,  0.3254E+04,  0.2449E+06, -0.7821E+05
1039      + ,  0.6456E+01,  0.1211E+04,  0.4583E+06, -0.1684E+07
1040      + , -0.4025E+01,  0.2031E+04,  0.3795E+07, -0.6651E+08
1041      + ,  0.2386E+01,  0.4497E+04,  0.2558E+07,  0.1403E+09,  4*0.0 /
1042       DATA ((TMP(I,J, 69),I=1,MAXPOW),J=1,MAXINT) /
1043      +   0.1729E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1044      + ,  0.5004E+04, -0.3383E+03,  0.1358E+02, -0.1932E+00
1045      + , -0.1721E+04,  0.6821E+04, -0.2191E+04,  0.2162E+03
1046      + , -0.5074E+03,  0.8420E+04,  0.1161E+05, -0.2277E+04
1047      + ,  0.4504E+03, -0.4717E+04,  0.9421E+05, -0.1104E+06
1048      + , -0.5189E+02,  0.3077E+04,  0.1808E+06, -0.6499E+05
1049      + , -0.5311E+02,  0.3348E+04,  0.2598E+06, -0.8657E+05
1050      + , -0.3887E+01,  0.2132E+04,  0.4637E+06, -0.1723E+07
1051      + ,  0.4124E+01,  0.2299E+04,  0.3956E+07, -0.7043E+08
1052      + ,  0.2512E+01,  0.4774E+04,  0.2690E+07,  0.1454E+09,  12*0.0 /
1053       DATA ((TMP(I,J, 70),I=1,MAXPOW),J=1,MAXINT) /
1054      +   0.2440E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1055      + ,  0.5923E+04, -0.5235E+03,  0.2450E+02, -0.3644E+00
1056      + , -0.1402E+04,  0.6524E+04, -0.1996E+04,  0.1849E+03
1057      + , -0.4649E+03,  0.1046E+05,  0.8110E+04, -0.1188E+04
1058      + ,  0.3077E+03, -0.2145E+04,  0.8374E+05, -0.9504E+05
1059      + , -0.5328E+02,  0.3291E+04,  0.1869E+06, -0.7332E+05
1060      + , -0.7311E+02,  0.4929E+04,  0.2602E+06, -0.1239E+06
1061      + , -0.2434E+02,  0.4269E+04,  0.4243E+06, -0.1379E+07
1062      + ,  0.4575E+01,  0.2312E+04,  0.4086E+07, -0.7562E+08
1063      + ,  0.2606E+01,  0.4980E+04,  0.2895E+07,  0.9923E+08,  12*0.0 /
1064       DATA ((TMP(I,J, 71),I=1,MAXPOW),J=1,MAXINT) /
1065      +   0.2123E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1066      + ,  0.8368E+04, -0.1105E+04,  0.5971E+02, -0.9746E+00
1067      + , -0.1536E+04,  0.7063E+04, -0.2291E+04,  0.2250E+03
1068      + , -0.4936E+03,  0.9095E+04,  0.1255E+05, -0.2717E+04
1069      + ,  0.2304E+03, -0.1379E+04,  0.8594E+05, -0.1012E+06
1070      + , -0.5213E+02,  0.3333E+04,  0.1974E+06, -0.7977E+05
1071      + , -0.5333E+02,  0.3627E+04,  0.2858E+06, -0.1071E+06
1072      + ,  0.3726E+01,  0.1486E+04,  0.5352E+06, -0.2222E+07
1073      + ,  0.5266E+01,  0.2137E+04,  0.4329E+07, -0.8327E+08
1074      + ,  0.2736E+01,  0.5209E+04,  0.3153E+07,  0.5239E+08,  12*0.0 /
1075       DATA ((TMP(I,J, 72),I=1,MAXPOW),J=1,MAXINT) /
1076      +   0.5463E+04, -0.2212E+03,  0.3566E+01, -0.1819E-01
1077      + ,  0.9488E+04, -0.1416E+04,  0.8076E+02, -0.1382E+01
1078      + , -0.1796E+04,  0.7903E+04, -0.2748E+04,  0.2932E+03
1079      + , -0.4835E+03,  0.9549E+04,  0.1284E+05, -0.2957E+04
1080      + ,  0.3480E+03, -0.3031E+04,  0.9762E+05, -0.1203E+06
1081      + , -0.6090E+02,  0.4112E+04,  0.2015E+06, -0.9934E+05
1082      + , -0.7386E+02,  0.5387E+04,  0.2859E+06, -0.1536E+06
1083      + , -0.1873E+02,  0.4102E+04,  0.4698E+06, -0.1572E+07
1084      + ,  0.4753E+01,  0.2673E+04,  0.4462E+07, -0.8960E+08
1085      + ,  0.2845E+01,  0.5554E+04,  0.3013E+07,  0.1535E+09,  12*0.0 /
1086       DATA ((TMP(I,J, 73),I=1,MAXPOW),J=1,MAXINT) /
1087      +  -0.9358E+04,  0.4382E+03, -0.5661E+01,  0.2301E-01
1088      + ,  0.3537E+05, -0.2591E+04,  0.6585E+02, -0.5643E+00
1089      + ,  0.9000E+04, -0.1376E+04,  0.8354E+02, -0.1487E+01
1090      + , -0.1944E+04,  0.8490E+04, -0.3040E+04,  0.3330E+03
1091      + , -0.4771E+03,  0.1003E+05,  0.1338E+05, -0.3254E+04
1092      + ,  0.1047E+03,  0.6524E+03,  0.8380E+05, -0.1025E+06
1093      + , -0.5521E+02,  0.3828E+04,  0.2145E+06, -0.1019E+06
1094      + , -0.4987E+02,  0.3641E+04,  0.3160E+06, -0.1250E+06
1095      + , -0.2176E+01,  0.2172E+04,  0.5708E+06, -0.2482E+07
1096      + ,  0.5194E+01,  0.2794E+04,  0.4645E+07, -0.9648E+08
1097      + ,  0.2973E+01,  0.5785E+04,  0.3320E+07,  0.9608E+08,  8*0.0 /
1098       DATA ((TMP(I,J, 74),I=1,MAXPOW),J=1,MAXINT) /
1099      +  -0.2084E+04,  0.1834E+03, -0.2686E+01,  0.1131E-01
1100      + , -0.1398E+05,  0.2451E+04, -0.1057E+03,  0.1381E+01
1101      + ,  0.9985E+04, -0.1816E+04,  0.1301E+03, -0.2747E+01
1102      + , -0.1826E+04,  0.8623E+04, -0.3114E+04,  0.3432E+03
1103      + , -0.4666E+03,  0.1057E+05,  0.1355E+05, -0.3503E+04
1104      + ,  0.1412E+03,  0.2762E+03,  0.8940E+05, -0.1129E+06
1105      + , -0.6206E+02,  0.4525E+04,  0.2194E+06, -0.1223E+06
1106      + , -0.1211E+03,  0.1058E+05,  0.2701E+06, -0.2556E+06
1107      + , -0.2105E+02,  0.4623E+04,  0.5109E+06, -0.1184E+07
1108      + ,  0.5276E+01,  0.2999E+04,  0.4828E+07, -0.1030E+09
1109      + ,  0.3099E+01,  0.6112E+04,  0.3348E+07,  0.1310E+09,  8*0.0 /
1110       DATA ((TMP(I,J, 75),I=1,MAXPOW),J=1,MAXINT) /
1111      +  -0.5759E+04,  0.4019E+03, -0.6479E+01,  0.3113E-01
1112      + ,  0.1235E+05, -0.2571E+04,  0.1879E+03, -0.4072E+01
1113      + , -0.1846E+04,  0.9196E+04, -0.3553E+04,  0.4261E+03
1114      + , -0.4438E+03,  0.1148E+05,  0.1253E+05, -0.3520E+04
1115      + ,  0.1961E+03, -0.3585E+03,  0.9550E+05, -0.1180E+06
1116      + , -0.5721E+02,  0.4284E+04,  0.2331E+06, -0.1270E+06
1117      + , -0.1192E+03,  0.1081E+05,  0.2828E+06, -0.2788E+06
1118      + , -0.8548E+01,  0.3122E+04,  0.5972E+06, -0.2571E+07
1119      + ,  0.4901E+01,  0.3621E+04,  0.4969E+07, -0.1085E+09
1120      + ,  0.3238E+01,  0.6600E+04,  0.3237E+07,  0.2204E+09,  12*0.0 /
1121       DATA ((TMP(I,J, 76),I=1,MAXPOW),J=1,MAXINT) /
1122      +   0.5189E+03,  0.1085E+03, -0.2362E+01,  0.1285E-01
1123      + ,  0.9943E+04, -0.1964E+04,  0.1499E+03, -0.3533E+01
1124      + , -0.1915E+04,  0.9625E+04, -0.3819E+04,  0.4641E+03
1125      + , -0.4312E+03,  0.1174E+05,  0.1298E+05, -0.3813E+04
1126      + , -0.7620E+02,  0.3809E+04,  0.8189E+05, -0.1090E+06
1127      + , -0.5946E+02,  0.4643E+04,  0.2401E+06, -0.1434E+06
1128      + , -0.1116E+03,  0.1035E+05,  0.3012E+06, -0.2961E+06
1129      + , -0.1803E+02,  0.4548E+04,  0.5700E+06, -0.2246E+07
1130      + ,  0.5735E+01,  0.3318E+04,  0.5232E+07, -0.1209E+09
1131      + ,  0.3351E+01,  0.6829E+04,  0.3509E+07,  0.1607E+09,  12*0.0 /
1132       DATA ((TMP(I,J, 77),I=1,MAXPOW),J=1,MAXINT) /
1133      +  -0.9441E+02,  0.1284E+03, -0.2563E+01,  0.1347E-01
1134      + ,  0.9926E+04, -0.2049E+04,  0.1630E+03, -0.4025E+01
1135      + , -0.2391E+04,  0.1087E+05, -0.4456E+04,  0.5612E+03
1136      + , -0.4371E+03,  0.1150E+05,  0.1539E+05, -0.4560E+04
1137      + ,  0.1540E+03, -0.7160E+03,  0.1147E+06, -0.1673E+06
1138      + , -0.5153E+02,  0.4097E+04,  0.2582E+06, -0.1419E+06
1139      + , -0.1366E+03,  0.1456E+05,  0.2707E+06, -0.3461E+06
1140      + , -0.1995E+02,  0.4906E+04,  0.6064E+06, -0.2648E+07
1141      + ,  0.6213E+01,  0.3458E+04,  0.5442E+07, -0.1286E+09
1142      + ,  0.3505E+01,  0.7126E+04,  0.3805E+07,  0.1114E+09,  12*0.0 /
1143       DATA ((TMP(I,J, 78),I=1,MAXPOW),J=1,MAXINT) /
1144      +   0.1220E+05, -0.3768E+03,  0.4167E+01, -0.1564E-01
1145      + ,  0.1396E+05, -0.1184E+04,  0.3650E+02, -0.3548E+00
1146      + ,  0.1172E+05, -0.2752E+04,  0.2212E+03, -0.4579E+01
1147      + , -0.2373E+04,  0.1129E+05, -0.4739E+04,  0.6166E+03
1148      + , -0.4333E+03,  0.1147E+05,  0.1710E+05, -0.5181E+04
1149      + , -0.7755E+02,  0.2847E+04,  0.1047E+06, -0.1599E+06
1150      + , -0.6719E+02,  0.5737E+04,  0.2565E+06, -0.1880E+06
1151      + , -0.1399E+03,  0.1605E+05,  0.2671E+06, -0.3742E+06
1152      + ,  0.1076E+02,  0.9617E+03,  0.7776E+06, -0.4207E+07
1153      + ,  0.6088E+01,  0.4001E+04,  0.5552E+07, -0.1300E+09
1154      + ,  0.3646E+01,  0.7453E+04,  0.3931E+07,  0.1226E+09,  8*0.0 /
1155       DATA ((TMP(I,J, 79),I=1,MAXPOW),J=1,MAXINT) /
1156      +   0.4105E+04,  0.2410E+02, -0.1837E+01,  0.1252E-01
1157      + ,  0.1581E+05, -0.4931E+04,  0.5197E+03, -0.1622E+02
1158      + , -0.2231E+04,  0.1142E+05, -0.4781E+04,  0.6012E+03
1159      + , -0.4203E+03,  0.1212E+05,  0.1690E+05, -0.5448E+04
1160      + ,  0.1443E+03, -0.4096E+03,  0.1249E+06, -0.1944E+06
1161      + , -0.8194E+02,  0.7574E+04,  0.2506E+06, -0.2318E+06
1162      + , -0.1134E+03,  0.1194E+05,  0.3355E+06, -0.3921E+06
1163      + , -0.3392E+01,  0.3070E+04,  0.7247E+06, -0.3499E+07
1164      + ,  0.7917E+01,  0.2918E+04,  0.6095E+07, -0.1596E+09
1165      + ,  0.3811E+01,  0.7846E+04,  0.4248E+07,  0.4571E+08,  12*0.0 /
1166       DATA ((TMP(I,J, 80),I=1,MAXPOW),J=1,MAXINT) /
1167      +   0.5142E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1168      + ,  0.1780E+05, -0.5970E+04,  0.6765E+03, -0.2176E+02
1169      + , -0.2317E+04,  0.1206E+05, -0.5219E+04,  0.7004E+03
1170      + , -0.3961E+03,  0.1286E+05,  0.1582E+05, -0.5482E+04
1171      + ,  0.2241E+03, -0.2195E+04,  0.1429E+06, -0.2325E+06
1172      + , -0.7533E+02,  0.7084E+04,  0.2675E+06, -0.2412E+06
1173      + , -0.1288E+03,  0.1524E+05,  0.3090E+06, -0.4397E+06
1174      + ,  0.4681E+00,  0.2405E+04,  0.7958E+06, -0.4232E+07
1175      + ,  0.6983E+01,  0.4309E+04,  0.5980E+07, -0.1518E+09
1176      + ,  0.3946E+01,  0.8311E+04,  0.4256E+07,  0.8942E+08,  12*0.0 /
1177       DATA ((TMP(I,J, 81),I=1,MAXPOW),J=1,MAXINT) /
1178      +   0.1823E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1179      + , -0.1010E+05,  0.1911E+04, -0.7721E+02,  0.9111E+00
1180      + ,  0.1488E+05, -0.4878E+04,  0.5536E+03, -0.1923E+02
1181      + , -0.2399E+04,  0.1271E+05, -0.5745E+04,  0.7873E+03
1182      + , -0.4087E+03,  0.1156E+05,  0.2169E+05, -0.7188E+04
1183      + ,  0.1551E+03, -0.9040E+03,  0.1404E+06, -0.2314E+06
1184      + , -0.6647E+02,  0.6324E+04,  0.2871E+06, -0.2439E+06
1185      + , -0.1405E+03,  0.2137E+05,  0.2340E+06, -0.4290E+06
1186      + , -0.1063E+02,  0.4195E+04,  0.7543E+06, -0.3747E+07
1187      + ,  0.7684E+01,  0.4048E+04,  0.6314E+07, -0.1714E+09
1188      + ,  0.4082E+01,  0.8543E+04,  0.4521E+07,  0.5178E+08,  8*0.0 /
1189       DATA ((TMP(I,J, 82),I=1,MAXPOW),J=1,MAXINT) /
1190      +   0.1569E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1191      + , -0.9651E+04,  0.1870E+04, -0.7761E+02,  0.9425E+00
1192      + ,  0.1609E+05, -0.5855E+04,  0.7149E+03, -0.2672E+02
1193      + , -0.2194E+04,  0.1239E+05, -0.5418E+04,  0.6794E+03
1194      + , -0.4016E+03,  0.1173E+05,  0.2316E+05, -0.7955E+04
1195      + ,  0.6749E+02,  0.6789E+03,  0.1401E+06, -0.2437E+06
1196      + , -0.7125E+02,  0.7126E+04,  0.2931E+06, -0.2779E+06
1197      + , -0.1363E+03,  0.1992E+05,  0.2721E+06, -0.4887E+06
1198      + , -0.8206E+01,  0.4058E+04,  0.7950E+06, -0.4115E+07
1199      + ,  0.7595E+01,  0.4692E+04,  0.6429E+07, -0.1747E+09
1200      + ,  0.4240E+01,  0.8883E+04,  0.4754E+07,  0.3448E+08,  8*0.0 /
1201       DATA ((TMP(I,J, 83),I=1,MAXPOW),J=1,MAXINT) /
1202      +  -0.4395E+03,  0.1185E+02,  0.2207E+00, -0.1865E-02
1203      + ,  0.4586E+04, -0.1052E+03, -0.2706E+01,  0.7877E-01
1204      + ,  0.1118E+05, -0.3159E+03, -0.9250E+01,  0.2918E+00
1205      + , -0.1578E+05,  0.3450E+04, -0.1741E+03,  0.2596E+01
1206      + ,  0.2124E+05, -0.8884E+04,  0.1240E+04, -0.5312E+02
1207      + , -0.1835E+04,  0.1188E+05, -0.5123E+04,  0.6229E+03
1208      + , -0.3943E+03,  0.1129E+05,  0.2671E+05, -0.9164E+04
1209      + ,  0.2402E+03, -0.2657E+04,  0.1637E+06, -0.2867E+06
1210      + , -0.6824E+02,  0.7005E+04,  0.3089E+06, -0.2951E+06
1211      + , -0.1337E+03,  0.1999E+05,  0.2903E+06, -0.5342E+06
1212      + , -0.1411E+01,  0.3094E+04,  0.8901E+06, -0.5203E+07
1213      + ,  0.9935E+01,  0.3289E+04,  0.6983E+07, -0.2020E+09
1214      + ,  0.4425E+01,  0.9396E+04,  0.5008E+07, -0.5280E+08
1215      +/
1216       DATA ((TMP(I,J, 84),I=1,MAXPOW),J=1,MAXINT) /
1217      +   0.1109E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1218      + ,  0.1939E+05, -0.1380E+04,  0.3526E+02, -0.3106E+00
1219      + ,  0.9152E+04, -0.3829E+04,  0.5535E+03, -0.2106E+02
1220      + , -0.1131E+04,  0.1041E+05, -0.4126E+04,  0.4235E+03
1221      + , -0.3899E+03,  0.1134E+05,  0.2963E+05, -0.1038E+05
1222      + ,  0.4918E+03, -0.1053E+05,  0.2383E+06, -0.4647E+06
1223      + , -0.7471E+02,  0.8106E+04,  0.3143E+06, -0.3397E+06
1224      + , -0.1335E+03,  0.2199E+05,  0.2786E+06, -0.5569E+06
1225      + ,  0.7789E+01,  0.1767E+04,  0.1003E+07, -0.6443E+07
1226      + ,  0.8525E+01,  0.5102E+04,  0.6990E+07, -0.2009E+09
1227      + ,  0.4633E+01,  0.9868E+04,  0.5142E+07,  0.1775E+08,  8*0.0 /
1228       DATA ((TMP(I,J, 85),I=1,MAXPOW),J=1,MAXINT) /
1229      +   0.5224E+03,  0.0000E+00,  0.0000E+00,  0.0000E+00
1230      + ,  0.2744E+05, -0.2310E+04,  0.6658E+02, -0.6458E+00
1231      + ,  0.6422E+04, -0.2954E+04,  0.4933E+03, -0.2019E+02
1232      + , -0.1933E+04,  0.1303E+05, -0.5806E+04,  0.6870E+03
1233      + , -0.3889E+03,  0.1171E+05,  0.3195E+05, -0.1160E+05
1234      + ,  0.3129E+03, -0.6378E+04,  0.2158E+06, -0.4073E+06
1235      + , -0.6880E+02,  0.7565E+04,  0.3411E+06, -0.3547E+06
1236      + , -0.1280E+03,  0.1941E+05,  0.3508E+06, -0.6558E+06
1237      + , -0.1695E+01,  0.3621E+04,  0.9665E+06, -0.5746E+07
1238      + ,  0.9913E+01,  0.4822E+04,  0.7370E+07, -0.2163E+09
1239      + ,  0.4871E+01,  0.1062E+05,  0.5352E+07, -0.2702E+08,  8*0.0 /
1240       DATA ((TMP(I,J, 86),I=1,MAXPOW),J=1,MAXINT) /
1241      +   0.1284E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1242      + , -0.7912E+03,  0.2111E+04, -0.1709E+03,  0.3589E+01
1243      + ,  0.1243E+05, -0.6150E+04,  0.1021E+04, -0.4644E+02
1244      + , -0.1804E+04,  0.1279E+05, -0.5624E+04,  0.6228E+03
1245      + , -0.3616E+03,  0.1088E+05,  0.3423E+05, -0.1248E+05
1246      + ,  0.1691E+03, -0.2424E+04,  0.1954E+06, -0.4008E+06
1247      + , -0.6756E+02,  0.7739E+04,  0.3363E+06, -0.3748E+06
1248      + , -0.1220E+03,  0.1976E+05,  0.3335E+06, -0.6657E+06
1249      + , -0.5899E+00,  0.3312E+04,  0.9833E+06, -0.6258E+07
1250      + ,  0.8633E+01,  0.5888E+04,  0.7122E+07, -0.2120E+09
1251      + ,  0.4841E+01,  0.1057E+05,  0.5286E+07,  0.1022E+08,  8*0.0 /
1252       DATA ((TMP(I,J, 87),I=1,MAXPOW),J=1,MAXINT) /
1253      +   0.1287E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1254      + , -0.9129E+04,  0.3064E+04, -0.1840E+03,  0.2869E+01
1255      + ,  0.1217E+05, -0.1865E+04, -0.1031E+04,  0.1980E+03
1256      + , -0.2201E+04,  0.1461E+05, -0.7082E+04,  0.9132E+03
1257      + , -0.3615E+03,  0.1150E+05,  0.3541E+05, -0.1363E+05
1258      + , -0.5769E+03,  0.1660E+05,  0.5526E+05, -0.7044E+05
1259      + , -0.6236E+02,  0.7249E+04,  0.3612E+06, -0.3893E+06
1260      + , -0.1213E+03,  0.2052E+05,  0.3455E+06, -0.7195E+06
1261      + , -0.8453E+01,  0.4660E+04,  0.9757E+06, -0.6010E+07
1262      + ,  0.1104E+02,  0.4630E+04,  0.7706E+07, -0.2445E+09
1263      + ,  0.5063E+01,  0.1115E+05,  0.5666E+07, -0.1082E+09,  8*0.0 /
1264       DATA ((TMP(I,J, 88),I=1,MAXPOW),J=1,MAXINT) /
1265      +   0.1415E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1266      + , -0.3003E+05,  0.8875E+04, -0.6798E+03,  0.1571E+02
1267      + ,  0.2941E+05, -0.1786E+05,  0.3591E+04, -0.2242E+03
1268      + , -0.1852E+04,  0.1361E+05, -0.6123E+04,  0.6294E+03
1269      + , -0.3547E+03,  0.1164E+05,  0.3753E+05, -0.1492E+05
1270      + ,  0.2933E+03, -0.6829E+04,  0.2528E+06, -0.5565E+06
1271      + , -0.6085E+02,  0.7287E+04,  0.3762E+06, -0.4161E+06
1272      + , -0.1192E+03,  0.2119E+05,  0.3500E+06, -0.7644E+06
1273      + , -0.1024E+02,  0.5310E+04,  0.9756E+06, -0.5773E+07
1274      + ,  0.1241E+02,  0.3727E+04,  0.8307E+07, -0.2833E+09
1275      + ,  0.5211E+01,  0.1168E+05,  0.5739E+07, -0.9046E+08,  8*0.0 /
1276       DATA ((TMP(I,J, 89),I=1,MAXPOW),J=1,MAXINT) /
1277      +   0.1258E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1278      + , -0.2747E+04,  0.3975E+03,  0.1941E+03, -0.1439E+02
1279      + ,  0.1958E+05, -0.9878E+04,  0.1388E+04, -0.1778E+02
1280      + , -0.2098E+04,  0.1506E+05, -0.7581E+04,  0.9825E+03
1281      + , -0.3474E+03,  0.1153E+05,  0.4113E+05, -0.1657E+05
1282      + ,  0.5789E+03, -0.1507E+05,  0.3328E+06, -0.7718E+06
1283      + , -0.6339E+02,  0.7915E+04,  0.3884E+06, -0.4637E+06
1284      + , -0.1182E+03,  0.2194E+05,  0.3621E+06, -0.8244E+06
1285      + ,  0.8142E+01,  0.2287E+04,  0.1185E+07, -0.8439E+07
1286      + ,  0.1003E+02,  0.6690E+04,  0.7961E+07, -0.2545E+09
1287      + ,  0.5479E+01,  0.1217E+05,  0.5897E+07, -0.2287E+08,  8*0.0 /
1288       DATA ((TMP(I,J, 90),I=1,MAXPOW),J=1,MAXINT) /
1289      +   0.1134E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1290      + ,  0.1398E+05, -0.6113E+04,  0.8572E+03, -0.2963E+02
1291      + , -0.1466E+04,  0.1349E+05, -0.6743E+04,  0.1013E+04
1292      + , -0.3308E+03,  0.1094E+05,  0.4517E+05, -0.1812E+05
1293      + ,  0.2265E+03, -0.3796E+04,  0.2424E+06, -0.5526E+06
1294      + , -0.5998E+02,  0.7669E+04,  0.4039E+06, -0.4843E+06
1295      + , -0.1152E+03,  0.2229E+05,  0.3694E+06, -0.8743E+06
1296      + , -0.1965E+02,  0.7238E+04,  0.9916E+06, -0.5869E+07
1297      + ,  0.1150E+02,  0.6085E+04,  0.8292E+07, -0.2758E+09
1298      + ,  0.5620E+01,  0.1245E+05,  0.6338E+07, -0.1549E+09,  12*0.0 /
1299       DATA ((TMP(I,J, 91),I=1,MAXPOW),J=1,MAXINT) /
1300      +   0.2194E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1301      + ,  0.8794E+04, -0.2841E+04,  0.2038E+03,  0.1370E+02
1302      + , -0.2042E+04,  0.1581E+05, -0.8451E+04,  0.1211E+04
1303      + , -0.3240E+03,  0.1078E+05,  0.5022E+05, -0.2015E+05
1304      + ,  0.1921E+04, -0.5267E+05,  0.6573E+06, -0.1562E+07
1305      + , -0.6436E+02,  0.8617E+04,  0.4180E+06, -0.5495E+06
1306      + , -0.1129E+03,  0.2138E+05,  0.4231E+06, -0.9875E+06
1307      + ,  0.1526E+02,  0.1356E+04,  0.1351E+07, -0.1082E+08
1308      + ,  0.1254E+02,  0.5894E+04,  0.8877E+07, -0.3120E+09
1309      + ,  0.5918E+01,  0.1306E+05,  0.6683E+07, -0.1778E+09,  12*0.0 /
1310       DATA ((TMP(I,J, 92),I=1,MAXPOW),J=1,MAXINT) /
1311      +   0.4400E+04, -0.1564E+03,  0.2241E+01, -0.1062E-01
1312      + , -0.2175E+04,  0.4327E+03, -0.1164E+02,  0.8742E-01
1313      + ,  0.1528E+05, -0.7106E+04,  0.1060E+04, -0.4036E+02
1314      + , -0.1807E+04,  0.1483E+05, -0.7405E+04,  0.8426E+03
1315      + , -0.3148E+03,  0.1094E+05,  0.5041E+05, -0.2121E+05
1316      + , -0.1100E+03,  0.4889E+04,  0.1908E+06, -0.4810E+06
1317      + , -0.6720E+02,  0.9489E+04,  0.4107E+06, -0.6000E+06
1318      + , -0.1033E+03,  0.1864E+05,  0.4878E+06, -0.1076E+07
1319      + ,  0.6281E+00,  0.3792E+04,  0.1281E+07, -0.1003E+08
1320      + ,  0.1146E+02,  0.7568E+04,  0.8699E+07, -0.3106E+09
1321      + ,  0.6017E+01,  0.1327E+05,  0.7079E+07, -0.2456E+09,  8*0.0 /
1322       DATA ((TMP(I,J, 93),I=1,MAXPOW),J=1,MAXINT) /
1323      +   0.2522E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1324      + ,  0.1668E+05, -0.8606E+04,  0.1441E+04, -0.6401E+02
1325      + , -0.2561E+04,  0.1858E+05, -0.1101E+05,  0.1791E+04
1326      + , -0.3232E+03,  0.1200E+05,  0.6096E+05, -0.2302E+05
1327      + , -0.4467E+02,  0.2778E+04,  0.2204E+06, -0.5257E+06
1328      + , -0.6355E+02,  0.9093E+04,  0.4439E+06, -0.6341E+06
1329      + , -0.1103E+03,  0.2501E+05,  0.3841E+06, -0.1044E+07
1330      + ,  0.5579E+01,  0.3489E+04,  0.1354E+07, -0.1078E+08
1331      + ,  0.1272E+02,  0.7422E+04,  0.9191E+07, -0.3336E+09
1332      + ,  0.6330E+01,  0.1392E+05,  0.7484E+07, -0.3064E+09,  12*0.0 /
1333       DATA ((TMP(I,J, 94),I=1,MAXPOW),J=1,MAXINT) /
1334      +   0.3508E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1335      + , -0.2939E+04,  0.8224E+03, -0.3576E+02,  0.5087E+00
1336      + ,  0.1270E+05, -0.6607E+04,  0.1093E+04, -0.3674E+02
1337      + , -0.2063E+04,  0.1760E+05, -0.1057E+05,  0.1695E+04
1338      + , -0.3181E+03,  0.1245E+05,  0.5096E+05, -0.2427E+05
1339      + ,  0.3925E+03, -0.1022E+05,  0.3515E+06, -0.9309E+06
1340      + , -0.6367E+02,  0.9444E+04,  0.4539E+06, -0.6824E+06
1341      + , -0.1065E+03,  0.2299E+05,  0.4461E+06, -0.1181E+07
1342      + , -0.1248E+02,  0.6703E+04,  0.1250E+07, -0.9411E+07
1343      + ,  0.1781E+02,  0.3158E+04,  0.1075E+08, -0.4693E+09
1344      + ,  0.6490E+01,  0.1422E+05,  0.8179E+07, -0.5340E+09,  8*0.0 /
1345       DATA ((TMP(I,J, 95),I=1,MAXPOW),J=1,MAXINT) /
1346      +   0.9600E+03,  0.0000E+00,  0.0000E+00,  0.0000E+00
1347      + ,  0.3669E+04, -0.1927E+04,  0.3556E+03, -0.1775E+02
1348      + , -0.8806E+05,  0.1824E+06, -0.1119E+06,  0.2150E+05
1349      + , -0.8786E+03,  0.1063E+05, -0.5553E+03, -0.1750E+04
1350      + , -0.3188E+03,  0.1325E+05,  0.5127E+05, -0.2586E+05
1351      + ,  0.3232E+03, -0.9027E+04,  0.3584E+06, -0.9931E+06
1352      + , -0.6201E+02,  0.9427E+04,  0.4778E+06, -0.7281E+06
1353      + , -0.9199E+02,  0.1716E+05,  0.6101E+06, -0.1350E+07
1354      + , -0.3730E+01,  0.5425E+04,  0.1376E+07, -0.1109E+08
1355      + ,  0.1966E+02,  0.2217E+04,  0.1154E+08, -0.5298E+09
1356      + ,  0.6763E+01,  0.1479E+05,  0.8647E+07, -0.6116E+09,  8*0.0 /
1357       DATA ((TMP(I,J, 96),I=1,MAXPOW),J=1,MAXINT) /
1358      +   0.1000E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1359      + ,  0.3148E+04, -0.1605E+04,  0.2954E+03, -0.1411E+02
1360      + , -0.7979E+05,  0.1716E+06, -0.1079E+06,  0.2118E+05
1361      + , -0.8139E+03,  0.1040E+05,  0.4547E+03, -0.2530E+04
1362      + , -0.3144E+03,  0.1378E+05,  0.5160E+05, -0.2733E+05
1363      + , -0.2912E+03,  0.1117E+05,  0.1679E+06, -0.4469E+06
1364      + , -0.5879E+02,  0.9126E+04,  0.4981E+06, -0.7597E+06
1365      + , -0.1010E+03,  0.2212E+05,  0.5193E+06, -0.1402E+07
1366      + , -0.1147E+02,  0.7190E+04,  0.1327E+07, -0.1016E+08
1367      + ,  0.2200E+02,  0.1083E+04,  0.1208E+08, -0.5684E+09
1368      + ,  0.7017E+01,  0.1439E+05,  0.9471E+07, -0.8527E+09,  8*0.0 /
1369       DATA ((TMP(I,J, 97),I=1,MAXPOW),J=1,MAXINT) /
1370      +   0.1100E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1371      + ,  0.3725E+04, -0.1998E+04,  0.3728E+03, -0.1783E+02
1372      + , -0.1348E+06,  0.2867E+06, -0.1865E+06,  0.3853E+05
1373      + , -0.1092E+04,  0.1257E+05, -0.3090E+04, -0.6861E+03
1374      + , -0.3108E+03,  0.1396E+05,  0.5425E+05, -0.2964E+05
1375      + , -0.8536E+02,  0.5200E+04,  0.2278E+06, -0.6067E+06
1376      + , -0.6035E+02,  0.9752E+04,  0.5095E+06, -0.8199E+06
1377      + , -0.1018E+03,  0.2538E+05,  0.4681E+06, -0.1141E+07
1378      + , -0.1308E+02,  0.7916E+04,  0.1354E+07, -0.1049E+08
1379      + ,  0.2374E+02,  0.1750E+03,  0.1285E+08, -0.6352E+09
1380      + ,  0.7219E+01,  0.1592E+05,  0.9491E+07, -0.8283E+09,  8*0.0 /
1381       DATA ((TMP(I,J, 98),I=1,MAXPOW),J=1,MAXINT) /
1382      +   0.1220E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1383      + ,  0.3846E+04, -0.2092E+04,  0.3854E+03, -0.1758E+02
1384      + , -0.5378E+05,  0.1380E+06, -0.9898E+05,  0.2185E+05
1385      + , -0.9355E+03,  0.1168E+05, -0.4663E+03, -0.2379E+04
1386      + , -0.3044E+03,  0.1361E+05,  0.6013E+05, -0.3273E+05
1387      + ,  0.3450E+03, -0.1057E+05,  0.4170E+06, -0.1248E+07
1388      + , -0.5797E+02,  0.9571E+04,  0.5352E+06, -0.8755E+06
1389      + , -0.9148E+02,  0.1942E+05,  0.6490E+06, -0.1673E+07
1390      + , -0.5416E+02,  0.1706E+05,  0.8630E+06, -0.1792E+07
1391      + ,  0.2771E+02, -0.2903E+04,  0.1419E+08, -0.7456E+09
1392      + ,  0.7544E+01,  0.1589E+05,  0.1052E+08, -0.1089E+10,  8*0.0 /
1393       DATA ((TMP(I,J, 99),I=1,MAXPOW),J=1,MAXINT) /
1394      +   0.1330E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1395      + , -0.8936E+03,  0.1145E+05,  0.1094E+04, -0.3692E+04
1396      + , -0.2152E+03,  0.7720E+04,  0.5372E+05, -0.2276E+05
1397      + , -0.2914E+03,  0.1268E+05,  0.6828E+05, -0.3620E+05
1398      + , -0.2568E+03,  0.1341E+05,  0.8685E+05, -0.5387E+05
1399      + , -0.2104E+03,  0.1155E+05,  0.1107E+06, -0.7279E+05
1400      + ,  0.1088E+03,  0.1536E+04,  0.2458E+06, -0.6005E+06
1401      + , -0.5691E+02,  0.9660E+04,  0.5556E+06, -0.9320E+06
1402      + , -0.9167E+02,  0.2050E+05,  0.6520E+06, -0.1784E+07
1403      + ,  0.9687E+00,  0.5331E+04,  0.1646E+07, -0.1518E+08
1404      + ,  0.3315E+02, -0.7298E+04,  0.1587E+08, -0.8866E+09
1405      + ,  0.7728E+01,  0.1675E+05,  0.1090E+08, -0.1231E+10,  4*0.0 /
1406       DATA ((TMP(I,J,100),I=1,MAXPOW),J=1,MAXINT) /
1407      +   0.1460E+04,  0.0000E+00,  0.0000E+00,  0.0000E+00
1408      + , -0.1082E+04,  0.1308E+05, -0.9100E+03, -0.2846E+04
1409      + , -0.1923E+03,  0.6709E+04,  0.6231E+05, -0.2410E+05
1410      + , -0.2559E+03,  0.1024E+05,  0.8602E+05, -0.3977E+05
1411      + , -0.2466E+03,  0.1276E+05,  0.9751E+05, -0.5925E+05
1412      + , -0.2304E+03,  0.1518E+05,  0.9220E+05, -0.7454E+05
1413      + , -0.8675E+02,  0.2132E+04,  0.2568E+06, -0.6538E+06
1414      + , -0.5579E+02,  0.1018E+05,  0.5793E+06, -0.1018E+07
1415      + , -0.8637E+02,  0.1890E+05,  0.7445E+06, -0.1949E+07
1416      + ,  0.5556E+00,  0.5648E+04,  0.1738E+07, -0.1662E+08
1417      + ,  0.2968E+02, -0.2026E+04,  0.1484E+08, -0.7825E+09
1418      + ,  0.8135E+01,  0.1726E+05,  0.1171E+08, -0.1346E+10,  4*0.0 /
1419 C Ionization energies
1420       DATA (TMP1(I), I = 1,56) /
1421      +     13.60, 24.59,
1422      +      5.39,  9.32,  8.30, 11.26, 14.53, 13.62, 17.42, 21.56,
1423      +      5.14,  7.65,  5.99,  8.15, 10.49, 10.36, 12.97, 15.76,
1424      +      4.34,  6.11,
1425      +      6.56,  6.83,  6.75,  6.77,  7.43,  7.90,  7.88,  7.64,
1426      +      7.73,  9.39,
1427      +                    6.00,  7.90,  9.82,  9.75, 11.81, 14.00,
1428      +      4.18,  5.69,
1429      +      6.22,  6.63,  6.76,  7.09,  7.28,  7.36,  7.46,  8.34,
1430      +      7.58,  8.99,  5.79,  7.34,  8.64,  9.01, 10.45, 12.13,
1431      +      3.89,  5.21 /
1432       DATA (TMP1(I), I = 57,MAXELZ) /
1433      +      5.58,  5.54,  5.46,  5.52,  5.55,  5.64,  5.67,  6.15,
1434      +      5.86,  5.94,  6.02,  6.11,  6.18,  6.25,
1435      +      5.43,  6.83,  7.89,  7.98,  7.88,  8.70,  9.10,  9.00,
1436      +      9.23, 10.44,
1437      +                    6.11,  7.42,  7.29,  8.42,  9.65, 10.75,
1438      +      3.97,  5.28,
1439      +      5.17,  6.08,  5.89,  6.19,  6.27,  6.06,  5.99,  6.02,
1440      +      6.23,  6.30,  6.42,  6.50 /
1441       DO 10 J = 1,MAXELZ
1442          DO 10 I = 1,MAXINT
1443             CRNGUP(I,J) = CTMP(I,J)
1444    10 CONTINUE
1445       DO 20 K = 1,MAXELZ
1446          DO 20 J = 1,MAXINT
1447             DO 20 I = 1,MAXPOW
1448                COFS(I,J,K) = TMP(I,J,K)
1449    20 CONTINUE
1450       DO 30 I = 1,MAXELZ
1451          GPOMIN(I) = MAX(TENEV,TMP1(I)*1.E-3)
1452    30 CONTINUE
1453       END