This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / GEANT321 / gphys / gphrin.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1995/10/24 10:21:29  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 GPHRIN
13 C.
14 C.    ******************************************************************
15 C.    *                                                                *
16 C     *       Initialize decay mode rates' constants for               *
17 C.    *       the photoelectric effect                                 *
18 C.    *                                                                *
19 C.    *                                                                *
20 C.    *    ==>Called by : GPHINI                                       *
21 C.    *       Author    J. Chwastowski                                 *
22 C.    *                                                                *
23 C.    ******************************************************************
24 C.
25       IMPLICIT NONE
26 #include "geant321/gcpmxz.inc"
27 #include "geant321/gcphpr.inc"
28 #include "geant321/gcphnr.inc"
29 #include "geant321/gcphrd.inc"
30       INTEGER I,J,ITRN
31       REAL TMP,TMP1
32       DIMENSION TMP(ISHLS,MAXELZ),TMP1(4,MAXELZ)
33       DIMENSION ITRN(29)
34       DATA ITRN / 3,4,6, 7,11,12,
35      +            3,4,6, 7,11,12,18,19,
36      +            4,5,8,10,13,17,20,
37      +            5,8,9,10,13,14,17,20 /
38       DATA (TMP(J,1),J=1,ISHLS) / ISHLS*0.0 /
39       DATA (TMP(J,2),J=1,ISHLS) / ISHLS*0.0 /
40       DATA (TMP(J,3),J=1,ISHLS) / ISHLS*0.0 /
41       DATA (TMP(J,4),J=1,ISHLS) / ISHLS*0.0 /
42       DATA (TMP(J,  5),J=1,ISHLS) /
43      +  0.333E+00, 0.667E+00, 0.000E+00, 0.000E+00
44      +, 0.000E+00, 0.000E+00, 0.333E+00, 0.667E+00
45      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
46      +, 0.000E+00, 0.000E+00, 0.100E+01, 14*0.0 /
47       DATA (TMP(J,  6),J=1,ISHLS) /
48      +  0.334E+00, 0.666E+00, 0.000E+00, 0.000E+00
49      +, 0.000E+00, 0.000E+00, 0.332E+00, 0.668E+00
50      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
51      +, 0.000E+00, 0.000E+00, 0.100E+01, 14*0.0 /
52       DATA (TMP(J,  7),J=1,ISHLS) /
53      +  0.334E+00, 0.666E+00, 0.000E+00, 0.000E+00
54      +, 0.000E+00, 0.000E+00, 0.332E+00, 0.668E+00
55      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
56      +, 0.000E+00, 0.000E+00, 0.100E+01, 14*0.0 /
57       DATA (TMP(J,  8),J=1,ISHLS) /
58      +  0.334E+00, 0.666E+00, 0.000E+00, 0.000E+00
59      +, 0.000E+00, 0.000E+00, 0.331E+00, 0.669E+00
60      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
61      +, 0.000E+00, 0.000E+00, 0.100E+01, 14*0.0 /
62       DATA (TMP(J,  9),J=1,ISHLS) /
63      +  0.334E+00, 0.666E+00, 0.000E+00, 0.000E+00
64      +, 0.000E+00, 0.000E+00, 0.331E+00, 0.669E+00
65      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
66      +, 0.000E+00, 0.000E+00, 0.100E+01, 14*0.0 /
67       DATA (TMP(J, 10),J=1,ISHLS) /
68      +  0.335E+00, 0.665E+00, 0.000E+00, 0.000E+00
69      +, 0.000E+00, 0.000E+00, 0.329E+00, 0.671E+00
70      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
71      +, 0.000E+00, 0.000E+00, 0.100E+01, 14*0.0 /
72       DATA (TMP(J, 11),J=1,ISHLS) /
73      +  0.335E+00, 0.665E+00, 0.000E+00, 0.000E+00
74      +, 0.000E+00, 0.000E+00, 0.328E+00, 0.672E+00
75      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
76      +, 0.000E+00, 0.000E+00, 0.233E-09, 0.100E+01
77      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
78      +, 0.000E+00, 0.100E+01,  7*0.0 /
79       DATA (TMP(J, 12),J=1,ISHLS) /
80      +  0.335E+00, 0.665E+00, 0.000E+00, 0.000E+00
81      +, 0.000E+00, 0.000E+00, 0.326E+00, 0.674E+00
82      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
83      +, 0.000E+00, 0.000E+00, 0.159E-09, 0.100E+01
84      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
85      +, 0.000E+00, 0.100E+01,  7*0.0 /
86       DATA (TMP(J, 13),J=1,ISHLS) /
87      +  0.333E+00, 0.661E+00, 0.203E-02, 0.404E-02
88      +, 0.000E+00, 0.000E+00, 0.291E+00, 0.608E+00
89      +, 0.342E-01, 0.663E-01, 0.000E+00, 0.000E+00
90      +, 0.000E+00, 0.000E+00, 0.239E-09, 0.100E+01
91      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
92      +, 0.000E+00, 0.100E+01,  7*0.0 /
93       DATA (TMP(J, 14),J=1,ISHLS) /
94      +  0.330E+00, 0.654E+00, 0.562E-02, 0.111E-01
95      +, 0.000E+00, 0.000E+00, 0.216E+00, 0.456E+00
96      +, 0.111E+00, 0.216E+00, 0.000E+00, 0.000E+00
97      +, 0.000E+00, 0.000E+00, 0.408E-09, 0.100E+01
98      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
99      +, 0.000E+00, 0.100E+01,  7*0.0 /
100       DATA (TMP(J, 15),J=1,ISHLS) /
101      +  0.325E+00, 0.645E+00, 0.100E-01, 0.199E-01
102      +, 0.000E+00, 0.000E+00, 0.141E+00, 0.300E+00
103      +, 0.190E+00, 0.369E+00, 0.000E+00, 0.000E+00
104      +, 0.000E+00, 0.000E+00, 0.703E-09, 0.100E+01
105      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
106      +, 0.000E+00, 0.100E+01,  7*0.0 /
107       DATA (TMP(J, 16),J=1,ISHLS) /
108      +  0.321E+00, 0.635E+00, 0.149E-01, 0.296E-01
109      +, 0.000E+00, 0.000E+00, 0.881E-01, 0.191E+00
110      +, 0.246E+00, 0.475E+00, 0.000E+00, 0.000E+00
111      +, 0.000E+00, 0.000E+00, 0.120E-08, 0.100E+01
112      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
113      +, 0.000E+00, 0.100E+01,  7*0.0 /
114       DATA (TMP(J, 17),J=1,ISHLS) /
115      +  0.315E+00, 0.625E+00, 0.201E-01, 0.398E-01
116      +, 0.000E+00, 0.000E+00, 0.559E-01, 0.124E+00
117      +, 0.281E+00, 0.540E+00, 0.000E+00, 0.000E+00
118      +, 0.000E+00, 0.000E+00, 0.199E-08, 0.100E+01
119      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
120      +, 0.000E+00, 0.100E+01,  7*0.0 /
121       DATA (TMP(J, 18),J=1,ISHLS) /
122      +  0.314E+00, 0.611E+00, 0.253E-01, 0.500E-01
123      +, 0.000E+00, 0.000E+00, 0.368E-01, 0.828E-01
124      +, 0.302E+00, 0.579E+00, 0.000E+00, 0.000E+00
125      +, 0.000E+00, 0.000E+00, 0.324E-08, 0.100E+01
126      +, 0.000E+00, 0.000E+00, 0.000E+00, 0.000E+00
127      +, 0.000E+00, 0.100E+01,  7*0.0 /
128       DATA (TMP(J, 19),J=1,ISHLS) /
129      +  0.307E+00, 0.606E+00, 0.293E-01, 0.580E-01
130      +, 0.000E+00, 0.000E+00, 0.260E-01, 0.597E-01
131      +, 0.313E+00, 0.601E+00, 0.000E+00, 0.000E+00
132      +, 0.000E+00, 0.000E+00, 0.480E-08, 0.976E+00
133      +, 0.000E+00, 0.238E-01, 0.000E+00, 0.000E+00
134      +, 0.000E+00, 0.976E+00, 0.000E+00, 0.000E+00
135      +, 0.238E-01,  4*0.0 /
136       DATA (TMP(J, 20),J=1,ISHLS) /
137      +  0.304E+00, 0.600E+00, 0.323E-01, 0.639E-01
138      +, 0.000E+00, 0.000E+00, 0.193E-01, 0.455E-01
139      +, 0.321E+00, 0.614E+00, 0.000E+00, 0.000E+00
140      +, 0.000E+00, 0.000E+00, 0.678E-08, 0.929E+00
141      +, 0.000E+00, 0.709E-01, 0.000E+00, 0.000E+00
142      +, 0.000E+00, 0.929E+00, 0.000E+00, 0.000E+00
143      +, 0.709E-01,  4*0.0 /
144       DATA (TMP(J, 21),J=1,ISHLS) /
145      +  0.303E+00, 0.597E+00, 0.335E-01, 0.661E-01
146      +, 0.000E+00, 0.000E+00, 0.155E-01, 0.376E-01
147      +, 0.326E+00, 0.621E+00, 0.000E+00, 0.000E+00
148      +, 0.000E+00, 0.000E+00, 0.504E-08, 0.452E+00
149      +, 0.513E+00, 0.347E-01, 0.000E+00, 0.000E+00
150      +, 0.000E+00, 0.461E+00, 0.506E-01, 0.453E+00
151      +, 0.353E-01,  4*0.0 /
152       DATA (TMP(J, 22),J=1,ISHLS) /
153      +  0.303E+00, 0.595E+00, 0.343E-01, 0.677E-01
154      +, 0.000E+00, 0.000E+00, 0.128E-01, 0.317E-01
155      +, 0.330E+00, 0.626E+00, 0.000E+00, 0.000E+00
156      +, 0.000E+00, 0.000E+00, 0.463E-08, 0.276E+00
157      +, 0.703E+00, 0.208E-01, 0.000E+00, 0.000E+00
158      +, 0.000E+00, 0.284E+00, 0.699E-01, 0.624E+00
159      +, 0.214E-01,  4*0.0 /
160       DATA (TMP(J, 23),J=1,ISHLS) /
161      +  0.302E+00, 0.594E+00, 0.350E-01, 0.690E-01
162      +, 0.000E+00, 0.000E+00, 0.107E-01, 0.273E-01
163      +, 0.333E+00, 0.629E+00, 0.000E+00, 0.000E+00
164      +, 0.000E+00, 0.000E+00, 0.471E-08, 0.191E+00
165      +, 0.795E+00, 0.139E-01, 0.000E+00, 0.000E+00
166      +, 0.000E+00, 0.198E+00, 0.793E-01, 0.708E+00
167      +, 0.145E-01,  4*0.0 /
168       DATA (TMP(J, 24),J=1,ISHLS) /
169      +  0.303E+00, 0.594E+00, 0.348E-01, 0.685E-01
170      +, 0.000E+00, 0.000E+00, 0.926E-02, 0.245E-01
171      +, 0.335E+00, 0.631E+00, 0.000E+00, 0.000E+00
172      +, 0.000E+00, 0.000E+00, 0.475E-08, 0.130E+00
173      +, 0.861E+00, 0.888E-02, 0.000E+00, 0.000E+00
174      +, 0.000E+00, 0.136E+00, 0.863E-01, 0.769E+00
175      +, 0.910E-02,  4*0.0 /
176       DATA (TMP(J, 25),J=1,ISHLS) /
177      +  0.302E+00, 0.592E+00, 0.359E-01, 0.708E-01
178      +, 0.000E+00, 0.000E+00, 0.778E-02, 0.213E-01
179      +, 0.337E+00, 0.633E+00, 0.000E+00, 0.000E+00
180      +, 0.000E+00, 0.000E+00, 0.573E-08, 0.111E+00
181      +, 0.881E+00, 0.759E-02, 0.000E+00, 0.000E+00
182      +, 0.000E+00, 0.117E+00, 0.882E-01, 0.787E+00
183      +, 0.803E-02,  4*0.0 /
184       DATA (TMP(J, 26),J=1,ISHLS) /
185      +  0.302E+00, 0.591E+00, 0.363E-01, 0.714E-01
186      +, 0.000E+00, 0.000E+00, 0.675E-02, 0.192E-01
187      +, 0.339E+00, 0.635E+00, 0.000E+00, 0.000E+00
188      +, 0.000E+00, 0.000E+00, 0.657E-08, 0.900E-01
189      +, 0.904E+00, 0.597E-02, 0.000E+00, 0.000E+00
190      +, 0.000E+00, 0.955E-01, 0.906E-01, 0.807E+00
191      +, 0.634E-02,  4*0.0 /
192       DATA (TMP(J, 27),J=1,ISHLS) /
193      +  0.302E+00, 0.590E+00, 0.366E-01, 0.719E-01
194      +, 0.000E+00, 0.000E+00, 0.589E-02, 0.176E-01
195      +, 0.341E+00, 0.635E+00, 0.000E+00, 0.000E+00
196      +, 0.000E+00, 0.000E+00, 0.767E-08, 0.751E-01
197      +, 0.920E+00, 0.483E-02, 0.000E+00, 0.000E+00
198      +, 0.000E+00, 0.803E-01, 0.924E-01, 0.822E+00
199      +, 0.515E-02,  4*0.0 /
200       DATA (TMP(J, 28),J=1,ISHLS) /
201      +  0.301E+00, 0.589E+00, 0.369E-01, 0.724E-01
202      +, 0.000E+00, 0.000E+00, 0.521E-02, 0.162E-01
203      +, 0.343E+00, 0.635E+00, 0.000E+00, 0.000E+00
204      +, 0.000E+00, 0.000E+00, 0.902E-08, 0.640E-01
205      +, 0.932E+00, 0.398E-02, 0.000E+00, 0.000E+00
206      +, 0.000E+00, 0.690E-01, 0.938E-01, 0.833E+00
207      +, 0.428E-02,  4*0.0 /
208       DATA (TMP(J, 29),J=1,ISHLS) /
209      +  0.302E+00, 0.590E+00, 0.366E-01, 0.717E-01
210      +, 0.000E+00, 0.000E+00, 0.470E-02, 0.153E-01
211      +, 0.345E+00, 0.635E+00, 0.000E+00, 0.000E+00
212      +, 0.000E+00, 0.000E+00, 0.104E-07, 0.537E-01
213      +, 0.943E+00, 0.325E-02, 0.000E+00, 0.000E+00
214      +, 0.000E+00, 0.585E-01, 0.954E-01, 0.845E+00
215      +, 0.134E-02,  4*0.0 /
216       DATA (TMP(J, 30),J=1,ISHLS) /
217      +  0.302E+00, 0.588E+00, 0.373E-01, 0.730E-01
218      +, 0.000E+00, 0.000E+00, 0.416E-02, 0.141E-01
219      +, 0.347E+00, 0.635E+00, 0.000E+00, 0.000E+00
220      +, 0.000E+00, 0.000E+00, 0.127E-07, 0.488E-01
221      +, 0.948E+00, 0.285E-02, 0.000E+00, 0.000E+00
222      +, 0.000E+00, 0.532E-01, 0.956E-01, 0.848E+00
223      +, 0.311E-02,  4*0.0 /
224       DATA (TMP(J, 31),J=1,ISHLS) /
225      +  0.301E+00, 0.586E+00, 0.381E-01, 0.744E-01
226      +, 0.217E-03, 0.414E-03, 0.365E-02, 0.131E-01
227      +, 0.346E+00, 0.631E+00, 0.208E-02, 0.373E-02
228      +, 0.000E+00, 0.000E+00, 0.153E-07, 0.449E-01
229      +, 0.952E+00, 0.345E-02, 0.000E+00, 0.000E+00
230      +, 0.000E+00, 0.494E-01, 0.958E-01, 0.851E+00
231      +, 0.378E-02,  4*0.0 /
232       DATA (TMP(J, 32),J=1,ISHLS) /
233      +  0.301E+00, 0.583E+00, 0.388E-01, 0.759E-01
234      +, 0.633E-03, 0.122E-02, 0.322E-02, 0.121E-01
235      +, 0.344E+00, 0.624E+00, 0.595E-02, 0.106E-01
236      +, 0.000E+00, 0.000E+00, 0.187E-07, 0.421E-01
237      +, 0.954E+00, 0.376E-02, 0.000E+00, 0.000E+00
238      +, 0.000E+00, 0.465E-01, 0.961E-01, 0.853E+00
239      +, 0.414E-02,  4*0.0 /
240       DATA (TMP(J, 33),J=1,ISHLS) /
241      +  0.299E+00, 0.580E+00, 0.396E-01, 0.772E-01
242      +, 0.121E-02, 0.232E-02, 0.282E-02, 0.112E-01
243      +, 0.341E+00, 0.614E+00, 0.110E-01, 0.197E-01
244      +, 0.000E+00, 0.000E+00, 0.229E-07, 0.398E-01
245      +, 0.956E+00, 0.398E-02, 0.000E+00, 0.000E+00
246      +, 0.000E+00, 0.442E-01, 0.962E-01, 0.855E+00
247      +, 0.441E-02,  4*0.0 /
248       DATA (TMP(J, 34),J=1,ISHLS) /
249      +  0.298E+00, 0.577E+00, 0.403E-01, 0.787E-01
250      +, 0.191E-02, 0.368E-02, 0.249E-02, 0.105E-01
251      +, 0.336E+00, 0.604E+00, 0.170E-01, 0.302E-01
252      +, 0.000E+00, 0.000E+00, 0.279E-07, 0.379E-01
253      +, 0.958E+00, 0.414E-02, 0.000E+00, 0.000E+00
254      +, 0.000E+00, 0.425E-01, 0.964E-01, 0.857E+00
255      +, 0.462E-02,  4*0.0 /
256       DATA (TMP(J, 35),J=1,ISHLS) /
257      +  0.297E+00, 0.574E+00, 0.410E-01, 0.799E-01
258      +, 0.273E-02, 0.526E-02, 0.220E-02, 0.979E-02
259      +, 0.331E+00, 0.592E+00, 0.235E-01, 0.417E-01
260      +, 0.000E+00, 0.000E+00, 0.340E-07, 0.364E-01
261      +, 0.959E+00, 0.427E-02, 0.000E+00, 0.000E+00
262      +, 0.000E+00, 0.411E-01, 0.965E-01, 0.858E+00
263      +, 0.480E-02,  4*0.0 /
264       DATA (TMP(J, 36),J=1,ISHLS) /
265      +  0.296E+00, 0.571E+00, 0.416E-01, 0.812E-01
266      +, 0.366E-02, 0.704E-02, 0.194E-02, 0.921E-02
267      +, 0.325E+00, 0.579E+00, 0.305E-01, 0.539E-01
268      +, 0.000E+00, 0.000E+00, 0.413E-07, 0.352E-01
269      +, 0.960E+00, 0.438E-02, 0.000E+00, 0.000E+00
270      +, 0.000E+00, 0.400E-01, 0.966E-01, 0.858E+00
271      +, 0.495E-02,  4*0.0 /
272       DATA (TMP(J, 37),J=1,ISHLS) /
273      +  0.295E+00, 0.567E+00, 0.423E-01, 0.824E-01
274      +, 0.440E-02, 0.855E-02, 0.173E-02, 0.873E-02
275      +, 0.321E+00, 0.569E+00, 0.358E-01, 0.636E-01
276      +, 0.000E+00, 0.000E+00, 0.502E-07, 0.341E-01
277      +, 0.961E+00, 0.463E-02, 0.000E+00, 0.130E-03
278      +, 0.000E+00, 0.390E-01, 0.966E-01, 0.859E+00
279      +, 0.529E-02, 0.000E+00, 0.000E+00, 0.148E-03, 0.0 /
280       DATA (TMP(J, 38),J=1,ISHLS) /
281      +  0.294E+00, 0.565E+00, 0.429E-01, 0.837E-01
282      +, 0.506E-02, 0.985E-02, 0.155E-02, 0.836E-02
283      +, 0.318E+00, 0.560E+00, 0.402E-01, 0.712E-01
284      +, 0.000E+00, 0.000E+00, 0.607E-07, 0.331E-01
285      +, 0.962E+00, 0.489E-02, 0.000E+00, 0.421E-03
286      +, 0.000E+00, 0.382E-01, 0.967E-01, 0.859E+00
287      +, 0.561E-02, 0.000E+00, 0.000E+00, 0.484E-03, 0.0 /
288       DATA (TMP(J, 39),J=1,ISHLS) /
289      +  0.293E+00, 0.562E+00, 0.436E-01, 0.848E-01
290      +, 0.554E-02, 0.108E-01, 0.140E-02, 0.807E-02
291      +, 0.317E+00, 0.554E+00, 0.432E-01, 0.762E-01
292      +, 0.000E+00, 0.000E+00, 0.730E-07, 0.322E-01
293      +, 0.957E+00, 0.501E-02, 0.496E-02, 0.465E-03
294      +, 0.000E+00, 0.374E-01, 0.962E-01, 0.855E+00
295      +, 0.579E-02, 0.494E-03, 0.437E-02, 0.535E-03, 0.0 /
296       DATA (TMP(J, 40),J=1,ISHLS) /
297      +  0.293E+00, 0.560E+00, 0.441E-01, 0.860E-01
298      +, 0.595E-02, 0.116E-01, 0.127E-02, 0.785E-02
299      +, 0.316E+00, 0.549E+00, 0.457E-01, 0.802E-01
300      +, 0.000E+00, 0.000E+00, 0.871E-07, 0.312E-01
301      +, 0.952E+00, 0.508E-02, 0.117E-01, 0.480E-03
302      +, 0.000E+00, 0.367E-01, 0.957E-01, 0.850E+00
303      +, 0.593E-02, 0.116E-02, 0.103E-01, 0.559E-03, 0.0 /
304       DATA (TMP(J, 41),J=1,ISHLS) /
305      +  0.292E+00, 0.558E+00, 0.447E-01, 0.871E-01
306      +, 0.623E-02, 0.121E-01, 0.116E-02, 0.770E-02
307      +, 0.316E+00, 0.546E+00, 0.472E-01, 0.823E-01
308      +, 0.000E+00, 0.000E+00, 0.103E-06, 0.303E-01
309      +, 0.941E+00, 0.506E-02, 0.232E-01, 0.198E-03
310      +, 0.000E+00, 0.359E-01, 0.946E-01, 0.841E+00
311      +, 0.595E-02, 0.232E-02, 0.204E-01, 0.233E-03, 0.0 /
312       DATA (TMP(J, 42),J=1,ISHLS) /
313      +  0.292E+00, 0.556E+00, 0.453E-01, 0.881E-01
314      +, 0.656E-02, 0.127E-01, 0.107E-02, 0.758E-02
315      +, 0.315E+00, 0.542E+00, 0.491E-01, 0.853E-01
316      +, 0.000E+00, 0.000E+00, 0.123E-06, 0.295E-01
317      +, 0.933E+00, 0.510E-02, 0.321E-01, 0.194E-03
318      +, 0.000E+00, 0.352E-01, 0.939E-01, 0.833E+00
319      +, 0.604E-02, 0.319E-02, 0.282E-01, 0.231E-03, 0.0 /
320       DATA (TMP(J, 43),J=1,ISHLS) /
321      +  0.291E+00, 0.554E+00, 0.458E-01, 0.890E-01
322      +, 0.688E-02, 0.133E-01, 0.978E-03, 0.749E-02
323      +, 0.315E+00, 0.538E+00, 0.508E-01, 0.878E-01
324      +, 0.000E+00, 0.000E+00, 0.145E-06, 0.288E-01
325      +, 0.924E+00, 0.510E-02, 0.416E-01, 0.189E-03
326      +, 0.000E+00, 0.347E-01, 0.930E-01, 0.825E+00
327      +, 0.610E-02, 0.413E-02, 0.364E-01, 0.226E-03, 0.0 /
328       DATA (TMP(J, 44),J=1,ISHLS) /
329      +  0.291E+00, 0.552E+00, 0.463E-01, 0.901E-01
330      +, 0.716E-02, 0.139E-01, 0.902E-03, 0.745E-02
331      +, 0.315E+00, 0.535E+00, 0.524E-01, 0.899E-01
332      +, 0.000E+00, 0.000E+00, 0.172E-06, 0.281E-01
333      +, 0.915E+00, 0.509E-02, 0.515E-01, 0.184E-03
334      +, 0.000E+00, 0.341E-01, 0.921E-01, 0.817E+00
335      +, 0.614E-02, 0.511E-02, 0.450E-01, 0.222E-03, 0.0 /
336       DATA (TMP(J, 45),J=1,ISHLS) /
337      +  0.291E+00, 0.550E+00, 0.468E-01, 0.909E-01
338      +, 0.742E-02, 0.144E-01, 0.834E-03, 0.744E-02
339      +, 0.315E+00, 0.531E+00, 0.538E-01, 0.918E-01
340      +, 0.000E+00, 0.000E+00, 0.202E-06, 0.274E-01
341      +, 0.906E+00, 0.508E-02, 0.617E-01, 0.177E-03
342      +, 0.000E+00, 0.336E-01, 0.911E-01, 0.809E+00
343      +, 0.618E-02, 0.612E-02, 0.540E-01, 0.216E-03, 0.0 /
344       DATA (TMP(J, 46),J=1,ISHLS) /
345      +  0.290E+00, 0.548E+00, 0.472E-01, 0.918E-01
346      +, 0.762E-02, 0.147E-01, 0.773E-03, 0.747E-02
347      +, 0.315E+00, 0.529E+00, 0.547E-01, 0.925E-01
348      +, 0.000E+00, 0.000E+00, 0.236E-06, 0.267E-01
349      +, 0.894E+00, 0.502E-02, 0.744E-01, 0.000E+00
350      +, 0.000E+00, 0.332E-01, 0.900E-01, 0.799E+00
351      +, 0.617E-02, 0.737E-02, 0.646E-01,  2*0.0 /
352       DATA (TMP(J, 47),J=1,ISHLS) /
353      +  0.290E+00, 0.547E+00, 0.477E-01, 0.927E-01
354      +, 0.790E-02, 0.153E-01, 0.719E-03, 0.751E-02
355      +, 0.315E+00, 0.525E+00, 0.564E-01, 0.947E-01
356      +, 0.000E+00, 0.000E+00, 0.276E-06, 0.262E-01
357      +, 0.886E+00, 0.503E-02, 0.828E-01, 0.164E-03
358      +, 0.000E+00, 0.329E-01, 0.891E-01, 0.791E+00
359      +, 0.624E-02, 0.820E-02, 0.723E-01, 0.202E-03, 0.0 /
360       DATA (TMP(J, 48),J=1,ISHLS) /
361      +  0.290E+00, 0.545E+00, 0.481E-01, 0.934E-01
362      +, 0.823E-02, 0.159E-01, 0.668E-03, 0.758E-02
363      +, 0.315E+00, 0.521E+00, 0.582E-01, 0.971E-01
364      +, 0.000E+00, 0.000E+00, 0.323E-06, 0.257E-01
365      +, 0.878E+00, 0.505E-02, 0.905E-01, 0.400E-03
366      +, 0.000E+00, 0.326E-01, 0.884E-01, 0.784E+00
367      +, 0.632E-02, 0.894E-02, 0.790E-01, 0.500E-03, 0.0 /
368       DATA (TMP(J, 49),J=1,ISHLS) /
369      +  0.290E+00, 0.543E+00, 0.484E-01, 0.941E-01
370      +, 0.853E-02, 0.165E-01, 0.623E-03, 0.766E-02
371      +, 0.315E+00, 0.516E+00, 0.598E-01, 0.995E-01
372      +, 0.507E-03, 0.778E-03, 0.378E-06, 0.253E-01
373      +, 0.871E+00, 0.507E-02, 0.978E-01, 0.503E-03
374      +, 0.000E+00, 0.324E-01, 0.877E-01, 0.778E+00
375      +, 0.641E-02, 0.965E-02, 0.855E-01, 0.634E-03, 0.0 /
376       DATA (TMP(J, 50),J=1,ISHLS) /
377      +  0.289E+00, 0.541E+00, 0.489E-01, 0.948E-01
378      +, 0.884E-02, 0.171E-01, 0.580E-03, 0.776E-02
379      +, 0.314E+00, 0.511E+00, 0.614E-01, 0.101E+00
380      +, 0.140E-02, 0.218E-02, 0.439E-06, 0.249E-01
381      +, 0.865E+00, 0.510E-02, 0.104E+00, 0.574E-03
382      +, 0.000E+00, 0.322E-01, 0.870E-01, 0.772E+00
383      +, 0.651E-02, 0.103E-01, 0.914E-01, 0.732E-03, 0.0 /
384       DATA (TMP(J, 51),J=1,ISHLS) /
385      +  0.289E+00, 0.540E+00, 0.492E-01, 0.953E-01
386      +, 0.915E-02, 0.177E-01, 0.541E-03, 0.787E-02
387      +, 0.313E+00, 0.505E+00, 0.629E-01, 0.103E+00
388      +, 0.260E-02, 0.406E-02, 0.511E-06, 0.245E-01
389      +, 0.859E+00, 0.513E-02, 0.111E+00, 0.634E-03
390      +, 0.000E+00, 0.321E-01, 0.864E-01, 0.766E+00
391      +, 0.661E-02, 0.109E-01, 0.967E-01, 0.816E-03, 0.0 /
392       DATA (TMP(J, 52),J=1,ISHLS) /
393      +  0.289E+00, 0.537E+00, 0.496E-01, 0.961E-01
394      +, 0.945E-02, 0.183E-01, 0.505E-03, 0.800E-02
395      +, 0.312E+00, 0.500E+00, 0.642E-01, 0.105E+00
396      +, 0.404E-02, 0.631E-02, 0.593E-06, 0.242E-01
397      +, 0.854E+00, 0.516E-02, 0.116E+00, 0.688E-03
398      +, 0.000E+00, 0.320E-01, 0.859E-01, 0.762E+00
399      +, 0.672E-02, 0.114E-01, 0.102E+00, 0.894E-03, 0.0 /
400       DATA (TMP(J, 53),J=1,ISHLS) /
401      +  0.289E+00, 0.536E+00, 0.499E-01, 0.967E-01
402      +, 0.974E-02, 0.189E-01, 0.473E-03, 0.816E-02
403      +, 0.311E+00, 0.494E+00, 0.656E-01, 0.107E+00
404      +, 0.571E-02, 0.891E-02, 0.688E-06, 0.240E-01
405      +, 0.849E+00, 0.519E-02, 0.121E+00, 0.738E-03
406      +, 0.000E+00, 0.319E-01, 0.854E-01, 0.757E+00
407      +, 0.683E-02, 0.119E-01, 0.106E+00, 0.968E-03, 0.0 /
408       DATA (TMP(J, 54),J=1,ISHLS) /
409      +  0.289E+00, 0.534E+00, 0.502E-01, 0.973E-01
410      +, 0.100E-01, 0.195E-01, 0.442E-03, 0.833E-02
411      +, 0.309E+00, 0.488E+00, 0.667E-01, 0.108E+00
412      +, 0.759E-02, 0.118E-01, 0.792E-06, 0.237E-01
413      +, 0.844E+00, 0.523E-02, 0.126E+00, 0.784E-03
414      +, 0.000E+00, 0.319E-01, 0.849E-01, 0.753E+00
415      +, 0.694E-02, 0.123E-01, 0.110E+00, 0.104E-02, 0.0 /
416       DATA (TMP(J, 55),J=1,ISHLS) /
417      +  0.288E+00, 0.533E+00, 0.505E-01, 0.978E-01
418      +, 0.103E-01, 0.200E-01, 0.415E-03, 0.853E-02
419      +, 0.309E+00, 0.482E+00, 0.679E-01, 0.109E+00
420      +, 0.900E-02, 0.142E-01, 0.915E-06, 0.235E-01
421      +, 0.840E+00, 0.526E-02, 0.130E+00, 0.857E-03
422      +, 0.000E+00, 0.319E-01, 0.845E-01, 0.749E+00
423      +, 0.706E-02, 0.128E-01, 0.114E+00, 0.115E-02, 0.0 /
424       DATA (TMP(J, 56),J=1,ISHLS) /
425      +  0.289E+00, 0.531E+00, 0.508E-01, 0.984E-01
426      +, 0.106E-01, 0.206E-01, 0.391E-03, 0.876E-02
427      +, 0.308E+00, 0.477E+00, 0.692E-01, 0.110E+00
428      +, 0.103E-01, 0.162E-01, 0.105E-05, 0.232E-01
429      +, 0.836E+00, 0.529E-02, 0.135E+00, 0.928E-03
430      +, 0.000E+00, 0.320E-01, 0.841E-01, 0.745E+00
431      +, 0.718E-02, 0.132E-01, 0.118E+00, 0.126E-02, 0.0 /
432       DATA (TMP(J, 57),J=1,ISHLS) /
433      +  0.289E+00, 0.530E+00, 0.510E-01, 0.989E-01
434      +, 0.108E-01, 0.211E-01, 0.368E-03, 0.903E-02
435      +, 0.308E+00, 0.472E+00, 0.705E-01, 0.112E+00
436      +, 0.112E-01, 0.175E-01, 0.121E-05, 0.229E-01
437      +, 0.831E+00, 0.531E-02, 0.139E+00, 0.980E-03
438      +, 0.108E-02, 0.320E-01, 0.836E-01, 0.740E+00
439      +, 0.730E-02, 0.135E-01, 0.121E+00, 0.134E-02, 0.103E-02 /
440       DATA (TMP(J, 58),J=1,ISHLS) /
441      +  0.288E+00, 0.529E+00, 0.514E-01, 0.995E-01
442      +, 0.109E-01, 0.212E-01, 0.348E-03, 0.938E-02
443      +, 0.310E+00, 0.471E+00, 0.710E-01, 0.111E+00
444      +, 0.104E-01, 0.161E-01, 0.140E-05, 0.229E-01
445      +, 0.833E+00, 0.530E-02, 0.138E+00, 0.927E-03
446      +, 0.000E+00, 0.323E-01, 0.839E-01, 0.742E+00
447      +, 0.736E-02, 0.135E-01, 0.120E+00, 0.128E-02, 0.0 /
448       DATA (TMP(J, 59),J=1,ISHLS) /
449      +  0.289E+00, 0.527E+00, 0.517E-01, 0.100E+00
450      +, 0.110E-01, 0.214E-01, 0.330E-03, 0.973E-02
451      +, 0.311E+00, 0.469E+00, 0.718E-01, 0.112E+00
452      +, 0.105E-01, 0.160E-01, 0.161E-05, 0.228E-01
453      +, 0.831E+00, 0.532E-02, 0.140E+00, 0.922E-03
454      +, 0.000E+00, 0.325E-01, 0.838E-01, 0.740E+00
455      +, 0.745E-02, 0.136E-01, 0.121E+00, 0.129E-02, 0.0 /
456       DATA (TMP(J, 60),J=1,ISHLS) /
457      +  0.289E+00, 0.526E+00, 0.520E-01, 0.101E+00
458      +, 0.112E-01, 0.217E-01, 0.313E-03, 0.101E-01
459      +, 0.313E+00, 0.466E+00, 0.725E-01, 0.112E+00
460      +, 0.105E-01, 0.158E-01, 0.185E-05, 0.227E-01
461      +, 0.830E+00, 0.532E-02, 0.141E+00, 0.918E-03
462      +, 0.000E+00, 0.327E-01, 0.834E-01, 0.740E+00
463      +, 0.753E-02, 0.136E-01, 0.122E+00, 0.129E-02, 0.0 /
464       DATA (TMP(J, 61),J=1,ISHLS) /
465      +  0.289E+00, 0.524E+00, 0.522E-01, 0.101E+00
466      +, 0.113E-01, 0.219E-01, 0.297E-03, 0.105E-01
467      +, 0.314E+00, 0.463E+00, 0.732E-01, 0.112E+00
468      +, 0.105E-01, 0.156E-01, 0.213E-05, 0.226E-01
469      +, 0.830E+00, 0.532E-02, 0.142E+00, 0.912E-03
470      +, 0.000E+00, 0.330E-01, 0.834E-01, 0.739E+00
471      +, 0.761E-02, 0.137E-01, 0.122E+00, 0.130E-02, 0.0 /
472       DATA (TMP(J, 62),J=1,ISHLS) /
473      +  0.289E+00, 0.523E+00, 0.525E-01, 0.101E+00
474      +, 0.114E-01, 0.221E-01, 0.283E-03, 0.110E-01
475      +, 0.316E+00, 0.461E+00, 0.740E-01, 0.112E+00
476      +, 0.105E-01, 0.155E-01, 0.243E-05, 0.225E-01
477      +, 0.829E+00, 0.531E-02, 0.142E+00, 0.905E-03
478      +, 0.000E+00, 0.333E-01, 0.835E-01, 0.737E+00
479      +, 0.772E-02, 0.138E-01, 0.123E+00, 0.131E-02, 0.0 /
480       DATA (TMP(J, 63),J=1,ISHLS) /
481      +  0.289E+00, 0.522E+00, 0.527E-01, 0.102E+00
482      +, 0.115E-01, 0.223E-01, 0.269E-03, 0.115E-01
483      +, 0.318E+00, 0.458E+00, 0.747E-01, 0.112E+00
484      +, 0.104E-01, 0.153E-01, 0.277E-05, 0.224E-01
485      +, 0.828E+00, 0.532E-02, 0.143E+00, 0.898E-03
486      +, 0.000E+00, 0.335E-01, 0.832E-01, 0.737E+00
487      +, 0.779E-02, 0.138E-01, 0.123E+00, 0.131E-02, 0.0 /
488       DATA (TMP(J, 64),J=1,ISHLS) /
489      +  0.290E+00, 0.521E+00, 0.529E-01, 0.102E+00
490      +, 0.116E-01, 0.226E-01, 0.257E-03, 0.120E-01
491      +, 0.318E+00, 0.454E+00, 0.755E-01, 0.112E+00
492      +, 0.111E-01, 0.161E-01, 0.315E-05, 0.223E-01
493      +, 0.826E+00, 0.533E-02, 0.145E+00, 0.933E-03
494      +, 0.865E-03, 0.338E-01, 0.830E-01, 0.734E+00
495      +, 0.790E-02, 0.139E-01, 0.125E+00, 0.138E-02, 0.798E-03 /
496       DATA (TMP(J, 65),J=1,ISHLS) /
497      +  0.290E+00, 0.520E+00, 0.532E-01, 0.103E+00
498      +, 0.117E-01, 0.226E-01, 0.245E-03, 0.126E-01
499      +, 0.321E+00, 0.453E+00, 0.760E-01, 0.112E+00
500      +, 0.104E-01, 0.148E-01, 0.360E-05, 0.223E-01
501      +, 0.827E+00, 0.532E-02, 0.144E+00, 0.887E-03
502      +, 0.000E+00, 0.342E-01, 0.832E-01, 0.735E+00
503      +, 0.799E-02, 0.139E-01, 0.124E+00, 0.133E-02, 0.0 /
504       DATA (TMP(J, 66),J=1,ISHLS) /
505      +  0.290E+00, 0.519E+00, 0.534E-01, 0.103E+00
506      +, 0.118E-01, 0.228E-01, 0.234E-03, 0.133E-01
507      +, 0.323E+00, 0.450E+00, 0.768E-01, 0.112E+00
508      +, 0.104E-01, 0.146E-01, 0.409E-05, 0.223E-01
509      +, 0.827E+00, 0.534E-02, 0.145E+00, 0.880E-03
510      +, 0.000E+00, 0.346E-01, 0.832E-01, 0.735E+00
511      +, 0.809E-02, 0.139E-01, 0.124E+00, 0.133E-02, 0.0 /
512       DATA (TMP(J, 67),J=1,ISHLS) /
513      +  0.291E+00, 0.517E+00, 0.536E-01, 0.104E+00
514      +, 0.118E-01, 0.229E-01, 0.224E-03, 0.139E-01
515      +, 0.325E+00, 0.447E+00, 0.774E-01, 0.111E+00
516      +, 0.104E-01, 0.144E-01, 0.465E-05, 0.222E-01
517      +, 0.826E+00, 0.534E-02, 0.146E+00, 0.876E-03
518      +, 0.000E+00, 0.349E-01, 0.831E-01, 0.734E+00
519      +, 0.819E-02, 0.139E-01, 0.125E+00, 0.134E-02, 0.0 /
520       DATA (TMP(J, 68),J=1,ISHLS) /
521      +  0.291E+00, 0.516E+00, 0.538E-01, 0.104E+00
522      +, 0.119E-01, 0.230E-01, 0.214E-03, 0.147E-01
523      +, 0.327E+00, 0.445E+00, 0.780E-01, 0.111E+00
524      +, 0.104E-01, 0.142E-01, 0.527E-05, 0.222E-01
525      +, 0.826E+00, 0.534E-02, 0.146E+00, 0.871E-03
526      +, 0.000E+00, 0.353E-01, 0.830E-01, 0.734E+00
527      +, 0.830E-02, 0.139E-01, 0.125E+00, 0.134E-02, 0.0 /
528       DATA (TMP(J, 69),J=1,ISHLS) /
529      +  0.291E+00, 0.515E+00, 0.540E-01, 0.104E+00
530      +, 0.120E-01, 0.232E-01, 0.205E-03, 0.155E-01
531      +, 0.329E+00, 0.442E+00, 0.788E-01, 0.111E+00
532      +, 0.104E-01, 0.140E-01, 0.598E-05, 0.222E-01
533      +, 0.825E+00, 0.535E-02, 0.146E+00, 0.866E-03
534      +, 0.000E+00, 0.357E-01, 0.830E-01, 0.733E+00
535      +, 0.840E-02, 0.140E-01, 0.125E+00, 0.135E-02, 0.0 /
536       DATA (TMP(J, 70),J=1,ISHLS) /
537      +  0.292E+00, 0.514E+00, 0.542E-01, 0.105E+00
538      +, 0.121E-01, 0.233E-01, 0.197E-03, 0.163E-01
539      +, 0.331E+00, 0.439E+00, 0.794E-01, 0.110E+00
540      +, 0.104E-01, 0.138E-01, 0.677E-05, 0.221E-01
541      +, 0.825E+00, 0.536E-02, 0.147E+00, 0.861E-03
542      +, 0.000E+00, 0.361E-01, 0.828E-01, 0.732E+00
543      +, 0.851E-02, 0.140E-01, 0.125E+00, 0.136E-02, 0.0 /
544       DATA (TMP(J, 71),J=1,ISHLS) /
545      +  0.292E+00, 0.513E+00, 0.543E-01, 0.105E+00
546      +, 0.121E-01, 0.235E-01, 0.189E-03, 0.172E-01
547      +, 0.332E+00, 0.435E+00, 0.801E-01, 0.110E+00
548      +, 0.110E-01, 0.145E-01, 0.765E-05, 0.221E-01
549      +, 0.823E+00, 0.537E-02, 0.148E+00, 0.890E-03
550      +, 0.649E-03, 0.365E-01, 0.827E-01, 0.730E+00
551      +, 0.863E-02, 0.140E-01, 0.126E+00, 0.142E-02, 0.572E-03 /
552       DATA (TMP(J, 72),J=1,ISHLS) /
553      +  0.293E+00, 0.512E+00, 0.545E-01, 0.105E+00
554      +, 0.122E-01, 0.237E-01, 0.181E-03, 0.181E-01
555      +, 0.334E+00, 0.430E+00, 0.810E-01, 0.110E+00
556      +, 0.116E-01, 0.152E-01, 0.858E-05, 0.220E-01
557      +, 0.821E+00, 0.539E-02, 0.149E+00, 0.921E-03
558      +, 0.160E-02, 0.369E-01, 0.825E-01, 0.728E+00
559      +, 0.875E-02, 0.141E-01, 0.127E+00, 0.149E-02, 0.142E-02 /
560       DATA (TMP(J, 73),J=1,ISHLS) /
561      +  0.293E+00, 0.511E+00, 0.546E-01, 0.106E+00
562      +, 0.124E-01, 0.239E-01, 0.174E-03, 0.191E-01
563      +, 0.335E+00, 0.426E+00, 0.818E-01, 0.110E+00
564      +, 0.122E-01, 0.158E-01, 0.968E-05, 0.220E-01
565      +, 0.818E+00, 0.539E-02, 0.151E+00, 0.950E-03
566      +, 0.278E-02, 0.373E-01, 0.822E-01, 0.726E+00
567      +, 0.888E-02, 0.142E-01, 0.127E+00, 0.156E-02, 0.250E-02 /
568       DATA (TMP(J, 74),J=1,ISHLS) /
569      +  0.293E+00, 0.509E+00, 0.548E-01, 0.106E+00
570      +, 0.125E-01, 0.242E-01, 0.167E-03, 0.202E-01
571      +, 0.337E+00, 0.421E+00, 0.826E-01, 0.110E+00
572      +, 0.128E-01, 0.165E-01, 0.109E-04, 0.219E-01
573      +, 0.816E+00, 0.541E-02, 0.152E+00, 0.980E-03
574      +, 0.420E-02, 0.377E-01, 0.820E-01, 0.723E+00
575      +, 0.901E-02, 0.143E-01, 0.128E+00, 0.162E-02, 0.377E-02 /
576       DATA (TMP(J, 75),J=1,ISHLS) /
577      +  0.294E+00, 0.508E+00, 0.549E-01, 0.106E+00
578      +, 0.126E-01, 0.244E-01, 0.161E-03, 0.213E-01
579      +, 0.338E+00, 0.417E+00, 0.835E-01, 0.109E+00
580      +, 0.134E-01, 0.170E-01, 0.122E-04, 0.219E-01
581      +, 0.813E+00, 0.543E-02, 0.153E+00, 0.101E-02
582      +, 0.581E-02, 0.381E-01, 0.817E-01, 0.721E+00
583      +, 0.914E-02, 0.144E-01, 0.129E+00, 0.169E-02, 0.522E-02 /
584       DATA (TMP(J, 76),J=1,ISHLS) /
585      +  0.294E+00, 0.507E+00, 0.550E-01, 0.106E+00
586      +, 0.127E-01, 0.246E-01, 0.155E-03, 0.225E-01
587      +, 0.340E+00, 0.412E+00, 0.844E-01, 0.109E+00
588      +, 0.140E-01, 0.176E-01, 0.137E-04, 0.218E-01
589      +, 0.810E+00, 0.544E-02, 0.154E+00, 0.104E-02
590      +, 0.762E-02, 0.384E-01, 0.814E-01, 0.718E+00
591      +, 0.928E-02, 0.145E-01, 0.130E+00, 0.176E-02, 0.684E-02 /
592       DATA (TMP(J, 77),J=1,ISHLS) /
593      +  0.295E+00, 0.506E+00, 0.552E-01, 0.107E+00
594      +, 0.128E-01, 0.249E-01, 0.149E-03, 0.239E-01
595      +, 0.342E+00, 0.408E+00, 0.854E-01, 0.109E+00
596      +, 0.144E-01, 0.178E-01, 0.153E-04, 0.217E-01
597      +, 0.806E+00, 0.545E-02, 0.155E+00, 0.105E-02
598      +, 0.109E-01, 0.388E-01, 0.810E-01, 0.714E+00
599      +, 0.940E-02, 0.145E-01, 0.131E+00, 0.181E-02, 0.958E-02 /
600       DATA (TMP(J, 78),J=1,ISHLS) /
601      +  0.295E+00, 0.505E+00, 0.553E-01, 0.107E+00
602      +, 0.129E-01, 0.251E-01, 0.144E-03, 0.253E-01
603      +, 0.343E+00, 0.403E+00, 0.863E-01, 0.109E+00
604      +, 0.151E-01, 0.184E-01, 0.172E-04, 0.217E-01
605      +, 0.803E+00, 0.547E-02, 0.157E+00, 0.109E-02
606      +, 0.125E-01, 0.393E-01, 0.807E-01, 0.711E+00
607      +, 0.955E-02, 0.146E-01, 0.132E+00, 0.188E-02, 0.111E-01 /
608       DATA (TMP(J, 79),J=1,ISHLS) /
609      +  0.296E+00, 0.503E+00, 0.554E-01, 0.107E+00
610      +, 0.130E-01, 0.254E-01, 0.140E-03, 0.270E-01
611      +, 0.348E+00, 0.392E+00, 0.881E-01, 0.110E+00
612      +, 0.158E-01, 0.190E-01, 0.192E-04, 0.216E-01
613      +, 0.799E+00, 0.549E-02, 0.158E+00, 0.111E-02
614      +, 0.148E-01, 0.397E-01, 0.803E-01, 0.708E+00
615      +, 0.969E-02, 0.147E-01, 0.132E+00, 0.195E-02, 0.131E-01 /
616       DATA (TMP(J, 80),J=1,ISHLS) /
617      +  0.296E+00, 0.502E+00, 0.555E-01, 0.107E+00
618      +, 0.131E-01, 0.256E-01, 0.133E-03, 0.283E-01
619      +, 0.347E+00, 0.393E+00, 0.882E-01, 0.108E+00
620      +, 0.163E-01, 0.194E-01, 0.214E-04, 0.216E-01
621      +, 0.796E+00, 0.551E-02, 0.159E+00, 0.114E-02
622      +, 0.165E-01, 0.401E-01, 0.801E-01, 0.705E+00
623      +, 0.984E-02, 0.147E-01, 0.133E+00, 0.203E-02, 0.147E-01 /
624       DATA (TMP(J, 81),J=1,ISHLS) /
625      +  0.297E+00, 0.501E+00, 0.556E-01, 0.108E+00
626      +, 0.132E-01, 0.259E-01, 0.129E-03, 0.301E-01
627      +, 0.348E+00, 0.388E+00, 0.891E-01, 0.108E+00
628      +, 0.170E-01, 0.200E-01, 0.240E-04, 0.215E-01
629      +, 0.793E+00, 0.553E-02, 0.160E+00, 0.117E-02
630      +, 0.183E-01, 0.406E-01, 0.798E-01, 0.703E+00
631      +, 0.999E-02, 0.148E-01, 0.134E+00, 0.210E-02, 0.163E-01 /
632       DATA (TMP(J, 82),J=1,ISHLS) /
633      +  0.297E+00, 0.500E+00, 0.557E-01, 0.108E+00
634      +, 0.133E-01, 0.261E-01, 0.124E-03, 0.319E-01
635      +, 0.350E+00, 0.382E+00, 0.901E-01, 0.107E+00
636      +, 0.176E-01, 0.206E-01, 0.267E-04, 0.215E-01
637      +, 0.790E+00, 0.555E-02, 0.161E+00, 0.120E-02
638      +, 0.200E-01, 0.411E-01, 0.795E-01, 0.700E+00
639      +, 0.101E-01, 0.149E-01, 0.135E+00, 0.218E-02, 0.179E-01 /
640       DATA (TMP(J, 83),J=1,ISHLS) /
641      +  0.298E+00, 0.498E+00, 0.557E-01, 0.108E+00
642      +, 0.134E-01, 0.263E-01, 0.120E-03, 0.338E-01
643      +, 0.352E+00, 0.377E+00, 0.910E-01, 0.107E+00
644      +, 0.183E-01, 0.212E-01, 0.298E-04, 0.215E-01
645      +, 0.787E+00, 0.557E-02, 0.163E+00, 0.123E-02
646      +, 0.216E-01, 0.415E-01, 0.792E-01, 0.697E+00
647      +, 0.103E-01, 0.150E-01, 0.135E+00, 0.226E-02, 0.194E-01 /
648       DATA (TMP(J, 84),J=1,ISHLS) /
649      +  0.299E+00, 0.497E+00, 0.558E-01, 0.108E+00
650      +, 0.135E-01, 0.266E-01, 0.116E-03, 0.358E-01
651      +, 0.354E+00, 0.371E+00, 0.921E-01, 0.106E+00
652      +, 0.190E-01, 0.218E-01, 0.332E-04, 0.214E-01
653      +, 0.785E+00, 0.560E-02, 0.164E+00, 0.126E-02
654      +, 0.232E-01, 0.420E-01, 0.789E-01, 0.694E+00
655      +, 0.105E-01, 0.150E-01, 0.136E+00, 0.234E-02, 0.208E-01 /
656       DATA (TMP(J, 85),J=1,ISHLS) /
657      +  0.299E+00, 0.496E+00, 0.559E-01, 0.108E+00
658      +, 0.135E-01, 0.268E-01, 0.112E-03, 0.381E-01
659      +, 0.356E+00, 0.365E+00, 0.930E-01, 0.106E+00
660      +, 0.197E-01, 0.223E-01, 0.369E-04, 0.214E-01
661      +, 0.782E+00, 0.562E-02, 0.165E+00, 0.129E-02
662      +, 0.247E-01, 0.426E-01, 0.786E-01, 0.692E+00
663      +, 0.106E-01, 0.151E-01, 0.137E+00, 0.243E-02, 0.222E-01 /
664       DATA (TMP(J, 86),J=1,ISHLS) /
665      +  0.300E+00, 0.495E+00, 0.560E-01, 0.109E+00
666      +, 0.136E-01, 0.270E-01, 0.108E-03, 0.405E-01
667      +, 0.358E+00, 0.360E+00, 0.940E-01, 0.105E+00
668      +, 0.204E-01, 0.227E-01, 0.411E-04, 0.214E-01
669      +, 0.779E+00, 0.565E-02, 0.166E+00, 0.133E-02
670      +, 0.262E-01, 0.431E-01, 0.783E-01, 0.689E+00
671      +, 0.108E-01, 0.151E-01, 0.138E+00, 0.251E-02, 0.235E-01 /
672       DATA (TMP(J, 87),J=1,ISHLS) /
673      +  0.301E+00, 0.493E+00, 0.560E-01, 0.109E+00
674      +, 0.137E-01, 0.273E-01, 0.105E-03, 0.430E-01
675      +, 0.360E+00, 0.354E+00, 0.951E-01, 0.104E+00
676      +, 0.211E-01, 0.232E-01, 0.457E-04, 0.214E-01
677      +, 0.776E+00, 0.568E-02, 0.167E+00, 0.136E-02
678      +, 0.276E-01, 0.437E-01, 0.781E-01, 0.687E+00
679      +, 0.110E-01, 0.152E-01, 0.138E+00, 0.260E-02, 0.248E-01 /
680       DATA (TMP(J, 88),J=1,ISHLS) /
681      +  0.301E+00, 0.492E+00, 0.561E-01, 0.109E+00
682      +, 0.138E-01, 0.275E-01, 0.101E-03, 0.456E-01
683      +, 0.362E+00, 0.348E+00, 0.960E-01, 0.103E+00
684      +, 0.218E-01, 0.236E-01, 0.508E-04, 0.214E-01
685      +, 0.774E+00, 0.571E-02, 0.168E+00, 0.139E-02
686      +, 0.291E-01, 0.425E-01, 0.779E-01, 0.685E+00
687      +, 0.112E-01, 0.153E-01, 0.139E+00, 0.269E-02, 0.261E-01 /
688       DATA (TMP(J, 89),J=1,ISHLS) /
689      +  0.302E+00, 0.491E+00, 0.561E-01, 0.109E+00
690      +, 0.139E-01, 0.277E-01, 0.978E-04, 0.485E-01
691      +, 0.364E+00, 0.341E+00, 0.971E-01, 0.103E+00
692      +, 0.224E-01, 0.240E-01, 0.564E-04, 0.215E-01
693      +, 0.771E+00, 0.574E-02, 0.170E+00, 0.142E-02
694      +, 0.305E-01, 0.448E-01, 0.775E-01, 0.681E+00
695      +, 0.113E-01, 0.153E-01, 0.140E+00, 0.277E-02, 0.272E-01 /
696       DATA (TMP(J, 90),J=1,ISHLS) /
697      +  0.303E+00, 0.490E+00, 0.561E-01, 0.109E+00
698      +, 0.139E-01, 0.280E-01, 0.946E-04, 0.515E-01
699      +, 0.366E+00, 0.335E+00, 0.981E-01, 0.102E+00
700      +, 0.231E-01, 0.243E-01, 0.626E-04, 0.215E-01
701      +, 0.769E+00, 0.577E-02, 0.171E+00, 0.145E-02
702      +, 0.319E-01, 0.454E-01, 0.773E-01, 0.679E+00
703      +, 0.115E-01, 0.153E-01, 0.140E+00, 0.286E-02, 0.284E-01 /
704       DATA (TMP(J, 91),J=1,ISHLS) /
705      +  0.304E+00, 0.488E+00, 0.562E-01, 0.110E+00
706      +, 0.140E-01, 0.282E-01, 0.913E-04, 0.547E-01
707      +, 0.368E+00, 0.329E+00, 0.991E-01, 0.101E+00
708      +, 0.237E-01, 0.244E-01, 0.695E-04, 0.215E-01
709      +, 0.766E+00, 0.580E-02, 0.172E+00, 0.147E-02
710      +, 0.326E-01, 0.460E-01, 0.770E-01, 0.677E+00
711      +, 0.117E-01, 0.154E-01, 0.141E+00, 0.294E-02, 0.290E-01 /
712       DATA (TMP(J, 92),J=1,ISHLS) /
713      +  0.304E+00, 0.487E+00, 0.562E-01, 0.110E+00
714      +, 0.140E-01, 0.284E-01, 0.882E-04, 0.582E-01
715      +, 0.371E+00, 0.322E+00, 0.100E+00, 0.997E-01
716      +, 0.243E-01, 0.245E-01, 0.772E-04, 0.215E-01
717      +, 0.764E+00, 0.584E-02, 0.173E+00, 0.150E-02
718      +, 0.335E-01, 0.466E-01, 0.768E-01, 0.675E+00
719      +, 0.119E-01, 0.154E-01, 0.142E+00, 0.302E-02, 0.299E-01 /
720       DATA (TMP(J, 93),J=1,ISHLS) /
721      +  0.305E+00, 0.486E+00, 0.562E-01, 0.110E+00
722      +, 0.141E-01, 0.286E-01, 0.852E-04, 0.619E-01
723      +, 0.373E+00, 0.316E+00, 0.101E+00, 0.987E-01
724      +, 0.249E-01, 0.246E-01, 0.855E-04, 0.216E-01
725      +, 0.762E+00, 0.587E-02, 0.174E+00, 0.153E-02
726      +, 0.344E-01, 0.473E-01, 0.766E-01, 0.673E+00
727      +, 0.121E-01, 0.155E-01, 0.142E+00, 0.310E-02, 0.304E-01 /
728       DATA (TMP(J, 94),J=1,ISHLS) /
729      +  0.306E+00, 0.484E+00, 0.561E-01, 0.110E+00
730      +, 0.141E-01, 0.298E-01, 0.822E-04, 0.658E-01
731      +, 0.376E+00, 0.309E+00, 0.102E+00, 0.974E-01
732      +, 0.254E-01, 0.246E-01, 0.949E-04, 0.216E-01
733      +, 0.760E+00, 0.591E-02, 0.175E+00, 0.155E-02
734      +, 0.351E-01, 0.480E-01, 0.764E-01, 0.671E+00
735      +, 0.123E-01, 0.156E-01, 0.143E+00, 0.318E-02, 0.308E-01 /
736       DATA (TMP(J, 95),J=1,ISHLS) /
737      +  0.307E+00, 0.483E+00, 0.562E-01, 0.110E+00
738      +, 0.142E-01, 0.291E-01, 0.793E-04, 0.700E-01
739      +, 0.378E+00, 0.302E+00, 0.103E+00, 0.962E-01
740      +, 0.259E-01, 0.245E-01, 0.105E-03, 0.216E-01
741      +, 0.758E+00, 0.595E-02, 0.176E+00, 0.158E-02
742      +, 0.359E-01, 0.486E-01, 0.763E-01, 0.669E+00
743      +, 0.125E-01, 0.156E-01, 0.143E+00, 0.326E-02, 0.314E-01 /
744       DATA (TMP(J, 96),J=1,ISHLS) /
745      +  0.308E+00, 0.482E+00, 0.562E-01, 0.111E+00
746      +, 0.142E-01, 0.293E-01, 0.764E-04, 0.745E-01
747      +, 0.381E+00, 0.295E+00, 0.104E+00, 0.948E-01
748      +, 0.265E-01, 0.245E-01, 0.116E-03, 0.217E-01
749      +, 0.756E+00, 0.598E-02, 0.178E+00, 0.160E-02
750      +, 0.369E-01, 0.493E-01, 0.761E-01, 0.667E+00
751      +, 0.127E-01, 0.156E-01, 0.144E+00, 0.335E-02, 0.322E-01 /
752       DATA (TMP(J, 97),J=1,ISHLS) /
753      +  0.309E+00, 0.481E+00, 0.561E-01, 0.111E+00
754      +, 0.143E-01, 0.295E-01, 0.736E-04, 0.792E-01
755      +, 0.383E+00, 0.287E+00, 0.106E+00, 0.934E-01
756      +, 0.271E-01, 0.244E-01, 0.129E-03, 0.218E-01
757      +, 0.754E+00, 0.603E-02, 0.179E+00, 0.163E-02
758      +, 0.376E-01, 0.500E-01, 0.757E-01, 0.665E+00
759      +, 0.129E-01, 0.157E-01, 0.144E+00, 0.343E-02, 0.327E-01 /
760       DATA (TMP(J, 98),J=1,ISHLS) /
761      +  0.310E+00, 0.479E+00, 0.561E-01, 0.111E+00
762      +, 0.143E-01, 0.297E-01, 0.708E-04, 0.843E-01
763      +, 0.386E+00, 0.280E+00, 0.106E+00, 0.919E-01
764      +, 0.276E-01, 0.242E-01, 0.143E-03, 0.218E-01
765      +, 0.752E+00, 0.607E-02, 0.180E+00, 0.165E-02
766      +, 0.381E-01, 0.507E-01, 0.756E-01, 0.664E+00
767      +, 0.131E-01, 0.157E-01, 0.145E+00, 0.351E-02, 0.330E-01 /
768       DATA (TMP(J, 99),J=1,ISHLS) /
769      +  0.311E+00, 0.478E+00, 0.561E-01, 0.111E+00
770      +, 0.143E-01, 0.299E-01, 0.679E-04, 0.897E-01
771      +, 0.388E+00, 0.272E+00, 0.107E+00, 0.903E-01
772      +, 0.281E-01, 0.240E-01, 0.158E-03, 0.219E-01
773      +, 0.751E+00, 0.611E-02, 0.181E+00, 0.168E-02
774      +, 0.388E-01, 0.515E-01, 0.754E-01, 0.662E+00
775      +, 0.133E-01, 0.157E-01, 0.145E+00, 0.359E-02, 0.335E-01 /
776       DATA (TMP(J,100),J=1,ISHLS) /
777      +  0.312E+00, 0.477E+00, 0.560E-01, 0.111E+00
778      +, 0.144E-01, 0.301E-01, 0.650E-04, 0.954E-01
779      +, 0.391E+00, 0.264E+00, 0.108E+00, 0.885E-01
780      +, 0.286E-01, 0.238E-01, 0.175E-03, 0.220E-01
781      +, 0.749E+00, 0.616E-02, 0.182E+00, 0.170E-02
782      +, 0.395E-01, 0.522E-01, 0.752E-01, 0.660E+00
783      +, 0.135E-01, 0.158E-01, 0.146E+00, 0.367E-02, 0.339E-01 /
784 C
785       DATA (TMP1(I,  1),I=1,4) / 4*0.0 /
786       DATA (TMP1(I,  2),I=1,4) / 4*0.0 /
787       DATA (TMP1(I,  3),I=1,4) / 4*0.0 /
788       DATA (TMP1(I,  4),I=1,4) / 4*0.0 /
789       DATA (TMP1(I,  5),I=1,4) /  0.170E-02, 3*0.0 /
790       DATA (TMP1(I,  6),I=1,4) /  0.280E-02, 3*0.0 /
791       DATA (TMP1(I,  7),I=1,4) /  0.520E-02, 3*0.0 /
792       DATA (TMP1(I,  8),I=1,4) /  0.830E-02, 3*0.0 /
793       DATA (TMP1(I,  9),I=1,4) /  0.130E-01, 3*0.0 /
794       DATA (TMP1(I, 10),I=1,4) /  0.180E-01, 3*0.0 /
795       DATA (TMP1(I, 11),I=1,4) /  0.230E-01, 3*0.0 /
796       DATA (TMP1(I, 12),I=1,4) /
797      +  0.300E-01, 0.290E-04, 0.120E-02, 0.120E-02 /
798       DATA (TMP1(I, 13),I=1,4) /
799      +  0.390E-01, 0.260E-04, 0.750E-03, 0.750E-03 /
800       DATA (TMP1(I, 14),I=1,4) /
801      +  0.500E-01, 0.300E-04, 0.370E-03, 0.380E-03 /
802       DATA (TMP1(I, 15),I=1,4) /
803      +  0.630E-01, 0.390E-04, 0.310E-03, 0.310E-03 /
804       DATA (TMP1(I, 16),I=1,4) /
805      +  0.780E-01, 0.740E-04, 0.260E-03, 0.260E-03 /
806       DATA (TMP1(I, 17),I=1,4) /
807      +  0.970E-01, 0.120E-03, 0.240E-03, 0.240E-03 /
808       DATA (TMP1(I, 18),I=1,4) /
809      +  0.118E+00, 0.180E-03, 0.220E-03, 0.220E-03 /
810       DATA (TMP1(I, 19),I=1,4) /
811      +  0.140E+00, 0.240E-03, 0.270E-03, 0.270E-03 /
812       DATA (TMP1(I, 20),I=1,4) /
813      +  0.163E+00, 0.310E-03, 0.330E-03, 0.330E-03 /
814       DATA (TMP1(I, 21),I=1,4) /
815      +  0.188E+00, 0.390E-03, 0.840E-03, 0.840E-03 /
816       DATA (TMP1(I, 22),I=1,4) /
817      +  0.214E+00, 0.470E-03, 0.150E-02, 0.150E-02 /
818       DATA (TMP1(I, 23),I=1,4) /
819      +  0.243E+00, 0.580E-03, 0.260E-02, 0.260E-02 /
820       DATA (TMP1(I, 24),I=1,4) /
821      +  0.275E+00, 0.710E-03, 0.370E-02, 0.370E-02 /
822       DATA (TMP1(I, 25),I=1,4) /
823      +  0.308E+00, 0.840E-03, 0.500E-02, 0.500E-02 /
824       DATA (TMP1(I, 26),I=1,4) /
825      +  0.340E+00, 0.100E-02, 0.630E-02, 0.630E-02 /
826       DATA (TMP1(I, 27),I=1,4) /
827      +  0.373E+00, 0.120E-02, 0.770E-02, 0.770E-02 /
828       DATA (TMP1(I, 28),I=1,4) /
829      +  0.406E+00, 0.140E-02, 0.860E-02, 0.930E-02 /
830       DATA (TMP1(I, 29),I=1,4) /
831      +  0.440E+00, 0.160E-02, 0.100E-01, 0.110E-01 /
832       DATA (TMP1(I, 30),I=1,4) /
833      +  0.474E+00, 0.180E-02, 0.110E-01, 0.120E-01 /
834       DATA (TMP1(I, 31),I=1,4) /
835      +  0.507E+00, 0.210E-02, 0.120E-01, 0.130E-01 /
836       DATA (TMP1(I, 32),I=1,4) /
837      +  0.535E+00, 0.240E-02, 0.130E-01, 0.150E-01 /
838       DATA (TMP1(I, 33),I=1,4) /
839      +  0.562E+00, 0.280E-02, 0.140E-01, 0.160E-01 /
840       DATA (TMP1(I, 34),I=1,4) /
841      +  0.589E+00, 0.320E-02, 0.160E-01, 0.180E-01 /
842       DATA (TMP1(I, 35),I=1,4) /
843      +  0.618E+00, 0.360E-02, 0.180E-01, 0.200E-01 /
844       DATA (TMP1(I, 36),I=1,4) /
845      +  0.643E+00, 0.410E-02, 0.200E-01, 0.220E-01 /
846       DATA (TMP1(I, 37),I=1,4) /
847      +  0.667E+00, 0.460E-02, 0.220E-01, 0.240E-01 /
848       DATA (TMP1(I, 38),I=1,4) /
849      +  0.690E+00, 0.510E-02, 0.240E-01, 0.260E-01 /
850       DATA (TMP1(I, 39),I=1,4) /
851      +  0.710E+00, 0.590E-02, 0.260E-01, 0.280E-01 /
852       DATA (TMP1(I, 40),I=1,4) /
853      +  0.730E+00, 0.680E-02, 0.280E-01, 0.310E-01 /
854       DATA (TMP1(I, 41),I=1,4) /
855      +  0.747E+00, 0.940E-02, 0.310E-01, 0.340E-01 /
856       DATA (TMP1(I, 42),I=1,4) /
857      +  0.765E+00, 0.100E-01, 0.340E-01, 0.370E-01 /
858       DATA (TMP1(I, 43),I=1,4) /
859      +  0.780E+00, 0.110E-01, 0.370E-01, 0.400E-01 /
860       DATA (TMP1(I, 44),I=1,4) /
861      +  0.794E+00, 0.120E-01, 0.400E-01, 0.430E-01 /
862       DATA (TMP1(I, 45),I=1,4) /
863      +  0.808E+00, 0.130E-01, 0.430E-01, 0.460E-01 /
864       DATA (TMP1(I, 46),I=1,4) /
865      +  0.820E+00, 0.140E-01, 0.470E-01, 0.490E-01 /
866       DATA (TMP1(I, 47),I=1,4) /
867      +  0.831E+00, 0.160E-01, 0.510E-01, 0.520E-01 /
868       DATA (TMP1(I, 48),I=1,4) /
869      +  0.843E+00, 0.180E-01, 0.560E-01, 0.560E-01 /
870       DATA (TMP1(I, 49),I=1,4) /
871      +  0.853E+00, 0.200E-01, 0.610E-01, 0.600E-01 /
872       DATA (TMP1(I, 50),I=1,4) /
873      +  0.862E+00, 0.370E-01, 0.650E-01, 0.640E-01 /
874       DATA (TMP1(I, 51),I=1,4) /
875      +  0.870E+00, 0.390E-01, 0.690E-01, 0.690E-01 /
876       DATA (TMP1(I, 52),I=1,4) /
877      +  0.877E+00, 0.410E-01, 0.740E-01, 0.740E-01 /
878       DATA (TMP1(I, 53),I=1,4) /
879      +  0.884E+00, 0.440E-01, 0.790E-01, 0.790E-01 /
880       DATA (TMP1(I, 54),I=1,4) /
881      +  0.891E+00, 0.460E-01, 0.830E-01, 0.850E-01 /
882       DATA (TMP1(I, 55),I=1,4) /
883      +  0.897E+00, 0.490E-01, 0.900E-01, 0.910E-01 /
884       DATA (TMP1(I, 56),I=1,4) /
885      +  0.902E+00, 0.520E-01, 0.960E-01, 0.970E-01 /
886       DATA (TMP1(I, 57),I=1,4) /
887      +  0.907E+00, 0.550E-01, 0.103E+00, 0.104E+00 /
888       DATA (TMP1(I, 58),I=1,4) /
889      +  0.912E+00, 0.580E-01, 0.110E+00, 0.111E+00 /
890       DATA (TMP1(I, 59),I=1,4) /
891      +  0.917E+00, 0.610E-01, 0.117E+00, 0.118E+00 /
892       DATA (TMP1(I, 60),I=1,4) /
893      +  0.921E+00, 0.640E-01, 0.124E+00, 0.125E+00 /
894       DATA (TMP1(I, 61),I=1,4) /
895      +  0.925E+00, 0.660E-01, 0.132E+00, 0.132E+00 /
896       DATA (TMP1(I, 62),I=1,4) /
897      +  0.929E+00, 0.710E-01, 0.140E+00, 0.139E+00 /
898       DATA (TMP1(I, 63),I=1,4) /
899      +  0.932E+00, 0.750E-01, 0.149E+00, 0.147E+00 /
900       DATA (TMP1(I, 64),I=1,4) /
901      +  0.935E+00, 0.790E-01, 0.158E+00, 0.155E+00 /
902       DATA (TMP1(I, 65),I=1,4) /
903      +  0.938E+00, 0.830E-01, 0.167E+00, 0.164E+00 /
904       DATA (TMP1(I, 66),I=1,4) /
905      +  0.941E+00, 0.890E-01, 0.178E+00, 0.174E+00 /
906       DATA (TMP1(I, 67),I=1,4) /
907      +  0.944E+00, 0.940E-01, 0.189E+00, 0.182E+00 /
908       DATA (TMP1(I, 68),I=1,4) /
909      +  0.947E+00, 0.100E+00, 0.200E+00, 0.192E+00 /
910       DATA (TMP1(I, 69),I=1,4) /
911      +  0.949E+00, 0.106E+00, 0.211E+00, 0.201E+00 /
912       DATA (TMP1(I, 70),I=1,4) /
913      +  0.951E+00, 0.112E+00, 0.222E+00, 0.210E+00 /
914       DATA (TMP1(I, 71),I=1,4) /
915      +  0.953E+00, 0.120E+00, 0.234E+00, 0.220E+00 /
916       DATA (TMP1(I, 72),I=1,4) /
917      +  0.955E+00, 0.128E+00, 0.246E+00, 0.231E+00 /
918       DATA (TMP1(I, 73),I=1,4) /
919      +  0.957E+00, 0.137E+00, 0.258E+00, 0.243E+00 /
920       DATA (TMP1(I, 74),I=1,4) /
921      +  0.958E+00, 0.147E+00, 0.270E+00, 0.255E+00 /
922       DATA (TMP1(I, 75),I=1,4) /
923      +  0.959E+00, 0.144E+00, 0.283E+00, 0.268E+00 /
924       DATA (TMP1(I, 76),I=1,4) /
925      +  0.961E+00, 0.130E+00, 0.295E+00, 0.281E+00 /
926       DATA (TMP1(I, 77),I=1,4) /
927      +  0.962E+00, 0.120E+00, 0.308E+00, 0.294E+00 /
928       DATA (TMP1(I, 78),I=1,4) /
929      +  0.963E+00, 0.114E+00, 0.321E+00, 0.306E+00 /
930       DATA (TMP1(I, 79),I=1,4) /
931      +  0.964E+00, 0.107E+00, 0.334E+00, 0.320E+00 /
932       DATA (TMP1(I, 80),I=1,4) /
933      +  0.965E+00, 0.107E+00, 0.347E+00, 0.333E+00 /
934       DATA (TMP1(I, 81),I=1,4) /
935      +  0.966E+00, 0.107E+00, 0.360E+00, 0.347E+00 /
936       DATA (TMP1(I, 82),I=1,4) /
937      +  0.967E+00, 0.112E+00, 0.373E+00, 0.360E+00 /
938       DATA (TMP1(I, 83),I=1,4) /
939      +  0.968E+00, 0.117E+00, 0.387E+00, 0.373E+00 /
940       DATA (TMP1(I, 84),I=1,4) /
941      +  0.968E+00, 0.122E+00, 0.401E+00, 0.386E+00 /
942       DATA (TMP1(I, 85),I=1,4) /
943      +  0.969E+00, 0.128E+00, 0.415E+00, 0.399E+00 /
944       DATA (TMP1(I, 86),I=1,4) /
945      +  0.969E+00, 0.134E+00, 0.429E+00, 0.411E+00 /
946       DATA (TMP1(I, 87),I=1,4) /
947      +  0.970E+00, 0.139E+00, 0.443E+00, 0.424E+00 /
948       DATA (TMP1(I, 88),I=1,4) /
949      +  0.970E+00, 0.146E+00, 0.456E+00, 0.437E+00 /
950       DATA (TMP1(I, 89),I=1,4) /
951      +  0.971E+00, 0.153E+00, 0.468E+00, 0.450E+00 /
952       DATA (TMP1(I, 90),I=1,4) /
953      +  0.971E+00, 0.161E+00, 0.479E+00, 0.463E+00 /
954       DATA (TMP1(I, 91),I=1,4) /
955      +  0.972E+00, 0.162E+00, 0.472E+00, 0.476E+00 /
956       DATA (TMP1(I, 92),I=1,4) /
957      +  0.972E+00, 0.176E+00, 0.467E+00, 0.489E+00 /
958       DATA (TMP1(I, 93),I=1,4) /
959      +  0.973E+00, 0.187E+00, 0.466E+00, 0.502E+00 /
960       DATA (TMP1(I, 94),I=1,4) /
961      +  0.973E+00, 0.205E+00, 0.464E+00, 0.514E+00 /
962       DATA (TMP1(I, 95),I=1,4) /
963      +  0.974E+00, 0.218E+00, 0.471E+00, 0.526E+00 /
964       DATA (TMP1(I, 96),I=1,4) /
965      +  0.974E+00, 0.228E+00, 0.479E+00, 0.539E+00 /
966       DATA (TMP1(I, 97),I=1,4) /
967      +  0.975E+00, 0.236E+00, 0.485E+00, 0.550E+00 /
968       DATA (TMP1(I, 98),I=1,4) /
969      +  0.975E+00, 0.244E+00, 0.490E+00, 0.560E+00 /
970       DATA (TMP1(I, 99),I=1,4) /
971      +  0.975E+00, 0.253E+00, 0.497E+00, 0.570E+00 /
972       DATA (TMP1(I,100),I=1,4) /
973      +  0.976E+00, 0.263E+00, 0.506E+00, 0.579E+00 /
974 C
975       DO 20 J = 1,MAXELZ
976          DO 10 I = 1,ISHLS
977             GPHRAT(I,J) = TMP(I,J)
978    10    CONTINUE
979    20 CONTINUE
980 C
981       DO 30 I = 1,29
982          ISHLTR(I) = ITRN(I)
983    30 CONTINUE
984 C
985       DO 40 I = 1,MAXELZ
986          GFLUPR(1,I) = TMP1(1,I)
987          GFLUPR(2,I) = TMP1(2,I)
988          GFLUPR(3,I) = TMP1(3,I)
989          GFLUPR(4,I) = TMP1(4,I)
990    40 CONTINUE
991 C
992 C K shell decay mode
993       DO 50 I = 5,MAXELZ
994          IF(I.LT.8) THEN
995             IGNRFN(1,I) = 2
996             IGNRFN(2,I) = 2
997          ELSEIF(I.LT.81) THEN
998             IGNRFN(1,I) = 3
999             IGNRFN(2,I) = 4
1000          ELSE
1001             IGNRFN(1,I) = 3
1002             IGNRFN(2,I) = 3
1003          ENDIF
1004    50 CONTINUE
1005 C L1 shell decay mode
1006       DO 60 I = 11,MAXELZ
1007          IGNRFN(3,I) = 4
1008          IF(I.LT.24) THEN
1009             IGNRFN(4,I) = 5
1010          ELSEIF(I.LT.49) THEN
1011             IGNRFN(4,I) = 9
1012          ELSEIF(I.LT.53) THEN
1013             IGNRFN(4,I) = 10
1014          ELSE
1015             IGNRFN(4,I) = 14
1016          ENDIF
1017    60 CONTINUE
1018 C L2 shell decay mode
1019       DO 70 I = 12,MAXELZ
1020          IF(I.LT.13) THEN
1021             IGNRFN(5,I) = 5
1022             IGNRFN(6,I) = 5
1023          ELSEIF(I.LT.16) THEN
1024             IGNRFN(5,I) = 5
1025             IGNRFN(6,I) = 6
1026          ELSEIF(I.LT.29) THEN
1027             IGNRFN(5,I) = 6
1028             IGNRFN(6,I) = 7
1029          ELSE
1030             IGNRFN(5,I) = 8
1031             IGNRFN(6,I) = 9
1032          ENDIF
1033    70 CONTINUE
1034 C L3 shell decay mode
1035       DO 80 I = 12,MAXELZ
1036          IF(I.LT.13) THEN
1037             IGNRFN(7,I) = 5
1038             IGNRFN(8,I) = 5
1039          ELSEIF(I.LT.16) THEN
1040             IGNRFN(7,I) = 5
1041             IGNRFN(8,I) = 7
1042          ELSEIF(I.LT.30) THEN
1043             IGNRFN(7,I) = 7
1044             IGNRFN(8,I) = 7
1045          ELSE
1046             IGNRFN(7,I) = 8
1047             IGNRFN(8,I) = 9
1048          ENDIF
1049    80 CONTINUE
1050       END