C 19 SEP 05 - MWS - ADD TRUE NUCLEAR CHARGE ARRAY TO INFOA COMMON
C  5 JUL 05 - MWS - SELECT NEW ATOM,BASIS,EFP,PCM,DAF DIMENSIONS
C 19 NOV 00 - MWS - BASHUZ: RETURN NORMALIZED CXINP COEFS
C 10 AUG 94 - MWS - INCREASE NUMBER OF DAF RECORDS
C 18 MAY 92 - MWS - BASHUZ: INITIALIZE ALL CONTRACTION COEFS
C 12 MAR 92 - MWS - REDIMENSION TO 500 ATOMS
C  4 JAN 92 - TLW - MAKE WRITES PARALLEL;ADD COMMON PAR
C  2 JUL 91 - MWS - PASS NSHELL,KSTART,ETC AS CALLING ARGUMENTS.
C 12 SEP 90 - MWS - INTRODUCE MXATM
C  7 AUG 90 - TLW - ADD CF AND CG TO COMMON NSHEL
C 10 MAY 90 - MWS - CLEAN UP SCALE FACTOR USAGE
C 10 AUG 88 - MWS - MXSH,MXGSH,MXGTOT FROM 120,10,440 TO 1000,30,5000
C 30 MAY 88 - MWS - USE PARAMETERS TO DIMENSION COMMON
C  7 JUL 86 - JAB - SANTIZE FLOATING POINT CONSTANTS
C 14 OCT 85 - STE - USE GENERIC SQRT,ABS
C 30 AUG 85 - MWS - COMPLETE TYPING OF TRANSITION METALS
C  9 AUG 85 - MWS - PROVIDE DEFAULT SCALE FACTORS FOR MAIN GROUP
C 13 JUN 85 - MWS - ALLOW INPUT OF SCALE FACTORS
C 10 JUN 85 - MWS - NEW SECTION, HUZINAGA MINI/MIDI BASES
C*MODULE BASHUZ  *DECK BASHUZ
      SUBROUTINE BASHUZ(NUCZ,MIDI,CSINP,CPINP,CDINP,CFINP,
     *                  SCFAC,IERR1,IERR2,INTYP,NANGM,NBFS,MINF,MAXF,
     *                  LOC,NGAUSS,NS,EX,CS,CP,CD,CF,KSTART,KATOM,KTYPE,
     *                  KNG,KLOC,KMIN,KMAX,NSHELL,MXGTOT,MXSH)
C
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
C
      LOGICAL MIDI,DONE,GOPARR,DSKWRK,MASWRK
C
      DIMENSION CSINP(MXGTOT),CPINP(MXGTOT),CDINP(MXGTOT),CFINP(MXGTOT),
     *          INTYP(*),NANGM(*),NBFS(*),MINF(*),MAXF(*),NS(*),
     *          EX(MXGTOT),CS(MXGTOT),CP(MXGTOT),CD(MXGTOT),CF(MXGTOT),
     *          KSTART(MXSH),KATOM(MXSH),KTYPE(MXSH),KNG(MXSH),
     *          KLOC(MXSH),KMIN(MXSH),KMAX(MXSH),SCFAC(4)
      DIMENSION EEX(45),COEF(45)
C
      PARAMETER (MXATM=2000)
C
      COMMON /IOFILE/ IR,IW,IP,IS,IPK,IDAF,NAV,IODA(950)
      COMMON /INFOA / NAT,ICH,MUL,NUM,NQMT,NE,NA,NB,
     *                ZAN(MXATM),C(3,MXATM),IAN(MXATM)
      COMMON /OUTPUT/ NPRINT,ITOL,ICUT,NORMF,NORMP,NOPK
      COMMON /PAR   / ME,MASTER,NPROC,IBTYP,IPTIM,GOPARR,DSKWRK,MASWRK
C
      DATA ZERO,ONE,PT5,PT75,PT1875/0.0D+00,1.0D+00,0.5D+00,0.75D+00,
     *            1.875D+00/
      DATA PI32,TM6,TM10/5.56832799683170D+00,1.0D-06,1.0D-10/
C
C     ----- GET THE HUZINAGA MINI OR MIDI BASIS SET -----
C     THESE BASIS SETS ARE TAKEN FROM "THE GREEN BOOK",
C     "GAUSSIAN BASIS SETS FOR MOLECULAR CALCULATIONS"
C     EDITED BY S. HUZINAGA, WITH COAUTHORS J. ANDZELM,
C     M. KLOBUKOWSKI, E.RADZIO-ANDZELM, Y.SAKAI, H.TATEWAKI.
C     ELSEVIER, AMSTERDAM, 1984.
C
C     THE MINI BASIS SET IS A MINIMAL BASIS SET, 3 GAUSSIANS
C     PER ATOMIC ORBITAL (THE SAME AS STO-3G). HOWEVER THE EXPONENTS
C     ARE OPTIMAL FOR EACH ATOM, AND S AND P EXPONENTS ARE DIFFERENT.
C     MIDI BASIS SETS ARE DERIVED FROM THE MINI BY FLOATING THE
C     OUTERMOST PRIMITIVE OF THE VALENCE ORBITALS, AND THEN
C     RENORMALIZING THE REMAINING 2 GAUSSIANS.
C
      SCALE = ZERO
      NG = -2**20
      IGAUSS = NG
      ITYP = NG
      IERR3=0
      DONE = .FALSE.
      DO 100 I = 1,45
         EEX(I) = ZERO
         COEF(I) = ZERO
  100 CONTINUE
C
C     ----- HYDROGEN TO HELIUM -----
C
      IF (NUCZ .GT. 2) GO TO 120
      CALL HONE(EEX,COEF,NUCZ)
      GO TO 200
C
C     ----- LITHIUM TO NEON -----
C
  120 IF (NUCZ .GT. 10) GO TO 130
      CALL HTWO(EEX,COEF,NUCZ)
      GO TO 200
C
C     ----- SODIUM TO ARGON -----
C
  130 IF (NUCZ .GT. 18) GO TO 140
      CALL HTHREE(EEX,COEF,NUCZ)
      GO TO 200
C
C     ----- POTASSIUM TO KRYPTON -----
C
  140 IF(NUCZ.GT.36) GO TO 150
      CALL HFOUR(EEX,COEF,NUCZ)
      GO TO 200
C
C     ----- RUBIDIUM TO XENON -----
C
  150 IF(NUCZ.GT.54) GO TO 160
      CALL HFIVE(EEX,COEF,NUCZ)
      GO TO 200
C
C     ----- CESIUM TO RADON -----
C
  160 IF(NUCZ.GT.86) GO TO 190
      CALL HSIX(EEX,COEF,NUCZ)
      GO TO 200
C
C     ----- PAST RADON DOES NOT EXIST -----
C
  190 CONTINUE
      CALL BERROR(3)
      STOP
C
C     ----- LOOP OVER EACH SHELL -----
C
  200 CONTINUE
      IPASS = 0
  210 IPASS = IPASS+1
      SCS = SCFAC(1)
      SCP = SCFAC(2)
      SCD = SCFAC(3)
      CALL HSHELL(IW,NUCZ,IPASS,MIDI,ITYP,IGAUSS,NG,SCALE,
     *            SCS,SCP,SCD,DONE)
      IF(DONE. AND. IERR3.NE.0) CALL ABRT
      IF(DONE) RETURN
C
C     ----- DEFINE THE CURRENT SHELL -----
C
      NSHELL = NSHELL+1
      IF(NSHELL.GT.MXSH) THEN
         IERR1=1
         RETURN
      END IF
      NS(NAT) = NS(NAT)+1
      KMIN(NSHELL) = MINF(ITYP)
      KMAX(NSHELL) = MAXF(ITYP)
      KSTART(NSHELL) = NGAUSS+1
      KATOM(NSHELL) = NAT
      KTYPE(NSHELL) = NANGM(ITYP)
      INTYP(NSHELL) = ITYP
      KNG(NSHELL) = IGAUSS
      KLOC(NSHELL) = LOC+1
      NGAUSS = NGAUSS+IGAUSS
      IF(NGAUSS.GT.MXGTOT) THEN
         IERR2=1
         RETURN
      END IF
      LOC = LOC+NBFS(ITYP)
      K1 = KSTART(NSHELL)
      K2 = K1+KNG(NSHELL)-1
      SCALE = SCALE*SCALE
      DO 440 I = 1,IGAUSS
         K = K1+I-1
         EX(K) = EEX(NG+I) * SCALE
         IF(ITYP.EQ.1) CSINP(K) = COEF(NG+I)
         IF(ITYP.EQ.2) CPINP(K) = COEF(NG+I)
         IF(ITYP.EQ.3) CDINP(K) = COEF(NG+I)
         IF(ITYP.EQ.4) CFINP(K) = COEF(NG+I)
         CS(K) = ZERO
         CP(K) = ZERO
         CD(K) = ZERO
         CF(K) = ZERO
         IF(ITYP.EQ.1) CS(K) = CSINP(K)
         IF(ITYP.EQ.2) CP(K) = CPINP(K)
         IF(ITYP.EQ.3) CD(K) = CDINP(K)
         IF(ITYP.EQ.4) CF(K) = CFINP(K)
  440 CONTINUE
C
C     ----- ALWAYS UNNORMALIZE PRIMITIVES -----
C
      DO 460 K = K1,K2
         EE = EX(K)+EX(K)
         FACS = PI32/(EE*SQRT(EE))
         FACP = PT5*FACS/EE
         FACD = PT75*FACS/(EE*EE)
         FACF = PT1875*FACS/(EE**3)
         IF(ITYP.EQ.1) CS(K) = CS(K)/SQRT(FACS)
         IF(ITYP.EQ.2) CP(K) = CP(K)/SQRT(FACP)
         IF(ITYP.EQ.3) CD(K) = CD(K)/SQRT(FACD)
         IF(ITYP.EQ.4) CF(K) = CF(K)/SQRT(FACF)
  460 CONTINUE
C
C     ----- IF(NORMF.EQ.0) NORMALIZE BASIS FUNCTIONS -----
C
      IF (NORMF .EQ. 1) GO TO 210
      FACS = ZERO
      FACP = ZERO
      FACD = ZERO
      FACF = ZERO
      DO 510 IG = K1,K2
         DO 500 JG = K1,IG
            EE = EX(IG)+EX(JG)
            FAC = EE*SQRT(EE)
            DUMS = CS(IG)*CS(JG)/FAC
            DUMP = PT5*CP(IG)*CP(JG)/(EE*FAC)
            DUMD = PT75*CD(IG)*CD(JG)/(EE*EE*FAC)
            DUMF = PT1875*CF(IG)*CF(JG)/(EE*EE*EE*FAC)
            IF (IG .EQ. JG) GO TO 480
               DUMS = DUMS+DUMS
               DUMP = DUMP+DUMP
               DUMD = DUMD+DUMD
               DUMF = DUMF+DUMF
  480       CONTINUE
            FACS = FACS+DUMS
            FACP = FACP+DUMP
            FACD = FACD+DUMD
            FACF = FACF+DUMF
  500    CONTINUE
  510 CONTINUE
C
      FAC=ZERO
      IF(ITYP.EQ.1 .AND. FACS.GT.TM10) FAC=ONE/SQRT(FACS*PI32)
      IF(ITYP.EQ.2 .AND. FACP.GT.TM10) FAC=ONE/SQRT(FACP*PI32)
      IF(ITYP.EQ.3 .AND. FACD.GT.TM10) FAC=ONE/SQRT(FACD*PI32)
      IF(ITYP.EQ.4 .AND. FACF.GT.TM10) FAC=ONE/SQRT(FACF*PI32)
C                                 VERIFY NORMALIZATION
      TNORM = ABS(ONE-FAC)
      IF(TNORM.LT.TM6) GO TO 520
      IF(MIDI  .AND.  IGAUSS.LT.3) GO TO 520
      IF(MIDI  .AND.  ITYP.EQ.4  .AND.  IGAUSS.EQ.3) GO TO 520
C
      IF (MASWRK) THEN
         WRITE(IW,9000) FAC,NUCZ,IPASS,ITYP
         WRITE(IW,9010) (EEX(NG+K),K=1,IGAUSS)
         WRITE(IW,9010) (COEF(NG+K),K=1,IGAUSS)
      ENDIF
      IERR3=IERR3+1
C
  520 CONTINUE
      DO 550 IG = K1,K2
         IF(ITYP.EQ.1) CS(IG) = FAC*CS(IG)
         IF(ITYP.EQ.2) CP(IG) = FAC*CP(IG)
         IF(ITYP.EQ.3) CD(IG) = FAC*CD(IG)
         IF(ITYP.EQ.4) CF(IG) = FAC*CF(IG)
         IF(ITYP.EQ.1) CSINP(IG) = FAC*CSINP(IG)
         IF(ITYP.EQ.2) CPINP(IG) = FAC*CPINP(IG)
         IF(ITYP.EQ.3) CDINP(IG) = FAC*CDINP(IG)
         IF(ITYP.EQ.4) CFINP(IG) = FAC*CFINP(IG)
  550 CONTINUE
      GO TO 210
C
 9000 FORMAT(/1X,'ERROR!!! NORMALIZATION FACTOR=',E16.8/
     *     1X,'FOR ATOM Z=',I3,' SHELL NO.',I3,' ITYP=',I4/
     *     1X,'(ITYP OF 1,2,3,4 MEANS AN S,P,D,F SHELL)'/
     *     1X,'CHECK BUILT IN HUZINAGA EXPONENTS AND CONT. COEFS')
 9010 FORMAT(3F16.10)
      END
C*MODULE BASHUZ  *DECK HSCALE
      DOUBLE PRECISION FUNCTION HSCALE(NUCZ)
      IMPLICIT DOUBLE PRECISION(A-H,O-Z)
      DIMENSION ROW1(2),ROW2(8),ROW3(8),ROW4(8),ROW5(8)
      DATA ROW1/1.25D+00,1.00D+00/
      DATA ROW2/1.091D+00,1.105D+00,1.118D+00,1.075D+00,
     *          1.035D+00,1.014D+00,1.007D+00,1.000D+00/
      DATA ROW3/1.104D+00,1.091D+00,1.126D+00,1.104D+00,
     *          1.036D+00,1.014D+00,1.004D+00,1.000D+00/
      DATA ROW4/1.131D+00,1.086D+00,1.158D+00,1.116D+00,
     *          1.040D+00,1.017D+00,1.006D+00,1.000D+00/
      DATA ROW5/1.134D+00,1.078D+00,1.200D+00,1.131D+00,
     *          1.053D+00,1.024D+00,1.011D+00,1.000D+00/
      DATA ONE/1.0D+00/
C
C     PROVIDE DEFAULT SCALE FACTORS FOR HUZINAGA MINI BASIS SETS.
C     THE MAIN GROUP SCALE FACTORS ARE DERIVED FROM RHF
C     GEOMETRY OPTIMIZATIONS ON SAMPLE MOLECULES, SUCH AS
C     HYDRIDES, CARRIED OUT AT NORTH DAKOTA STATE UNIVERSITY
C     IN THE SUMMER OF 1985 BY JOHN DEISZ.
C     S AND P SCALE FACTORS ARE THE SAME, WHICH DOES NOT SEEM
C     TO BE A SEVERE CONSTAINT. FOR A GIVEN ELEMENT, THE S SCALE
C     IS SOMETIMES SMALLER THAN THE P IN ONE MOLECULE, LARGER IN
C     ANOTHER, AND ABOUT EQUAL IN A THIRD.
C     OTHER ELEMENTS (TRANSITION METALS, Z>54) ARE UNSCALED.
C
      IF(NUCZ.GT.2) GO TO 020
         HSCALE=ROW1(NUCZ)
         RETURN
  020 IF(NUCZ.GT.10) GO TO 030
         NUC=NUCZ-2
         HSCALE=ROW2(NUC)
         RETURN
  030 IF(NUCZ.GT.18) GO TO 040
         NUC=NUCZ-10
         HSCALE=ROW3(NUC)
         RETURN
  040 IF(NUCZ.GT.20) GO TO 050
         NUC=NUCZ-18
         HSCALE=ROW4(NUC)
         RETURN
  050 IF(NUCZ.GT.30) GO TO 060
         HSCALE=ONE
         RETURN
  060 IF(NUCZ.GT.36) GO TO 070
         NUC=NUCZ-28
         HSCALE=ROW4(NUC)
         RETURN
  070 IF(NUCZ.GT.38) GO TO 080
         NUC=NUCZ-36
         HSCALE=ROW5(NUC)
         RETURN
  080 IF(NUCZ.GT.48) GO TO 090
         HSCALE=ONE
         RETURN
  090 IF(NUCZ.GT.54) GO TO 100
         NUC=NUCZ-46
         HSCALE=ROW5(NUC)
         RETURN
  100 CONTINUE
      HSCALE=ONE
      RETURN
      END
C*MODULE BASHUZ  *DECK HSHELL
      SUBROUTINE HSHELL(IW,NUCZ,IPASS,MIDI,ITYP,IGAUSS,NG,
     *                  SCALE,SCS,SCP,SCD,DONE)
      IMPLICIT DOUBLE PRECISION(A-H,O-Z)
      LOGICAL MIDI,DONE,GOPARR,DSKWRK,MASWRK
      DIMENSION IOFF(15,15),JOFF(17,15),KTYP(15,3)
C
      COMMON /PAR   / ME,MASTER,NPROC,IBTYP,IPTIM,GOPARR,DSKWRK,MASWRK
C
      DATA KTYP/1, 1,2, 1,2,3, 1,2,3,   1,2,3,  1,2*0,
     *          1, 1,2, 1,2,3, 1,2,3,   1,2, 4, 1,2*0,
     *          1, 1,2, 1,2,3, 1,2,3,4, 1,2,3,  1,2/
C                         POINTERS TO THE MINI BASIS
      DATA IOFF/0,  14*0,
     2          0, 3,  13*0,
     *          0, 3, 6,  12*0,
     3          0, 3, 9, 6,  11*0,
     *          0, 3, 9, 6,12,  10*0,
     4          0, 3,12, 6,15, 9,  9*0,
     *          0, 3,12, 6,15,18, 9,  8*0,
     *          0, 3,12, 6,15,21, 9,18,  7*0,
     5          0, 3,15, 6,18,24, 9,21,12,  6*0,
     *          0, 3,15, 6,18,24, 9,21,27,12,  5*0,
     *          0, 3,15, 6,18,27, 9,21,30,12,24,  4*0,
     6          0, 3, 6, 9,12,15,18,21,24,27,30,33,  3*0,
     *          0, 3, 6, 9,12,15,18,21,24,27,30,33,37, 2*0,
     *          0, 3, 6, 9,12,15,18,21,24,27,30,33,36,39, 0,
     *          0, 3, 6, 9,12,15,18,21,24,27,30,33,36,39,42/
C                         POINTERS TO THE MIDI BASIS
      DATA JOFF/0, 2, 15*0,
     2          0, 3, 5,  14*0,
     *          0, 3, 5, 6, 8,  12*0,
     3          0, 3, 9, 6, 8,  12*0,
     *          0, 3, 9, 6, 8,12,14,  10*0,
     4          0, 3,12, 6,15, 9,11,  10*0,
     *          0, 3,12, 6,15,18,20, 9,11,  8*0,
     *          0, 3,12, 6,15,21, 9,11,18,20,  7*0,
     5          0, 3,15, 6,18,24, 9,21,12,14,  7*0,
     *          0, 3,15, 6,18,24, 9,21,27,29,12,14,  5*0,
     *          0, 3,15, 6,18,27, 9,21,30,12,14,24,26,  4*0,
     6          0, 3, 6, 9,12,15,18,21,24,27,30,33,35,  4*0,
     *          0, 3, 6, 9,12,15,18,21,24,27,30,33,36,37,39,  2*0,
     *          0, 3, 6, 9,12,15,18,21,24,27,30,33,36,38,39,41, 0,
     *          0, 3, 6, 9,12,15,18,21,24,27,30,33,36,39,41,42,44/
      DATA ZERO,ONE/0.0D+00,1.0D+00/
C
C     SET VALUES FOR THE CURRENT HUZINAGA MINI/MIDI SHELL
C     RETURN IGAUSS = NUMBER OF GAUSSIANS IN CURRENT SHELL
C            ITYP   = 1,2,3,4 FOR S,P,D,F SHELL
C            NG     = OFFSET IN EEX,COEF ARRAYS FOR CURRENT SHELL
C            SCALE  = SCALING FACTOR FOR CURRENT SHELL
C
      KIND=1
      IF(NUCZ.GT. 2) KIND=2
      IF(NUCZ.GT. 4) KIND=3
      IF(NUCZ.GT.10) KIND=4
      IF(NUCZ.GT.12) KIND=5
      IF(NUCZ.GT.18) KIND=6
      IF(NUCZ.GT.20) KIND=7
      IF(NUCZ.GT.30) KIND=8
      IF(NUCZ.GT.36) KIND=9
      IF(NUCZ.GT.38) KIND=10
      IF(NUCZ.GT.48) KIND=11
      IF(NUCZ.GT.54) KIND=12
      IF(NUCZ.GT.56) KIND=13
      IF(NUCZ.GT.70) KIND=14
      IF(NUCZ.GT.80) KIND=15
      IF(NUCZ.GT.86) THEN
         IF (MASWRK) WRITE(IW,*)
     *      'HUZINAGA BASIS SETS ONLY GO TO RADON.'
         CALL ABRT
      END IF
C
      KIND2 = 1
      IF(NUCZ.GT.56) KIND2=2
      IF(NUCZ.GT.70) KIND2=3
C
      MXPASS=KIND
C                 ALL ATOMS HAVE TWO VALENCE SHELLS, EXCEPT ALKALIS
      NSPLIT=2
      IF(KIND.LE.2  .OR.  KIND.EQ.4  .OR.  KIND.EQ.6  .OR.
     *   KIND.EQ.9  .OR.  KIND.EQ.12) NSPLIT=1
      IF(MIDI) MXPASS=MXPASS+NSPLIT
C
      IF(IPASS.GT.MXPASS) DONE=.TRUE.
      IF(DONE) RETURN
C              ******
C
      IGAUSS = 3
      IF(MIDI) THEN
         MI = MXPASS-IPASS
         IF(MI.LT.2*NSPLIT) THEN
            IF(MOD(MI,2).EQ.1) IGAUSS=2
            IF(MOD(MI,2).EQ.0) IGAUSS=1
            IF(MI.GE.2) THEN
               ITYP=KTYP(KIND-1,KIND2)
            ELSE
               ITYP=KTYP(KIND,KIND2)
            END IF
         ELSE
            ITYP=KTYP(IPASS,KIND2)
         END IF
         NG  =JOFF(IPASS,KIND)
      ELSE
         ITYP=KTYP(IPASS,KIND2)
         NG  =IOFF(IPASS,KIND)
      END IF
C                LANTHANIDES HAVE 4 (OR 3,1) GAUSSIANS IN VALENCE F
      IF(KIND2.EQ.2  .AND.  ITYP.EQ.4  .AND.  IGAUSS.NE.1)
     *   IGAUSS=IGAUSS+1
C                ALKALI'S DON'T YET FILL THE PREVIOUS D LEVEL
      IF(KIND.EQ.6  .AND.  IPASS.GE.KIND) ITYP=1
      IF(KIND.EQ.9  .AND.  IPASS.GE.KIND) ITYP=1
      IF(KIND.EQ.12 .AND.  IPASS.GE.KIND) ITYP=1
C
C          DO NOT SCALE MIDI, OR MINI INNER SHELLS
C
      SCALE=ONE
      IF(MIDI) RETURN
      IF(IPASS.LE.(KIND-NSPLIT)) RETURN
      IF(ITYP.EQ.1) SCALE=SCS
      IF(ITYP.EQ.2) SCALE=SCP
      IF(ITYP.EQ.3) SCALE=SCD
      IF(SCALE.EQ.ZERO) SCALE=HSCALE(NUCZ)
      RETURN
      END
C*MODULE BASHUZ  *DECK HONE
      SUBROUTINE HONE(E,C,NUCZ)
      IMPLICIT DOUBLE PRECISION(A-H,O-Z)
      DIMENSION E(3),C(3)
C
C     FIRST ROW BASIS SETS TAKEN FROM F. B. VAN DUIJNEVELDT,
C     IBM TECHNICAL REPORT RJ-945.
C
      IF(NUCZ.EQ.2) GO TO 200
C
C           HYDROGEN
C
      E(1)=4.501800D+00
      E(2)=0.681444D+00
      E(3)=0.151398D+00
      C(1)=0.070452D+00
      C(2)=0.407826D+00
      C(3)=0.647752D+00
      RETURN
C
C           HELIUM
C
  200 CONTINUE
      E(1)=13.626736D+00
      E(2)=1.999349D+00
      E(3)=0.382993D+00
      C(1)=0.080241D+00
      C(2)=0.409143D+00
      C(3)=0.657278D+00
      RETURN
      END
C*MODULE BASHUZ  *DECK HTWO
      SUBROUTINE HTWO(E,C,NUCZ)
      IMPLICIT DOUBLE PRECISION(A-H,O-Z)
      DIMENSION E(9),C(9)
C
      IBR=NUCZ-2
      GO TO (100,200,300,400,500,600,700,800), IBR
C
C           LITHIUM
C
  100 CONTINUE
      E(1)=35.04615D+00
      E(2)=5.20169D+00
      E(3)=1.05624D+00
      C(1)=0.0737603D+00
      C(2)=0.3974715D+00
      C(3)=0.6650924D+00
      E(4)=0.71517D+00
      E(5)=0.07053D+00
      E(6)=0.02735D+00
      C(4)=-0.09397D+00
      C(5)=0.5701001D+00
      C(6)=0.4997501D+00
      RETURN
C
C           BERYLLIUM
C
  200 CONTINUE
      E(1)=66.953540D+00
      E(2)=9.9392900D+00
      E(3)=2.05713D+00
      C(1)=0.0702000D+00
      C(2)=0.3919101D+00
      C(3)=0.6699701D+00
      E(4)=1.9122100D+00
      E(5)=0.16132D+00
      E(6)=0.05524D+00
      C(4)=-0.0828204D+00
      C(5)=0.5575528D+00
      C(6)=0.5160426D+00
      RETURN
C
C           BORON
C
  300 CONTINUE
      E(1)=108.43704D+00
      E(2)=16.120560D+00
      E(3)=3.3734300D+00
      C(1)=0.0686505D+00
      C(2)=0.3899329D+00
      C(3)=0.6713951D+00
      E(4)=3.5665D+00
      E(5)=0.29547D+00
      E(6)=0.09805D+00
      C(4)=-0.0824192D+00
      C(5)=0.5590644D+00
      C(6)=0.5167948D+00
      E(7)=2.57207D+00
      E(8)=0.51694D+00
      E(9)=0.12314D+00
      C(7)=0.1059001D+00
      C(8)=0.4571805D+00
      C(9)=0.6318607D+00
      RETURN
C
C           CARBON
C
  400 CONTINUE
      E(1)=153.17226D+00
      E(2)=23.073030D+00
      E(3)=4.9232900D+00
      C(1)=0.0707401D+00
      C(2)=0.3953804D+00
      C(3)=0.6633107D+00
      E(4)=5.7255700D+00
      E(5)=0.45504000D+00
      E(6)=0.14707D+00
      C(4)=-0.0813805D+00
      C(5)=0.5748532D+00
      C(6)=0.5024128D+00
      E(7)=4.25131D+00
      E(8)=0.86327D+00
      E(9)=0.20135D+00
      C(7)=0.1099306D+00
      C(8)=0.4627127D+00
      C(9)=0.6275137D+00
      RETURN
C
C           NITROGEN
C
  500 CONTINUE
      E(1)=218.36449D+00
      E(2)=32.59889D+00
      E(3)=6.91739D+00
      C(1)=0.0678703D+00
      C(2)=0.3902018D+00
      C(3)=0.6700832D+00
      E(4)=8.32638D+00
      E(5)=0.65919D+00
      E(6)=0.21009D+00
      C(4)=-0.0808903D+00
      C(5)=0.5672020D+00
      C(6)=0.5110918D+00
      E(7)=6.12035D+00
      E(8)=1.25938D+00
      E(9)=0.29145D+00
      C(7)=0.1159195D+00
      C(8)=0.4699582D+00
      C(9)=0.6184476D+00
      RETURN
C
C           OXYGEN
C
  600 CONTINUE
      E(1)=281.86658D+00
      E(2)=42.416D+00
      E(3)=9.09562D+00
      C(1)=0.0690599D+00
      C(2)=0.3931595D+00
      C(3)=0.6656691D+00
      E(4)=11.46603D+00
      E(5)=0.88786D+00
      E(6)=0.2788D+00
      C(4)=-0.0808199D+00
      C(5)=0.5820895D+00
      C(6)=0.4971596D+00
      E(7)=8.04724D+00
      E(8)=1.66842D+00
      E(9)=0.37251D+00
      C(7)=0.1242709D+00
      C(8)=0.4765935D+00
      C(9)=0.6130445D+00
      RETURN
C
C           FLUORINE
C
  700 CONTINUE
      E(1)=368.37112D+00
      E(2)=55.06106D+00
      E(3)=11.74767D+00
      C(1)=0.0670398D+00
      C(2)=0.3892491D+00
      C(3)=0.6707884D+00
      E(4)=15.15184D+00
      E(5)=1.15137D+00
      E(6)=0.35811D+00
      C(4)=-0.0805499D+00
      C(5)=0.587729D+00
      C(6)=0.4919792D+00
      E(7)=10.57707D+00
      E(8)=2.19498D+00
      E(9)=0.47937D+00
      C(7)=0.1262695D+00
      C(8)=0.4779482D+00
      C(9)=0.6140077D+00
      RETURN
C
C           NEON
C
  800 CONTINUE
      E(1)=456.95285D+00
      E(2)=68.365430D+00
      E(3)=14.619760D+00
      C(1)=0.0669098D+00
      C(2)=0.3893491D+00
      C(3)=0.6705184D+00
      E(4)=19.327190D+00
      E(5)=1.4418200D+00
      E(6)=0.44408000D+00
      C(4)=-0.0802497D+00
      C(5)=0.5952977D+00
      C(6)=0.4848681D+00
      E(7)=13.352520D+00
      E(8)=2.77947D+00
      E(9)=0.60097D+00
      C(7)=0.1288403D+00
      C(8)=0.4804412D+00
      C(9)=0.6116716D+00
      RETURN
      END
C*MODULE BASHUZ  *DECK HTHREE
      SUBROUTINE HTHREE(E,C,NUCZ)
      IMPLICIT DOUBLE PRECISION(A-H,O-Z)
      DIMENSION E(15),C(15)
C
      IBR=NUCZ-10
      GO TO (100,200,300,400,500,600,700,800), IBR
C
C           SODIUM
C
  100 CONTINUE
      E(1)=542.76053D+00
      E(2)=81.959470D+00
      E(3)=17.723770D+00
      C(1)=0.0684104D+00
      C(2)=0.3920923D+00
      C(3)=0.6660838D+00
      E(4)=23.280420D+00
      E(5)=1.8683400D+00
      E(6)=0.6232500D+00
      C(4)=-0.0838006D+00
      C(5)=0.5827943D+00
      C(6)=0.4924737D+00
      E(7)=0.5067300D+00
      E(8)=0.0535100D+00
      E(9)=0.0208000D+00
      C(7)=-0.1157622D+00
      C(8)=0.6958634D+00
      C(9)=0.3810473D+00
      E(10)=17.836360D+00
      E(11)=3.7956900D+00
      E(12)=0.8775100D+00
      C(10)=0.1257103D+00
      C(11)=0.4804611D+00
      C(12)=0.6022814D+00
      RETURN
C
C           MAGNESIUM
C
  200 CONTINUE
      E(1)=650.643367D+00
      E(2)=98.370780D+00
      E(3)=21.322490D+00
      C(1)=0.0680297D+00
      C(2)=0.3907384D+00
      C(3)=0.6672673D+00
      E(4)=27.977380D+00
      E(5)=2.3265200D+00
      E(6)=0.8180800D+00
      C(4)=-0.0867195D+00
      C(5)=0.5856969D+00
      C(6)=0.4864974D+00
      E(7)=0.9113400D+00
      E(8)=0.0994000D+00
      E(9)=0.0362300D+00
      C(7)=-0.1276505D+00
      C(8)=0.6507728D+00
      C(9)=0.4362719D+00
      E(10)=23.216620D+00
      E(11)=5.0022200D+00
      E(12)=1.2046500D+00
      C(10)=0.1214603D+00
      C(11)=0.4792914D+00
      C(12)=0.5989417D+00
      RETURN
C
C           ALUMINUM
C
  300 CONTINUE
      E(1)=777.44334D+00
      E(2)=117.23153D+00
      E(3)=25.376297D+00
      C(1)=0.0668865D+00
      C(2)=0.3877681D+00
      C(3)=0.6707027D+00
      E(4)=33.356253D+00
      E(5)=2.8013153D+00
      E(6)=1.0227334D+00
      C(4)=-0.0889564D+00
      C(5)=0.6010610D+00
      C(6)=0.4687856D+00
      E(7)=1.3291880D+00
      E(8)=0.1616630D+00
      E(9)=0.05898876D+00
      C(7)=-0.1513892D+00
      C(8)=0.6593862D+00
      C(9)=0.4386928D+00
      E(10)=30.569580D+00
      E(11)=6.6447007D+00
      E(12)=1.6543954D+00
      C(10)=0.1123529D+00
      C(11)=0.4674668D+00
      C(12)=0.6097823D+00
      E(13)=0.0446518D+00
      E(14)=0.29965103D+00
      E(15)=0.11584816D+00
      C(13)=0.3964698D+00
      C(14)=0.2264326D+00
      C(15)=0.5005859D+00
      RETURN
C
C           SILICON
C
  400 CONTINUE
      E(1)=909.23487D+00
      E(2)=137.12456D+00
      E(3)=29.714810D+00
      C(1)=0.0664053D+00
      C(2)=0.3862219D+00
      C(3)=0.6722400D+00
      E(4)=39.129423D+00
      E(5)=3.3359807D+00
      E(6)=1.2512594D+00
      C(4)=-0.0909994D+00
      C(5)=0.6116149D+00
      C(6)=0.4568602D+00
      E(7)=1.8031012D+00
      E(8)=0.22638923D+00
      E(9)=0.08239593D+00
      C(7)=-0.1687328D+00
      C(8)=0.6754527D+00
      C(9)=0.4294192D+00
      E(10)=37.881761D+00
      E(11)=8.3045980D+00
      E(12)=2.1207924D+00
      C(10)=0.1087525D+00
      C(11)=0.4635152D+00
      C(12)=0.6113343D+00
      E(13)=0.44780246D+00
      E(14)=0.06236175D+00
      E(15)=0.17083752D+00
      C(13)=0.2389130D+00
      C(14)=0.3454528D+00
      C(15)=0.5422949D+00
      RETURN
C
C           PHOSPHORUS
C
  500 CONTINUE
      E(1)=1053.2658D+00
      E(2)=158.79044D+00
      E(3)=34.424407D+00
      C(1)=0.0658649D+00
      C(2)=0.3845776D+00
      C(3)=0.6739634D+00
      E(4)=45.450377D+00
      E(5)=3.8999265D+00
      E(6)=1.4885071D+00
      C(4)=-0.0926548D+00
      C(5)=0.6265125D+00
      C(6)=0.4410392D+00
      E(7)=2.3008408D+00
      E(8)=0.29895928D+00
      E(9)=0.10885358D+00
      C(7)=-0.1805490D+00
      C(8)=0.6809516D+00
      C(9)=0.4291418D+00
      E(10)=46.100019D+00
      E(11)=10.165057D+00
      E(12)=2.6447940D+00
      C(10)=0.1053880D+00
      C(11)=0.4597120D+00
      C(12)=0.6137136D+00
      E(13)=0.63268608D+00
      E(14)=0.24021895D+00
      E(15)=0.08644653D+00
      C(13)=0.2358855D+00
      C(14)=0.5541603D+00
      C(15)=0.3365297D+00
      RETURN
C
C           SULFUR
C
  600 CONTINUE
      E(1)=1201.4584D+00
      E(2)=181.39212D+00
      E(3)=39.404795D+00
      C(1)=0.0657647D+00
      C(2)=0.3839476D+00
      C(3)=0.6743724D+00
      E(4)=52.139030D+00
      E(5)=4.5287990D+00
      E(6)=1.7549380D+00
      C(4)=-0.0942319D+00
      C(5)=0.6354680D+00
      C(6)=0.4315059D+00
      E(7)=2.8404372D+00
      E(8)=0.38143256D+00
      E(9)=0.13878583D+00
      C(7)=0.1900423D+00
      C(8)=-0.6855267D+00
      C(9)=-0.4292724D+00
      E(10)=54.644071D+00
      E(11)=12.122902D+00
      E(12)=3.2065037D+00
      C(10)=0.1036728D+00
      C(11)=0.4581901D+00
      C(12)=0.6134001D+00
      E(13)=0.86327394D+00
      E(14)=0.10867870D+00
      E(15)=0.31812973D+00
      C(13)=0.2294364D+00
      C(14)=0.3537001D+00
      C(15)=0.5529600D+00
      RETURN
C
C           CHLORINE
C
  700 CONTINUE
      E(1)=1362.0220D+00
      E(2)=205.81110D+00
      E(3)=44.772167D+00
      C(1)=0.0655436D+00
      C(2)=0.3829867D+00
      C(3)=0.6752103D+00
      E(4)=59.225732D+00
      E(5)=5.2139024D+00
      E(6)=2.0473460D+00
      C(4)=-0.0956202D+00
      C(5)=0.6414256D+00
      C(6)=0.4251534D+00
      E(7)=3.4197116D+00
      E(8)=0.47001739D+00
      E(9)=0.16995836D+00
      C(7)=0.1964008D+00
      C(8)=-0.6923599D+00
      C(9)=-0.4261933D+00
      E(10)=64.099958D+00
      E(11)=14.287139D+00
      E(12)=3.8281352D+00
      C(10)=0.1017893D+00
      C(11)=0.4561069D+00
      C(12)=0.6142825D+00
      E(13)=1.0951257D+00
      E(14)=0.13217647D+00
      E(15)=0.39600382D+00
      C(13)=0.2359029D+00
      C(14)=0.3465996D+00
      C(15)=0.5580658D+00
      RETURN
C
C           ARGON
C
  800 CONTINUE
      E(1)=1536.9325D+00
      E(2)=232.17765D+00
      E(3)=50.521685D+00
      C(1)=0.0651589D+00
      C(2)=0.3818067D+00
      C(3)=0.6764462D+00
      E(4)=66.933949D+00
      E(5)=5.9185516D+00
      E(6)=2.3393420D+00
      C(4)=-0.0967405D+00
      C(5)=0.6527492D+00
      C(6)=0.4135728D+00
      E(7)=4.0453072D+00
      E(8)=0.56570104D+00
      E(9)=0.20406479D+00
      C(7)=0.2007356D+00
      C(8)=-0.6962680D+00
      C(9)=-0.4248435D+00
      E(10)=74.352915D+00
      E(11)=16.631346D+00
      E(12)=4.5039275D+00
      C(10)=0.1000791D+00
      C(11)=0.4542258D+00
      C(12)=0.6152592D+00
      E(13)=1.3570912D+00
      E(14)=0.48811284D+00
      E(15)=0.16212589D+00
      C(13)=0.2372757D+00
      C(14)=0.5583599D+00
      C(15)=0.3461648D+00
      RETURN
      END
C*MODULE BASHUZ  *DECK HFOUR
      SUBROUTINE HFOUR(E,C,NUCZ)
      IMPLICIT DOUBLE PRECISION(A-H,O-Z)
      DIMENSION E(24),C(24)
C
      IF(NUCZ.EQ.19) GO TO 100
      IF(NUCZ.EQ.20) GO TO 200
      IF(NUCZ.LE.30) CALL HTMONE(E,C,NUCZ)
      IF(NUCZ.LE.30) RETURN
      IBR=NUCZ-30
      GO TO (300,400,500,600,700,800), IBR
C
C           POTASSIUM
C
  100 CONTINUE
      E(1)=1721.1755D+00
      E(2)=260.01633D+00
      E(3)=56.624554D+00
      C(1)=0.0648747D+00
      C(2)=0.3808593D+00
      C(3)=0.6773681D+00
      E(4)=75.055600D+00
      E(5)=6.6911626D+00
      E(6)=2.6671665D+00
      C(4)=-0.0978731D+00
      C(5)=0.6595598D+00
      C(6)=0.4065295D+00
      E(7)=4.6689367D+00
      E(8)=0.70001332D+00
      E(9)=0.27533364D+00
      C(7)=-0.2129302D+00
      C(8)=0.6892491D+00
      C(9)=0.4339254D+00
      E(10)=0.25231153D+00
      E(11)=0.03763114D+00
      E(12)=0.01621835D+00
      C(10)=-0.1529800D+00
      C(11)=0.6817284D+00
      C(12)=0.4144829D+00
      E(13)=85.789846D+00
      E(14)=19.254794D+00
      E(15)=5.2686239D+00
      C(13)=0.0980354D+00
      C(14)=0.4510408D+00
      C(15)=0.6177076D+00
      E(16)=1.6831452D+00
      E(17)=0.62580945D+00
      E(18)=0.22389802D+00
      C(16)=0.2386360D+00
      C(17)=0.5710568D+00
      C(18)=0.3165357D+00
      RETURN
C
C           CALCIUM
C
  200 CONTINUE
      E(1)=1915.4348D+00
      E(2)=289.53324D+00
      E(3)=63.106352D+00
      C(1)=0.0646237D+00
      C(2)=0.3798376D+00
      C(3)=0.6783294D+00
      E(4)=83.633281D+00
      E(5)=7.5118396D+00
      E(6)=3.0145846D+00
      C(4)=-0.0988822D+00
      C(5)=0.6660273D+00
      C(6)=0.3999119D+00
      E(7)=5.3707540D+00
      E(8)=0.83837962D+00
      E(9)=0.34622613D+00
      C(7)=-0.2235309D+00
      C(8)=0.7028863D+00
      C(9)=0.4235525D+00
      E(10)=0.42974055D+00
      E(11)=0.06217342D+00
      E(12)=0.02494803D+00
      C(10)=-0.1793098D+00
      C(11)=0.6728792D+00
      C(12)=0.4356664D+00
      E(13)=97.974592D+00
      E(14)=22.067384D+00
      E(15)=6.0938756D+00
      C(13)=0.0963163D+00
      C(14)=0.4481076D+00
      C(15)=0.6199212D+00
      E(16)=2.0178863D+00
      E(17)=0.76665059D+00
      E(18)=0.28431923D+00
      C(16)=0.2450293D+00
      C(17)=0.5843849D+00
      C(18)=0.2876662D+00
      RETURN
C
C           GALLIUM
C
  300 CONTINUE
      E(1)=4751.8979D+00
      E(2)=718.92054D+00
      E(3)=157.44592D+00
      C(1)=0.0628396D+00
      C(2)=0.3736112D+00
      C(3)=0.6843626D+00
      E(4)=210.53833D+00
      E(5)=19.794805D+00
      E(6)=8.0497965D+00
      C(4)=-0.1054027D+00
      C(5)=0.7070778D+00
      C(6)=0.3599497D+00
      E(7)=16.424933D+00
      E(8)=2.6576527D+00
      E(9)=1.0609873D+00
      C(7)=0.2270732D+00
      C(8)=-0.7111408D+00
      C(9)=-0.4213859D+00
      E(10)=1.3533154D+00
      E(11)=0.17615808D+00
      E(12)=0.06371947D+00
      C(10)=-0.1507084D+00
      C(11)=0.6413441D+00
      C(12)=0.4600643D+00
      E(13)=273.64584D+00
      E(14)=62.815151D+00
      E(15)=18.134221D+00
      C(13)=0.0908956D+00
      C(14)=0.4456493D+00
      C(15)=0.6179026D+00
      E(16)=6.8387139D+00
      E(17)=2.5958668D+00
      E(18)=0.98389370D+00
      C(16)=0.2739015D+00
      C(17)=0.5689512D+00
      C(18)=0.2756075D+00
      E(19)=0.31568846D+00
      E(20)=0.10961938D+00
      E(21)=0.03994387D+00
      C(19)=0.2062964D+00
      C(20)=0.5563565D+00
      C(21)=0.3670988D+00
      E(22)=21.292530D+00
      E(23)=5.3931662D+00
      E(24)=1.3338828D+00
      C(22)=0.1619895D+00
      C(23)=0.5116739D+00
      C(24)=0.5898732D+00
      RETURN
C
C           GERMANIUM
C
  400 CONTINUE
      E(1)=5073.7499D+00
      E(2)=767.72417D+00
      E(3)=168.18881D+00
      C(1)=0.0627249D+00
      C(2)=0.3731671D+00
      C(3)=0.6847867D+00
      E(4)=224.98366D+00
      E(5)=21.187659D+00
      E(6)=8.6256370D+00
      C(4)=-0.1057096D+00
      C(5)=0.7094592D+00
      C(6)=0.3575586D+00
      E(7)=17.764305D+00
      E(8)=2.8551503D+00
      E(9)=1.1585827D+00
      C(7)=-0.2280577D+00
      C(8)=0.7250544D+00
      C(9)=0.4064506D+00
      E(10)=1.5859493D+00
      E(11)=0.22692409D+00
      E(12)=0.08372494D+00
      C(10)=-0.1769108D+00
      C(11)=0.6604420D+00
      C(12)=0.4525469D+00
      E(13)=293.73174D+00
      E(14)=67.552546D+00
      E(15)=19.568631D+00
      C(13)=0.0906975D+00
      C(14)=0.4451940D+00
      C(15)=0.6178940D+00
      E(16)=7.4326168D+00
      E(17)=2.8426842D+00
      E(18)=1.1041168D+00
      C(16)=0.2751638D+00
      C(17)=0.5703353D+00
      C(18)=0.2690321D+00
      E(19)=0.38224145D+00
      E(20)=0.14147296D+00
      E(21)=0.05262765D+00
      C(19)=0.2627754D+00
      C(20)=0.5701236D+00
      C(21)=0.2901737D+00
      E(22)=24.321421D+00
      E(23)=6.2238135D+00
      E(24)=1.5887375D+00
      C(22)=0.1577985D+00
      C(23)=0.5114922D+00
      C(24)=0.5857703D+00
      RETURN
C
C           ARSENIC
C
  500 CONTINUE
      E(1)=5407.6138D+00
      E(2)=818.17436D+00
      E(3)=179.26569D+00
      C(1)=0.0626011D+00
      C(2)=0.3727790D+00
      C(3)=0.6851842D+00
      E(4)=239.85464D+00
      E(5)=22.64906D+00
      E(6)=9.2280039D+00
      C(4)=-0.1060438D+00
      C(5)=0.7107219D+00
      C(6)=0.3563753D+00
      E(7)=19.058432D+00
      E(8)=3.0854133D+00
      E(9)=1.2727348D+00
      C(7)=-0.2308017D+00
      C(8)=0.7325742D+00
      C(9)=0.3994769D+00
      E(10)=1.8198631D+00
      E(11)=0.27883046D+00
      E(12)=0.10372930D+00
      C(10)=-0.1983584D+00
      C(11)=0.6759512D+00
      C(12)=0.4476004D+00
      E(13)=315.02252D+00
      E(14)=72.557901D+00
      E(15)=21.079855D+00
      C(13)=0.0903237D+00
      C(14)=0.4443841D+00
      C(15)=0.6183953D+00
      E(16)=8.1044166D+00
      E(17)=3.1282025D+00
      E(18)=1.2391193D+00
      C(16)=0.2728475D+00
      C(17)=0.5717222D+00
      C(18)=0.2664488D+00
      E(19)=0.45260113D+00
      E(20)=0.17204402D+00
      E(21)=0.06573731D+00
      C(19)=0.3086721D+00
      C(20)=0.5671752D+00
      C(21)=0.2410732D+00
      E(22)=27.437209D+00
      E(23)=7.0840440D+00
      E(24)=1.8558226D+00
      C(22)=0.1544952D+00
      C(23)=0.5114318D+00
      C(24)=0.5821935D+00
      RETURN
C
C           SELENIUM
C
  600 CONTINUE
      E(1)=5751.3215D+00
      E(2)=870.25721D+00
      E(3)=190.72949D+00
      C(1)=0.0624934D+00
      C(2)=0.3723683D+00
      C(3)=0.6855799D+00
      E(4)=255.18802D+00
      E(5)=24.155395D+00
      E(6)=9.8525390D+00
      C(4)=-0.1063403D+00
      C(5)=0.7120577D+00
      C(6)=0.3550792D+00
      E(7)=20.463481D+00
      E(8)=3.3125921D+00
      E(9)=1.3791545D+00
      C(7)=-0.2326525D+00
      C(8)=0.7465053D+00
      C(9)=0.3857492D+00
      E(10)=2.0492750D+00
      E(11)=0.33641155D+00
      E(12)=0.12612688D+00
      C(10)=-0.2169517D+00
      C(11)=0.6852937D+00
      C(12)=0.4486601D+00
      E(13)=337.08753D+00
      E(14)=77.748281D+00
      E(15)=22.650400D+00
      C(13)=0.0899627D+00
      C(14)=0.4435785D+00
      C(15)=0.6189101D+00
      E(16)=8.6914000D+00
      E(17)=3.3668794D+00
      E(18)=1.3546532D+00
      C(16)=0.2789803D+00
      C(17)=0.5760239D+00
      C(18)=0.2529785D+00
      E(19)=0.54258450D+00
      E(20)=0.20662960D+00
      E(21)=0.07760208D+00
      C(19)=0.3286253D+00
      C(20)=0.5489677D+00
      C(21)=0.2444426D+00
      E(22)=30.627464D+00
      E(23)=7.9712764D+00
      E(24)=2.1348097D+00
      C(22)=0.1519858D+00
      C(23)=0.5116403D+00
      C(24)=0.5786936D+00
      RETURN
C
C           BROMINE
C
  700 CONTINUE
      E(1)=6103.2899D+00
      E(2)=923.69743D+00
      E(3)=202.52031D+00
      C(1)=0.0624175D+00
      C(2)=0.3720414D+00
      C(3)=0.6858728D+00
      E(4)=270.99649D+00
      E(5)=25.706395D+00
      E(6)=10.501175D+00
      C(4)=-0.1066146D+00
      C(5)=0.7133792D+00
      C(6)=0.3537553D+00
      E(7)=21.859513D+00
      E(8)=3.5681587D+00
      E(9)=1.5090061D+00
      C(7)=-0.2354458D+00
      C(8)=0.7528340D+00
      C(9)=0.3801608D+00
      E(10)=2.3404284D+00
      E(11)=0.38825413D+00
      E(12)=0.14598986D+00
      C(10)=-0.2363212D+00
      C(11)=0.7085531D+00
      C(12)=0.4329347D+00
      E(13)=360.49973D+00
      E(14)=83.244591D+00
      E(15)=24.311838D+00
      C(13)=0.0894281D+00
      C(14)=0.4423065D+00
      C(15)=0.6200053D+00
      E(16)=9.4221830D+00
      E(17)=3.6771841D+00
      E(18)=1.4935778D+00
      C(16)=0.2767522D+00
      C(17)=0.5795347D+00
      C(18)=0.2491162D+00
      E(19)=0.63454669D+00
      E(20)=0.24252803D+00
      E(21)=0.09053005D+00
      C(19)=0.3466285D+00
      C(20)=0.5413298D+00
      C(21)=0.2347667D+00
      E(22)=33.965097D+00
      E(23)=8.9008312D+00
      E(24)=2.4284360D+00
      C(22)=0.1496666D+00
      C(23)=0.5117475D+00
      C(24)=0.5759148D+00
      RETURN
C
C           KRYPTON
C
  800 CONTINUE
      E(1)=6446.6307D+00
      E(2)=976.87570D+00
      E(3)=214.47955D+00
      C(1)=0.0625398D+00
      C(2)=0.3721075D+00
      C(3)=0.6856107D+00
      E(4)=287.01158D+00
      E(5)=27.328672D+00
      E(6)=11.174440D+00
      C(4)=-0.1068933D+00
      C(5)=0.7139493D+00
      C(6)=0.3532914D+00
      E(7)=23.331340D+00
      E(8)=3.8133080D+00
      E(9)=1.6134940D+00
      C(7)=-0.2376841D+00
      C(8)=0.7698488D+00
      C(9)=0.3641119D+00
      E(10)=2.5964926D+00
      E(11)=0.4490900D+00
      E(12)=0.17012362D+00
      C(10)=-0.2498594D+00
      C(11)=0.7128783D+00
      C(12)=0.4361011D+00
      E(13)=384.37950D+00
      E(14)=88.868156D+00
      E(15)=26.019696D+00
      C(13)=0.0890288D+00
      C(14)=0.4413316D+00
      C(15)=0.6207416D+00
      E(16)=10.184792D+00
      E(17)=4.0094356D+00
      E(18)=1.6459659D+00
      C(16)=0.2745449D+00
      C(17)=0.5813334D+00
      C(18)=0.2470820D+00
      E(19)=0.72962570D+00
      E(20)=0.28060560D+00
      E(21)=0.10521386D+00
      C(19)=0.3620055D+00
      C(20)=0.5352020D+00
      C(21)=0.2242007D+00
      E(22)=37.368103D+00
      E(23)=9.8543131D+00
      E(24)=2.7327955D+00
      C(22)=0.1479466D+00
      C(23)=0.5121719D+00
      C(24)=0.5729498D+00
      RETURN
      END
C*MODULE BASHUZ  *DECK HTMONE
      SUBROUTINE HTMONE(E,C,NUCZ)
      IMPLICIT DOUBLE PRECISION(A-H,O-Z)
      DIMENSION E(21),C(21)
C
      IBR=NUCZ-20
      GO TO (100,200,300,400,500,600,700,800,900,1000), IBR
C
C     SCANDIUM (4-F)
C
  100 CONTINUE
      E(1)=2121.1153D+00
      E(2)=320.52014D+00
      E(3)=69.895734D+00
      C(1)=0.0643811D+00
      C(2)=0.3791356D+00
      C(3)=0.6790384D+00
      E(4)=92.827788D+00
      E(5)=8.3661954D+00
      E(6)=3.3619693D+00
      C(4)=-0.0997668D+00
      C(5)=0.6751087D+00
      C(6)=0.3908713D+00
      E(7)=6.1732210D+00
      E(8)=0.95574656D+00
      E(9)=0.37934186D+00
      C(7)=-0.2224250D+00
      C(8)=0.71394D+00
      C(9)=0.4142398D+00
      E(10)=0.45254574D+00
      E(11)=0.05968487D+00
      E(12)=0.02352408D+00
      C(10)=-0.1502541D+00
      C(11)=0.651097D+00
      C(12)=0.4444372D+00
      E(13)=110.16022D+00
      E(14)=24.886185D+00
      E(15)=6.9221188D+00
      C(13)=0.0957608D+00
      C(14)=0.4482512D+00
      C(15)=0.6189104D+00
      E(16)=2.3903679D+00
      E(17)=0.90248233D+00
      E(18)=0.32475569D+00
      C(16)=0.2403846D+00
      C(17)=0.5797456D+00
      C(18)=0.3020693D+00
      E(19)=4.526837D+00
      E(20)=1.0394154D+00
      E(21)=0.21850767D+00
      C(19)=0.1644737D+00
      C(20)=0.4787843D+00
      C(21)=0.6556389D+00
      RETURN
C
C     TITANIUM (5-F)
C
  200 CONTINUE
      E( 1)=2336.7153D+00
      E( 2)=353.15331D+00
      E( 3)=77.059812D+00
      C( 1)=0.0641708D+00
      C( 2)=0.3783819D+00
      C( 3)=0.6797627D+00
      E( 4)=102.45710D+00
      E( 5)=9.2769361D+00
      E( 6)=3.7368558D+00
      C( 4)=-0.1005637D+00
      C( 5)=0.6815896D+00
      C( 6)=0.3844102D+00
      E( 7)=7.0008762D+00
      E( 8)=1.0875122D+00
      E( 9)=0.42697745D+00
      C( 7)=-0.2231174D+00
      C( 8)=0.7199595D+00
      C( 9)=0.4094673D+00
      E(10)=0.54148021D+00
      E(11)=0.06648139D+00
      E(12)=0.02576476D+00
      C(10)=-0.1434300D+00
      C(11)=0.6030041D+00
      C(12)=0.4881755D+00
      E(13)=123.33273D+00
      E(14)=27.926218D+00
      E(15)=7.8132163D+00
      C(13)=0.0949569D+00
      C(14)=0.4477326D+00
      C(15)=0.6188790D+00
      E(16)=2.75228040D+00
      E(17)=1.0394584D+00
      E(18)=0.37298279D+00
      C(16)=0.2438232D+00
      C(17)=0.5781427D+00
      C(18)=0.3010109D+00
      E(19)=5.9418882D+00
      E(20)=1.3963895D+00
      E(21)=0.30602356D+00
      C(19)=0.1638319D+00
      C(20)=0.4905779D+00
      C(21)=0.6362293D+00
      RETURN
C
C     VANADIUM (6-D)
C
  300 CONTINUE
      E( 1)=2563.4389D+00
      E( 2)=387.46176D+00
      E( 3)=84.577619D+00
      C( 1)=0.0639578D+00
      C( 2)=0.3776589D+00
      C( 3)=0.6804836D+00
      E( 4)=112.57688D+00
      E( 5)=10.242910D+00
      E( 6)=4.1346585D+00
      C( 4)=-0.1013124D+00
      C( 5)=0.6865452D+00
      C( 6)=0.3795191D+00
      E( 7)=7.8770045D+00
      E( 8)=1.2266238D+00
      E( 9)=0.47757702D+00
      C( 7)=-0.2233367D+00
      C( 8)=0.7226660D+00
      C( 9)=0.4075932D+00
      E(10)=0.61583188D+00
      E(11)=0.07247550D+00
      E(12)=0.02753023D+00
      C(10)=-0.1329102D+00
      C(11)=0.5794076D+00
      C(12)=0.5081949D+00
      E(13)=136.91997D+00
      E(14)=31.081206D+00
      E(15)=8.7448393D+00
      C(13)=0.0945211D+00
      C(14)=0.4477583D+00
      C(15)=0.6181758D+00
      E(16)=3.1310711D+00
      E(17)=1.1822484D+00
      E(18)=0.42352383D+00
      C(16)=0.2473272D+00
      C(17)=0.5762748D+00
      C(18)=0.3001113D+00
      E(19)=7.2698474D+00
      E(20)=1.7314426D+00
      E(21)=0.38675247D+00
      C(19)=0.1642244D+00
      C(20)=0.4967508D+00
      C(21)=0.6258979D+00
      RETURN
C
C     CHROMIUM (7-S)
C
  400 CONTINUE
      E( 1)=2799.7104D+00
      E( 2)=423.21821D+00
      E( 3)=92.434349D+00
      C( 1)=0.0637916D+00
      C( 2)=0.3770772D+00
      C( 3)=0.6810364D+00
      E( 4)=123.22634D+00
      E( 5)=11.246630D+00
      E( 6)=4.5372546D+00
      C( 4)=-0.1019598D+00
      C( 5)=0.6925012D+00
      C( 6)=0.3737091D+00
      E( 7)=8.7804852D+00
      E( 8)=1.3793596D+00
      E( 9)=0.53502318D+00
      C( 7)=-0.2235841D+00
      C( 8)=0.7175520D+00
      C( 9)=0.4136269D+00
      E(10)=0.72655885D+00
      E(11)=0.07639020D+00
      E(12)=0.02860476D+00
      C(10)=-0.1294804D+00
      C(11)=0.5711774D+00
      C(12)=0.5126118D+00
      E(13)=151.28977D+00
      E(14)=34.410311D+00
      E(15)=9.7258583D+00
      C(13)=0.0940601D+00
      C(14)=0.4477828D+00
      C(15)=0.6176532D+00
      E(16)=3.4776582D+00
      E(17)=1.3079405D+00
      E(18)=0.46753229D+00
      C(16)=0.2581111D+00
      C(17)=0.5757324D+00
      C(18)=0.2908906D+00
      E(19)=8.6262099D+00
      E(20)=2.0754716D+00
      E(21)=0.47058738D+00
      C(19)=0.1647178D+00
      C(20)=0.5011371D+00
      C(21)=0.6180587D+00
      RETURN
C
C     MANGANESE (6-D)
C
  500 CONTINUE
      E( 1)=3044.3379D+00
      E( 2)=460.32835D+00
      E( 3)=100.60952D+00
      C( 1)=0.0636857D+00
      C( 2)=0.3766602D+00
      C( 3)=0.6813959D+00
      E( 4)=134.04812D+00
      E( 5)=12.355650D+00
      E( 6)=5.0118328D+00
      C( 4)=-0.1026805D+00
      C( 5)=0.6915740D+00
      C( 6)=0.3747452D+00
      E( 7)=9.7428339D+00
      E( 8)=1.5375154D+00
      E( 9)=0.59408091D+00
      C( 7)=-0.2239616D+00
      C( 8)=0.7171970D+00
      C( 9)=0.4146924D+00
      E(10)=0.76261459D+00
      E(11)=0.07836474D+00
      E(12)=0.02930169D+00
      C(10)=-0.1159459D+00
      C(11)=0.5638312D+00
      C(12)=0.5157105D+00
      E(13)=166.73180D+00
      E(14)=37.977416D+00
      E(15)=10.774333D+00
      C(13)=0.0933798D+00
      C(14)=0.4472796D+00
      C(15)=0.6178799D+00
      E(16)=3.9731007D+00
      E(17)=1.5014460D+00
      E(18)=0.53680586D+00
      C(16)=0.2500751D+00
      C(17)=0.5725861D+00
      C(18)=0.3020932D+00
      E(19)=9.8656162D+00
      E(20)=2.3866048D+00
      E(21)=0.53817952D+00
      C(19)=0.1674847D+00
      C(20)=0.5041187D+00
      C(21)=0.6143492D+00
      RETURN
C
C     IRON (5-F)
C
  600 CONTINUE
      E( 1)=3302.0297D+00
      E( 2)=499.32535D+00
      E( 3)=109.17229D+00
      C( 1)=0.0635283D+00
      C( 2)=0.3761252D+00
      C( 3)=0.6819206D+00
      E( 4)=145.61613D+00
      E( 5)=13.468079D+00
      E( 6)=5.4624842D+00
      C( 4)=-0.1032386D+00
      C( 5)=0.6955826D+00
      C( 6)=0.3708869D+00
      E( 7)=10.736268D+00
      E( 8)=1.7009576D+00
      E( 9)=0.65364599D+00
      C( 7)=-0.2241245D+00
      C( 8)=0.7177415D+00
      C( 9)=0.4148646D+00
      E(10)=0.81194929D+00
      E(11)=0.07967335D+00
      E(12)=0.02982454D+00
      C(10)=-0.1083327D+00
      C(11)=0.5566047D+00
      C(12)=0.5194924D+00
      E(13)=182.55775D+00
      E(14)=41.647822D+00
      E(15)=11.857806D+00
      C(13)=0.0929956D+00
      C(14)=0.4473347D+00
      C(15)=0.6174492D+00
      E(16)=4.4051684D+00
      E(17)=1.6609361D+00
      E(18)=0.59330383D+00
      C(16)=0.2543695D+00
      C(17)=0.5722388D+00
      C(18)=0.2986741D+00
      E(19)=11.198379D+00
      E(20)=2.7235552D+00
      E(21)=0.61395523D+00
      C(19)=0.1694177D+00
      C(20)=0.5066518D+00
      C(21)=0.6105528D+00
      RETURN
C
C     COBALT (4-F)
C
  700 CONTINUE
      E( 1)=3569.1023D+00
      E( 2)=539.75057D+00
      E( 3)=118.06735D+00
      C( 1)=0.0634072D+00
      C( 2)=0.3757053D+00
      C( 3)=0.6823127D+00
      E( 4)=157.60903D+00
      E( 5)=14.633143D+00
      E( 6)=5.9354094D+00
      C( 4)=-0.1037551D+00
      C( 5)=0.6988574D+00
      C( 6)=0.3677667D+00
      E( 7)=11.790403D+00
      E( 8)=1.8690733D+00
      E( 9)=0.71451064D+00
      C( 7)=-0.2238506D+00
      C( 8)=0.7186408D+00
      C( 9)=0.4142824D+00
      E(10)=0.83085499D+00
      E(11)=0.08093542D+00
      E(12)=0.03021553D+00
      C(10)=-0.0999744D+00
      C(11)=0.5545757D+00
      C(12)=0.5193358D+00
      E(13)=198.85110D+00
      E(14)=45.443908D+00
      E(15)=12.984531D+00
      C(13)=0.0927924D+00
      C(14)=0.4476071D+00
      C(15)=0.6167211D+00
      E(16)=4.8835867D+00
      E(17)=1.8407375D+00
      E(18)=0.65687534D+00
      C(16)=0.2550456D+00
      C(17)=0.5708333D+00
      C(18)=0.2998891D+00
      E(19)=12.634196D+00
      E(20)=3.0847767D+00
      E(21)=0.69473162D+00
      C(19)=0.1704471D+00
      C(20)=0.5085159D+00
      C(21)=0.6080700D+00
      RETURN
C
C     NICKEL (3-D)
C
  800 CONTINUE
      E( 1)=3849.2417D+00
      E( 2)=582.08183D+00
      E( 3)=127.35676D+00
      C( 1)=0.0632462D+00
      C( 2)=0.3751737D+00
      C( 3)=0.6828494D+00
      E( 4)=170.11363D+00
      E( 5)=15.852695D+00
      E( 6)=6.4342950D+00
      C( 4)=-0.1042343D+00
      C( 5)=0.7012328D+00
      C( 6)=0.3655163D+00
      E( 7)=12.889647D+00
      E( 8)=2.0461387D+00
      E( 9)=0.77923389D+00
      C( 7)=-0.2236128D+00
      C( 8)=0.7180189D+00
      C( 9)=0.4151786D+00
      E(10)=0.84339749D+00
      E(11)=0.08247463D+00
      E(12)=0.03073497D+00
      C(10)=-0.0938590D+00
      C(11)=0.5512875D+00
      C(12)=0.5211700D+00
      E(13)=216.01321D+00
      E(14)=49.433819D+00
      E(15)=14.166185D+00
      C(13)=0.0925038D+00
      C(14)=0.4477120D+00
      C(15)=0.6162927D+00
      E(16)=5.3861988D+00
      E(17)=2.0323051D+00
      E(18)=0.72556764D+00
      C(16)=0.2552417D+00
      C(17)=0.5684974D+00
      C(18)=0.3023861D+00
      E(19)=14.087823D+00
      E(20)=3.4535186D+00
      E(21)=0.77782386D+00
      C(19)=0.1719093D+00
      C(20)=0.5102939D+00
      C(21)=0.6052430D+00
      RETURN
C
C     COPPER (2-S)
C
  900 CONTINUE
      E(1)=4137.5632D+00
      E(2)=625.73294D+00
      E(3)=136.95969D+00
      C(1)=0.0631365D+00
      C(2)=0.3747929D+00
      C(3)=0.6832095D+00
      E(4)=183.1051D+00
      E(5)=17.122418D+00
      E(6)=6.9491752D+00
      C(4)=-0.1046861D+00
      C(5)=0.7034840D+00
      C(6)=0.3634281D+00
      E(7)=13.999937D+00
      E(8)=2.2431225D+00
      E(9)=0.85310594D+00
      C(7)=-0.2239634D+00
      C(8)=0.7113033D+00
      C(9)=0.4227015D+00
      E(10)=0.83040542D+00
      E(11)=0.08540467D+00
      E(12)=0.03181982D+00
      C(10)=-0.0890494D+00
      C(11)=0.5338648D+00
      C(12)=0.5380713D+00
      E(13)=233.62353D+00
      E(14)=53.545569D+00
      E(15)=15.388041D+00
      C(13)=0.0923643D+00
      C(14)=0.4480449D+00
      C(15)=0.6155840D+00
      E(16)=5.8732783D+00
      E(17)=2.2084129D+00
      E(18)=0.78648338D+00
      C(16)=0.2593135D+00
      C(17)=0.5687950D+00
      C(18)=0.2987106D+00
      E(19)=15.586486D+00
      E(20)=3.8355780D+00
      E(21)=0.86504248D+00
      C(19)=0.1734535D+00
      C(20)=0.5121846D+00
      C(21)=0.6020146D+00
      RETURN
C
C     ZINC (1-S)
C
 1000 CONTINUE
      E( 1)=4432.2885D+00
      E( 2)=670.66012D+00
      E( 3)=146.90245D+00
      C( 1)=0.0630928D+00
      C( 2)=0.3745038D+00
      C( 3)=0.6834160D+00
      E( 4)=196.48911D+00
      E( 5)=18.441663D+00
      E( 6)=7.4929688D+00
      C( 4)=-0.1050863D+00
      C( 5)=0.7049994D+00
      C( 6)=0.3620077D+00
      E( 7)=15.180864D+00
      E( 8)=2.4467928D+00
      E( 9)=0.95293505D+00
      C( 7)=-0.2254217D+00
      C( 8)=0.7096746D+00
      C( 9)=0.4236169D+00
      E(10)=1.1210370D+00
      E(11)=0.1251838D+00
      E(12)=0.04406728D+00
      C(10)=-0.1166693D+00
      C(11)=0.5950219D+00
      C(12)=0.4923302D+00
      E(13)=252.47779D+00
      E(14)=57.934825D+00
      E(15)=16.689667D+00
      C(13)=0.0919491D+00
      C(14)=0.4476288D+00
      C(15)=0.6158034D+00
      E(16)=6.3410797D+00
      E(17)=2.3880503D+00
      E(18)=0.87307830D+00
      C(16)=0.2675452D+00
      C(17)=0.5701029D+00
      C(18)=0.2856214D+00
      E(19)=18.368202D+00
      E(20)=4.5913041D+00
      E(21)=1.0902026D+00
      C(19)=0.1672532D+00
      C(20)=0.5122790D+00
      C(21)=0.5948177D+00
      RETURN
      END
C*MODULE BASHUZ  *DECK HFIVE
      SUBROUTINE HFIVE(E,C,NUCZ)
      IMPLICIT DOUBLE PRECISION(A-H,O-Z)
      DIMENSION E(33),C(33)
C
      IF(NUCZ.EQ.37) GO TO 100
      IF(NUCZ.EQ.38) GO TO 200
      IF(NUCZ.LE.48) CALL HTMTWO(E,C,NUCZ)
      IF(NUCZ.LE.48) RETURN
      IBR=NUCZ-48
      GO TO (300,400,500,600,700,800), IBR
C
C           RUBIDIUM
C
  100 CONTINUE
      E(1)=6816.7225D+00
      E(2)=1033.0007D+00
      E(3)=226.90861D+00
      C(1)=0.0624962D+00
      C(2)=0.3719500D+00
      C(3)=0.6857293D+00
      E(4)=302.83628D+00
      E(5)=29.201467D+00
      E(6)=12.082845D+00
      C(4)=-0.1074004D+00
      C(5)=0.7049410D+00
      C(6)=0.3622692D+00
      E(7)=24.692611D+00
      E(8)=4.2016508D+00
      E(9)=1.8522603D+00
      C(7)=-0.2430518D+00
      C(8)=0.7375963D+00
      C(9)=0.3988455D+00
      E(10)=2.9359876D+00
      E(11)=0.52542640D+00
      E(12)=0.21916888D+00
      C(10)=-0.2806358D+00
      C(11)=0.7159840D+00
      C(12)=0.4402934D+00
      E(13)=0.19166785D+00
      E(14)=0.03343655D+00
      E(15)=0.01473352D+00
      C(13)=0.1580968D+00
      C(14)=-0.6855144D+00
      C(15)=-0.4191493D+00
      E(16)=407.86841D+00
      E(17)=94.455492D+00
      E(18)=27.745510D+00
      C(16)=0.0890053D+00
      C(17)=0.4411370D+00
      C(18)=0.6204330D+00
      E(19)=10.715052D+00
      E(20)=4.2311626D+00
      E(21)=1.7434701D+00
      C(19)=0.2867814D+00
      C(20)=0.5863660D+00
      C(21)=0.2277478D+00
      E(22)=0.84991897D+00
      E(23)=0.34353368D+00
      E(24)=0.13941053D+00
      C(22)=0.3796538D+00
      C(23)=0.5370593D+00
      C(24)=0.1869077D+00
      E(25)=40.866031D+00
      E(26)=10.840885D+00
      E(27)=3.0508341D+00
      C(25)=0.1466037D+00
      C(26)=0.5127252D+00
      C(27)=0.5699804D+00
      RETURN
C
C           STRONTIUM
C
  200 CONTINUE
      E(1)=7215.4735D+00
      E(2)=1092.8519D+00
      E(3)=239.98182D+00
      C(1)=0.0622818D+00
      C(2)=0.3713101D+00
      C(3)=0.6864439D+00
      E(4)=320.46089D+00
      E(5)=30.742778D+00
      E(6)=12.612079D+00
      C(4)=-0.1074604D+00
      C(5)=0.7140333D+00
      C(6)=0.3533614D+00
      E(7)=26.234117D+00
      E(8)=4.4599505D+00
      E(9)=1.9839250D+00
      C(7)=-0.2455863D+00
      C(8)=0.7556468D+00
      C(9)=0.3814041D+00
      E(10)=3.1962046D+00
      E(11)=0.60836651D+00
      E(12)=0.26573519D+00
      C(10)=-0.2974230D+00
      C(11)=0.7364992D+00
      C(12)=0.4298103D+00
      E(13)=0.33420177D+00
      E(14)=0.05169948D+00
      E(15)=0.02174977D+00
      C(13)=0.1961392D+00
      C(14)=-0.6888851D+00
      C(15)=-0.4263861D+00
      E(16)=434.56915D+00
      E(17)=100.70360D+00
      E(18)=29.625034D+00
      C(16)=0.0882772D+00
      C(17)=0.4395651D+00
      C(18)=0.6220621D+00
      E(19)=11.572090D+00
      E(20)=4.5937454D+00
      E(21)=1.9060808D+00
      C(19)=0.2829866D+00
      C(20)=0.5896899D+00
      C(21)=0.2263017D+00
      E(22)=0.95506174D+00
      E(23)=0.39453919D+00
      E(24)=0.15759333D+00
      C(22)=0.4130687D+00
      C(23)=0.5407330D+00
      C(24)=0.1403073D+00
      E(25)=44.566115D+00
      E(26)=11.881489D+00
      E(27)=3.3875579D+00
      C(25)=0.1451271D+00
      C(26)=0.5130677D+00
      C(27)=0.5676640D+00
      RETURN
C
C           INDIUM
C
  300 CONTINUE
      E(1)=12214.547D+00
      E(2)=1848.9136D+00
      E(3)=406.36833D+00
      C(1)=0.0612476D+00
      C(2)=0.3676754D+00
      C(3)=0.6901359D+00
      E(4)=541.34239D+00
      E(5)=53.018379D+00
      E(6)=22.190283D+00
      C(4)=-0.1096575D+00
      C(5)=0.7156786D+00
      C(6)=0.3516058D+00
      E(7)=46.414454D+00
      E(8)=8.3531661D+00
      E(9)=3.8649088D+00
      C(7)=-0.2690027D+00
      C(8)=0.8333355D+00
      C(9)=0.3162081D+00
      E(10)=7.1775944D+00
      E(11)=1.621446D+00
      E(12)=0.71996916D+00
      C(10)=-0.3439307D+00
      C(11)=0.7769627D+00
      C(12)=0.4273271D+00
      E(13)=1.0654870D+00
      E(14)=0.14431502D+00
      E(15)=0.05458735D+00
      C(13)=-0.2008092D+00
      C(14)=0.6599925D+00
      C(15)=0.4553545D+00
      E(16)=781.86921D+00
      E(17)=182.79818D+00
      E(18)=54.211142D+00
      C(16)=0.0821886D+00
      C(17)=0.4267795D+00
      C(18)=0.6357629D+00
      E(19)=114.03634D+00
      E(20)=18.713126D+00
      E(21)=7.1660478D+00
      C(19)=-0.0311407D+00
      C(20)=0.4133453D+00
      C(21)=0.6578615D+00
      E(22)=2.79637815D+00
      E(23)=1.1841215D+00
      E(24)=0.48773463D+00
      C(22)=0.4068683D+00
      C(23)=0.5479905D+00
      C(24)=0.1323252D+00
      E(25)=0.22993404D+00
      E(26)=0.08355539D+00
      E(27)=0.03249889D+00
      C(25)=0.2630070D+00
      C(26)=0.5733118D+00
      C(27)=0.2829178D+00
      E(28)=102.17356D+00
      E(29)=28.394632D+00
      E(30)=8.9248045D+00
      C(28)=0.1205559D+00
      C(29)=0.4884976D+00
      C(30)=0.5850190D+00
      E(31)=4.5353637D+00
      E(32)=1.5371481D+00
      E(33)=0.49949226D+00
      C(31)=0.2508574D+00
      C(32)=0.5693113D+00
      C(33)=0.3840635D+00
      RETURN
C
C           TIN
C
  400 CONTINUE
      E(1)=12741.674D+00
      E(2)=1928.4692D+00
      E(3)=423.80797D+00
      C(1)=0.0611353D+00
      C(2)=0.3672929D+00
      C(3)=0.6905447D+00
      E(4)=564.68220D+00
      E(5)=55.241831D+00
      E(6)=23.121421D+00
      C(4)=-0.1097488D+00
      C(5)=0.7181120D+00
      C(6)=0.3490973D+00
      E(7)=48.493098D+00
      E(8)=8.7703153D+00
      E(9)=4.0378104D+00
      C(7)=-0.2708984D+00
      C(8)=0.8409064D+00
      C(9)=0.3101030D+00
      E(10)=7.6154955D+00
      E(11)=1.7298929D+00
      E(12)=0.77188815D+00
      C(10)=-0.3475656D+00
      C(11)=0.7937583D+00
      C(12)=0.4132006D+00
      E(13)=1.1976547D+00
      E(14)=0.17303676D+00
      E(15)=0.06695953D+00
      C(13)=-0.2385589D+00
      C(14)=0.7059861D+00
      C(15)=0.4236124D+00
      E(16)=816.77323D+00
      E(17)=191.12772D+00
      E(18)=56.771595D+00
      C(16)=0.0820729D+00
      C(17)=0.4264331D+00
      C(18)=0.6359016D+00
      E(19)=121.24487D+00
      E(20)=19.577904D+00
      E(21)=7.5717131D+00
      C(19)=-0.0315643D+00
      C(20)=0.4151667D+00
      C(21)=0.6549931D+00
      E(22)=3.0474031D+00
      E(23)=1.3022184D+00
      E(24)=0.54462907D+00
      C(22)=0.3982670D+00
      C(23)=0.5557792D+00
      C(24)=0.1301395D+00
      E(25)=0.26752856D+00
      E(26)=0.10418254D+00
      E(27)=0.04162604D+00
      C(25)=0.3210906D+00
      C(26)=0.5659043D+00
      C(27)=0.2218115D+00
      E(28)=108.05630D+00
      E(29)=30.131576D+00
      E(30)=9.5300359D+00
      C(28)=0.1198237D+00
      C(29)=0.4875910D+00
      C(30)=0.5849874D+00
      E(31)=4.9626098D+00
      E(32)=1.7120829D+00
      E(33)=0.57719451D+00
      C(31)=0.2529487D+00
      C(32)=0.5727612D+00
      C(33)=0.3690387D+00
      RETURN
C
C           ANTIMONY
C
  500 CONTINUE
      E(1)=13289.383D+00
      E(2)=2010.5218D+00
      E(3)=441.69815D+00
      C(1)=0.0609843D+00
      C(2)=0.3668487D+00
      C(3)=0.6910501D+00
      E(4)=588.38790D+00
      E(5)=57.523192D+00
      E(6)=24.036693D+00
      C(4)=-0.1098365D+00
      C(5)=0.7203891D+00
      C(6)=0.3468674D+00
      E(7)=50.787743D+00
      E(8)=9.1473306D+00
      E(9)=4.1723695D+00
      C(7)=-0.2714342D+00
      C(8)=0.8547971D+00
      C(9)=0.2964847D+00
      E(10)=7.9995999D+00
      E(11)=1.8381194D+00
      E(12)=0.83126332D+00
      C(10)=-0.3512204D+00
      C(11)=0.8121005D+00
      C(12)=0.3980049D+00
      E(13)=1.3445260D+00
      E(14)=0.21146841D+00
      E(15)=0.08689142D+00
      C(13)=-0.2714728D+00
      C(14)=0.6682413D+00
      C(15)=0.4725305D+00
      E(16)=860.08072D+00
      E(17)=201.01421D+00
      E(18)=59.657596D+00
      C(16)=0.0809758D+00
      C(17)=0.4241487D+00
      C(18)=0.6387257D+00
      E(19)=122.27109D+00
      E(20)=20.929665D+00
      E(21)=8.0980350D+00
      C(19)=-0.0318112D+00
      C(20)=0.4045588D+00
      C(21)=0.6653749D+00
      E(22)=3.3003687D+00
      E(23)=1.4287737D+00
      E(24)=0.59505134D+00
      C(22)=0.3889519D+00
      C(23)=0.5640495D+00
      C(24)=0.1289287D+00
      E(25)=0.32801846D+00
      E(26)=0.13661293D+00
      E(27)=0.05411368D+00
      C(25)=0.3110910D+00
      C(26)=0.5652580D+00
      C(27)=0.2257358D+00
      E(28)=115.80955D+00
      E(29)=32.305835D+00
      E(30)=10.250328D+00
      C(28)=0.1166279D+00
      C(29)=0.4834363D+00
      C(30)=0.5901395D+00
      E(31)=5.4862102D+00
      E(32)=1.9216196D+00
      E(33)=0.66606265D+00
      C(31)=0.2483656D+00
      C(32)=0.5743154D+00
      C(33)=0.3643044D+00
      RETURN
C
C           TELLURIUM
C
  600 CONTINUE
      E(1)=13796.560D+00
      E(2)=2088.8798D+00
      E(3)=459.39319D+00
      C(1)=0.0610862D+00
      C(2)=0.3669629D+00
      C(3)=0.6907944D+00
      E(4)=611.58438D+00
      E(5)=60.221194D+00
      E(6)=25.491369D+00
      C(4)=-0.1101441D+00
      C(5)=0.7134762D+00
      C(6)=0.3534463D+00
      E(7)=53.042328D+00
      E(8)=9.6632025D+00
      E(9)=4.4188576D+00
      C(7)=-0.2728363D+00
      C(8)=0.8475640D+00
      C(9)=0.3052385D+00
      E(10)=8.4582805D+00
      E(11)=1.9753929D+00
      E(12)=0.90818331D+00
      C(10)=-0.3587862D+00
      C(11)=0.8112858D+00
      C(12)=0.4033116D+00
      E(13)=1.4270256D+00
      E(14)=0.23880222D+00
      E(15)=0.09319910D+00
      C(13)=-0.2806370D+00
      C(14)=0.7509597D+00
      C(15)=0.4038957D+00
      E(16)=895.07497D+00
      E(17)=209.57862D+00
      E(18)=62.405210D+00
      C(16)=0.0811099D+00
      C(17)=0.4236824D+00
      C(18)=0.6386030D+00
      E(19)=131.44799D+00
      E(20)=21.663180D+00
      E(21)=8.4667304D+00
      C(19)=-0.0320210D+00
      C(20)=0.4119187D+00
      C(21)=0.6571079D+00
      E(22)=3.4281D+00
      E(23)=1.4943404D+00
      E(24)=0.63848516D+00
      C(22)=0.4159312D+00
      C(23)=0.5489605D+00
      C(24)=0.1133174D+00
      E(25)=0.36614463D+00
      E(26)=0.14887871D+00
      E(27)=0.05976329D+00
      C(25)=0.3690254D+00
      C(26)=0.5355049D+00
      C(27)=0.2011813D+00
      E(28)=121.40830D+00
      E(29)=34.015217D+00
      E(30)=10.869138D+00
      C(28)=0.1169136D+00
      C(29)=0.4835555D+00
      C(30)=0.5883864D+00
      E(31)=5.8031113D+00
      E(32)=2.0580658D+00
      E(33)=0.73283021D+00
      C(31)=0.2601939D+00
      C(32)=0.5797758D+00
      C(33)=0.3405580D+00
      RETURN
C
C           IODINE
C
  700 CONTINUE
      E(1)=14351.186D+00
      E(2)=2173.0741D+00
      E(3)=477.87205D+00
      C(1)=0.0610028D+00
      C(2)=0.3666398D+00
      C(3)=0.6911306D+00
      E(4)=636.28295D+00
      E(5)=62.368944D+00
      E(6)=26.035562D+00
      C(4)=-0.1100819D+00
      C(5)=0.7219464D+00
      C(6)=0.3454775D+00
      E(7)=55.135577D+00
      E(8)=10.126452D+00
      E(9)=4.5928632D+00
      C(7)=-0.2754113D+00
      C(8)=0.8543504D+00
      C(9)=0.3006473D+00
      E(10)=8.9049791D+00
      E(11)=2.0925689D+00
      E(12)=0.9735899D+00
      C(10)=-0.3630477D+00
      C(11)=0.8279763D+00
      C(12)=0.3894209D+00
      E(13)=1.5539994D+00
      E(14)=0.27488087D+00
      E(15)=0.10826672D+00
      C(13)=0.2967653D+00
      C(14)=-0.7583335D+00
      C(15)=-0.4067087D+00
      E(16)=933.48281D+00
      E(17)=218.50559D+00
      E(18)=65.325776D+00
      C(16)=0.0810928D+00
      C(17)=0.4231510D+00
      C(18)=0.6385674D+00
      E(19)=158.48285D+00
      E(20)=21.215135D+00
      E(21)=8.4623149D+00
      C(19)=-0.0283465D+00
      C(20)=0.4622270D+00
      C(21)=0.6046503D+00
      E(22)=3.6295595D+00
      E(23)=1.6133740D+00
      E(24)=0.70733879D+00
      C(22)=0.4210620D+00
      C(23)=0.5427157D+00
      C(24)=0.1108784D+00
      E(25)=0.42307578D+00
      E(26)=0.17128391D+00
      E(27)=0.06847909D+00
      C(25)=0.3814674D+00
      C(26)=0.5321926D+00
      C(27)=0.1922032D+00
      E(28)=128.09026D+00
      E(29)=35.982378D+00
      E(30)=11.551116D+00
      C(28)=0.1158636D+00
      C(29)=0.4820494D+00
      C(30)=0.5894448D+00
      E(31)=6.1461523D+00
      E(32)=2.2209370D+00
      E(33)=0.80991202D+00
      C(31)=0.2681817D+00
      C(32)=0.5800614D+00
      C(33)=0.3262263D+00
      RETURN
C
C           XENON
C
  800 CONTINUE
      E(1)=14902.236D+00
      E(2)=2256.5383D+00
      E(3)=496.37317D+00
      C(1)=0.0609969D+00
      C(2)=0.3666290D+00
      C(3)=0.6911155D+00
      E(4)=660.94947D+00
      E(5)=64.863039D+00
      E(6)=27.111403D+00
      C(4)=-0.1102247D+00
      C(5)=0.7225832D+00
      C(6)=0.3448171D+00
      E(7)=57.515261D+00
      E(8)=10.536085D+00
      E(9)=4.7042710D+00
      C(7)=-0.2761463D+00
      C(8)=0.8678121D+00
      C(9)=0.2879474D+00
      E(10)=9.3786082D+00
      E(11)=2.2401181D+00
      E(12)=1.0615968D+00
      C(10)=-0.3673742D+00
      C(11)=0.8242752D+00
      C(12)=0.3963011D+00
      E(13)=1.6958079D+00
      E(14)=0.30651180D+00
      E(15)=0.12049354D+00
      C(13)=-0.3127449D+00
      C(14)=0.7841594D+00
      C(15)=0.3890762D+00
      E(16)=971.98851D+00
      E(17)=227.72243D+00
      E(18)=68.162563D+00
      C(16)=0.0809562D+00
      C(17)=0.4227290D+00
      C(18)=0.6388428D+00
      E(19)=164.53820D+00
      E(20)=22.210478D+00
      E(21)=8.90000556D+00
      C(19)=-0.0285155D+00
      C(20)=0.4626758D+00
      C(21)=0.6037699D+00
      E(22)=3.7649228D+00
      E(23)=1.6522375D+00
      E(24)=0.67649625D+00
      C(22)=0.4503453D+00
      C(23)=0.5400257D+00
      C(24)=0.0828350D+00
      E(25)=0.47842852D+00
      E(26)=0.19781294D+00
      E(27)=0.07915387D+00
      C(25)=-0.3900001D+00
      C(26)=-0.5233464D+00
      C(27)=-0.1906046D+00
      E(28)=134.91331D+00
      E(29)=37.956387D+00
      E(30)=12.227475D+00
      C(28)=0.1150105D+00
      C(29)=0.4815952D+00
      C(30)=0.5896133D+00
      E(31)=6.6004928D+00
      E(32)=2.3980513D+00
      E(33)=0.88648239D+00
      C(31)=0.2718844D+00
      C(32)=0.5855569D+00
      C(33)=0.3127456D+00
      RETURN
      END
C*MODULE BASHUZ  *DECK HTMTWO
      SUBROUTINE HTMTWO(E,C,NUCZ)
      IMPLICIT DOUBLE PRECISION(A-H,O-Z)
      DIMENSION E(30),C(30)
C
      IBR=NUCZ-38
      GO TO (100,200,300,400,500,600,700,800,900,1000), IBR
C
C     YTTRIUM (4-F)
C
  100 CONTINUE
      E( 1)=7647.6932D+00
      E( 2)=1157.0020D+00
      E( 3)=253.75700D+00
      C( 1)=0.0618818D+00
      C( 2)=0.3701570D+00
      C( 3)=0.6878017D+00
      E( 4)=339.10364D+00
      E( 5)=32.403264D+00
      E( 6)=13.281887D+00
      C( 4)=-0.1075786D+00
      C( 5)=0.7181663D+00
      C( 6)=0.3490842D+00
      E( 7)=27.996257D+00
      E( 8)=4.6829212D+00
      E( 9)=2.0566163D+00
      C( 7)=-0.2457049D+00
      C( 8)=0.7878323D+00
      C( 9)=0.3488907D+00
      E(10)=3.4864304D+00
      E(11)=0.66744514D+00
      E(12)=0.27959697D+00
      C(10)=-0.3000605D+00
      C(11)=0.7858715D+00
      C(12)=0.3857544D+00
      E(13)=0.46867231D+00
      E(14)=0.04978055D+00
      E(15)=0.02040820D+00
      C(13)=-0.2438980D+00
      C(14)=0.7576396D+00
      C(15)=0.3475036D+00
      E(16)=464.62389D+00
      E(17)=107.54806D+00
      E(18)=31.650408D+00
      C(16)=0.0870002D+00
      C(17)=0.4370123D+00
      C(18)=0.6250957D+00
      E(19)=12.443905D+00
      E(20)=4.9934558D+00
      E(21)=2.1033868D+00
      C(19)=0.2788493D+00
      C(20)=0.5885743D+00
      C(21)=0.2292208D+00
      E(22)=1.0947920D+00
      E(23)=0.46072626D+00
      E(24)=0.18428969D+00
      C(22)=0.4060140D+00
      C(23)=0.5320897D+00
      C(24)=0.1575146D+00
      E(25)=50.348619D+00
      E(26)=13.497089D+00
      E(27)=3.9166237D+00
      C(25)=0.1367063D+00
      C(26)=0.5029592D+00
      C(27)=0.5788425D+00
      E(28)=1.4046174D+00
      E(29)=0.45202594D+00
      E(30)=0.12497283D+00
      C(28)=0.1778019D+00
      C(29)=0.4652198D+00
      C(30)=0.5916440D+00
      RETURN
C
C     ZIRCONIUM (5-F)
C
  200 CONTINUE
      E( 1)=8074.5440D+00
      E( 2)=1220.5300D+00
      E( 3)=267.54751D+00
      C( 1)=0.0616606D+00
      C( 2)=0.3696140D+00
      C( 3)=0.6884549D+00
      E( 4)=357.68109D+00
      E( 5)=34.142568D+00
      E( 6)=13.947566D+00
      C( 4)=-0.1077339D+00
      C( 5)=0.7216851D+00
      C( 6)=0.3456921D+00
      E( 7)=29.610522D+00
      E( 8)=4.9809391D+00
      E( 9)=2.1850378D+00
      C( 7)=-0.2483862D+00
      C( 8)=0.8005982D+00
      C( 9)=0.3376774D+00
      E(10)=3.8399904D+00
      E(11)=0.73978700D+00
      E(12)=0.30961916D+00
      C(10)=-0.3104305D+00
      C(11)=0.8022937D+00
      C(12)=0.3741656D+00
      E(13)=0.55502132D+00
      E(14)=0.05532887D+00
      E(15)=0.02137166D+00
      C(13)=-0.2288857D+00
      C(14)=0.7817074D+00
      C(15)=0.3195193D+00
      E(16)=491.62198D+00
      E(17)=113.93040D+00
      E(18)=33.605515D+00
      C(16)=0.0867409D+00
      C(17)=0.4362728D+00
      C(18)=0.6255724D+00
      E(19)=13.247292D+00
      E(20)=5.3389914D+00
      E(21)=2.2525019D+00
      C(19)=0.2803844D+00
      C(20)=0.5922702D+00
      C(21)=0.2225147D+00
      E(22)=1.2314752D+00
      E(23)=0.52978846D+00
      E(24)=0.21408758D+00
      C(22)=0.4058518D+00
      C(23)=0.5246687D+00
      C(24)=0.1633953D+00
      E(25)=54.716068D+00
      E(26)=14.743888D+00
      E(27)=4.3333204D+00
      C(25)=0.1347642D+00
      C(26)=0.5014995D+00
      C(27)=0.5786696D+00
      E(28)=1.6859057D+00
      E(29)=0.53743517D+00
      E(30)=0.15022104D+00
      C(28)=0.1925337D+00
      C(29)=0.5052917D+00
      C(30)=0.5385693D+00
      RETURN
C
C     NIOBIUM (6-S)
C
  300 CONTINUE
      E(1)=8465.0874D+00
      E(2)=1281.1163D+00
      E(3)=281.21159D+00
      C(1)=0.0618144D+00
      C(2)=0.3698298D+00
      C(3)=0.6880474D+00
      E(4)=375.85025D+00
      E(5)=36.053370D+00
      E(6)=14.816642D+00
      C(4)=-0.1080131D+00
      C(5)=0.7199938D+00
      C(6)=0.3472657D+00
      E(7)=31.365510D+00
      E(8)=5.2718932D+00
      E(9)=2.3070864D+00
      C(7)=-0.2499774D+00
      C(8)=0.8156328D+00
      C(9)=0.3232937D+00
      E(10)=4.1677904D+00
      E(11)=0.8215584D+00
      E(12)=0.34324945D+00
      C(10)=-0.3163261D+00
      C(11)=0.8073422D+00
      C(12)=0.3740729D+00
      E(13)=0.58716698D+00
      E(14)=0.06616704D+00
      E(15)=0.02593466D+00
      C(13)=-0.2224673D+00
      C(14)=0.7048521D+00
      C(15)=0.4026700D+00
      E(16)=519.28752D+00
      E(17)=120.45899D+00
      E(18)=35.608105D+00
      C(16)=0.0865238D+00
      C(17)=0.4357148D+00
      C(18)=0.6258786D+00
      E(19)=13.976575D+00
      E(20)=5.6388612D+00
      E(21)=2.3712274D+00
      C(19)=0.2867436D+00
      C(20)=0.5975733D+00
      C(21)=0.2096306D+00
      E(22)=1.3455223D+00
      E(23)=0.56846749D+00
      E(24)=0.22696230D+00
      C(22)=0.4308749D+00
      C(23)=0.5215640D+00
      C(24)=0.1404843D+00
      E(25)=59.035804D+00
      E(26)=15.997120D+00
      E(27)=4.7574796D+00
      C(25)=0.1335596D+00
      C(26)=0.5006254D+00
      C(27)=0.5776791D+00
      E(28)=1.9792676D+00
      E(29)=0.62926194D+00
      E(30)=0.17752293D+00
      C(28)=0.2015902D+00
      C(29)=0.5271681D+00
      C(30)=0.5065745D+00
      RETURN
C
C     MOLYBDENUM (7-S)
C
  400 CONTINUE
      E(1)=8898.9752D+00
      E(2)=1346.6930D+00
      E(3)=295.62876D+00
      C(1)=0.0617106D+00
      C(2)=0.3694639D+00
      C(3)=0.6884201D+00
      E(4)=395.08162D+00
      E(5)=37.964578D+00
      E(6)=15.608898D+00
      C(4)=-0.1082198D+00
      C(5)=0.7205080D+00
      C(6)=0.3468176D+00
      E(7)=33.085228D+00
      E(8)=5.6180351D+00
      E(9)=2.4825318D+00
      C(7)=-0.2526233D+00
      C(8)=0.8177164D+00
      C(9)=0.3225675D+00
      E(10)=4.5000086D+00
      E(11)=0.90593940D+00
      E(12)=0.37777344D+00
      C(10)=-0.3206422D+00
      C(11)=0.8092539D+00
      C(12)=0.3763588D+00
      E(13)=0.65419456D+00
      E(14)=0.07258102D+00
      E(15)=0.02797177D+00
      C(13)=-0.2116815D+00
      C(14)=0.6864946D+00
      C(15)=0.4187549D+00
      E(16)=548.86388D+00
      E(17)=127.39268D+00
      E(18)=37.712842D+00
      C(16)=0.0860665D+00
      C(17)=0.4346698D+00
      C(18)=0.6268737D+00
      E(19)=14.744588D+00
      E(20)=5.9600351D+00
      E(21)=2.4937534D+00
      C(19)=0.2917026D+00
      C(20)=0.6020321D+00
      C(21)=0.1991546D+00
      E(22)=1.4621554D+00
      E(23)=0.60827575D+00
      E(24)=0.23912471D+00
      C(22)=0.4518315D+00
      C(23)=0.5172613D+00
      C(24)=0.1232182D+00
      E(25)=63.792046D+00
      E(26)=17.356419D+00
      E(27)=5.2109516D+00
      C(25)=0.1316317D+00
      C(26)=0.4990747D+00
      C(27)=0.5782636D+00
      E(28)=2.2861428D+00
      E(29)=0.72931951D+00
      E(30)=0.2094156D+00
      C(28)=0.2072519D+00
      C(29)=0.5396099D+00
      C(30)=0.4846239D+00
      RETURN
C
C     TECHNETIUM (6-D)
C
  500 CONTINUE
      E(1)=9331.8660D+00
      E(2)=1412.7444D+00
      E(3)=310.28294D+00
      C(1)=0.0616947D+00
      C(2)=0.3693039D+00
      C(3)=0.6885214D+00
      E(4)=414.56635D+00
      E(5)=39.944899D+00
      E(6)=16.451697D+00
      C(4)=-0.1084438D+00
      C(5)=0.7206125D+00
      C(6)=0.3467395D+00
      E(7)=34.887843D+00
      E(8)=5.9470507D+00
      E(9)=2.6240416D+00
      C(7)=-0.2545651D+00
      C(8)=0.8276066D+00
      C(9)=0.3138042D+00
      E(10)=4.8541288D+00
      E(11)=0.98792524D+00
      E(12)=0.41093096D+00
      C(10)=-0.3241571D+00
      C(11)=0.8151655D+00
      C(12)=0.3733329D+00
      E(13)=0.71023291D+00
      E(14)=0.07580316D+00
      E(15)=0.02868756D+00
      C(13)=-0.1967391D+00
      C(14)=0.6579450D+00
      C(15)=0.4429358D+00
      E(16)=580.94360D+00
      E(17)=135.36036D+00
      E(18)=39.889938D+00
      C(16)=0.0845326D+00
      C(17)=0.4321914D+00
      C(18)=0.6304900D+00
      E(19)=87.347195D+00
      E(20)=13.287886D+00
      E(21)=4.9090562D+00
      C(19)=-0.0295523D+00
      C(20)=0.4201419D+00
      C(21)=0.6547726D+00
      E(22)=1.9258524D+00
      E(23)=0.84385557D+00
      E(24)=0.33618541D+00
      C(22)=0.3011680D+00
      C(23)=0.5688683D+00
      C(24)=0.2252629D+00
      E(25)=69.016818D+00
      E(26)=18.830107D+00
      E(27)=5.6925288D+00
      C(25)=0.1290297D+00
      C(26)=0.4966878D+00
      C(27)=0.5805887D+00
      E(28)=2.6360931D+00
      E(29)=0.84019534D+00
      E(30)=0.23965135D+00
      C(28)=0.2078727D+00
      C(29)=0.5439701D+00
      C(30)=0.4804963D+00
      RETURN
C
C     RUTHENIUM (5-F)
C
  600 CONTINUE
      E(1)=9786.1429D+00
      E(2)=1481.4714D+00
      E(3)=325.40260D+00
      C(1)=0.0616049D+00
      C(2)=0.3689840D+00
      C(3)=0.6888455D+00
      E(4)=434.79530D+00
      E(5)=41.918783D+00
      E(6)=17.257795D+00
      C(4)=-0.1086056D+00
      C(5)=0.7222292D+00
      C(6)=0.3451630D+00
      E(7)=36.735183D+00
      E(8)=6.2780444D+00
      E(9)=2.7519759D+00
      C(7)=-0.2564513D+00
      C(8)=0.8399185D+00
      C(9)=0.3027171D+00
      E(10)=5.2173878D+00
      E(11)=1.0759512D+00
      E(12)=0.44710855D+00
      C(10)=-0.3265667D+00
      C(11)=0.8168014D+00
      C(12)=0.3742870D+00
      E(13)=0.79343625D+00
      E(14)=0.07742209D+00
      E(15)=0.02934842D+00
      C(13)=-0.1864361D+00
      C(14)=0.6339363D+00
      C(15)=0.4599697D+00
      E(16)=613.23511D+00
      E(17)=142.89492D+00
      E(18)=42.136712D+00
      C(16)=0.0839306D+00
      C(17)=0.4310037D+00
      C(18)=0.6318141D+00
      E(19)=91.243693D+00
      E(20)=14.138100D+00
      E(21)=5.2568481D+00
      C(19)=-0.0298490D+00
      C(20)=0.4184301D+00
      C(21)=0.6558081D+00
      E(22)=2.0542379D+00
      E(23)=0.88489893D+00
      E(24)=0.35005918D+00
      C(22)=0.3259420D+00
      C(23)=0.5693077D+00
      C(24)=0.2001975D+00
      E(25)=74.088526D+00
      E(26)=20.283721D+00
      E(27)=6.1778559D+00
      C(25)=0.1274352D+00
      C(26)=0.4953767D+00
      C(27)=0.5811668D+00
      E(28)=2.8983401D+00
      E(29)=0.92735408D+00
      E(30)=0.26618991D+00
      C(28)=0.2189550D+00
      C(29)=0.5502884D+00
      C(30)=0.4631699D+00
      RETURN
C
C     RHODIUM (4-F)
C
  700 CONTINUE
      E(1)=10225.249D+00
      E(2)=1549.1390D+00
      E(3)=340.56851D+00
      C(1)=0.0616846D+00
      C(2)=0.3690559D+00
      C(3)=0.6886519D+00
      E(4)=455.17966D+00
      E(5)=43.986950D+00
      E(6)=18.132961D+00
      C(4)=-0.1088006D+00
      C(5)=0.7224507D+00
      C(6)=0.3449743D+00
      E(7)=38.751877D+00
      E(8)=6.6037791D+00
      E(9)=2.8952705D+00
      C(7)=-0.2573715D+00
      C(8)=0.8517085D+00
      C(9)=0.2909376D+00
      E(10)=5.6152907D+00
      E(11)=1.1613908D+00
      E(12)=0.48020132D+00
      C(10)=-0.3266525D+00
      C(11)=0.8214561D+00
      C(12)=0.3704174D+00
      E(13)=0.85451427D+00
      E(14)=0.07947300D+00
      E(15)=0.03009295D+00
      C(13)=-0.1758253D+00
      C(14)=0.6050953D+00
      C(15)=0.4842742D+00
      E(16)=644.23213D+00
      E(17)=150.28013D+00
      E(18)=44.384507D+00
      C(16)=0.0837350D+00
      C(17)=0.4304508D+00
      C(18)=0.6321828D+00
      E(19)=95.833303D+00
      E(20)=14.963302D+00
      E(21)=5.6013156D+00
      C(19)=-0.0300917D+00
      C(20)=0.4189458D+00
      C(21)=0.6546474D+00
      E(22)=2.1513076D+00
      E(23)=0.90932095D+00
      E(24)=0.35848243D+00
      C(22)=0.3612858D+00
      C(23)=0.5622628D+00
      C(24)=0.1716521D+00
      E(25)=79.284881D+00
      E(26)=21.777632D+00
      E(27)=6.6786285D+00
      C(25)=0.1260477D+00
      C(26)=0.4942743D+00
      C(27)=0.5815451D+00
      E(28)=3.2225210D+00
      E(29)=1.0349052D+00
      E(30)=0.29821448D+00
      C(28)=0.2230469D+00
      C(29)=0.5533923D+00
      C(30)=0.4551932D+00
      RETURN
C
C     PALLADIUM (3-D)
C
  800 CONTINUE
      E(1)=10732.268D+00
      E(2)=1624.2557D+00
      E(3)=356.77443D+00
      C(1)=0.0614091D+00
      C(2)=0.3683372D+00
      C(3)=0.6895216D+00
      E(4)=476.45981D+00
      E(5)=46.089079D+00
      E(6)=19.011457D+00
      C(4)=-0.1089646D+00
      C(5)=0.7229184D+00
      C(6)=0.3445285D+00
      E(7)=40.501816D+00
      E(8)=7.0317180D+00
      E(9)=3.1102881D+00
      C(7)=-0.2610873D+00
      C(8)=0.8473020D+00
      C(9)=0.2980309D+00
      E(10)=5.9535538D+00
      E(11)=1.2737777D+00
      E(12)=0.52932056D+00
      C(10)=-0.3315604D+00
      C(11)=0.8065730D+00
      C(12)=0.3906202D+00
      E(13)=0.91680596D+00
      E(14)=0.08014658D+00
      E(15)=0.03033468D+00
      C(13)=-0.1665563D+00
      C(14)=0.5960441D+00
      C(15)=0.4889716D+00
      E(16)=676.16554D+00
      E(17)=157.85238D+00
      E(18)=46.699134D+00
      C(16)=0.0835399D+00
      C(17)=0.4299341D+00
      C(18)=0.6325113D+00
      E(19)=101.25606D+00
      E(20)=15.760330D+00
      E(21)=5.9416372D+00
      C(19)=-0.0303593D+00
      C(20)=0.4214283D+00
      C(21)=0.6515363D+00
      E(22)=2.3082437D+00
      E(23)=0.97075110D+00
      E(24)=0.38294812D+00
      C(22)=0.3727721D+00
      C(23)=0.5581097D+00
      C(24)=0.1642278D+00
      E(25)=84.456904D+00
      E(26)=23.277977D+00
      E(27)=7.1869656D+00
      C(25)=0.1251206D+00
      C(26)=0.4936143D+00
      C(27)=0.5812637D+00
      E(28)=3.4886548D+00
      E(29)=1.1266213D+00
      E(30)=0.32670438D+00
      C(28)=0.2326503D+00
      C(29)=0.5566155D+00
      C(30)=0.4419046D+00
      RETURN
C
C     SILVER (2-S)
C
  900 CONTINUE
      E(1)=11195.393D+00
      E(2)=1695.6114D+00
      E(3)=372.71564D+00
      C(1)=0.0614652D+00
      C(2)=0.3683407D+00
      C(3)=0.6894170D+00
      E(4)=498.08789D+00
      E(5)=48.189273D+00
      E(6)=19.879233D+00
      C(4)=-0.1090984D+00
      C(5)=0.7248356D+00
      C(6)=0.3425905D+00
      E(7)=42.617304D+00
      E(8)=7.3549843D+00
      E(9)=3.2193663D+00
      C(7)=-0.2616127D+00
      C(8)=0.8638930D+00
      C(9)=0.2814949D+00
      E(10)=6.4471843D+00
      E(11)=1.3406340D+00
      E(12)=0.55082491D+00
      C(10)=-0.3264484D+00
      C(11)=0.8244194D+00
      C(12)=0.3684570D+00
      E(13)=0.896024D+00
      E(14)=0.07873945D+00
      E(15)=0.02942569D+00
      C(13)=-0.1464338D+00
      C(14)=0.6078123D+00
      C(15)=0.474578D+00
      E(16)=709.79348D+00
      E(17)=165.75065D+00
      E(18)=49.079374D+00
      C(16)=0.0832016D+00
      C(17)=0.4292951D+00
      C(18)=0.6331456D+00
      E(19)=105.90189D+00
      E(20)=16.746072D+00
      E(21)=6.3433500D+00
      C(19)=-0.0305549D+00
      C(20)=0.4170359D+00
      C(21)=0.6552444D+00
      E(22)=2.511517D+00
      E(23)=1.0560803D+00
      E(24)=0.41595731D+00
      C(22)=0.3678284D+00
      C(23)=0.5598336D+00
      C(24)=0.1679555D+00
      E(25)=90.022978D+00
      E(26)=24.875706D+00
      E(27)=7.7224910D+00
      C(25)=0.1237838D+00
      C(26)=0.4925467D+00
      C(27)=0.5818767D+00
      E(28)=3.7854315D+00
      E(29)=1.2278562D+00
      E(30)=0.35902458D+00
      C(28)=0.2395459D+00
      C(29)=0.5591501D+00
      C(30)=0.4311417D+00
      RETURN
C
C     CADMIUM (1-S)
C
 1000 CONTINUE
      E(1)=11686.086D+00
      E(2)=1770.1114D+00
      E(3)=389.20899D+00
      C(1)=0.0614265D+00
      C(2)=0.3681567D+00
      C(3)=0.6895722D+00
      E(4)=519.91106D+00
      E(5)=50.397993D+00
      E(6)=20.834269D+00
      C(4)=-0.1092676D+00
      C(5)=0.7250389D+00
      C(6)=0.3423508D+00
      E(7)=44.652532D+00
      E(8)=7.7506717D+00
      E(9)=3.3818076D+00
      C(7)=-0.2634441D+00
      C(8)=0.8694815D+00
      C(9)=0.2772414D+00
      E(10)=6.7897003D+00
      E(11)=1.4664135D+00
      E(12)=0.61806063D+00
      C(10)=-0.3343146D+00
      C(11)=0.8177517D+00
      C(12)=0.3812846D+00
      E(13)=0.91773482D+00
      E(14)=0.11202343D+00
      E(15)=0.04082857D+00
      C(13)=-0.1509803D+00
      C(14)=0.6086784D+00
      C(15)=0.4891503D+00
      E(16)=745.10479D+00
      E(17)=174.04696D+00
      E(18)=51.552438D+00
      C(16)=0.0827347D+00
      C(17)=0.4283685D+00
      C(18)=0.6341921D+00
      E(19)=106.27841D+00
      E(20)=17.924992D+00
      E(21)=6.7904099D+00
      C(19)=-0.0305639D+00
      C(20)=0.4096720D+00
      C(21)=0.6626091D+00
      E(22)=2.5792466D+00
      E(23)=1.0573897D+00
      E(24)=0.40668464D+00
      C(22)=0.415516D+00
      C(23)=0.5559935D+00
      C(24)=0.1202681D+00
      E(25)=95.472743D+00
      E(26)=26.481959D+00
      E(27)=8.2828858D+00
      C(25)=0.1230828D+00
      C(26)=0.4916768D+00
      C(27)=0.5815408D+00
      E(28)=4.0821413D+00
      E(29)=1.3572792D+00
      E(30)=0.42083076D+00
      C(28)=0.2516043D+00
      C(29)=0.5652168D+00
      C(30)=0.3995299D+00
      RETURN
      END
