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