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 13 FEB 05 - MWS - PAD COMMON BLOCK NSHEL
C  4 FEB 05 - FRJ - ADD PC-N BASIS SETS FOR H,C,N,O,F,SI,P,S,CL.
C              NOTE THAT THE PC-4 BASIS SET CONTAINS H-FUNCTIONS,
C              WHICH GAMESS CURRENTLY CANNOT HANDLE. THESE FUNCTIONS
C              HAVE BEEN COMMENTED OUT AND THE HYDROGEN G-FUNCTIONS
C              HAS BEEN REMOVED AS WELL FOR CONSISTENCY
C*MODULE BASPCN  *DECK PCNBAS
      SUBROUTINE PCNBAS(NUCZ,CSINP,CPINP,CDINP,CFINP,CGINP,IERR1,IERR2,
     *                 INTYP,NANGM,NBFS,MINF,MAXF,LOC,NGAUSS,NS,NPOL)
C
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
C
      LOGICAL GOPARR,DSKWRK,MASWRK
      LOGICAL LDIFF
C
      DIMENSION CSINP(*),CPINP(*),CDINP(*),CFINP(*),CGINP(*),
     *          INTYP(*),NANGM(*),NBFS(*),MINF(*),MAXF(*),NS(*)
      DIMENSION EEX(90),CCS(90),CCP(90),CCD(90),CCF(90),CCG(90),
     *          ITYPES(32),NGBFS(32)
C
      PARAMETER (MXSH=5000, MXGTOT=20000, MXATM=2000)
C
      COMMON /INFOA / NAT,ICH,MUL,NUM,NQMT,NE,NA,NB,
     *                ZAN(MXATM),C(3,MXATM),IAN(MXATM)
      COMMON /NSHEL / EX(MXGTOT),CS(MXGTOT),CP(MXGTOT),CD(MXGTOT),
     *                CF(MXGTOT),CG(MXGTOT),CH(MXGTOT),CI(MXGTOT),
     *                KSTART(MXSH),KATOM(MXSH),KTYPE(MXSH),KNG(MXSH),
     *                KLOC(MXSH),KMIN(MXSH),KMAX(MXSH),NSHELL
      COMMON /OUTPUT/ NPRINT,ITOL,ICUT,NORMF,NORMP,NOPK
      COMMON /PAR   / ME,MASTER,NPROC,IBTYP,IPTIM,GOPARR,DSKWRK,MASWRK
C
      PARAMETER (ZERO=0.0D+00, ONE=1.0D+00, PT5=0.5D+00, PT75=0.75D+00)
      PARAMETER (PI32=5.56832799683170D+00, TOL=1.0D-10)
      PARAMETER (PT1875=1.875D+00, PT65625=6.5625D+00)
C
C     ----- SET UP THE PC-N BASIS SETS -----
C
      DO 100 I = 1,90
         EEX(I) = ZERO
         CCS(I) = ZERO
         CCP(I) = ZERO
         CCD(I) = ZERO
         CCF(I) = ZERO
         CCG(I) = ZERO
  100 CONTINUE
      DO 101 I = 1,32
         ITYPES(I) = 0
         NGBFS(I)  = 0
  101 CONTINUE
      LDIFF = .FALSE.
      NTMP = NPOL
      IF (NTMP.GE.10) THEN
         LDIFF = .TRUE.
         NTMP = NTMP - 10
      ENDIF
      ILAST = 0
C
C     ----- HYDROGEN TO HELIUM -----
C
      IF (NUCZ .GT. 2) GO TO 220
      CALL PCN001(EEX,CCS,CCP,CCD,CCF,NUCZ,NTMP,ILAST)
      IF (LDIFF) CALL APCN01(EEX,CCS,CCP,CCD,CCF,NUCZ,NTMP,ILAST)
      GO TO 300
C
C     ----- LITHIUM TO NEON -----
C
  220 CONTINUE
      IF (NUCZ.GT.10) GO TO 230
      CALL PCN002(EEX,CCS,CCP,CCD,CCF,CCG,NUCZ,NTMP,ILAST)
      IF (LDIFF) CALL APCN02(EEX,CCS,CCP,CCD,CCF,CCG,NUCZ,NTMP,ILAST)
      GO TO 300
C
C     ----- SODIUM TO ARGON -----
C
  230 CONTINUE
      IF(NUCZ.GT.18) GO TO 250
      CALL PCN003(EEX,CCS,CCP,CCD,CCF,CCG,NUCZ,NTMP,ILAST)
      IF (LDIFF) CALL APCN03(EEX,CCS,CCP,CCD,CCF,CCG,NUCZ,NTMP,ILAST)
      GO TO 300
C
  250 CONTINUE
      IF (MASWRK) WRITE(6,*)'PC-N BASIS SETS EXIST ONLY UP TO CL.'
      CALL BERROR(3)
      STOP
C
C     ----- NOW SET UP EACH SHELL -----
C
  300 CONTINUE
      CALL PCNSHL(NUCZ,MXPASS,ITYPES,NGBFS,NTMP,LDIFF)
C
      IGX = 0
      IPASS = 0
  310 CONTINUE
      IPASS = IPASS + 1
      IF(IPASS.GT.MXPASS) RETURN
C                         ******
      IGAUSS = NGBFS(IPASS)
      ITYP = ITYPES(IPASS)
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
      DO 720 I = 1,IGAUSS
         K = K1+I-1
         EX(K)    = EEX(IGX+I)
         CSINP(K) = CCS(IGX+I)
         CPINP(K) = CCP(IGX+I)
         CDINP(K) = CCD(IGX+I)
         CFINP(K) = CCF(IGX+I)
         CGINP(K) = CCG(IGX+I)
         CS(K) = CSINP(K)
         CP(K) = CPINP(K)
         CD(K) = CDINP(K)
         CF(K) = CFINP(K)
         CG(K) = CGINP(K)
  720 CONTINUE
      IGX = IGX + IGAUSS
C
C     ----- ALWAYS UNNORMALIZE PRIMITIVES -----
C
      DO 740 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*EE*EE)
         FACG = PT65625*FACS/(EE*EE*EE*EE)
         CS(K) = CS(K)/SQRT(FACS)
         CP(K) = CP(K)/SQRT(FACP)
         CD(K) = CD(K)/SQRT(FACD)
         CF(K) = CF(K)/SQRT(FACF)
         CG(K) = CG(K)/SQRT(FACG)
  740 CONTINUE
C
C     ----- IF(NORMF.EQ.0) NORMALIZE BASIS FUNCTIONS. -----
C
      IF (NORMF .EQ. 1) GO TO 820
      FACS = ZERO
      FACP = ZERO
      FACD = ZERO
      FACF = ZERO
      FACG = ZERO
      DO 780 IG = K1,K2
         DO 770 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)
            DUMG = PT65625*CG(IG)*CG(JG)/(EE*EE*EE*EE*FAC)
            IF (IG .NE. JG) THEN
               DUMS = DUMS+DUMS
               DUMP = DUMP+DUMP
               DUMD = DUMD+DUMD
               DUMF = DUMF+DUMF
               DUMG = DUMG+DUMG
            END IF
            FACS = FACS+DUMS
            FACP = FACP+DUMP
            FACD = FACD+DUMD
            FACF = FACF+DUMF
            FACG = FACG+DUMG
  770    CONTINUE
  780 CONTINUE
      IF (FACS .GT. TOL) FACS = ONE/SQRT(FACS*PI32)
      IF (FACP .GT. TOL) FACP = ONE/SQRT(FACP*PI32)
      IF (FACD .GT. TOL) FACD = ONE/SQRT(FACD*PI32)
CFRJ
      IF (FACF .GT. TOL) FACF = ONE/SQRT(FACF*PI32)
      IF (FACG .GT. TOL) FACG = ONE/SQRT(FACG*PI32)
      DO 800 IG = K1,K2
         CS(IG) = CS(IG) * FACS
         CP(IG) = CP(IG) * FACP
         CD(IG) = CD(IG) * FACD
         CF(IG) = CF(IG) * FACF
         CG(IG) = CG(IG) * FACG
         CSINP(IG) = CSINP(IG) * FACS
         CPINP(IG) = CPINP(IG) * FACP
         CDINP(IG) = CDINP(IG) * FACD
         CFINP(IG) = CFINP(IG) * FACF
         CGINP(IG) = CGINP(IG) * FACG
  800 CONTINUE
C
C        AND BRANCH BACK FOR THE NEXT SHELL
C
  820 CONTINUE
      GO TO 310
      END
C*MODULE BASPCN  *DECK PCNSHL
      SUBROUTINE PCNSHL(NUCZ,MXPASS,ITYPES,NGBFS,NTMP,LDIFF)
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      DIMENSION ITYPES(32),NGBFS(32)
      DIMENSION NPASS(5,3)
      DIMENSION ITYPS(32,5,3),NGS(32,5,3)
      LOGICAL LDIFF
C
CFRJ  NOTE THAT THE LAST VALUE IN EACH NPASS SEQUENCE SHOULD BE
C     INCREMENTED BY ONE WHEN H-FUNCTIONS ARE IMPLEMENTED, AND
C     CORRESPONDING CHANGES IN ITYPS AND NGS
C
C     NPASS: H ; C-F ; SI-CL
      DATA NPASS/2,3,6,12,18,   5,6,10,18,26,   7,8,12,18,24/
C
C     ITYPS DEFINES THE TYPE OF SHELL, 1=S,2=P,3=D,4=F,5=G
C     BEWARE 6 DOES NOT EQUAL H !!
C
      DATA ITYPS/2*1,30*0,   2*1,1*2,29*0,   3*1,2*2,1*3,26*0,
     *           5*1,4*2,2*3,1*4,20*0,  7*1,6*2,3*3,2*4,14*0,
     *           3*1,2*2,27*0,   3*1,2*2,1*3,26*0,
     *           4*1,3*2,2*3,1*4,22*0,
     *           6*1,5*2,4*3,2*4,1*5,14*0,  8*1,7*2,6*3,3*4,2*5,6*0,
     *           4*1,3*2,25*0,   4*1,3*2,1*3,24*0,
     *           5*1,4*2,2*3,1*4,20*0,
     *           6*1,5*2,4*3,2*4,1*5,14*0,  7*1,6*2,6*3,3*4,2*5,8*0/
C
C       NGS TELLS THE NUMBER OF GAUSSIANS IN EACH SHELL
C
      DATA NGS/  2,1,30*0,   3,1,1,29*0,   4,5*1,26*0,
     *           5,11*1,20*0,  5,17*1,14*0,
     *           2*3,1,2,1,27*0,   2*5,1,3,2*1,26*0,
     *           2*7,2*1,4,5*1,22*0,
     *           9,7,4*1,5,11*1,14*0,  11,8,6*1,5,17*1,6*0,
     *           3*5,1,2*4,1,25*0,   3*8,1,2*6,2*1,24*0,
     *           9,7,8,2*1,2*7,5*1,20*0,
     *           12,10,9,3*1,2*9,10*1,14*0,
     *           14,10,11,4*1,11,9,15*1,8*0/
C
      NTMP = NTMP + 1
      NROW = 99
      IF (NUCZ.LE.18) NROW = 3
      IF (NUCZ.LE.10) NROW = 2
      IF (NUCZ.LE.2)  NROW = 1
      MXPASS = NPASS(NTMP,NROW)
      DO 100 I=1,MXPASS
         ITYPES(I) = ITYPS(I,NTMP,NROW)
         NGBFS(I)  = NGS(I,NTMP,NROW)
C        A FEW UGLY HACKS FOR ADJUSTING NGBFS ACCORDING TO NUCZ.
C        THE DAVIDSON PURIFICATION TO PARTLY SEGMENT THE GENERAL
C        CONTRACTION MEANS THAT THE NUMBER OF FUNCTIONS IS NOT
C        NECESSARILY THE SAME FOR ALL ELEMENTS IN A GIVEN ROW
         IF (NUCZ.EQ.14 .AND. NTMP.EQ.2 .AND. I.EQ.3) NGBFS(I) = 7
         IF (NUCZ.EQ.14 .AND. NTMP.EQ.3 .AND. I.EQ.2) NGBFS(I) = 8
         IF (NUCZ.EQ.14 .AND. NTMP.EQ.3 .AND. I.EQ.3) NGBFS(I) = 9
         IF (NUCZ.EQ.15 .AND. NTMP.EQ.5 .AND. I.EQ.9) NGBFS(I) = 8
         IF (NUCZ.EQ.17 .AND. NTMP.EQ.2 .AND. I.EQ.6) NGBFS(I) = 5
         IF (NUCZ.EQ.17 .AND. NTMP.EQ.4 .AND. I.EQ.3) NGBFS(I) = 10
  100 CONTINUE
C     POSSIBLY ADJUST FOR DIFFUSE FUNCTIONS
      IF (LDIFF) THEN
         NEND = NTMP
         IF (NROW.GT.1) NEND = NEND + 1
C        ADJUST FOR THE MISSING H-FUNCTIONS
         IF (NTMP.EQ.5) NEND = NEND - 1
         DO I = 1,NEND
            ITYPES(MXPASS+1) = I
            NGBFS(MXPASS+1)  = 1
            MXPASS = MXPASS + 1
         ENDDO
      ENDIF
      RETURN
      END
C*MODULE BASPCN  *DECK PCN001
      SUBROUTINE PCN001(E,S,P,D,F,N,NPOL,ILAST)
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      DIMENSION E(90),S(90),P(90),D(90),F(90)
C
      LOGICAL GOPARR,DSKWRK,MASWRK
C
      COMMON /PAR   / ME,MASTER,NPROC,IBTYP,IPTIM,GOPARR,DSKWRK,MASWRK
C
      GO TO (100,120),N
C
C     ----- HYDROGEN -----
C
  100 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(1) = 0.43448D+01
         S(1) = 0.79299D-01
         E(2) = 0.66049D+00
         S(2) = 0.42422D+00
         E(3) = 0.13669D+00
         S(3) = 0.10000D+01
         ILAST = 3
      ELSE IF (NPOL.EQ.1) THEN
         E(1) = 0.12252D+02
         S(1) = 0.22822D-01
         E(2) = 0.18687D+01
         S(2) = 0.15564D+00
         E(3) = 0.41821D+00
         S(3) = 0.48898D+00
         E(4) = 0.10610D+00
         S(4) = 0.10000D+01
         E(5) = 0.10000D+01
         P(5) = 0.10000D+01
         ILAST = 5
      ELSE IF (NPOL.EQ.2) THEN
         E(1) = 0.75423D+02
         S(1) = 0.24065D-02
         E(2) = 0.11350D+02
         S(2) = 0.18487D-01
         E(3) = 0.25993D+01
         S(3) = 0.89742D-01
         E(4) = 0.73513D+00
         S(4) = 0.28111D+00
         E(5) = 0.23167D+00
         S(5) = 0.10000D+01
         E(6) = 0.74147D-01
         S(6) = 0.10000D+01
         E(7) = 0.16000D+01
         P(7) = 0.10000D+01
         E(8) = 0.45000D+00
         P(8) = 0.10000D+01
         E(9) = 0.12500D+01
         D(9) = 0.10000D+01
         ILAST = 9
      ELSE IF (NPOL.EQ.3) THEN
         E(1) = 0.77469D+03
         S(1) = 0.13163D-03
         E(2) = 0.11616D+03
         S(2) = 0.10226D-02
         E(3) = 0.26449D+02
         S(3) = 0.53629D-02
         E(4) = 0.75015D+01
         S(4) = 0.22285D-01
         E(5) = 0.24568D+01
         S(5) = 0.75374D-01
         E(6) = 0.88774D+00
         S(6) = 0.10000D+01
         E(7) = 0.34085D+00
         S(7) = 0.10000D+01
         E(8) = 0.13477D+00
         S(8) = 0.10000D+01
         E(9) = 0.52263D-01
         S(9) = 0.10000D+01
         E(10) = 0.65000D+01
         P(10) = 0.10000D+01
         E(11) = 0.15500D+01
         P(11) = 0.10000D+01
         E(12) = 0.55000D+00
         P(12) = 0.10000D+01
         E(13) = 0.20000D+00
         P(13) = 0.10000D+01
         E(14) = 0.19500D+01
         D(14) = 0.10000D+01
         E(15) = 0.65000D+00
         D(15) = 0.10000D+01
         E(16) = 0.22000D+01
         F(16) = 0.10000D+01
         ILAST = 16
      ELSE IF (NPOL.EQ.4) THEN
         E(1) = 0.30797D+04
         S(1) = 0.23473D-04
         E(2) = 0.46152D+03
         S(2) = 0.18245D-03
         E(3) = 0.10506D+03
         S(3) = 0.95933D-03
         E(4) = 0.29764D+02
         S(4) = 0.40461D-02
         E(5) = 0.97183D+01
         S(5) = 0.14607D-01
         E(6) = 0.35179D+01
         S(6) = 0.10000D+01
         E(7) = 0.13770D+01
         S(7) = 0.10000D+01
         E(8) = 0.56831D+00
         S(8) = 0.10000D+01
         E(9) = 0.24341D+00
         S(9) = 0.10000D+01
         E(10) = 0.10505D+00
         S(10) = 0.10000D+01
         E(11) = 0.44198D-01
         S(11) = 0.10000D+01
         E(12) = 0.30000D+02
         P(12) = 0.10000D+01
         E(13) = 0.66000D+01
         P(13) = 0.10000D+01
         E(14) = 0.21000D+01
         P(14) = 0.10000D+01
         E(15) = 0.95000D+00
         P(15) = 0.10000D+01
         E(16) = 0.45000D+00
         P(16) = 0.10000D+01
         E(17) = 0.15000D+00
         P(17) = 0.10000D+01
         E(18) = 0.21000D+01
         D(18) = 0.10000D+01
         E(19) = 0.70000D+00
         D(19) = 0.10000D+01
         E(20) = 0.30000D+00
         D(20) = 0.10000D+01
         E(21) = 0.25000D+01
         F(21) = 0.10000D+01
         E(22) = 0.90000D+00
         F(22) = 0.10000D+01
CFRJ     E(23) = 0.32000D+01
CFRJ     G(23) = 0.10000D+01
         ILAST = 22
      ENDIF
      RETURN
C
C     ----- HE -----
C
  120 CONTINUE
      IF (MASWRK) WRITE(6,*) 'PC-N BASIS DOES NOT EXIST FOR HE'
      CALL BERROR(3)
      STOP
      END
C*MODULE BASPCN  *DECK PCN002
      SUBROUTINE PCN002(E,S,P,D,F,G,N,NPOL,ILAST)
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      DIMENSION E(90),S(90),P(90),D(90),F(90),G(90)
C
      LOGICAL GOPARR,DSKWRK,MASWRK
C
      COMMON /PAR   / ME,MASTER,NPROC,IBTYP,IPTIM,GOPARR,DSKWRK,MASWRK
C
      NN = N-2
      GO TO (100,120,140,160,180,200,220,240),NN
C
  100 CONTINUE
  120 CONTINUE
  140 CONTINUE
      IF (MASWRK) WRITE(6,*) 'PC-N BASIS DOES NOT EXIST FOR LI,BE,B'
      CALL BERROR(3)
      STOP
C
C     ----- C  -----
C
  160 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(1) = 0.44651D+03
         S(1) = 0.19777D-01
         E(2) = 0.67462D+02
         S(2) = 0.14304D+00
         E(3) = 0.15067D+02
         S(3) = 0.20390D+00
         E(4) = 0.67462D+02
         S(4) = -.80909D-02
         E(5) = 0.15067D+02
         S(5) = 0.27790D+00
         E(6) = 0.39388D+01
         S(6) = 0.58957D+00
         E(7) = 0.29991D+00
         S(7) = 0.10000D+01
         E(8) = 0.43677D+01
         P(8) = 0.11350D+00
         E(9) = 0.86982D+00
         P(9) = 0.47074D+00
         E(10) = 0.18893D+00
         P(10) = 0.10000D+01
         ILAST = 10
      ELSE IF (NPOL.EQ.1) THEN
         E(1) = 0.12526D+04
         S(1) = 0.55734D-02
         E(2) = 0.18857D+03
         S(2) = 0.41492D-01
         E(3) = 0.42839D+02
         S(3) = 0.18263D+00
         E(4) = 0.11818D+02
         S(4) = 0.46118D+00
         E(5) = 0.35567D+01
         S(5) = 0.44940D+00
         E(6) = 0.18857D+03
         S(6) = 0.27745D-03
         E(7) = 0.42839D+02
         S(7) = 0.25602D-02
         E(8) = 0.11818D+02
         S(8) = 0.33485D-01
         E(9) = 0.35567D+01
         S(9) = 0.87579D-01
         E(10) = 0.54258D+00
         S(10) = -.53739D+00
         E(11) = 0.16058D+00
         S(11) = 0.10000D+01
         E(12) = 0.91426D+01
         P(12) = 0.44464D-01
         E(13) = 0.19298D+01
         P(13) = 0.22886D+00
         E(14) = 0.52522D+00
         P(14) = 0.51223D+00
         E(15) = 0.13608D+00
         P(15) = 0.10000D+01
         E(16) = 0.80000D+00
         D(16) = 0.10000D+01
         ILAST = 16
      ELSE IF (NPOL.EQ.2) THEN
         E(1) = 0.78571D+04
         S(1) = 0.56822D-03
         E(2) = 0.11787D+04
         S(2) = 0.43912D-02
         E(3) = 0.26832D+03
         S(3) = 0.22503D-01
         E(4) = 0.75948D+02
         S(4) = 0.86643D-01
         E(5) = 0.24559D+02
         S(5) = 0.24406D+00
         E(6) = 0.86212D+01
         S(6) = 0.44133D+00
         E(7) = 0.31278D+01
         S(7) = 0.35342D+00
         E(8) = 0.11787D+04
         S(8) = 0.96183D-06
         E(9) = 0.26832D+03
         S(9) = 0.62126D-04
         E(10) = 0.75948D+02
         S(10) = 0.77125D-03
         E(11) = 0.24559D+02
         S(11) = 0.73217D-02
         E(12) = 0.86212D+01
         S(12) = 0.39106D-01
         E(13) = 0.31278D+01
         S(13) = 0.89222D-01
         E(14) = 0.82202D+00
         S(14) = -.21859D+00
         E(15) = 0.33017D+00
         S(15) = 0.10000D+01
         E(16) = 0.11463D+00
         S(16) = 0.10000D+01
         E(17) = 0.33775D+02
         P(17) = 0.59791D-02
         E(18) = 0.76766D+01
         P(18) = 0.43276D-01
         E(19) = 0.22357D+01
         P(19) = 0.16125D+00
         E(20) = 0.76447D+00
         P(20) = 0.36373D+00
         E(21) = 0.26232D+00
         P(21) = 0.10000D+01
         E(22) = 0.84638D-01
         P(22) = 0.10000D+01
         E(23) = 0.14000D+01
         D(23) = 0.10000D+01
         E(24) = 0.45000D+00
         D(24) = 0.10000D+01
         E(25) = 0.95000D+00
         F(25) = 0.10000D+01
         ILAST = 25
      ELSE IF (NPOL.EQ.3) THEN
         E(1) = 0.87141D+05
         S(1) = 0.28180D-04
         E(2) = 0.13052D+05
         S(2) = 0.21825D-03
         E(3) = 0.29706D+04
         S(3) = 0.11508D-02
         E(4) = 0.84146D+03
         S(4) = 0.47884D-02
         E(5) = 0.27459D+03
         S(5) = 0.17126D-01
         E(6) = 0.99215D+02
         S(6) = 0.51160D-01
         E(7) = 0.38694D+02
         S(7) = 0.12794D+00
         E(8) = 0.15911D+02
         S(8) = 0.22811D+00
         E(9) = 0.67668D+01
         S(9) = 0.24690D+00
         E(10) = 0.84146D+03
         S(10) = 0.11375D-03
         E(11) = 0.27459D+03
         S(11) = 0.27554D-03
         E(12) = 0.99215D+02
         S(12) = 0.50008D-02
         E(13) = 0.38694D+02
         S(13) = 0.22297D-01
         E(14) = 0.15911D+02
         S(14) = 0.14027D+00
         E(15) = 0.67668D+01
         S(15) = 0.41571D+00
         E(16) = 0.28984D+01
         S(16) = 0.82977D+00
         E(17) = 0.11246D+01
         S(17) = 0.10000D+01
         E(18) = 0.51446D+00
         S(18) = 0.10000D+01
         E(19) = 0.21176D+00
         S(19) = 0.10000D+01
         E(20) = 0.85017D-01
         S(20) = 0.10000D+01
         E(21) = 0.14824D+03
         P(21) = 0.50719D-03
         E(22) = 0.34912D+02
         P(22) = 0.41794D-02
         E(23) = 0.10894D+02
         P(23) = 0.21293D-01
         E(24) = 0.41934D+01
         P(24) = 0.62269D-01
         E(25) = 0.17134D+01
         P(25) = 0.16415D+00
         E(26) = 0.73566D+00
         P(26) = 0.10000D+01
         E(27) = 0.31637D+00
         P(27) = 0.10000D+01
         E(28) = 0.13242D+00
         P(28) = 0.10000D+01
         E(29) = 0.51901D-01
         P(29) = 0.10000D+01
         E(30) = 0.52000D+01
         D(30) = 0.10000D+01
         E(31) = 0.13500D+01
         D(31) = 0.10000D+01
         E(32) = 0.52000D+00
         D(32) = 0.10000D+01
         E(33) = 0.20000D+00
         D(33) = 0.10000D+01
         E(34) = 0.15000D+01
         F(34) = 0.10000D+01
         E(35) = 0.65000D+00
         F(35) = 0.10000D+01
         E(36) = 0.14000D+01
         G(36) = 0.10000D+01
         ILAST = 36
      ELSE IF (NPOL.EQ.4) THEN
         E(1) = 0.52778D+06
         S(1) = 0.29660D-05
         E(2) = 0.79105D+05
         S(2) = 0.23134D-04
         E(3) = 0.18006D+05
         S(3) = 0.12091D-03
         E(4) = 0.51008D+04
         S(4) = 0.51430D-03
         E(5) = 0.16647D+04
         S(5) = 0.18471D-02
         E(6) = 0.60166D+03
         S(6) = 0.60209D-02
         E(7) = 0.23540D+03
         S(7) = 0.17055D-01
         E(8) = 0.98263D+02
         S(8) = 0.44034D-01
         E(9) = 0.43269D+02
         S(9) = 0.89625D-01
         E(10) = 0.19858D+02
         S(10) = 0.14845D+00
         E(11) = 0.95692D+01
         S(11) = 0.69092D-01
         E(12) = 0.16647D+04
         S(12) = 0.16830D-04
         E(13) = 0.60166D+03
         S(13) = 0.98080D-05
         E(14) = 0.23540D+03
         S(14) = 0.68920D-03
         E(15) = 0.98263D+02
         S(15) = 0.26428D-02
         E(16) = 0.43269D+02
         S(16) = 0.19725D-01
         E(17) = 0.19858D+02
         S(17) = 0.66280D-01
         E(18) = 0.95692D+01
         S(18) = 0.26732D+00
         E(19) = 0.49373D+01
         S(19) = 0.33642D+00
         E(20) = 0.26687D+01
         S(20) = 0.10000D+01
         E(21) = 0.14674D+01
         S(21) = 0.10000D+01
         E(22) = 0.63863D+00
         S(22) = 0.10000D+01
         E(23) = 0.30240D+00
         S(23) = 0.10000D+01
         E(24) = 0.14116D+00
         S(24) = 0.10000D+01
         E(25) = 0.63749D-01
         S(25) = 0.10000D+01
         E(26) = 0.38221D+03
         P(26) = 0.97576D-04
         E(27) = 0.90489D+02
         P(27) = 0.85864D-03
         E(28) = 0.28871D+02
         P(28) = 0.47882D-02
         E(29) = 0.10385D+02
         P(29) = 0.20388D-01
         E(30) = 0.44972D+01
         P(30) = 0.50005D-01
         E(31) = 0.20274D+01
         P(31) = 0.10000D+01
         E(32) = 0.95154D+00
         P(32) = 0.10000D+01
         E(33) = 0.44815D+00
         P(33) = 0.10000D+01
         E(34) = 0.20721D+00
         P(34) = 0.10000D+01
         E(35) = 0.94870D-01
         P(35) = 0.10000D+01
         E(36) = 0.40843D-01
         P(36) = 0.10000D+01
         E(37) = 0.84000D+01
         D(37) = 0.10000D+01
         E(38) = 0.30000D+01
         D(38) = 0.10000D+01
         E(39) = 0.12000D+01
         D(39) = 0.10000D+01
         E(40) = 0.55000D+00
         D(40) = 0.10000D+01
         E(41) = 0.25000D+00
         D(41) = 0.10000D+01
         E(42) = 0.10000D+00
         D(42) = 0.10000D+01
         E(43) = 0.21000D+01
         F(43) = 0.10000D+01
         E(44) = 0.90000D+00
         F(44) = 0.10000D+01
         E(45) = 0.40000D+00
         F(45) = 0.10000D+01
         E(46) = 0.21000D+01
         G(46) = 0.10000D+01
         E(47) = 0.90000D+00
         G(47) = 0.10000D+01
CFRJ     E(48) = 0.18500D+01
CFRJ     H(48) = 0.10000D+01
         ILAST = 47
      ENDIF
      RETURN
C
C     ----- N -----
C
  180 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(1) = 0.61585D+03
         S(1) = 0.19588D-01
         E(2) = 0.93044D+02
         S(2) = 0.14694D+00
         E(3) = 0.20848D+02
         S(3) = 0.15565D+00
         E(4) = 0.93044D+02
         S(4) = -.12052D-01
         E(5) = 0.20848D+02
         S(5) = 0.30201D+00
         E(6) = 0.54959D+01
         S(6) = 0.54102D+00
         E(7) = 0.42624D+00
         S(7) = 0.10000D+01
         E(8) = 0.62680D+01
         P(8) = 0.12135D+00
         E(9) = 0.12690D+01
         P(9) = 0.48013D+00
         E(10) = 0.27390D+00
         P(10) = 0.10000D+01
         ILAST = 10
      ELSE IF (NPOL.EQ.1) THEN
         E(1) = 0.17422D+04
         S(1) = 0.54603D-02
         E(2) = 0.26221D+03
         S(2) = 0.40706D-01
         E(3) = 0.59585D+02
         S(3) = 0.18038D+00
         E(4) = 0.16489D+02
         S(4) = 0.46009D+00
         E(5) = 0.49945D+01
         S(5) = 0.45173D+00
         E(6) = 0.26221D+03
         S(6) = 0.30104D-03
         E(7) = 0.59585D+02
         S(7) = 0.25511D-02
         E(8) = 0.16489D+02
         S(8) = 0.34900D-01
         E(9) = 0.49945D+01
         S(9) = 0.87429D-01
         E(10) = 0.78673D+00
         S(10) = -.52660D+00
         E(11) = 0.22837D+00
         S(11) = 0.10000D+01
         E(12) = 0.13132D+02
         P(12) = 0.45515D-01
         E(13) = 0.28108D+01
         P(13) = 0.23740D+00
         E(14) = 0.76351D+00
         P(14) = 0.51392D+00
         E(15) = 0.19560D+00
         P(15) = 0.10000D+01
         E(16) = 0.90000D+00
         D(16) = 0.10000D+01
         ILAST = 16
      ELSE IF (NPOL.EQ.2) THEN
         E(1) = 0.11101D+05
         S(1) = 0.54570D-03
         E(2) = 0.16653D+04
         S(2) = 0.42181D-02
         E(3) = 0.37908D+03
         S(3) = 0.21645D-01
         E(4) = 0.10730D+03
         S(4) = 0.83776D-01
         E(5) = 0.34747D+02
         S(5) = 0.23881D+00
         E(6) = 0.12249D+02
         S(6) = 0.43994D+00
         E(7) = 0.44743D+01
         S(7) = 0.36364D+00
         E(8) = 0.16653D+04
         S(8) = 0.15006D-05
         E(9) = 0.37908D+03
         S(9) = 0.58753D-04
         E(10) = 0.10730D+03
         S(10) = 0.77453D-03
         E(11) = 0.34747D+02
         S(11) = 0.72161D-02
         E(12) = 0.12249D+02
         S(12) = 0.39338D-01
         E(13) = 0.44743D+01
         S(13) = 0.90611D-01
         E(14) = 0.12404D+01
         S(14) = -.19882D+00
         E(15) = 0.48743D+00
         S(15) = 0.10000D+01
         E(16) = 0.16484D+00
         S(16) = 0.10000D+01
         E(17) = 0.47934D+02
         P(17) = 0.61706D-02
         E(18) = 0.10999D+02
         P(18) = 0.44260D-01
         E(19) = 0.32472D+01
         P(19) = 0.16825D+00
         E(20) = 0.11114D+01
         P(20) = 0.36848D+00
         E(21) = 0.37988D+00
         P(21) = 0.10000D+01
         E(22) = 0.12068D+00
         P(22) = 0.10000D+01
         E(23) = 0.18000D+01
         D(23) = 0.10000D+01
         E(24) = 0.55000D+00
         D(24) = 0.10000D+01
         E(25) = 0.10200D+01
         F(25) = 0.10000D+01
         ILAST = 25
      ELSE IF (NPOL.EQ.3) THEN
         E(1) = 0.12455D+06
         S(1) = 0.26672D-04
         E(2) = 0.18656D+05
         S(2) = 0.20684D-03
         E(3) = 0.42451D+04
         S(3) = 0.10893D-02
         E(4) = 0.12021D+04
         S(4) = 0.45463D-02
         E(5) = 0.39214D+03
         S(5) = 0.16237D-01
         E(6) = 0.14158D+03
         S(6) = 0.48957D-01
         E(7) = 0.55136D+02
         S(7) = 0.12317D+00
         E(8) = 0.22643D+02
         S(8) = 0.22425D+00
         E(9) = 0.96566D+01
         S(9) = 0.23876D+00
         E(10) = 0.12021D+04
         S(10) = 0.81439D-04
         E(11) = 0.39214D+03
         S(11) = 0.32439D-03
         E(12) = 0.14158D+03
         S(12) = 0.43621D-02
         E(13) = 0.55136D+02
         S(13) = 0.22838D-01
         E(14) = 0.22643D+02
         S(14) = 0.13462D+00
         E(15) = 0.96566D+01
         S(15) = 0.42501D+00
         E(16) = 0.41685D+01
         S(16) = 0.82517D+00
         E(17) = 0.16678D+01
         S(17) = 0.10000D+01
         E(18) = 0.74574D+00
         S(18) = 0.10000D+01
         E(19) = 0.30516D+00
         S(19) = 0.10000D+01
         E(20) = 0.12089D+00
         S(20) = 0.10000D+01
         E(21) = 0.21829D+03
         P(21) = 0.48281D-03
         E(22) = 0.51488D+02
         P(22) = 0.40420D-02
         E(23) = 0.16130D+02
         P(23) = 0.20946D-01
         E(24) = 0.61307D+01
         P(24) = 0.64757D-01
         E(25) = 0.25052D+01
         P(25) = 0.17042D+00
         E(26) = 0.10730D+01
         P(26) = 0.10000D+01
         E(27) = 0.45838D+00
         P(27) = 0.10000D+01
         E(28) = 0.18948D+00
         P(28) = 0.10000D+01
         E(29) = 0.73359D-01
         P(29) = 0.10000D+01
         E(30) = 0.60000D+01
         D(30) = 0.10000D+01
         E(31) = 0.16500D+01
         D(31) = 0.10000D+01
         E(32) = 0.65000D+00
         D(32) = 0.10000D+01
         E(33) = 0.25000D+00
         D(33) = 0.10000D+01
         E(34) = 0.18000D+01
         F(34) = 0.10000D+01
         E(35) = 0.77000D+00
         F(35) = 0.10000D+01
         E(36) = 0.15000D+01
         G(36) = 0.10000D+01
         ILAST = 36
      ELSE IF (NPOL.EQ.4) THEN
         E(1) = 0.89476D+06
         S(1) = 0.22632D-05
         E(2) = 0.13483D+06
         S(2) = 0.17541D-04
         E(3) = 0.30808D+05
         S(3) = 0.91168D-04
         E(4) = 0.87515D+04
         S(4) = 0.38740D-03
         E(5) = 0.28620D+04
         S(5) = 0.13870D-02
         E(6) = 0.10353D+04
         S(6) = 0.45503D-02
         E(7) = 0.40473D+03
         S(7) = 0.13009D-01
         E(8) = 0.16833D+03
         S(8) = 0.34696D-01
         E(9) = 0.73695D+02
         S(9) = 0.74804D-01
         E(10) = 0.33753D+02
         S(10) = 0.13457D+00
         E(11) = 0.16157D+02
         S(11) = 0.11672D+00
         E(12) = 0.28620D+04
         S(12) = 0.13703D-04
         E(13) = 0.10353D+04
         S(13) = -.14862D-04
         E(14) = 0.40473D+03
         S(14) = 0.48102D-03
         E(15) = 0.16833D+03
         S(15) = 0.15384D-02
         E(16) = 0.73695D+02
         S(16) = 0.13555D-01
         E(17) = 0.33753D+02
         S(17) = 0.47311D-01
         E(18) = 0.16157D+02
         S(18) = 0.20160D+00
         E(19) = 0.79688D+01
         S(19) = 0.40713D+00
         E(20) = 0.39521D+01
         S(20) = 0.10000D+01
         E(21) = 0.19454D+01
         S(21) = 0.10000D+01
         E(22) = 0.92929D+00
         S(22) = 0.10000D+01
         E(23) = 0.44085D+00
         S(23) = 0.10000D+01
         E(24) = 0.20340D+00
         S(24) = 0.10000D+01
         E(25) = 0.91226D-01
         S(25) = 0.10000D+01
         E(26) = 0.56623D+03
         P(26) = 0.92433D-04
         E(27) = 0.13410D+03
         P(27) = 0.80898D-03
         E(28) = 0.43013D+02
         P(28) = 0.45463D-02
         E(29) = 0.15612D+02
         P(29) = 0.19622D-01
         E(30) = 0.66427D+01
         P(30) = 0.51754D-01
         E(31) = 0.29748D+01
         P(31) = 0.10000D+01
         E(32) = 0.13885D+01
         P(32) = 0.10000D+01
         E(33) = 0.65012D+00
         P(33) = 0.10000D+01
         E(34) = 0.30028D+00
         P(34) = 0.10000D+01
         E(35) = 0.13620D+00
         P(35) = 0.10000D+01
         E(36) = 0.57845D-01
         P(36) = 0.10000D+01
         E(37) = 0.12800D+02
         D(37) = 0.10000D+01
         E(38) = 0.44000D+01
         D(38) = 0.10000D+01
         E(39) = 0.17500D+01
         D(39) = 0.10000D+01
         E(40) = 0.80000D+00
         D(40) = 0.10000D+01
         E(41) = 0.37000D+00
         D(41) = 0.10000D+01
         E(42) = 0.15000D+00
         D(42) = 0.10000D+01
         E(43) = 0.23000D+01
         F(43) = 0.10000D+01
         E(44) = 0.10000D+01
         F(44) = 0.10000D+01
         E(45) = 0.45000D+00
         F(45) = 0.10000D+01
         E(46) = 0.25500D+01
         G(46) = 0.10000D+01
         E(47) = 0.11000D+01
         G(47) = 0.10000D+01
CFRJ     E(48) = 0.19500D+01
CFRJ     H(48) = 0.10000D+01
         ILAST = 47
      ENDIF
      RETURN
C
C     ----- O  ------
C
  200 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(1) = 0.81212D+03
         S(1) = 0.19448D-01
         E(2) = 0.12269D+03
         S(2) = 0.14795D+00
         E(3) = 0.27561D+02
         S(3) = 0.13420D+00
         E(4) = 0.12269D+03
         S(4) = -.13358D-01
         E(5) = 0.27561D+02
         S(5) = 0.31565D+00
         E(6) = 0.73173D+01
         S(6) = 0.52789D+00
         E(7) = 0.57924D+00
         S(7) = 0.10000D+01
         E(8) = 0.82178D+01
         P(8) = 0.12882D+00
         E(9) = 0.16754D+01
         P(9) = 0.48612D+00
         E(10) = 0.35067D+00
         P(10) = 0.10000D+01
         ILAST = 10
      ELSE IF (NPOL.EQ.1) THEN
         E(1) = 0.23067D+04
         S(1) = 0.53938D-02
         E(2) = 0.34715D+03
         S(2) = 0.40244D-01
         E(3) = 0.78890D+02
         S(3) = 0.17921D+00
         E(4) = 0.21876D+02
         S(4) = 0.45968D+00
         E(5) = 0.66646D+01
         S(5) = 0.45241D+00
         E(6) = 0.34715D+03
         S(6) = 0.31840D-03
         E(7) = 0.78890D+02
         S(7) = 0.25952D-02
         E(8) = 0.21876D+02
         S(8) = 0.36317D-01
         E(9) = 0.66646D+01
         S(9) = 0.88039D-01
         E(10) = 0.10669D+01
         S(10) = -.53589D+00
         E(11) = 0.30700D+00
         S(11) = 0.10000D+01
         E(12) = 0.17022D+02
         P(12) = 0.48906D-01
         E(13) = 0.36838D+01
         P(13) = 0.24878D+00
         E(14) = 0.99234D+00
         P(14) = 0.51295D+00
         E(15) = 0.24487D+00
         P(15) = 0.10000D+01
         E(16) = 0.10000D+01
         D(16) = 0.10000D+01
         ILAST = 16
      ELSE IF (NPOL.EQ.2) THEN
         E(1) = 0.14782D+05
         S(1) = 0.53516D-03
         E(2) = 0.22173D+04
         S(2) = 0.41373D-02
         E(3) = 0.50474D+03
         S(3) = 0.21244D-01
         E(4) = 0.14287D+03
         S(4) = 0.82447D-01
         E(5) = 0.46300D+02
         S(5) = 0.23670D+00
         E(6) = 0.16337D+02
         S(6) = 0.44028D+00
         E(7) = 0.59828D+01
         S(7) = 0.36473D+00
         E(8) = 0.22173D+04
         S(8) = 0.21421D-05
         E(9) = 0.50474D+03
         S(9) = 0.57894D-04
         E(10) = 0.14287D+03
         S(10) = 0.80314D-03
         E(11) = 0.46300D+02
         S(11) = 0.73144D-02
         E(12) = 0.16337D+02
         S(12) = 0.40696D-01
         E(13) = 0.59828D+01
         S(13) = 0.91877D-01
         E(14) = 0.16718D+01
         S(14) = -.21057D+00
         E(15) = 0.64662D+00
         S(15) = 0.10000D+01
         E(16) = 0.21669D+00
         S(16) = 0.10000D+01
         E(17) = 0.60424D+02
         P(17) = 0.68724D-02
         E(18) = 0.13935D+02
         P(18) = 0.49025D-01
         E(19) = 0.41531D+01
         P(19) = 0.18189D+00
         E(20) = 0.14158D+01
         P(20) = 0.37559D+00
         E(21) = 0.47549D+00
         P(21) = 0.10000D+01
         E(22) = 0.14529D+00
         P(22) = 0.10000D+01
         E(23) = 0.22000D+01
         D(23) = 0.10000D+01
         E(24) = 0.65000D+00
         D(24) = 0.10000D+01
         E(25) = 0.11000D+01
         F(25) = 0.10000D+01
         ILAST = 25
      ELSE IF (NPOL.EQ.3) THEN
         E(1) = 0.16403D+06
         S(1) = 0.26520D-04
         E(2) = 0.24567D+05
         S(2) = 0.20574D-03
         E(3) = 0.55910D+04
         S(3) = 0.10826D-02
         E(4) = 0.15836D+04
         S(4) = 0.45202D-02
         E(5) = 0.51667D+03
         S(5) = 0.16132D-01
         E(6) = 0.18656D+03
         S(6) = 0.48733D-01
         E(7) = 0.72696D+02
         S(7) = 0.12262D+00
         E(8) = 0.29992D+02
         S(8) = 0.22209D+00
         E(9) = 0.12899D+02
         S(9) = 0.23890D+00
         E(10) = 0.15836D+04
         S(10) = 0.75682D-04
         E(11) = 0.51667D+03
         S(11) = 0.35568D-03
         E(12) = 0.18656D+03
         S(12) = 0.43918D-02
         E(13) = 0.72696D+02
         S(13) = 0.23820D-01
         E(14) = 0.29992D+02
         S(14) = 0.13830D+00
         E(15) = 0.12899D+02
         S(15) = 0.43028D+00
         E(16) = 0.56135D+01
         S(16) = 0.85179D+00
         E(17) = 0.22986D+01
         S(17) = 0.10000D+01
         E(18) = 0.99980D+00
         S(18) = 0.10000D+01
         E(19) = 0.40746D+00
         S(19) = 0.10000D+01
         E(20) = 0.15922D+00
         S(20) = 0.10000D+01
         E(21) = 0.27202D+03
         P(21) = 0.54881D-03
         E(22) = 0.64243D+02
         P(22) = 0.45733D-02
         E(23) = 0.20249D+02
         P(23) = 0.23534D-01
         E(24) = 0.76438D+01
         P(24) = 0.74362D-01
         E(25) = 0.31232D+01
         P(25) = 0.18659D+00
         E(26) = 0.13251D+01
         P(26) = 0.10000D+01
         E(27) = 0.55551D+00
         P(27) = 0.10000D+01
         E(28) = 0.22229D+00
         P(28) = 0.10000D+01
         E(29) = 0.84010D-01
         P(29) = 0.10000D+01
         E(30) = 0.70000D+01
         D(30) = 0.10000D+01
         E(31) = 0.20000D+01
         D(31) = 0.10000D+01
         E(32) = 0.80000D+00
         D(32) = 0.10000D+01
         E(33) = 0.30000D+00
         D(33) = 0.10000D+01
         E(34) = 0.21000D+01
         F(34) = 0.10000D+01
         E(35) = 0.90000D+00
         F(35) = 0.10000D+01
         E(36) = 0.16000D+01
         G(36) = 0.10000D+01
         ILAST = 36
      ELSE IF (NPOL.EQ.4) THEN
         E(1) = 0.12597D+07
         S(1) = 0.20683D-05
         E(2) = 0.19011D+06
         S(2) = 0.15946D-04
         E(3) = 0.43449D+05
         S(3) = 0.83352D-04
         E(4) = 0.12342D+05
         S(4) = 0.35168D-03
         E(5) = 0.40350D+04
         S(5) = 0.12720D-02
         E(6) = 0.14597D+04
         S(6) = 0.41256D-02
         E(7) = 0.57066D+03
         S(7) = 0.11983D-01
         E(8) = 0.23750D+03
         S(8) = 0.31496D-01
         E(9) = 0.10411D+03
         S(9) = 0.70151D-01
         E(10) = 0.47581D+02
         S(10) = 0.12667D+00
         E(11) = 0.22546D+02
         S(11) = 0.12061D+00
         E(12) = 0.40350D+04
         S(12) = 0.62153D-05
         E(13) = 0.14597D+04
         S(13) = 0.27247D-04
         E(14) = 0.57066D+03
         S(14) = 0.33552D-03
         E(15) = 0.23750D+03
         S(15) = 0.17960D-02
         E(16) = 0.10411D+03
         S(16) = 0.11237D-01
         E(17) = 0.47581D+02
         S(17) = 0.46776D-01
         E(18) = 0.22546D+02
         S(18) = 0.19054D+00
         E(19) = 0.11093D+02
         S(19) = 0.40519D+00
         E(20) = 0.55409D+01
         S(20) = 0.10000D+01
         E(21) = 0.27794D+01
         S(21) = 0.10000D+01
         E(22) = 0.12837D+01
         S(22) = 0.10000D+01
         E(23) = 0.59649D+00
         S(23) = 0.10000D+01
         E(24) = 0.27057D+00
         S(24) = 0.10000D+01
         E(25) = 0.12028D+00
         S(25) = 0.10000D+01
         E(26) = 0.70861D+03
         P(26) = 0.10433D-03
         E(27) = 0.16786D+03
         P(27) = 0.91174D-03
         E(28) = 0.53967D+02
         P(28) = 0.50886D-02
         E(29) = 0.19807D+02
         P(29) = 0.21544D-01
         E(30) = 0.83832D+01
         P(30) = 0.58879D-01
         E(31) = 0.37447D+01
         P(31) = 0.10000D+01
         E(32) = 0.17322D+01
         P(32) = 0.10000D+01
         E(33) = 0.79977D+00
         P(33) = 0.10000D+01
         E(34) = 0.36323D+00
         P(34) = 0.10000D+01
         E(35) = 0.15963D+00
         P(35) = 0.10000D+01
         E(36) = 0.65664D-01
         P(36) = 0.10000D+01
         E(37) = 0.18200D+02
         D(37) = 0.10000D+01
         E(38) = 0.60000D+01
         D(38) = 0.10000D+01
         E(39) = 0.24000D+01
         D(39) = 0.10000D+01
         E(40) = 0.11000D+01
         D(40) = 0.10000D+01
         E(41) = 0.50000D+00
         D(41) = 0.10000D+01
         E(42) = 0.20000D+00
         D(42) = 0.10000D+01
         E(43) = 0.26000D+01
         F(43) = 0.10000D+01
         E(44) = 0.11500D+01
         F(44) = 0.10000D+01
         E(45) = 0.50000D+00
         F(45) = 0.10000D+01
         E(46) = 0.30500D+01
         G(46) = 0.10000D+01
         E(47) = 0.13000D+01
         G(47) = 0.10000D+01
CFRJ     E(48) = 0.20500D+01
CFRJ     H(48) = 0.10000D+01
         ILAST = 47
      ENDIF
      RETURN
C
C     ----- F  -----
C
  220 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(1) = 0.10351D+04
         S(1) = 0.19350D-01
         E(2) = 0.15638D+03
         S(2) = 0.14900D+00
         E(3) = 0.35196D+02
         S(3) = 0.11613D+00
         E(4) = 0.15638D+03
         S(4) = -.14455D-01
         E(5) = 0.35196D+02
         S(5) = 0.32585D+00
         E(6) = 0.93921D+01
         S(6) = 0.51561D+00
         E(7) = 0.75342D+00
         S(7) = 0.10000D+01
         E(8) = 0.10599D+02
         P(8) = 0.13258D+00
         E(9) = 0.21733D+01
         P(9) = 0.48902D+00
         E(10) = 0.44871D+00
         P(10) = 0.10000D+01
         ILAST = 10
      ELSE IF (NPOL.EQ.1) THEN
         E(1) = 0.29546D+04
         S(1) = 0.53336D-02
         E(2) = 0.44462D+03
         S(2) = 0.39824D-01
         E(3) = 0.10105D+03
         S(3) = 0.17802D+00
         E(4) = 0.28068D+02
         S(4) = 0.45910D+00
         E(5) = 0.85837D+01
         S(5) = 0.45358D+00
         E(6) = 0.44462D+03
         S(6) = 0.33106D-03
         E(7) = 0.10105D+03
         S(7) = 0.26085D-02
         E(8) = 0.28068D+02
         S(8) = 0.37216D-01
         E(9) = 0.85837D+01
         S(9) = 0.88336D-01
         E(10) = 0.13986D+01
         S(10) = -.53303D+00
         E(11) = 0.39969D+00
         S(11) = 0.10000D+01
         E(12) = 0.21883D+02
         P(12) = 0.50238D-01
         E(13) = 0.47769D+01
         P(13) = 0.25374D+00
         E(14) = 0.12823D+01
         P(14) = 0.51223D+00
         E(15) = 0.31127D+00
         P(15) = 0.10000D+01
         E(16) = 0.11000D+01
         D(16) = 0.10000D+01
         ILAST = 16
      ELSE IF (NPOL.EQ.2) THEN
         E(1) = 0.19127D+05
         S(1) = 0.52247D-03
         E(2) = 0.28689D+04
         S(2) = 0.40398D-02
         E(3) = 0.65304D+03
         S(3) = 0.20758D-01
         E(4) = 0.18485D+03
         S(4) = 0.80796D-01
         E(5) = 0.59949D+02
         S(5) = 0.23360D+00
         E(6) = 0.21190D+02
         S(6) = 0.43957D+00
         E(7) = 0.77846D+01
         S(7) = 0.37047D+00
         E(8) = 0.28689D+04
         S(8) = 0.23704D-05
         E(9) = 0.65304D+03
         S(9) = 0.56239D-04
         E(10) = 0.18485D+03
         S(10) = 0.80276D-03
         E(11) = 0.59949D+02
         S(11) = 0.72679D-02
         E(12) = 0.21190D+02
         S(12) = 0.40910D-01
         E(13) = 0.77846D+01
         S(13) = 0.93109D-01
         E(14) = 0.22305D+01
         S(14) = -.20301D+00
         E(15) = 0.85409D+00
         S(15) = 0.10000D+01
         E(16) = 0.28336D+00
         S(16) = 0.10000D+01
         E(17) = 0.77639D+02
         P(17) = 0.70197D-02
         E(18) = 0.17965D+02
         P(18) = 0.50148D-01
         E(19) = 0.54005D+01
         P(19) = 0.18513D+00
         E(20) = 0.18423D+01
         P(20) = 0.37517D+00
         E(21) = 0.61532D+00
         P(21) = 0.10000D+01
         E(22) = 0.18489D+00
         P(22) = 0.10000D+01
         E(23) = 0.28000D+01
         D(23) = 0.10000D+01
         E(24) = 0.80000D+00
         D(24) = 0.10000D+01
         E(25) = 0.12000D+01
         F(25) = 0.10000D+01
         ILAST = 25
      ELSE IF (NPOL.EQ.3) THEN
         E(1) = 0.21566D+06
         S(1) = 0.25377D-04
         E(2) = 0.32299D+05
         S(2) = 0.19695D-03
         E(3) = 0.73504D+04
         S(3) = 0.10359D-02
         E(4) = 0.20818D+04
         S(4) = 0.43300D-02
         E(5) = 0.67909D+03
         S(5) = 0.15459D-01
         E(6) = 0.24510D+03
         S(6) = 0.46915D-01
         E(7) = 0.95462D+02
         S(7) = 0.11877D+00
         E(8) = 0.39422D+02
         S(8) = 0.21727D+00
         E(9) = 0.17018D+02
         S(9) = 0.23687D+00
         E(10) = 0.20818D+04
         S(10) = 0.65619D-04
         E(11) = 0.67909D+03
         S(11) = 0.35364D-03
         E(12) = 0.24510D+03
         S(12) = 0.41021D-02
         E(13) = 0.95462D+02
         S(13) = 0.23212D-01
         E(14) = 0.39422D+02
         S(14) = 0.13371D+00
         E(15) = 0.17018D+02
         S(15) = 0.42345D+00
         E(16) = 0.74525D+01
         S(16) = 0.86202D+00
         E(17) = 0.31329D+01
         S(17) = 0.10000D+01
         E(18) = 0.13395D+01
         S(18) = 0.10000D+01
         E(19) = 0.54218D+00
         S(19) = 0.10000D+01
         E(20) = 0.20904D+00
         S(20) = 0.10000D+01
         E(21) = 0.35122D+03
         P(21) = 0.55257D-03
         E(22) = 0.82978D+02
         P(22) = 0.46211D-02
         E(23) = 0.26267D+02
         P(23) = 0.23794D-01
         E(24) = 0.99307D+01
         P(24) = 0.76241D-01
         E(25) = 0.40670D+01
         P(25) = 0.18912D+00
         E(26) = 0.17246D+01
         P(26) = 0.10000D+01
         E(27) = 0.71976D+00
         P(27) = 0.10000D+01
         E(28) = 0.28624D+00
         P(28) = 0.10000D+01
         E(29) = 0.10789D+00
         P(29) = 0.10000D+01
         E(30) = 0.84000D+01
         D(30) = 0.10000D+01
         E(31) = 0.26000D+01
         D(31) = 0.10000D+01
         E(32) = 0.10000D+01
         D(32) = 0.10000D+01
         E(33) = 0.35000D+00
         D(33) = 0.10000D+01
         E(34) = 0.26000D+01
         F(34) = 0.10000D+01
         E(35) = 0.10500D+01
         F(35) = 0.10000D+01
         E(36) = 0.17000D+01
         G(36) = 0.10000D+01
         ILAST = 36
      ELSE IF (NPOL.EQ.4) THEN
         E(1) = 0.18560D+07
         S(1) = 0.17171D-05
         E(2) = 0.28007D+06
         S(2) = 0.13237D-04
         E(3) = 0.64078D+05
         S(3) = 0.69031D-04
         E(4) = 0.18216D+05
         S(4) = 0.29141D-03
         E(5) = 0.59578D+04
         S(5) = 0.10528D-02
         E(6) = 0.21550D+04
         S(6) = 0.34267D-02
         E(7) = 0.84196D+03
         S(7) = 0.99773D-02
         E(8) = 0.34982D+03
         S(8) = 0.26593D-01
         E(9) = 0.15279D+03
         S(9) = 0.60475D-01
         E(10) = 0.69367D+02
         S(10) = 0.11449D+00
         E(11) = 0.32574D+02
         S(11) = 0.11554D+00
         E(12) = 0.59578D+04
         S(12) = 0.52091D-05
         E(13) = 0.21550D+04
         S(13) = 0.14902D-04
         E(14) = 0.84196D+03
         S(14) = 0.25857D-03
         E(15) = 0.34982D+03
         S(15) = 0.12838D-02
         E(16) = 0.15279D+03
         S(16) = 0.87364D-02
         E(17) = 0.69367D+02
         S(17) = 0.36990D-01
         E(18) = 0.32574D+02
         S(18) = 0.16180D+00
         E(19) = 0.15889D+02
         S(19) = 0.37507D+00
         E(20) = 0.78651D+01
         S(20) = 0.10000D+01
         E(21) = 0.39199D+01
         S(21) = 0.10000D+01
         E(22) = 0.17861D+01
         S(22) = 0.10000D+01
         E(23) = 0.82065D+00
         S(23) = 0.10000D+01
         E(24) = 0.36937D+00
         S(24) = 0.10000D+01
         E(25) = 0.16198D+00
         S(25) = 0.10000D+01
         E(26) = 0.91789D+03
         P(26) = 0.10430D-03
         E(27) = 0.21741D+03
         P(27) = 0.91436D-03
         E(28) = 0.69939D+02
         P(28) = 0.51142D-02
         E(29) = 0.25831D+02
         P(29) = 0.21545D-01
         E(30) = 0.10969D+02
         P(30) = 0.59682D-01
         E(31) = 0.49215D+01
         P(31) = 0.10000D+01
         E(32) = 0.22839D+01
         P(32) = 0.10000D+01
         E(33) = 0.10579D+01
         P(33) = 0.10000D+01
         E(34) = 0.48135D+00
         P(34) = 0.10000D+01
         E(35) = 0.20902D+00
         P(35) = 0.10000D+01
         E(36) = 0.85556D-01
         P(36) = 0.10000D+01
         E(37) = 0.25500D+02
         D(37) = 0.10000D+01
         E(38) = 0.80000D+01
         D(38) = 0.10000D+01
         E(39) = 0.32000D+01
         D(39) = 0.10000D+01
         E(40) = 0.14500D+01
         D(40) = 0.10000D+01
         E(41) = 0.65000D+00
         D(41) = 0.10000D+01
         E(42) = 0.25000D+00
         D(42) = 0.10000D+01
         E(43) = 0.31000D+01
         F(43) = 0.10000D+01
         E(44) = 0.13500D+01
         F(44) = 0.10000D+01
         E(45) = 0.60000D+00
         F(45) = 0.10000D+01
         E(46) = 0.36000D+01
         G(46) = 0.10000D+01
         E(47) = 0.15500D+01
         G(47) = 0.10000D+01
CFRJ     E(48) = 0.21500D+01
CFRJ     H(48) = 0.10000D+01
         ILAST = 47
      ENDIF
      RETURN
C
C     ----- NE -----
C
  240 CONTINUE
      IF (MASWRK) WRITE(6,*) 'PC-N BASIS DOES NOT EXIST FOR NE.'
      CALL BERROR(3)
      STOP
      END
C*MODULE BASPCN  *DECK PCN003
      SUBROUTINE PCN003(E,S,P,D,F,G,N,NPOL,ILAST)
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      DIMENSION E(90),S(90),P(90),D(90),F(90),G(90)
      LOGICAL GOPARR,DSKWRK,MASWRK
C
      COMMON /PAR   / ME,MASTER,NPROC,IBTYP,IPTIM,GOPARR,DSKWRK,MASWRK
C
      NN = N-10
      GO TO (100,120,140,160,180,200,220,240),NN
C
  100 CONTINUE
  120 CONTINUE
  140 CONTINUE
      IF (MASWRK) WRITE(6,*) 'PC-N BASIS DOES NOT EXIST FOR NA,MG,AL.'
      CALL BERROR(3)
      STOP
C
C     ----- SI -----
C
  160 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(1) = 0.77430D+04
         S(1) = 0.48825D-02
         E(2) = 0.11647D+04
         S(2) = 0.36741D-01
         E(3) = 0.26482D+03
         S(3) = 0.16695D+00
         E(4) = 0.73968D+02
         S(4) = 0.45657D+00
         E(5) = 0.22912D+02
         S(5) = 0.49043D+00
         E(6) = 0.11647D+04
         S(6) = -.29111D-03
         E(7) = 0.26482D+03
         S(7) = -.14122D-02
         E(8) = 0.73968D+02
         S(8) = -.26046D-01
         E(9) = 0.22912D+02
         S(9) = -.50933D-01
         E(10) = 0.43741D+01
         S(10) = 0.10810D+00
         E(11) = 0.26482D+03
         S(11) = 0.44437D-03
         E(12) = 0.73968D+02
         S(12) = 0.38067D-02
         E(13) = 0.22912D+02
         S(13) = 0.17078D-01
         E(14) = 0.43741D+01
         S(14) = -.10569D+00
         E(15) = 0.15391D+01
         S(15) = -.20429D+00
         E(16) = 0.16775D+00
         S(16) = 0.10000D+01
         E(17) = 0.96080D+02
         P(17) = 0.28035D-01
         E(18) = 0.22150D+02
         P(18) = 0.16819D+00
         E(19) = 0.65630D+01
         P(19) = 0.45197D+00
         E(20) = 0.20250D+01
         P(20) = 0.50132D+00
         E(21) = 0.22150D+02
         P(21) = -.72797D-03
         E(22) = 0.65630D+01
         P(22) = -.86644D-02
         E(23) = 0.20250D+01
         P(23) = -.35344D-02
         E(24) = 0.47780D+00
         P(24) = 0.31085D+00
         E(25) = 0.12824D+00
         P(25) = 0.10000D+01
         ILAST = 25
      ELSE IF (NPOL.EQ.1) THEN
         E(1) = 0.41263D+05
         S(1) = 0.60948D-03
         E(2) = 0.61896D+04
         S(2) = 0.47040D-02
         E(3) = 0.14087D+04
         S(3) = 0.24136D-01
         E(4) = 0.39877D+03
         S(4) = 0.92762D-01
         E(5) = 0.12969D+03
         S(5) = 0.26047D+00
         E(6) = 0.46425D+02
         S(6) = 0.44181D+00
         E(7) = 0.17455D+02
         S(7) = 0.30755D+00
         E(8) = 0.42222D+01
         S(8) = 0.29197D-01
         E(9) = 0.61896D+04
         S(9) = -.22459D-04
         E(10) = 0.14087D+04
         S(10) = -.41760D-04
         E(11) = 0.39877D+03
         S(11) = -.18570D-02
         E(12) = 0.12969D+03
         S(12) = -.95113D-02
         E(13) = 0.46425D+02
         S(13) = -.65800D-01
         E(14) = 0.17455D+02
         S(14) = -.77546D-01
         E(15) = 0.42222D+01
         S(15) = 0.50303D+00
         E(16) = 0.14933D+01
         S(16) = 0.63824D+00
         E(17) = 0.39877D+03
         S(17) = 0.74278D-04
         E(18) = 0.12969D+03
         S(18) = 0.74023D-03
         E(19) = 0.46425D+02
         S(19) = 0.40107D-02
         E(20) = 0.17455D+02
         S(20) = 0.80137D-02
         E(21) = 0.42222D+01
         S(21) = -.74091D-01
         E(22) = 0.14933D+01
         S(22) = -.21382D+00
         E(23) = 0.27174D+00
         S(23) = 0.61263D+00
         E(24) = 0.98502D-01
         S(24) = 0.10000D+01
         E(25) = 0.31073D+03
         P(25) = 0.41198D-02
         E(26) = 0.72823D+02
         P(26) = 0.31981D-01
         E(27) = 0.22985D+02
         P(27) = 0.13224D+00
         E(28) = 0.83213D+01
         P(28) = 0.32779D+00
         E(29) = 0.31353D+01
         P(29) = 0.45107D+00
         E(30) = 0.11693D+01
         P(30) = 0.26345D+00
         E(31) = 0.72823D+02
         P(31) = 0.12637D-03
         E(32) = 0.22985D+02
         P(32) = -.10252D-02
         E(33) = 0.83213D+01
         P(33) = -.21337D-02
         E(34) = 0.31353D+01
         P(34) = -.18774D-01
         E(35) = 0.11693D+01
         P(35) = 0.57528D-01
         E(36) = 0.29281D+00
         P(36) = 0.53571D+00
         E(37) = 0.82704D-01
         P(37) = 0.10000D+01
         E(38) = 0.45000D+00
         D(38) = 0.10000D+01
         ILAST = 38
      ELSE IF (NPOL.EQ.2) THEN
         E(1) = 0.12004D+06
         S(1) = 0.16075D-03
         E(2) = 0.17991D+05
         S(2) = 0.12472D-02
         E(3) = 0.40948D+04
         S(3) = 0.65040D-02
         E(4) = 0.11596D+04
         S(4) = 0.26665D-01
         E(5) = 0.37800D+03
         S(5) = 0.88816D-01
         E(6) = 0.13593D+03
         S(6) = 0.22932D+00
         E(7) = 0.52411D+02
         S(7) = 0.40025D+00
         E(8) = 0.20927D+02
         S(8) = 0.33828D+00
         E(9) = 0.77130D+01
         S(9) = 0.65512D-01
         E(10) = 0.40948D+04
         S(10) = -.21124D-05
         E(11) = 0.11596D+04
         S(11) = -.69157D-04
         E(12) = 0.37800D+03
         S(12) = -.48080D-03
         E(13) = 0.13593D+03
         S(13) = -.41694D-02
         E(14) = 0.52411D+02
         S(14) = -.17674D-01
         E(15) = 0.20927D+02
         S(15) = -.42891D-01
         E(16) = 0.77130D+01
         S(16) = 0.45368D-01
         E(17) = 0.31604D+01
         S(17) = 0.13904D+00
         E(18) = 0.40948D+04
         S(18) = 0.39897D-05
         E(19) = 0.11596D+04
         S(19) = 0.50203D-04
         E(20) = 0.37800D+03
         S(20) = 0.57553D-03
         E(21) = 0.13593D+03
         S(21) = 0.41524D-02
         E(22) = 0.52411D+02
         S(22) = 0.20437D-01
         E(23) = 0.20927D+02
         S(23) = 0.46998D-01
         E(24) = 0.77130D+01
         S(24) = -.43986D-01
         E(25) = 0.31604D+01
         S(25) = -.33957D+00
         E(26) = 0.12348D+01
         S(26) = -.35038D+00
         E(27) = 0.26775D+00
         S(27) = 0.10000D+01
         E(28) = 0.94067D-01
         S(28) = 0.10000D+01
         E(29) = 0.67713D+03
         P(29) = 0.10925D-02
         E(30) = 0.16067D+03
         P(30) = 0.89601D-02
         E(31) = 0.51585D+02
         P(31) = 0.44744D-01
         E(32) = 0.18948D+02
         P(32) = 0.14748D+00
         E(33) = 0.76163D+01
         P(33) = 0.31480D+00
         E(34) = 0.31317D+01
         P(34) = 0.41339D+00
         E(35) = 0.12703D+01
         P(35) = 0.26440D+00
         E(36) = 0.16067D+03
         P(36) = 0.79014D-05
         E(37) = 0.51585D+02
         P(37) = -.13613D-03
         E(38) = 0.18948D+02
         P(38) = -.90378D-03
         E(39) = 0.76163D+01
         P(39) = -.55863D-02
         E(40) = 0.31317D+01
         P(40) = -.10141D-01
         E(41) = 0.12703D+01
         P(41) = 0.18661D-01
         E(42) = 0.43332D+00
         P(42) = 0.28890D+00
         E(43) = 0.16088D+00
         P(43) = 0.10000D+01
         E(44) = 0.54883D-01
         P(44) = 0.10000D+01
         E(45) = 0.16800D+01
         D(45) = 0.10000D+01
         E(46) = 0.38000D+00
         D(46) = 0.10000D+01
         E(47) = 0.54000D+00
         F(47) = 0.10000D+01
         ILAST = 47
      ELSE IF (NPOL.EQ.3) THEN
         E(1) = 0.91848D+06
         S(1) = 0.12650D-04
         E(2) = 0.13759D+06
         S(2) = 0.98336D-04
         E(3) = 0.31289D+05
         S(3) = 0.51783D-03
         E(4) = 0.88507D+04
         S(4) = 0.21853D-02
         E(5) = 0.28819D+04
         S(5) = 0.79183D-02
         E(6) = 0.10381D+04
         S(6) = 0.25222D-01
         E(7) = 0.40397D+03
         S(7) = 0.70531D-01
         E(8) = 0.16697D+03
         S(8) = 0.16612D+00
         E(9) = 0.72202D+02
         S(9) = 0.30560D+00
         E(10) = 0.32334D+02
         S(10) = 0.35748D+00
         E(11) = 0.14771D+02
         S(11) = 0.19583D+00
         E(12) = 0.59240D+01
         S(12) = 0.23686D-01
         E(13) = 0.88507D+04
         S(13) = -.40074D-05
         E(14) = 0.28819D+04
         S(14) = -.55577D-05
         E(15) = 0.10381D+04
         S(15) = -.15833D-03
         E(16) = 0.40397D+03
         S(16) = -.68701D-03
         E(17) = 0.16697D+03
         S(17) = -.51741D-02
         E(18) = 0.72202D+02
         S(18) = -.19323D-01
         E(19) = 0.32334D+02
         S(19) = -.60309D-01
         E(20) = 0.14771D+02
         S(20) = -.51880D-01
         E(21) = 0.59240D+01
         S(21) = 0.17152D+00
         E(22) = 0.26285D+01
         S(22) = 0.30338D+00
         E(23) = 0.10381D+04
         S(23) = 0.93287D-05
         E(24) = 0.40397D+03
         S(24) = 0.67519D-04
         E(25) = 0.16697D+03
         S(25) = 0.40170D-03
         E(26) = 0.72202D+02
         S(26) = 0.18138D-02
         E(27) = 0.32334D+02
         S(27) = 0.55784D-02
         E(28) = 0.14771D+02
         S(28) = 0.69224D-02
         E(29) = 0.59240D+01
         S(29) = -.26192D-01
         E(30) = 0.26285D+01
         S(30) = -.12508D+00
         E(31) = 0.11360D+01
         S(31) = -.17545D+00
         E(32) = 0.35611D+00
         S(32) = 0.10000D+01
         E(33) = 0.16384D+00
         S(33) = 0.10000D+01
         E(34) = 0.68062D-01
         S(34) = 0.10000D+01
         E(35) = 0.25529D+04
         P(35) = 0.10925D-03
         E(36) = 0.60606D+03
         P(36) = 0.95808D-03
         E(37) = 0.19555D+03
         P(37) = 0.54477D-02
         E(38) = 0.72844D+02
         P(38) = 0.23607D-01
         E(39) = 0.30401D+02
         P(39) = 0.72280D-01
         E(40) = 0.13565D+02
         P(40) = 0.17561D+00
         E(41) = 0.62479D+01
         P(41) = 0.30405D+00
         E(42) = 0.28977D+01
         P(42) = 0.35739D+00
         E(43) = 0.13321D+01
         P(43) = 0.23612D+00
         E(44) = 0.60606D+03
         P(44) = 0.62893D-05
         E(45) = 0.19555D+03
         P(45) = -.38321D-05
         E(46) = 0.72844D+02
         P(46) = 0.10549D-03
         E(47) = 0.30401D+02
         P(47) = -.41779D-03
         E(48) = 0.13565D+02
         P(48) = -.51137D-03
         E(49) = 0.62479D+01
         P(49) = -.70295D-02
         E(50) = 0.28977D+01
         P(50) = -.43975D-02
         E(51) = 0.13321D+01
         P(51) = 0.51784D-02
         E(52) = 0.58716D+00
         P(52) = 0.13923D+00
         E(53) = 0.26362D+00
         P(53) = 0.10000D+01
         E(54) = 0.10828D+00
         P(54) = 0.10000D+01
         E(55) = 0.41189D-01
         P(55) = 0.10000D+01
         E(56) = 0.66600D+01
         D(56) = 0.10000D+01
         E(57) = 0.17500D+01
         D(57) = 0.10000D+01
         E(58) = 0.52000D+00
         D(58) = 0.10000D+01
         E(59) = 0.19000D+00
         D(59) = 0.10000D+01
         E(60) = 0.11500D+01
         F(60) = 0.10000D+01
         E(61) = 0.38000D+00
         F(61) = 0.10000D+01
         E(62) = 0.65000D+00
         G(62) = 0.10000D+01
         ILAST = 62
      ELSE IF (NPOL.EQ.4) THEN
         E(1) = 0.71901D+06
         S(1) = 0.12134D-04
         E(2) = 0.16999D+06
         S(2) = 0.60969D-04
         E(3) = 0.49400D+05
         S(3) = 0.25124D-03
         E(4) = 0.16385D+05
         S(4) = 0.90088D-03
         E(5) = 0.59854D+04
         S(5) = 0.29039D-02
         E(6) = 0.23541D+04
         S(6) = 0.85852D-02
         E(7) = 0.98140D+03
         S(7) = 0.23421D-01
         E(8) = 0.42890D+03
         S(8) = 0.58422D-01
         E(9) = 0.19454D+03
         S(9) = 0.12977D+00
         E(10) = 0.90777D+02
         S(10) = 0.24037D+00
         E(11) = 0.43715D+02
         S(11) = 0.32480D+00
         E(12) = 0.21491D+02
         S(12) = 0.26684D+00
         E(13) = 0.10669D+02
         S(13) = 0.91324D-01
         E(14) = 0.48882D+01
         S(14) = 0.53424D-02
         E(15) = 0.23541D+04
         S(15) = 0.42175D-05
         E(16) = 0.98140D+03
         S(16) = 0.20469D-04
         E(17) = 0.42890D+03
         S(17) = 0.14261D-03
         E(18) = 0.19454D+03
         S(18) = 0.63795D-03
         E(19) = 0.90777D+02
         S(19) = 0.28664D-02
         E(20) = 0.43715D+02
         S(20) = 0.75695D-02
         E(21) = 0.21491D+02
         S(21) = 0.15473D-01
         E(22) = 0.10669D+02
         S(22) = -.15913D-02
         E(23) = 0.48882D+01
         S(23) = -.46344D-01
         E(24) = 0.23361D+01
         S(24) = -.60218D-01
         E(25) = 0.23541D+04
         S(25) = -.57997D-04
         E(26) = 0.98140D+03
         S(26) = -.37736D-03
         E(27) = 0.42890D+03
         S(27) = -.21933D-02
         E(28) = 0.19454D+03
         S(28) = -.10675D-01
         E(29) = 0.90777D+02
         S(29) = -.45499D-01
         E(30) = 0.43715D+02
         S(30) = -.12587D+00
         E(31) = 0.21491D+02
         S(31) = -.25051D+00
         E(32) = 0.10669D+02
         S(32) = 0.10433D-01
         E(33) = 0.48882D+01
         S(33) = 0.86160D+00
         E(34) = 0.23361D+01
         S(34) = 0.14367D+01
         E(35) = 0.10967D+01
         S(35) = 0.76721D+00
         E(36) = 0.42312D+00
         S(36) = 0.10000D+01
         E(37) = 0.24276D+00
         S(37) = 0.10000D+01
         E(38) = 0.11703D+00
         S(38) = 0.10000D+01
         E(39) = 0.52861D-01
         S(39) = 0.10000D+01
         E(40) = 0.64533D+04
         P(40) = 0.21653D-04
         E(41) = 0.15339D+04
         P(41) = 0.19162D-03
         E(42) = 0.49838D+03
         P(42) = 0.11143D-02
         E(43) = 0.18898D+03
         P(43) = 0.50830D-02
         E(44) = 0.77887D+02
         P(44) = 0.19269D-01
         E(45) = 0.35514D+02
         P(45) = 0.53091D-01
         E(46) = 0.17016D+02
         P(46) = 0.12639D+00
         E(47) = 0.84977D+01
         P(47) = 0.22312D+00
         E(48) = 0.43544D+01
         P(48) = 0.29854D+00
         E(49) = 0.22575D+01
         P(49) = 0.28651D+00
         E(50) = 0.11597D+01
         P(50) = 0.17168D+00
         E(51) = 0.18898D+03
         P(51) = 0.19921D-04
         E(52) = 0.77887D+02
         P(52) = -.38970D-04
         E(53) = 0.35514D+02
         P(53) = 0.63327D-04
         E(54) = 0.17016D+02
         P(54) = -.14057D-02
         E(55) = 0.84977D+01
         P(55) = -.12080D-02
         E(56) = 0.43544D+01
         P(56) = -.11390D-01
         E(57) = 0.22575D+01
         P(57) = 0.46079D-02
         E(58) = 0.11597D+01
         P(58) = -.57688D-03
         E(59) = 0.59343D+00
         P(59) = 0.13174D+00
         E(60) = 0.29012D+00
         P(60) = 0.10000D+01
         E(61) = 0.13825D+00
         P(61) = 0.10000D+01
         E(62) = 0.64594D-01
         P(62) = 0.10000D+01
         E(63) = 0.28354D-01
         P(63) = 0.10000D+01
         E(64) = 0.19520D+02
         D(64) = 0.10000D+01
         E(65) = 0.55800D+01
         D(65) = 0.10000D+01
         E(66) = 0.19200D+01
         D(66) = 0.10000D+01
         E(67) = 0.69000D+00
         D(67) = 0.10000D+01
         E(68) = 0.29900D+00
         D(68) = 0.10000D+01
         E(69) = 0.12500D+00
         D(69) = 0.10000D+01
         E(70) = 0.26400D+01
         F(70) = 0.10000D+01
         E(71) = 0.62000D+00
         F(71) = 0.10000D+01
         E(72) = 0.26000D+00
         F(72) = 0.10000D+01
         E(73) = 0.13500D+01
         G(73) = 0.10000D+01
         E(74) = 0.45000D+00
         G(74) = 0.10000D+01
CFRJ     E(75) = 0.78000D+00
CFRJ     H(75) = 0.10000D+01
         ILAST = 74
      ENDIF
      RETURN
C
C     ----- P  -----
C
  180 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(1) = 0.89908D+04
         S(1) = 0.48206D-02
         E(2) = 0.13524D+04
         S(2) = 0.36392D-01
         E(3) = 0.30750D+03
         S(3) = 0.16563D+00
         E(4) = 0.85941D+02
         S(4) = 0.46217D+00
         E(5) = 0.26671D+02
         S(5) = 0.50582D+00
         E(6) = 0.13524D+04
         S(6) = -.28900D-03
         E(7) = 0.30750D+03
         S(7) = -.12638D-02
         E(8) = 0.85941D+02
         S(8) = -.24790D-01
         E(9) = 0.26671D+02
         S(9) = -.46165D-01
         E(10) = 0.51881D+01
         S(10) = 0.70180D-01
         E(11) = 0.30750D+03
         S(11) = 0.54512D-03
         E(12) = 0.85941D+02
         S(12) = 0.46838D-02
         E(13) = 0.26671D+02
         S(13) = 0.21160D-01
         E(14) = 0.51881D+01
         S(14) = -.12534D+00
         E(15) = 0.18778D+01
         S(15) = -.22427D+00
         E(16) = 0.22010D+00
         S(16) = 0.10000D+01
         E(17) = 0.11164D+03
         P(17) = 0.28663D-01
         E(18) = 0.25784D+02
         P(18) = 0.17265D+00
         E(19) = 0.76725D+01
         P(19) = 0.46181D+00
         E(20) = 0.23996D+01
         P(20) = 0.49174D+00
         E(21) = 0.25784D+02
         P(21) = -.68545D-03
         E(22) = 0.76725D+01
         P(22) = -.11296D-01
         E(23) = 0.23996D+01
         P(23) = 0.32962D-02
         E(24) = 0.54943D+00
         P(24) = 0.41613D+00
         E(25) = 0.14748D+00
         P(25) = 0.10000D+01
         ILAST = 25
      ELSE IF (NPOL.EQ.1) THEN
         E(1) = 0.48073D+05
         S(1) = 0.59919D-03
         E(2) = 0.72114D+04
         S(2) = 0.46241D-02
         E(3) = 0.16417D+04
         S(3) = 0.23724D-01
         E(4) = 0.46492D+03
         S(4) = 0.91303D-01
         E(5) = 0.15128D+03
         S(5) = 0.25750D+00
         E(6) = 0.54189D+02
         S(6) = 0.44080D+00
         E(7) = 0.20397D+02
         S(7) = 0.31175D+00
         E(8) = 0.50241D+01
         S(8) = 0.31118D-01
         E(9) = 0.72114D+04
         S(9) = -.21309D-04
         E(10) = 0.16417D+04
         S(10) = -.44215D-04
         E(11) = 0.46492D+03
         S(11) = -.18221D-02
         E(12) = 0.15128D+03
         S(12) = -.96366D-02
         E(13) = 0.54189D+02
         S(13) = -.66072D-01
         E(14) = 0.20397D+02
         S(14) = -.81824D-01
         E(15) = 0.50241D+01
         S(15) = 0.49850D+00
         E(16) = 0.18246D+01
         S(16) = 0.64492D+00
         E(17) = 0.16417D+04
         S(17) = 0.71340D-05
         E(18) = 0.46492D+03
         S(18) = 0.84500D-04
         E(19) = 0.15128D+03
         S(19) = 0.85919D-03
         E(20) = 0.54189D+02
         S(20) = 0.46478D-02
         E(21) = 0.20397D+02
         S(21) = 0.95811D-02
         E(22) = 0.50241D+01
         S(22) = -.84066D-01
         E(23) = 0.18246D+01
         S(23) = -.23124D+00
         E(24) = 0.35871D+00
         S(24) = 0.62350D+00
         E(25) = 0.12887D+00
         S(25) = 0.10000D+01
         E(26) = 0.36624D+03
         P(26) = 0.40977D-02
         E(27) = 0.85931D+02
         P(27) = 0.31956D-01
         E(28) = 0.27172D+02
         P(28) = 0.13326D+00
         E(29) = 0.98685D+01
         P(29) = 0.33228D+00
         E(30) = 0.37533D+01
         P(30) = 0.45335D+00
         E(31) = 0.14302D+01
         P(31) = 0.25058D+00
         E(32) = 0.85931D+02
         P(32) = 0.12506D-03
         E(33) = 0.27172D+02
         P(33) = -.12166D-02
         E(34) = 0.98685D+01
         P(34) = -.31813D-02
         E(35) = 0.37533D+01
         P(35) = -.21475D-01
         E(36) = 0.14302D+01
         P(36) = 0.69991D-01
         E(37) = 0.38653D+00
         P(37) = 0.56081D+00
         E(38) = 0.11146D+00
         P(38) = 0.10000D+01
         E(39) = 0.50000D+00
         D(39) = 0.10000D+01
         ILAST = 39
      ELSE IF (NPOL.EQ.2) THEN
         E(1) = 0.14568D+06
         S(1) = 0.15016D-03
         E(2) = 0.21836D+05
         S(2) = 0.11650D-02
         E(3) = 0.49699D+04
         S(3) = 0.60798D-02
         E(4) = 0.14073D+04
         S(4) = 0.24968D-01
         E(5) = 0.45878D+03
         S(5) = 0.83637D-01
         E(6) = 0.16508D+03
         S(6) = 0.21863D+00
         E(7) = 0.63744D+02
         S(7) = 0.39234D+00
         E(8) = 0.25598D+02
         S(8) = 0.35038D+00
         E(9) = 0.99307D+01
         S(9) = 0.77287D-01
         E(10) = 0.14073D+04
         S(10) = -.87601D-04
         E(11) = 0.45878D+03
         S(11) = -.52745D-03
         E(12) = 0.16508D+03
         S(12) = -.49575D-02
         E(13) = 0.63744D+02
         S(13) = -.20535D-01
         E(14) = 0.25598D+02
         S(14) = -.55030D-01
         E(15) = 0.99307D+01
         S(15) = 0.41583D-01
         E(16) = 0.39686D+01
         S(16) = 0.17029D+00
         E(17) = 0.14073D+04
         S(17) = 0.30955D-04
         E(18) = 0.45878D+03
         S(18) = 0.35776D-03
         E(19) = 0.16508D+03
         S(19) = 0.26123D-02
         E(20) = 0.63744D+02
         S(20) = 0.13204D-01
         E(21) = 0.25598D+02
         S(21) = 0.32848D-01
         E(22) = 0.99307D+01
         S(22) = -.18182D-01
         E(23) = 0.39686D+01
         S(23) = -.25722D+00
         E(24) = 0.15565D+01
         S(24) = -.30682D+00
         E(25) = 0.35404D+00
         S(25) = 0.10000D+01
         E(26) = 0.12394D+00
         S(26) = 0.10000D+01
         E(27) = 0.82623D+03
         P(27) = 0.10225D-02
         E(28) = 0.19602D+03
         P(28) = 0.84472D-02
         E(29) = 0.62951D+02
         P(29) = 0.42664D-01
         E(30) = 0.23222D+02
         P(30) = 0.14239D+00
         E(31) = 0.93885D+01
         P(31) = 0.31059D+00
         E(32) = 0.38997D+01
         P(32) = 0.41672D+00
         E(33) = 0.16048D+01
         P(33) = 0.26618D+00
         E(34) = 0.19602D+03
         P(34) = 0.41082D-05
         E(35) = 0.62951D+02
         P(35) = -.15057D-03
         E(36) = 0.23222D+02
         P(36) = -.11141D-02
         E(37) = 0.93885D+01
         P(37) = -.64942D-02
         E(38) = 0.38997D+01
         P(38) = -.12469D-01
         E(39) = 0.16048D+01
         P(39) = 0.23458D-01
         E(40) = 0.56053D+00
         P(40) = 0.32270D+00
         E(41) = 0.21175D+00
         P(41) = 0.10000D+01
         E(42) = 0.73486D-01
         P(42) = 0.10000D+01
         E(43) = 0.18900D+01
         D(43) = 0.10000D+01
         E(44) = 0.43000D+00
         D(44) = 0.10000D+01
         E(45) = 0.60000D+00
         F(45) = 0.10000D+01
         ILAST = 45
      ELSE IF (NPOL.EQ.3) THEN
         E(1) = 0.10637D+07
         S(1) = 0.12509D-04
         E(2) = 0.16000D+06
         S(2) = 0.96671D-04
         E(3) = 0.36513D+05
         S(3) = 0.50682D-03
         E(4) = 0.10359D+05
         S(4) = 0.21310D-02
         E(5) = 0.33828D+04
         S(5) = 0.76925D-02
         E(6) = 0.12223D+04
         S(6) = 0.24421D-01
         E(7) = 0.47730D+03
         S(7) = 0.68122D-01
         E(8) = 0.19809D+03
         S(8) = 0.16069D+00
         E(9) = 0.85994D+02
         S(9) = 0.29867D+00
         E(10) = 0.38627D+02
         S(10) = 0.35878D+00
         E(11) = 0.17688D+02
         S(11) = 0.20616D+00
         E(12) = 0.72551D+01
         S(12) = 0.27386D-01
         E(13) = 0.10359D+05
         S(13) = -.34502D-05
         E(14) = 0.33828D+04
         S(14) = -.53427D-05
         E(15) = 0.12223D+04
         S(15) = -.13814D-03
         E(16) = 0.47730D+03
         S(16) = -.60882D-03
         E(17) = 0.19809D+03
         S(17) = -.45064D-02
         E(18) = 0.85994D+02
         S(18) = -.17071D-01
         E(19) = 0.38627D+02
         S(19) = -.53954D-01
         E(20) = 0.17688D+02
         S(20) = -.51494D-01
         E(21) = 0.72551D+01
         S(21) = 0.13758D+00
         E(22) = 0.32699D+01
         S(22) = 0.27525D+00
         E(23) = 0.12223D+04
         S(23) = 0.15801D-04
         E(24) = 0.47730D+03
         S(24) = 0.11353D-03
         E(25) = 0.19809D+03
         S(25) = 0.67329D-03
         E(26) = 0.85994D+02
         S(26) = 0.30314D-02
         E(27) = 0.38627D+02
         S(27) = 0.92292D-02
         E(28) = 0.17688D+02
         S(28) = 0.11936D-01
         E(29) = 0.72551D+01
         S(29) = -.35725D-01
         E(30) = 0.32699D+01
         S(30) = -.15875D+00
         E(31) = 0.14313D+01
         S(31) = -.21236D+00
         E(32) = 0.47047D+00
         S(32) = 0.10000D+01
         E(33) = 0.21289D+00
         S(33) = 0.10000D+01
         E(34) = 0.87742D-01
         S(34) = 0.10000D+01
         E(35) = 0.27622D+04
         P(35) = 0.12618D-03
         E(36) = 0.65478D+03
         P(36) = 0.11048D-02
         E(37) = 0.21258D+03
         P(37) = 0.61483D-02
         E(38) = 0.79895D+02
         P(38) = 0.26307D-01
         E(39) = 0.33366D+02
         P(39) = 0.79866D-01
         E(40) = 0.14986D+02
         P(40) = 0.18840D+00
         E(41) = 0.69680D+01
         P(41) = 0.31362D+00
         E(42) = 0.32818D+01
         P(42) = 0.35326D+00
         E(43) = 0.15247D+01
         P(43) = 0.21366D+00
         E(44) = 0.65478D+03
         P(44) = 0.59592D-05
         E(45) = 0.21258D+03
         P(45) = -.68086D-05
         E(46) = 0.79895D+02
         P(46) = 0.74513D-04
         E(47) = 0.33366D+02
         P(47) = -.60664D-03
         E(48) = 0.14986D+02
         P(48) = -.12604D-02
         E(49) = 0.69680D+01
         P(49) = -.94719D-02
         E(50) = 0.32818D+01
         P(50) = -.44014D-02
         E(51) = 0.15247D+01
         P(51) = 0.12466D-01
         E(52) = 0.66924D+00
         P(52) = 0.21343D+00
         E(53) = 0.29929D+00
         P(53) = 0.10000D+01
         E(54) = 0.12811D+00
         P(54) = 0.10000D+01
         E(55) = 0.50992D-01
         P(55) = 0.10000D+01
         E(56) = 0.77100D+01
         D(56) = 0.10000D+01
         E(57) = 0.20300D+01
         D(57) = 0.10000D+01
         E(58) = 0.60000D+00
         D(58) = 0.10000D+01
         E(59) = 0.22000D+00
         D(59) = 0.10000D+01
         E(60) = 0.13400D+01
         F(60) = 0.10000D+01
         E(61) = 0.45000D+00
         F(61) = 0.10000D+01
         E(62) = 0.72000D+00
         G(62) = 0.10000D+01
         ILAST = 62
      ELSE IF (NPOL.EQ.4) THEN
         E(1) = 0.82213D+06
         S(1) = 0.12207D-04
         E(2) = 0.19425D+06
         S(2) = 0.61440D-04
         E(3) = 0.56424D+05
         S(3) = 0.25322D-03
         E(4) = 0.18716D+05
         S(4) = 0.90752D-03
         E(5) = 0.68393D+04
         S(5) = 0.29230D-02
         E(6) = 0.26931D+04
         S(6) = 0.86138D-02
         E(7) = 0.11265D+04
         S(7) = 0.23329D-01
         E(8) = 0.49521D+03
         S(8) = 0.57641D-01
         E(9) = 0.22638D+03
         S(9) = 0.12694D+00
         E(10) = 0.10662D+03
         S(10) = 0.23408D+00
         E(11) = 0.51905D+02
         S(11) = 0.31805D+00
         E(12) = 0.25851D+02
         S(12) = 0.27002D+00
         E(13) = 0.13022D+02
         S(13) = 0.10197D+00
         E(14) = 0.60568D+01
         S(14) = 0.80899D-02
         E(15) = 0.26931D+04
         S(15) = 0.36778D-05
         E(16) = 0.11265D+04
         S(16) = 0.16467D-04
         E(17) = 0.49521D+03
         S(17) = 0.11861D-03
         E(18) = 0.22638D+03
         S(18) = 0.50936D-03
         E(19) = 0.10662D+03
         S(19) = 0.22939D-02
         E(20) = 0.51905D+02
         S(20) = 0.59503D-02
         E(21) = 0.25851D+02
         S(21) = 0.12532D-01
         E(22) = 0.13022D+02
         S(22) = 0.34088D-03
         E(23) = 0.60568D+01
         S(23) = -.33461D-01
         E(24) = 0.29574D+01
         S(24) = -.50380D-01
         E(25) = 0.26931D+04
         S(25) = -.75718D-04
         E(26) = 0.11265D+04
         S(26) = -.48033D-03
         E(27) = 0.49521D+03
         S(27) = -.27731D-02
         E(28) = 0.22638D+03
         S(28) = -.13224D-01
         E(29) = 0.10662D+03
         S(29) = -.55695D-01
         E(30) = 0.51905D+02
         S(30) = -.15313D+00
         E(31) = 0.25851D+02
         S(31) = -.31055D+00
         E(32) = 0.13022D+02
         S(32) = -.34248D-01
         E(33) = 0.60568D+01
         S(33) = 0.97684D+00
         E(34) = 0.29574D+01
         S(34) = 0.18179D+01
         E(35) = 0.14078D+01
         S(35) = 0.10730D+01
         E(36) = 0.57187D+00
         S(36) = 0.10000D+01
         E(37) = 0.31078D+00
         S(37) = 0.10000D+01
         E(38) = 0.14947D+00
         S(38) = 0.10000D+01
         E(39) = 0.68313D-01
         S(39) = 0.10000D+01
         E(40) = 0.77486D+04
         P(40) = 0.20870D-04
         E(41) = 0.18336D+04
         P(41) = 0.18596D-03
         E(42) = 0.59555D+03
         P(42) = 0.10799D-02
         E(43) = 0.22661D+03
         P(43) = 0.48946D-02
         E(44) = 0.94191D+02
         P(44) = 0.18320D-01
         E(45) = 0.43332D+02
         P(45) = 0.50669D-01
         E(46) = 0.20907D+02
         P(46) = 0.12156D+00
         E(47) = 0.10502D+02
         P(47) = 0.21931D+00
         E(48) = 0.54255D+01
         P(48) = 0.29640D+00
         E(49) = 0.28818D+01
         P(49) = 0.28179D+00
         E(50) = 0.15284D+01
         P(50) = 0.17853D+00
         E(51) = 0.22661D+03
         P(51) = 0.15973D-04
         E(52) = 0.94191D+02
         P(52) = -.47396D-04
         E(53) = 0.20907D+02
         P(53) = -.15832D-02
         E(54) = 0.10502D+02
         P(54) = -.18117D-02
         E(55) = 0.54255D+01
         P(55) = -.13465D-01
         E(56) = 0.28818D+01
         P(56) = 0.61503D-02
         E(57) = 0.15284D+01
         P(57) = -.71444D-02
         E(58) = 0.81081D+00
         P(58) = 0.14315D+00
         E(59) = 0.39028D+00
         P(59) = 0.10000D+01
         E(60) = 0.18645D+00
         P(60) = 0.10000D+01
         E(61) = 0.87567D-01
         P(61) = 0.10000D+01
         E(62) = 0.38673D-01
         P(62) = 0.10000D+01
         E(63) = 0.21230D+02
         D(63) = 0.10000D+01
         E(64) = 0.60700D+01
         D(64) = 0.10000D+01
         E(65) = 0.20900D+01
         D(65) = 0.10000D+01
         E(66) = 0.75000D+00
         D(66) = 0.10000D+01
         E(67) = 0.32500D+00
         D(67) = 0.10000D+01
         E(68) = 0.13500D+00
         D(68) = 0.10000D+01
         E(69) = 0.31200D+01
         F(69) = 0.10000D+01
         E(70) = 0.73000D+00
         F(70) = 0.10000D+01
         E(71) = 0.30000D+00
         F(71) = 0.10000D+01
         E(72) = 0.15500D+01
         G(72) = 0.10000D+01
         E(73) = 0.52000D+00
         G(73) = 0.10000D+01
CFRJ     E(74) = 0.87000D+00
CFRJ     H(74) = 0.10000D+01
         ILAST = 73
      ENDIF
      RETURN
C
C     ----- S  -----
C
  200 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(1) = 0.10330D+05
         S(1) = 0.47685D-02
         E(2) = 0.15537D+04
         S(2) = 0.36209D-01
         E(3) = 0.35330D+03
         S(3) = 0.16491D+00
         E(4) = 0.98805D+02
         S(4) = 0.47548D+00
         E(5) = 0.30716D+02
         S(5) = 0.53319D+00
         E(6) = 0.15537D+04
         S(6) = -.28903D-03
         E(7) = 0.35330D+03
         S(7) = -.11575D-02
         E(8) = 0.98805D+02
         S(8) = -.23993D-01
         E(9) = 0.30716D+02
         S(9) = -.42830D-01
         E(10) = 0.60633D+01
         S(10) = 0.42063D-01
         E(11) = 0.35330D+03
         S(11) = 0.62501D-03
         E(12) = 0.98805D+02
         S(12) = 0.53881D-02
         E(13) = 0.30716D+02
         S(13) = 0.24521D-01
         E(14) = 0.60633D+01
         S(14) = -.14132D+00
         E(15) = 0.22475D+01
         S(15) = -.23948D+00
         E(16) = 0.27853D+00
         S(16) = 0.10000D+01
         E(17) = 0.12976D+03
         P(17) = 0.28706D-01
         E(18) = 0.30025D+02
         P(18) = 0.17420D+00
         E(19) = 0.89733D+01
         P(19) = 0.46709D+00
         E(20) = 0.28440D+01
         P(20) = 0.48571D+00
         E(21) = 0.30025D+02
         P(21) = -.11785D-03
         E(22) = 0.89733D+01
         P(22) = -.14483D-01
         E(23) = 0.28440D+01
         P(23) = 0.15939D-01
         E(24) = 0.63361D+00
         P(24) = 0.48257D+00
         E(25) = 0.16955D+00
         P(25) = 0.10000D+01
         ILAST = 25
      ELSE IF (NPOL.EQ.1) THEN
         E(1) = 0.55510D+05
         S(1) = 0.58899D-03
         E(2) = 0.83267D+04
         S(2) = 0.45460D-02
         E(3) = 0.18956D+04
         S(3) = 0.23335D-01
         E(4) = 0.53677D+03
         S(4) = 0.89965D-01
         E(5) = 0.17467D+03
         S(5) = 0.25486D+00
         E(6) = 0.62578D+02
         S(6) = 0.43998D+00
         E(7) = 0.23573D+02
         S(7) = 0.31549D+00
         E(8) = 0.58908D+01
         S(8) = 0.32826D-01
         E(9) = 0.83267D+04
         S(9) = -.20480D-04
         E(10) = 0.18956D+04
         S(10) = -.45429D-04
         E(11) = 0.53677D+03
         S(11) = -.17937D-02
         E(12) = 0.17467D+03
         S(12) = -.97252D-02
         E(13) = 0.62578D+02
         S(13) = -.66399D-01
         E(14) = 0.23573D+02
         S(14) = -.85502D-01
         E(15) = 0.58908D+01
         S(15) = 0.49588D+00
         E(16) = 0.21852D+01
         S(16) = 0.64988D+00
         E(17) = 0.18956D+04
         S(17) = 0.83461D-05
         E(18) = 0.53677D+03
         S(18) = 0.98306D-04
         E(19) = 0.17467D+03
         S(19) = 0.10040D-02
         E(20) = 0.62578D+02
         S(20) = 0.54409D-02
         E(21) = 0.23573D+02
         S(21) = 0.11400D-01
         E(22) = 0.58908D+01
         S(22) = -.95134D-01
         E(23) = 0.21852D+01
         S(23) = -.24735D+00
         E(24) = 0.45282D+00
         S(24) = 0.63973D+00
         E(25) = 0.16202D+00
         S(25) = 0.10000D+01
         E(26) = 0.39545D+03
         P(26) = 0.46266D-02
         E(27) = 0.92834D+02
         P(27) = 0.35898D-01
         E(28) = 0.29266D+02
         P(28) = 0.14854D+00
         E(29) = 0.10595D+02
         P(29) = 0.35878D+00
         E(30) = 0.40272D+01
         P(30) = 0.45417D+00
         E(31) = 0.15484D+01
         P(31) = 0.20545D+00
         E(32) = 0.92834D+02
         P(32) = 0.10934D-03
         E(33) = 0.29266D+02
         P(33) = -.16358D-02
         E(34) = 0.10595D+02
         P(34) = -.50632D-02
         E(35) = 0.40272D+01
         P(35) = -.23960D-01
         E(36) = 0.15484D+01
         P(36) = 0.10388D+00
         E(37) = 0.46134D+00
         P(37) = 0.58063D+00
         E(38) = 0.13131D+00
         P(38) = 0.10000D+01
         E(39) = 0.55000D+00
         D(39) = 0.10000D+01
         ILAST = 39
      ELSE IF (NPOL.EQ.2) THEN
         E(1) = 0.17203D+06
         S(1) = 0.14353D-03
         E(2) = 0.25783D+05
         S(2) = 0.11137D-02
         E(3) = 0.58682D+04
         S(3) = 0.58143D-02
         E(4) = 0.16617D+04
         S(4) = 0.23903D-01
         E(5) = 0.54175D+03
         S(5) = 0.80354D-01
         E(6) = 0.19501D+03
         S(6) = 0.21175D+00
         E(7) = 0.75396D+02
         S(7) = 0.38664D+00
         E(8) = 0.30407D+02
         S(8) = 0.35747D+00
         E(9) = 0.12121D+02
         S(9) = 0.86131D-01
         E(10) = 0.16617D+04
         S(10) = -.90480D-04
         E(11) = 0.54175D+03
         S(11) = -.50105D-03
         E(12) = 0.19501D+03
         S(12) = -.49510D-02
         E(13) = 0.75396D+02
         S(13) = -.20252D-01
         E(14) = 0.30407D+02
         S(14) = -.57388D-01
         E(15) = 0.12121D+02
         S(15) = 0.34847D-01
         E(16) = 0.47777D+01
         S(16) = 0.16833D+00
         E(17) = 0.16617D+04
         S(17) = 0.28367D-04
         E(18) = 0.54175D+03
         S(18) = 0.32591D-03
         E(19) = 0.19501D+03
         S(19) = 0.24131D-02
         E(20) = 0.75396D+02
         S(20) = 0.12352D-01
         E(21) = 0.30407D+02
         S(21) = 0.32119D-01
         E(22) = 0.12121D+02
         S(22) = -.11611D-01
         E(23) = 0.47777D+01
         S(23) = -.25542D+00
         E(24) = 0.18906D+01
         S(24) = -.31384D+00
         E(25) = 0.44686D+00
         S(25) = 0.10000D+01
         E(26) = 0.15621D+00
         S(26) = 0.10000D+01
         E(27) = 0.96812D+03
         P(27) = 0.10048D-02
         E(28) = 0.22970D+03
         P(28) = 0.83313D-02
         E(29) = 0.73833D+02
         P(29) = 0.42260D-01
         E(30) = 0.27312D+02
         P(30) = 0.14206D+00
         E(31) = 0.11074D+02
         P(31) = 0.31256D+00
         E(32) = 0.46307D+01
         P(32) = 0.41933D+00
         E(33) = 0.19299D+01
         P(33) = 0.25924D+00
         E(34) = 0.22970D+03
         P(34) = -.80693D-05
         E(35) = 0.73833D+02
         P(35) = -.16943D-03
         E(36) = 0.27312D+02
         P(36) = -.15243D-02
         E(37) = 0.11074D+02
         P(37) = -.71993D-02
         E(38) = 0.46307D+01
         P(38) = -.16113D-01
         E(39) = 0.19299D+01
         P(39) = 0.34234D-01
         E(40) = 0.68783D+00
         P(40) = 0.35518D+00
         E(41) = 0.25663D+00
         P(41) = 0.10000D+01
         E(42) = 0.87251D-01
         P(42) = 0.10000D+01
         E(43) = 0.21000D+01
         D(43) = 0.10000D+01
         E(44) = 0.48000D+00
         D(44) = 0.10000D+01
         E(45) = 0.66000D+00
         F(45) = 0.10000D+01
         ILAST = 45
      ELSE IF (NPOL.EQ.3) THEN
         E(1) = 0.12504D+07
         S(1) = 0.12023D-04
         E(2) = 0.18808D+06
         S(2) = 0.92935D-04
         E(3) = 0.42918D+05
         S(3) = 0.48720D-03
         E(4) = 0.12179D+05
         S(4) = 0.20479D-02
         E(5) = 0.39787D+04
         S(5) = 0.73897D-02
         E(6) = 0.14381D+04
         S(6) = 0.23481D-01
         E(7) = 0.56167D+03
         S(7) = 0.65685D-01
         E(8) = 0.23311D+03
         S(8) = 0.15597D+00
         E(9) = 0.10122D+03
         S(9) = 0.29322D+00
         E(10) = 0.45493D+02
         S(10) = 0.36011D+00
         E(11) = 0.20847D+02
         S(11) = 0.21506D+00
         E(12) = 0.86492D+01
         S(12) = 0.30782D-01
         E(13) = 0.12179D+05
         S(13) = -.30383D-05
         E(14) = 0.39787D+04
         S(14) = -.46456D-05
         E(15) = 0.14381D+04
         S(15) = -.12316D-03
         E(16) = 0.56167D+03
         S(16) = -.54723D-03
         E(17) = 0.23311D+03
         S(17) = -.40702D-02
         E(18) = 0.10122D+03
         S(18) = -.15652D-01
         E(19) = 0.45493D+02
         S(19) = -.50309D-01
         E(20) = 0.20847D+02
         S(20) = -.51517D-01
         E(21) = 0.86492D+01
         S(21) = 0.11823D+00
         E(22) = 0.39456D+01
         S(22) = 0.25628D+00
         E(23) = 0.14381D+04
         S(23) = 0.19556D-04
         E(24) = 0.56167D+03
         S(24) = 0.14386D-03
         E(25) = 0.23311D+03
         S(25) = 0.84419D-03
         E(26) = 0.10122D+03
         S(26) = 0.38621D-02
         E(27) = 0.45493D+02
         S(27) = 0.11813D-01
         E(28) = 0.20847D+02
         S(28) = 0.16114D-01
         E(29) = 0.86492D+01
         S(29) = -.41492D-01
         E(30) = 0.39456D+01
         S(30) = -.18517D+00
         E(31) = 0.17472D+01
         S(31) = -.23963D+00
         E(32) = 0.59953D+00
         S(32) = 0.10000D+01
         E(33) = 0.27127D+00
         S(33) = 0.10000D+01
         E(34) = 0.11080D+00
         S(34) = 0.10000D+01
         E(35) = 0.32858D+04
         P(35) = 0.12062D-03
         E(36) = 0.77893D+03
         P(36) = 0.10587D-02
         E(37) = 0.25276D+03
         P(37) = 0.59191D-02
         E(38) = 0.95100D+02
         P(38) = 0.25405D-01
         E(39) = 0.39820D+02
         P(39) = 0.77859D-01
         E(40) = 0.17925D+02
         P(40) = 0.18578D+00
         E(41) = 0.83699D+01
         P(41) = 0.31347D+00
         E(42) = 0.39655D+01
         P(42) = 0.35718D+00
         E(43) = 0.18540D+01
         P(43) = 0.21179D+00
         E(44) = 0.77893D+03
         P(44) = 0.42397D-05
         E(45) = 0.25276D+03
         P(45) = -.74333D-05
         E(46) = 0.95100D+02
         P(46) = 0.27726D-04
         E(47) = 0.39820D+02
         P(47) = -.65028D-03
         E(48) = 0.17925D+02
         P(48) = -.18933D-02
         E(49) = 0.83699D+01
         P(49) = -.10191D-01
         E(50) = 0.39655D+01
         P(50) = -.76291D-02
         E(51) = 0.18540D+01
         P(51) = 0.21032D-01
         E(52) = 0.81971D+00
         P(52) = 0.24347D+00
         E(53) = 0.36117D+00
         P(53) = 0.10000D+01
         E(54) = 0.15097D+00
         P(54) = 0.10000D+01
         E(55) = 0.59620D-01
         P(55) = 0.10000D+01
         E(56) = 0.87600D+01
         D(56) = 0.10000D+01
         E(57) = 0.23000D+01
         D(57) = 0.10000D+01
         E(58) = 0.68000D+00
         D(58) = 0.10000D+01
         E(59) = 0.25000D+00
         D(59) = 0.10000D+01
         E(60) = 0.15300D+01
         F(60) = 0.10000D+01
         E(61) = 0.51000D+00
         F(61) = 0.10000D+01
         E(62) = 0.79000D+00
         G(62) = 0.10000D+01
         ILAST = 62
      ELSE IF (NPOL.EQ.4) THEN
         E(1) = 0.98389D+06
         S(1) = 0.11563D-04
         E(2) = 0.22977D+06
         S(2) = 0.59017D-04
         E(3) = 0.66276D+05
         S(3) = 0.24474D-03
         E(4) = 0.21884D+05
         S(4) = 0.88159D-03
         E(5) = 0.79696D+04
         S(5) = 0.28494D-02
         E(6) = 0.31299D+04
         S(6) = 0.84273D-02
         E(7) = 0.13052D+04
         S(7) = 0.22951D-01
         E(8) = 0.57170D+03
         S(8) = 0.57038D-01
         E(9) = 0.26058D+03
         S(9) = 0.12616D+00
         E(10) = 0.12248D+03
         S(10) = 0.23373D+00
         E(11) = 0.59518D+02
         S(11) = 0.31883D+00
         E(12) = 0.29627D+02
         S(12) = 0.27098D+00
         E(13) = 0.14943D+02
         S(13) = 0.10217D+00
         E(14) = 0.70479D+01
         S(14) = 0.81326D-02
         E(15) = 0.31299D+04
         S(15) = 0.34121D-05
         E(16) = 0.13052D+04
         S(16) = 0.15187D-04
         E(17) = 0.57170D+03
         S(17) = 0.11093D-03
         E(18) = 0.26058D+03
         S(18) = 0.47967D-03
         E(19) = 0.12248D+03
         S(19) = 0.21682D-02
         E(20) = 0.59518D+02
         S(20) = 0.56908D-02
         E(21) = 0.29627D+02
         S(21) = 0.11915D-01
         E(22) = 0.14943D+02
         S(22) = 0.27168D-03
         E(23) = 0.70479D+01
         S(23) = -.30365D-01
         E(24) = 0.35040D+01
         S(24) = -.47588D-01
         E(25) = 0.31299D+04
         S(25) = -.84316D-04
         E(26) = 0.13052D+04
         S(26) = -.54837D-03
         E(27) = 0.57170D+03
         S(27) = -.31574D-02
         E(28) = 0.26058D+03
         S(28) = -.15292D-01
         E(29) = 0.12248D+03
         S(29) = -.64288D-01
         E(30) = 0.59518D+02
         S(30) = -.17969D+00
         E(31) = 0.29627D+02
         S(31) = -.36055D+00
         E(32) = 0.14943D+02
         S(32) = -.42524D-01
         E(33) = 0.70479D+01
         S(33) = 0.11075D+01
         E(34) = 0.35040D+01
         S(34) = 0.20784D+01
         E(35) = 0.16997D+01
         S(35) = 0.12396D+01
         E(36) = 0.73101D+00
         S(36) = 0.10000D+01
         E(37) = 0.37948D+00
         S(37) = 0.10000D+01
         E(38) = 0.18185D+00
         S(38) = 0.10000D+01
         E(39) = 0.83961D-01
         S(39) = 0.10000D+01
         E(40) = 0.93414D+04
         P(40) = 0.19493D-04
         E(41) = 0.22086D+04
         P(41) = 0.17432D-03
         E(42) = 0.71627D+03
         P(42) = 0.10161D-02
         E(43) = 0.27251D+03
         P(43) = 0.46172D-02
         E(44) = 0.11377D+03
         P(44) = 0.17163D-01
         E(45) = 0.52571D+02
         P(45) = 0.48132D-01
         E(46) = 0.25397D+02
         P(46) = 0.11704D+00
         E(47) = 0.12770D+02
         P(47) = 0.21602D+00
         E(48) = 0.66224D+01
         P(48) = 0.29528D+00
         E(49) = 0.35388D+01
         P(49) = 0.28825D+00
         E(50) = 0.18756D+01
         P(50) = 0.17982D+00
         E(51) = 0.27251D+03
         P(51) = 0.10626D-04
         E(52) = 0.11377D+03
         P(52) = -.50176D-04
         E(53) = 0.52571D+02
         P(53) = -.78047D-04
         E(54) = 0.25397D+02
         P(54) = -.16188D-02
         E(55) = 0.12770D+02
         P(55) = -.26551D-02
         E(56) = 0.66224D+01
         P(56) = -.13777D-01
         E(57) = 0.35388D+01
         P(57) = 0.24246D-02
         E(58) = 0.18756D+01
         P(58) = 0.58572D-03
         E(59) = 0.99018D+00
         P(59) = 0.16518D+00
         E(60) = 0.47723D+00
         P(60) = 0.10000D+01
         E(61) = 0.22692D+00
         P(61) = 0.10000D+01
         E(62) = 0.10473D+00
         P(62) = 0.10000D+01
         E(63) = 0.45303D-01
         P(63) = 0.10000D+01
         E(64) = 0.22930D+02
         D(64) = 0.10000D+01
         E(65) = 0.65500D+01
         D(65) = 0.10000D+01
         E(66) = 0.22600D+01
         D(66) = 0.10000D+01
         E(67) = 0.81000D+00
         D(67) = 0.10000D+01
         E(68) = 0.35100D+00
         D(68) = 0.10000D+01
         E(69) = 0.14600D+00
         D(69) = 0.10000D+01
         E(70) = 0.35900D+01
         F(70) = 0.10000D+01
         E(71) = 0.84000D+00
         F(71) = 0.10000D+01
         E(72) = 0.35000D+00
         F(72) = 0.10000D+01
         E(73) = 0.17600D+01
         G(73) = 0.10000D+01
         E(74) = 0.59000D+00
         G(74) = 0.10000D+01
CFRJ     E(75) = 0.95000D+00
CFRJ     H(75) = 0.10000D+01
         ILAST = 74
      ENDIF
      RETURN
C
C     ----- CL -----
C
  220 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(1) = 0.11759D+05
         S(1) = 0.47248D-02
         E(2) = 0.17686D+04
         S(2) = 0.36546D-01
         E(3) = 0.40219D+03
         S(3) = 0.16597D+00
         E(4) = 0.11254D+03
         S(4) = 0.52530D+00
         E(5) = 0.35039D+02
         S(5) = 0.61955D+00
         E(6) = 0.17686D+04
         S(6) = -.29053D-03
         E(7) = 0.40219D+03
         S(7) = -.10664D-02
         E(8) = 0.11254D+03
         S(8) = -.23406D-01
         E(9) = 0.35039D+02
         S(9) = -.39915D-01
         E(10) = 0.69956D+01
         S(10) = 0.17937D-01
         E(11) = 0.40219D+03
         S(11) = 0.69197D-03
         E(12) = 0.11254D+03
         S(12) = 0.59760D-02
         E(13) = 0.35039D+02
         S(13) = 0.27394D-01
         E(14) = 0.69956D+01
         S(14) = -.15516D+00
         E(15) = 0.26455D+01
         S(15) = -.25028D+00
         E(16) = 0.34157D+00
         S(16) = 0.10000D+01
         E(17) = 0.15066D+03
         P(17) = 0.28299D-01
         E(18) = 0.34933D+02
         P(18) = 0.17336D+00
         E(19) = 0.10492D+02
         P(19) = 0.46812D+00
         E(20) = 0.33726D+01
         P(20) = 0.48321D+00
         E(21) = 0.34933D+02
         P(21) = 0.13598D-03
         E(22) = 0.10492D+02
         P(22) = -.16611D-01
         E(23) = 0.33726D+01
         P(23) = 0.22351D-01
         E(24) = 0.75989D+00
         P(24) = 0.51192D+00
         E(25) = 0.20289D+00
         P(25) = 0.10000D+01
         ILAST = 25
      ELSE IF (NPOL.EQ.1) THEN
         E(1) = 0.63474D+05
         S(1) = 0.58028D-03
         E(2) = 0.95213D+04
         S(2) = 0.44792D-02
         E(3) = 0.21674D+04
         S(3) = 0.23000D-01
         E(4) = 0.61374D+03
         S(4) = 0.88815D-01
         E(5) = 0.19972D+03
         S(5) = 0.25258D+00
         E(6) = 0.71564D+02
         S(6) = 0.43921D+00
         E(7) = 0.26979D+02
         S(7) = 0.31868D+00
         E(8) = 0.68260D+01
         S(8) = 0.34348D-01
         E(9) = 0.95213D+04
         S(9) = -.19824D-04
         E(10) = 0.21674D+04
         S(10) = -.47199D-04
         E(11) = 0.61374D+03
         S(11) = -.17713D-02
         E(12) = 0.19972D+03
         S(12) = -.98068D-02
         E(13) = 0.71564D+02
         S(13) = -.66719D-01
         E(14) = 0.26979D+02
         S(14) = -.88711D-01
         E(15) = 0.68260D+01
         S(15) = 0.49396D+00
         E(16) = 0.25767D+01
         S(16) = 0.65392D+00
         E(17) = 0.21674D+04
         S(17) = 0.90584D-05
         E(18) = 0.61374D+03
         S(18) = 0.10695D-03
         E(19) = 0.19972D+03
         S(19) = 0.10981D-02
         E(20) = 0.71564D+02
         S(20) = 0.59898D-02
         E(21) = 0.26979D+02
         S(21) = 0.12762D-01
         E(22) = 0.68260D+01
         S(22) = -.10307D+00
         E(23) = 0.25767D+01
         S(23) = -.25941D+00
         E(24) = 0.55752D+00
         S(24) = 0.64899D+00
         E(25) = 0.19825D+00
         S(25) = 0.10000D+01
         E(26) = 0.41312D+03
         P(26) = 0.54251D-02
         E(27) = 0.97054D+02
         P(27) = 0.41644D-01
         E(28) = 0.30482D+02
         P(28) = 0.16990D+00
         E(29) = 0.10962D+02
         P(29) = 0.39275D+00
         E(30) = 0.41276D+01
         P(30) = 0.45045D+00
         E(31) = 0.15338D+01
         P(31) = 0.15308D+00
         E(32) = 0.30482D+02
         P(32) = -.21441D-02
         E(33) = 0.10962D+02
         P(33) = -.91298D-02
         E(34) = 0.41276D+01
         P(34) = -.19906D-01
         E(35) = 0.15338D+01
         P(35) = 0.15582D+00
         E(36) = 0.51439D+00
         P(36) = 0.58772D+00
         E(37) = 0.15026D+00
         P(37) = 0.10000D+01
         E(38) = 0.60000D+00
         D(38) = 0.10000D+01
         ILAST = 38
      ELSE IF (NPOL.EQ.2) THEN
         E(1) = 0.20115D+06
         S(1) = 0.13753D-03
         E(2) = 0.30146D+05
         S(2) = 0.10671D-02
         E(3) = 0.68610D+04
         S(3) = 0.55737D-02
         E(4) = 0.19428D+04
         S(4) = 0.22935D-01
         E(5) = 0.63342D+03
         S(5) = 0.77354D-01
         E(6) = 0.22810D+03
         S(6) = 0.20534D+00
         E(7) = 0.88288D+02
         S(7) = 0.38093D+00
         E(8) = 0.35750D+02
         S(8) = 0.36361D+00
         E(9) = 0.14571D+02
         S(9) = 0.95182D-01
         E(10) = 0.19428D+04
         S(10) = -.94653D-04
         E(11) = 0.63342D+03
         S(11) = -.48337D-03
         E(12) = 0.22810D+03
         S(12) = -.50009D-02
         E(13) = 0.88288D+02
         S(13) = -.20230D-01
         E(14) = 0.35750D+02
         S(14) = -.60209D-01
         E(15) = 0.14571D+02
         S(15) = 0.29195D-01
         E(16) = 0.56580D+01
         S(16) = 0.16993D+00
         E(17) = 0.19428D+04
         S(17) = 0.24868D-04
         E(18) = 0.63342D+03
         S(18) = 0.28967D-03
         E(19) = 0.22810D+03
         S(19) = 0.21513D-02
         E(20) = 0.88288D+02
         S(20) = 0.11195D-01
         E(21) = 0.35750D+02
         S(21) = 0.30180D-01
         E(22) = 0.14571D+02
         S(22) = -.59456D-02
         E(23) = 0.56580D+01
         S(23) = -.24702D+00
         E(24) = 0.22554D+01
         S(24) = -.31412D+00
         E(25) = 0.54947D+00
         S(25) = 0.10000D+01
         E(26) = 0.19126D+00
         S(26) = 0.10000D+01
         E(27) = 0.11299D+04
         P(27) = 0.97621D-03
         E(28) = 0.26806D+03
         P(28) = 0.81260D-02
         E(29) = 0.86222D+02
         P(29) = 0.41433D-01
         E(30) = 0.31979D+02
         P(30) = 0.14035D+00
         E(31) = 0.13008D+02
         P(31) = 0.31207D+00
         E(32) = 0.54766D+01
         P(32) = 0.42099D+00
         E(33) = 0.23090D+01
         P(33) = 0.25664D+00
         E(34) = 0.26806D+03
         P(34) = -.11596D-04
         E(35) = 0.86222D+02
         P(35) = -.18289D-03
         E(36) = 0.31979D+02
         P(36) = -.17170D-02
         E(37) = 0.13008D+02
         P(37) = -.78974D-02
         E(38) = 0.54766D+01
         P(38) = -.18237D-01
         E(39) = 0.23090D+01
         P(39) = 0.40035D-01
         E(40) = 0.84601D+00
         P(40) = 0.36717D+00
         E(41) = 0.31659D+00
         P(41) = 0.10000D+01
         E(42) = 0.10716D+00
         P(42) = 0.10000D+01
         E(43) = 0.23100D+01
         D(43) = 0.10000D+01
         E(44) = 0.52000D+00
         D(44) = 0.10000D+01
         E(45) = 0.72000D+00
         F(45) = 0.10000D+01
         ILAST = 45
      ELSE IF (NPOL.EQ.3) THEN
         E(1) = 0.14611D+07
         S(1) = 0.11534D-04
         E(2) = 0.21959D+06
         S(2) = 0.89235D-04
         E(3) = 0.50085D+05
         S(3) = 0.46811D-03
         E(4) = 0.14206D+05
         S(4) = 0.19692D-02
         E(5) = 0.46389D+04
         S(5) = 0.71112D-02
         E(6) = 0.16763D+04
         S(6) = 0.22625D-01
         E(7) = 0.65451D+03
         S(7) = 0.63488D-01
         E(8) = 0.27153D+03
         S(8) = 0.15173D+00
         E(9) = 0.11787D+03
         S(9) = 0.28826D+00
         E(10) = 0.52984D+02
         S(10) = 0.36113D+00
         E(11) = 0.24292D+02
         S(11) = 0.22303D+00
         E(12) = 0.10184D+02
         S(12) = 0.34031D-01
         E(13) = 0.14206D+05
         S(13) = -.27093D-05
         E(14) = 0.46389D+04
         S(14) = -.43191D-05
         E(15) = 0.16763D+04
         S(15) = -.10907D-03
         E(16) = 0.65451D+03
         S(16) = -.48853D-03
         E(17) = 0.27153D+03
         S(17) = -.36327D-02
         E(18) = 0.11787D+03
         S(18) = -.14142D-01
         E(19) = 0.52984D+02
         S(19) = -.46167D-01
         E(20) = 0.24292D+02
         S(20) = -.50107D-01
         E(21) = 0.10184D+02
         S(21) = 0.99999D-01
         E(22) = 0.46928D+01
         S(22) = 0.23750D+00
         E(23) = 0.46389D+04
         S(23) = 0.28740D-05
         E(24) = 0.16763D+04
         S(24) = 0.24949D-04
         E(25) = 0.65451D+03
         S(25) = 0.18272D-03
         E(26) = 0.27153D+03
         S(26) = 0.10821D-02
         E(27) = 0.11787D+03
         S(27) = 0.49846D-02
         E(28) = 0.52984D+02
         S(28) = 0.15414D-01
         E(29) = 0.24292D+02
         S(29) = 0.21774D-01
         E(30) = 0.10184D+02
         S(30) = -.48881D-01
         E(31) = 0.46928D+01
         S(31) = -.21808D+00
         E(32) = 0.20974D+01
         S(32) = -.27412D+00
         E(33) = 0.74867D+00
         S(33) = 0.10000D+01
         E(34) = 0.33613D+00
         S(34) = 0.10000D+01
         E(35) = 0.13573D+00
         S(35) = 0.10000D+01
         E(36) = 0.38457D+04
         P(36) = 0.11653D-03
         E(37) = 0.91123D+03
         P(37) = 0.10250D-02
         E(38) = 0.29564D+03
         P(38) = 0.57481D-02
         E(39) = 0.11140D+03
         P(39) = 0.24700D-01
         E(40) = 0.46767D+02
         P(40) = 0.76263D-01
         E(41) = 0.21106D+02
         P(41) = 0.18343D+00
         E(42) = 0.99058D+01
         P(42) = 0.31219D+00
         E(43) = 0.47347D+01
         P(43) = 0.35836D+00
         E(44) = 0.22364D+01
         P(44) = 0.21176D+00
         E(45) = 0.91123D+03
         P(45) = 0.35761D-05
         E(46) = 0.29564D+03
         P(46) = -.77768D-05
         E(47) = 0.11140D+03
         P(47) = 0.67666D-05
         E(48) = 0.46767D+02
         P(48) = -.69394D-03
         E(49) = 0.21106D+02
         P(49) = -.22563D-02
         E(50) = 0.99058D+01
         P(50) = -.11018D-01
         E(51) = 0.47347D+01
         P(51) = -.93663D-02
         E(52) = 0.22364D+01
         P(52) = 0.25171D-01
         E(53) = 0.10071D+01
         P(53) = 0.25574D+00
         E(54) = 0.44369D+00
         P(54) = 0.10000D+01
         E(55) = 0.18488D+00
         P(55) = 0.10000D+01
         E(56) = 0.73268D-01
         P(56) = 0.10000D+01
         E(57) = 0.98100D+01
         D(57) = 0.10000D+01
         E(58) = 0.25800D+01
         D(58) = 0.10000D+01
         E(59) = 0.76000D+00
         D(59) = 0.10000D+01
         E(60) = 0.28000D+00
         D(60) = 0.10000D+01
         E(61) = 0.17200D+01
         F(61) = 0.10000D+01
         E(62) = 0.57000D+00
         F(62) = 0.10000D+01
         E(63) = 0.86000D+00
         G(63) = 0.10000D+01
         ILAST = 63
      ELSE IF (NPOL.EQ.4) THEN
         E(1) = 0.11411D+07
         S(1) = 0.11243D-04
         E(2) = 0.26443D+06
         S(2) = 0.58000D-04
         E(3) = 0.75791D+05
         S(3) = 0.24215D-03
         E(4) = 0.24925D+05
         S(4) = 0.87550D-03
         E(5) = 0.90526D+04
         S(5) = 0.28376D-02
         E(6) = 0.35481D+04
         S(6) = 0.84086D-02
         E(7) = 0.14780D+04
         S(7) = 0.22906D-01
         E(8) = 0.64735D+03
         S(8) = 0.56895D-01
         E(9) = 0.29533D+03
         S(9) = 0.12569D+00
         E(10) = 0.13898D+03
         S(10) = 0.23327D+00
         E(11) = 0.67549D+02
         S(11) = 0.31894D+00
         E(12) = 0.33661D+02
         S(12) = 0.27117D+00
         E(13) = 0.17027D+02
         S(13) = 0.10247D+00
         E(14) = 0.81504D+01
         S(14) = 0.83067D-02
         E(15) = 0.35481D+04
         S(15) = 0.58169D-05
         E(16) = 0.14780D+04
         S(16) = 0.26634D-04
         E(17) = 0.64735D+03
         S(17) = 0.19172D-03
         E(18) = 0.29533D+03
         S(18) = 0.83597D-03
         E(19) = 0.13898D+03
         S(19) = 0.37502D-02
         E(20) = 0.67549D+02
         S(20) = 0.99495D-02
         E(21) = 0.33661D+02
         S(21) = 0.20640D-01
         E(22) = 0.17027D+02
         S(22) = 0.52063D-03
         E(23) = 0.81504D+01
         S(23) = -.50048D-01
         E(24) = 0.41319D+01
         S(24) = -.83332D-01
         E(25) = 0.35481D+04
         S(25) = -.58469D-04
         E(26) = 0.14780D+04
         S(26) = -.37803D-03
         E(27) = 0.64735D+03
         S(27) = -.21881D-02
         E(28) = 0.29533D+03
         S(28) = -.10582D-01
         E(29) = 0.13898D+03
         S(29) = -.44410D-01
         E(30) = 0.67549D+02
         S(30) = -.12491D+00
         E(31) = 0.33661D+02
         S(31) = -.24877D+00
         E(32) = 0.17027D+02
         S(32) = -.29947D-01
         E(33) = 0.81504D+01
         S(33) = 0.72881D+00
         E(34) = 0.41319D+01
         S(34) = 0.13743D+01
         E(35) = 0.20457D+01
         S(35) = 0.78685D+00
         E(36) = 0.94029D+00
         S(36) = 0.10000D+01
         E(37) = 0.46901D+00
         S(37) = 0.10000D+01
         E(38) = 0.22219D+00
         S(38) = 0.10000D+01
         E(39) = 0.10220D+00
         S(39) = 0.10000D+01
         E(40) = 0.11035D+05
         P(40) = 0.18514D-04
         E(41) = 0.26115D+04
         P(41) = 0.16527D-03
         E(42) = 0.84802D+03
         P(42) = 0.96233D-03
         E(43) = 0.32300D+03
         P(43) = 0.43777D-02
         E(44) = 0.13533D+03
         P(44) = 0.16218D-01
         E(45) = 0.62684D+02
         P(45) = 0.46104D-01
         E(46) = 0.30298D+02
         P(46) = 0.11348D+00
         E(47) = 0.15240D+02
         P(47) = 0.21353D+00
         E(48) = 0.79259D+01
         P(48) = 0.29438D+00
         E(49) = 0.42853D+01
         P(49) = 0.28778D+00
         E(50) = 0.23043D+01
         P(50) = 0.18269D+00
         E(51) = 0.32300D+03
         P(51) = 0.92007D-05
         E(52) = 0.13533D+03
         P(52) = -.51118D-04
         E(53) = 0.62684D+02
         P(53) = -.10597D-03
         E(54) = 0.30298D+02
         P(54) = -.16689D-02
         E(55) = 0.15240D+02
         P(55) = -.30998D-02
         E(56) = 0.79259D+01
         P(56) = -.14628D-01
         E(57) = 0.42853D+01
         P(57) = 0.12345D-02
         E(58) = 0.23043D+01
         P(58) = 0.15315D-02
         E(59) = 0.12357D+01
         P(59) = 0.16993D+00
         E(60) = 0.60021D+00
         P(60) = 0.10000D+01
         E(61) = 0.28739D+00
         P(61) = 0.10000D+01
         E(62) = 0.13185D+00
         P(62) = 0.10000D+01
         E(63) = 0.56970D-01
         P(63) = 0.10000D+01
         E(64) = 0.24630D+02
         D(64) = 0.10000D+01
         E(65) = 0.70400D+01
         D(65) = 0.10000D+01
         E(66) = 0.24300D+01
         D(66) = 0.10000D+01
         E(67) = 0.87000D+00
         D(67) = 0.10000D+01
         E(68) = 0.37700D+00
         D(68) = 0.10000D+01
         E(69) = 0.15700D+00
         D(69) = 0.10000D+01
         E(70) = 0.40600D+01
         F(70) = 0.10000D+01
         E(71) = 0.95000D+00
         F(71) = 0.10000D+01
         E(72) = 0.39000D+00
         F(72) = 0.10000D+01
         E(73) = 0.19700D+01
         G(73) = 0.10000D+01
         E(74) = 0.66000D+00
         G(74) = 0.10000D+01
CFRJ     E(75) = 0.10400D+01
CFRJ     H(75) = 0.10000D+01
         ILAST = 74
      ENDIF
      RETURN
C
C     ----- AR -----
C
  240 CONTINUE
      IF (MASWRK) WRITE(6,*) 'PC-N BASIS DOES NOT EXIST FOR AR.'
      CALL BERROR(3)
      STOP
      END
C*MODULE BASPCN  *DECK APCN01
      SUBROUTINE APCN01(E,S,P,D,F,N,NPOL,ILAST)
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      DIMENSION E(90),S(90),P(90),D(90),F(90)
C
      LOGICAL GOPARR,DSKWRK,MASWRK
C
      COMMON /PAR   / ME,MASTER,NPROC,IBTYP,IPTIM,GOPARR,DSKWRK,MASWRK
C
      GO TO (100,120),N
C
C     ----- HYDROGEN -----
C
  100 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(ILAST+1) = 0.27174D-01
         S(ILAST+1) = 0.10000D+01
      ELSE IF (NPOL.EQ.1) THEN
         E(ILAST+1) = 0.25628D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.76884D-01
         P(ILAST+2) = 0.10000D+01
      ELSE IF (NPOL.EQ.2) THEN
         E(ILAST+1) = 0.22334D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.67001D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.89334D-01
         D(ILAST+3) = 0.10000D+01
      ELSE IF (NPOL.EQ.3) THEN
         E(ILAST+1) = 0.18800D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.56399D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.75199D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.93999D-01
         F(ILAST+4) = 0.10000D+01
      ELSE IF (NPOL.EQ.4) THEN
         E(ILAST+1) = 0.17131D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.51393D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.68524D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.85656D-01
         F(ILAST+4) = 0.10000D+01
CFRJ     E(ILAST+5) = 0.10279D+00
CFRJ     G(ILAST+5) = 0.10000D+01
      ENDIF
      RETURN
C
C     ----- HE -----
C
  120 CONTINUE
      IF (MASWRK) WRITE(6,*) 'AUG-PC-N BASIS DOES NOT EXIST FOR HE'
      CALL BERROR(3)
      STOP
      END
C*MODULE BASPCN  *DECK APCN02
      SUBROUTINE APCN02(E,S,P,D,F,G,N,NPOL,ILAST)
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      DIMENSION E(90),S(90),P(90),D(90),F(90),G(90)
C
      LOGICAL GOPARR,DSKWRK,MASWRK
C
      COMMON /PAR   / ME,MASTER,NPROC,IBTYP,IPTIM,GOPARR,DSKWRK,MASWRK
C
      NN = N-2
      GO TO (100,120,140,160,180,200,220,240),NN
C
  100 CONTINUE
  120 CONTINUE
  140 CONTINUE
      IF (MASWRK) WRITE(6,*) 'AUG-PC-N BASIS DOES NOT EXIST FOR LI,BE,B'
      CALL BERROR(3)
      STOP
C
C     ----- C  -----
C
  160 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(ILAST+1) = 0.69748D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.39361D-01
         P(ILAST+2) = 0.10000D+01
      ELSE IF (NPOL.EQ.1) THEN
         E(ILAST+1) = 0.44856D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.33518D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.10055D+00
         D(ILAST+3) = 0.10000D+01
      ELSE IF (NPOL.EQ.2) THEN
         E(ILAST+1) = 0.37217D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.26449D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.79348D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.10580D+00
         F(ILAST+4) = 0.10000D+01
      ELSE IF (NPOL.EQ.3) THEN
         E(ILAST+1) = 0.31605D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.18873D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.56619D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.75492D-01
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.94365D-01
         G(ILAST+5) = 0.10000D+01
      ELSE IF (NPOL.EQ.4) THEN
         E(ILAST+1) = 0.26452D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.16207D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.48622D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.64830D-01
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.81037D-01
         G(ILAST+5) = 0.10000D+01
CFRJ     E(ILAST+6) = 0.97244D-01
CFRJ     H(ILAST+6) = 0.10000D+01
      ENDIF
      RETURN
C
C     ----- N -----
C
  180 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(ILAST+1) = 0.97986D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.56709D-01
         P(ILAST+2) = 0.10000D+01
      ELSE IF (NPOL.EQ.1) THEN
         E(ILAST+1) = 0.62739D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.47707D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.14312D+00
         D(ILAST+3) = 0.10000D+01
      ELSE IF (NPOL.EQ.2) THEN
         E(ILAST+1) = 0.52165D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.36023D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.10807D+00
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.14409D+00
         F(ILAST+4) = 0.10000D+01
      ELSE IF (NPOL.EQ.3) THEN
         E(ILAST+1) = 0.44444D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.26388D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.79164D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.10555D+00
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.13194D+00
         G(ILAST+5) = 0.10000D+01
      ELSE IF (NPOL.EQ.4) THEN
         E(ILAST+1) = 0.37542D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.22684D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.68053D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.90737D-01
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.11342D+00
         G(ILAST+5) = 0.10000D+01
CFRJ     E(ILAST+6) = 0.13611D+00
CFRJ     H(ILAST+6) = 0.10000D+01
      ENDIF
      RETURN
C
C     ----- O  ------
C
  200 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(ILAST+1) = 0.13165D+00
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.70416D-01
         P(ILAST+2) = 0.10000D+01
      ELSE IF (NPOL.EQ.1) THEN
         E(ILAST+1) = 0.83425D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.57616D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.17285D+00
         D(ILAST+3) = 0.10000D+01
      ELSE IF (NPOL.EQ.2) THEN
         E(ILAST+1) = 0.68141D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.41871D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.12561D+00
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.16748D+00
         F(ILAST+4) = 0.10000D+01
      ELSE IF (NPOL.EQ.3) THEN
         E(ILAST+1) = 0.57690D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.29477D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.88431D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.11791D+00
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.14739D+00
         G(ILAST+5) = 0.10000D+01
      ELSE IF (NPOL.EQ.4) THEN
         E(ILAST+1) = 0.49094D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.24967D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.74902D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.99870D-01
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.12484D+00
         G(ILAST+5) = 0.10000D+01
CFRJ     E(ILAST+6) = 0.14980D+00
CFRJ     H(ILAST+6) = 0.10000D+01
      ENDIF
      RETURN
C
C     ----- F  -----
C
  220 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(ILAST+1) = 0.16931D+00
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.89029D-01
         P(ILAST+2) = 0.10000D+01
      ELSE IF (NPOL.EQ.1) THEN
         E(ILAST+1) = 0.10802D+00
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.72054D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.21616D+00
         D(ILAST+3) = 0.10000D+01
      ELSE IF (NPOL.EQ.2) THEN
         E(ILAST+1) = 0.88276D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.52376D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.15713D+00
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.20951D+00
         F(ILAST+4) = 0.10000D+01
      ELSE IF (NPOL.EQ.3) THEN
         E(ILAST+1) = 0.74924D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.37856D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.11357D+00
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.15142D+00
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.18928D+00
         G(ILAST+5) = 0.10000D+01
      ELSE IF (NPOL.EQ.4) THEN
         E(ILAST+1) = 0.65313D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.32408D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.97223D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.12963D+00
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.16204D+00
         G(ILAST+5) = 0.10000D+01
CFRJ     E(ILAST+6) = 0.19445D+00
CFRJ     H(ILAST+6) = 0.10000D+01
      ENDIF
      RETURN
C
C     ----- NE -----
C
  240 CONTINUE
      IF (MASWRK) WRITE(6,*) 'AUG-PC-N BASIS DOES NOT EXIST FOR NE.'
      CALL BERROR(3)
      STOP
      END
C*MODULE BASPCN  *DECK APCN03
      SUBROUTINE APCN03(E,S,P,D,F,G,N,NPOL,ILAST)
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      DIMENSION E(90),S(90),P(90),D(90),F(90),G(90)
      LOGICAL GOPARR,DSKWRK,MASWRK
C
      COMMON /PAR   / ME,MASTER,NPROC,IBTYP,IPTIM,GOPARR,DSKWRK,MASWRK
C
      NN = N-10
      GO TO (100,120,140,160,180,200,220,240),NN
C
  100 CONTINUE
  120 CONTINUE
  140 CONTINUE
      IF(MASWRK)WRITE(6,*)'AUG-PC-N BASIS DOES NOT EXIST FOR NA,MG,AL.'
      CALL BERROR(3)
      STOP
C
C     ----- SI -----
C
  160 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(ILAST+1) = 0.42730D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.32666D-01
         P(ILAST+2) = 0.10000D+01
      ELSE IF (NPOL.EQ.1) THEN
         E(ILAST+1) = 0.33293D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.22111D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.66332D-01
         D(ILAST+3) = 0.10000D+01
      ELSE IF (NPOL.EQ.2) THEN
         E(ILAST+1) = 0.30878D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.17527D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.52580D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.70107D-01
         F(ILAST+4) = 0.10000D+01
      ELSE IF (NPOL.EQ.3) THEN
         E(ILAST+1) = 0.26105D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.14560D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.43680D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.58240D-01
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.72800D-01
         G(ILAST+5) = 0.10000D+01
      ELSE IF (NPOL.EQ.4) THEN
         E(ILAST+1) = 0.21897D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.11442D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.34326D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.45768D-01
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.57211D-01
         G(ILAST+5) = 0.10000D+01
CFRJ     E(ILAST+6) = 0.68653D-01
CFRJ     H(ILAST+6) = 0.10000D+01
      ENDIF
      RETURN
C
C     ----- P  -----
C
  180 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(ILAST+1) = 0.56069D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.37569D-01
         P(ILAST+2) = 0.10000D+01
      ELSE IF (NPOL.EQ.1) THEN
         E(ILAST+1) = 0.43194D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.30386D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.91157D-01
         D(ILAST+3) = 0.10000D+01
      ELSE IF (NPOL.EQ.2) THEN
         E(ILAST+1) = 0.40548D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.23848D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.71543D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.95390D-01
         F(ILAST+4) = 0.10000D+01
      ELSE IF (NPOL.EQ.3) THEN
         E(ILAST+1) = 0.33408D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.18800D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.56399D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.75199D-01
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.93999D-01
         G(ILAST+5) = 0.10000D+01
      ELSE IF (NPOL.EQ.4) THEN
         E(ILAST+1) = 0.28607D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.15693D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.47080D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.62774D-01
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.78467D-01
         G(ILAST+5) = 0.10000D+01
CFRJ     E(ILAST+6) = 0.94161D-01
CFRJ     H(ILAST+6) = 0.10000D+01
      ENDIF
      RETURN
C
C     ----- S  -----
C
  200 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(ILAST+1) = 0.70749D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.43067D-01
         P(ILAST+2) = 0.10000D+01
      ELSE IF (NPOL.EQ.1) THEN
         E(ILAST+1) = 0.54100D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.35364D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.10609D+00
         D(ILAST+3) = 0.10000D+01
      ELSE IF (NPOL.EQ.2) THEN
         E(ILAST+1) = 0.51037D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.27776D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.83327D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.11110D+00
         F(ILAST+4) = 0.10000D+01
      ELSE IF (NPOL.EQ.3) THEN
         E(ILAST+1) = 0.41837D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.21820D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.65461D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.87282D-01
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.10910D+00
         G(ILAST+5) = 0.10000D+01
      ELSE IF (NPOL.EQ.4) THEN
         E(ILAST+1) = 0.35487D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.18036D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.54107D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.72142D-01
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.90178D-01
         G(ILAST+5) = 0.10000D+01
CFRJ     E(ILAST+6) = 0.10821D+00
CFRJ     H(ILAST+6) = 0.10000D+01
      ENDIF
      RETURN
C
C     ----- CL -----
C
  220 CONTINUE
      IF (NPOL.EQ.0) THEN
         E(ILAST+1) = 0.86577D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.51427D-01
         P(ILAST+2) = 0.10000D+01
      ELSE IF (NPOL.EQ.1) THEN
         E(ILAST+1) = 0.65816D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.41467D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.12440D+00
         D(ILAST+3) = 0.10000D+01
      ELSE IF (NPOL.EQ.2) THEN
         E(ILAST+1) = 0.62241D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.33970D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.10191D+00
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.13588D+00
         F(ILAST+4) = 0.10000D+01
      ELSE IF (NPOL.EQ.3) THEN
         E(ILAST+1) = 0.50711D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.26904D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.80711D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.10761D+00
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.13452D+00
         G(ILAST+5) = 0.10000D+01
      ELSE IF (NPOL.EQ.4) THEN
         E(ILAST+1) = 0.43046D-01
         S(ILAST+1) = 0.10000D+01
         E(ILAST+2) = 0.22657D-01
         P(ILAST+2) = 0.10000D+01
         E(ILAST+3) = 0.67972D-01
         D(ILAST+3) = 0.10000D+01
         E(ILAST+4) = 0.90629D-01
         F(ILAST+4) = 0.10000D+01
         E(ILAST+5) = 0.11329D+00
         G(ILAST+5) = 0.10000D+01
CFRJ     E(ILAST+6) = 0.13594D+00
CFRJ     H(ILAST+6) = 0.10000D+01
      ENDIF
      RETURN
C
C     ----- AR -----
C
  240 CONTINUE
      IF (MASWRK) WRITE(6,*) 'AUG-PC-N BASIS DOES NOT EXIST FOR AR.'
      CALL BERROR(3)
      STOP
      END
