C 19 SEP 05 - MK  - INCLUDE MODEL CORE POTENTIAL LIBRARY
C
C*MODULE MCPLIB  *DECK TM2130
      SUBROUTINE TM2130(MMPMOD,NUCZ,MMPTYP,SHINFO,CDFCTR,NOSH,KTYP,KNG,
     2                  EX,C,ZCORE,NOAN0,NOAN1,AN0,ALPN0,AN1,
     3                  ALPN1,BPAR,ITYP,IERR)
C
C     --- MMPMOD: 0  READ VALENCE BASIS-SET
C     ---         1  READ MODEL-POTENTIAL-PARAMETERS
C     ---            AND CORE-SHELL BASIS-SET
C
C     --- NUCZ:    NUCLEAR CHARGE
C     --- MMPTYP:  1  NON-RELATIVISTIC   MOD.POT.
C     ---          2  QUASI-RELATIVISTIC MOD.POT.
C     ---          3  READ IN MOD.POT.FROM $MMP GROUP
C     ---          4  TYPE NOT SPECIFIED
C     ---    W*10**I  FOR I-TH SHELL (I>=1)
C     ---    W: 1,2,3,4,6  S,P,D,F,L - SHELL
C     --- CDFCTR(I) CODED DEFAULT CONTRACTION OF VALENCE-SHELL I
C     --- NOSH     NO.OF (VALENCE OR CORE) SHELLS
C     --- KTYP(I)  TYPE OF SHELL I
C     --- KNG(I)   NO.OF GTFS IN SHELL I
C     --- EX(J)    EXPONENT OF J-TH GTF
C     --- C(J)     CONTRACTION COEFF.OF J-TH GTF
C     --- ZCORE    CORE-CHARGE
C     --- NOAN0    NO.OF A-TERMS FOR N=0
C     --- NOAN1    NO.OF A-TERMS FOR N=1
C     --- AN0(I)   A-PARAMETERS FOR N=0
C     --- ALPN0(I) ALPHA-EXPONENTS FOR N=0
C     --- AN1(I)   A-PARAMETERS FOR N=1
C     --- ALPN1(I) ALPHA-EXPONENTS FOR N=1
C     --- BPAR(I)  B-PARAMETERS
C     --- IERR(I)  ERROR-STATUS   FOR MMPMOD=0:
C     ---          0  OK; 1  NOT MMP-ATOM; 2 MMPTYP WRONG
C
C
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      DIMENSION AN0(*),ALPN0(*),AN1(*),ALPN1(*)
      DIMENSION BPAR(*),KNG(*),KTYP(*)
      DIMENSION C(*),EX(*)
      CHARACTER*20 CDFCTR(*),SHINFO
C
      DIMENSION ITYP(10)
      COMMON /IOFILE/ IR,IW,IP,IJK,IPK,IDAF,NAV,IODA(950)
C
      IZSAVE=NUCZ
      IERR=0
C
C     --- ITYP(I) IS THE TYPE (MMPTYP) OF THE I-TH BUILT-IN
C     --- MODEL-POTENTIAL FOR A GIVEN ATOM (IZSAVE)
C
      DO 40 I=1,10
         ITYP(I)=0
  40  CONTINUE
C ... ../TOOLS_MCP/TM2130:V-ZN.DSP.MPPAR
C ... ../TOOLS_MCP/TM2130:V-ZN.DS.MPPAR
C
C     --- 
C
      IF(IZSAVE.EQ.21) THEN
        ITYP(1)=111
        ITYP(2)=112
C     --- 
        IF (MMPTYP.EQ.112) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : SC/NR             CONFIGURATION : S1P6D02      *
C   VALENCE : (4S/3P/3D)        BASIS SET     : ( 6/ 4/ 5)   *
C          EFFECTIVE CORE CHARGE :  9.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =     -29.589342 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.1925 A.U.
C  ***     ORBITAL ENERGY( 3P ) =   -1.4374 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.2153 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE -  9.0D+00
          SHINFO='4S/3P/3D'
          NOSH= 3
          CDFCTR(1)='51'
          CDFCTR(2)='4'
          CDFCTR(3)='41'
             KNG(1)= 6
            KTYP(1)= 1
             KNG(2)= 4
            KTYP(2)= 2
             KNG(3)= 5
            KTYP(3)= 3
          EX( 1)= 0.16097900D+03
           C( 1)=-0.24468195D-01
          EX( 2)= 0.61185900D+01
           C( 2)= 0.13295449D+00
          EX( 3)= 0.85294300D+00
           C( 3)=-0.25091526D+00
          EX( 4)= 0.37138700D+00
           C( 4)=-0.22218935D+00
          EX( 5)= 0.62370400D-01
           C( 5)= 0.65261211D+00
          EX( 6)= 0.24988500D-01
           C( 6)= 0.48876012D+00
          EX( 7)= 0.65552900D+02
           C( 7)=-0.81279696D-01
          EX( 8)= 0.13325100D+02
           C( 8)=-0.27231868D+00
          EX( 9)= 0.11931800D+01
           C( 9)= 0.68791025D+00
          EX(10)= 0.37973100D+00
           C(10)= 0.39609531D+00
          EX(11)= 0.12969800D+02
           C(11)= 0.44836600D-01
          EX(12)= 0.31060700D+01
           C(12)= 0.21104768D+00
          EX(13)= 0.92017600D+00
           C(13)= 0.41224799D+00
          EX(14)= 0.26976100D+00
           C(14)= 0.45887073D+00
          EX(15)= 0.75763100D-01
           C(15)= 0.28235395D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE -  9.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.70003394D+00
          ALPN0( 1)= 0.23007947D+03
            AN0( 2)= 0.30875239D+00
          ALPN0( 2)= 0.55847240D+01
            AN0( 3)= 0.73325900D-01
          ALPN0( 3)= 0.15574928D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.67433688D+02
          ALPN1( 1)= 0.16800002D+04
            AN1( 2)= 0.54557319D+01
          ALPN1( 2)= 0.38511091D+02
            AN1( 3)= 0.85561099D-01
          ALPN1( 3)= 0.88141537D+00
C         --- CORE-SHELL INFORMATION 
          NOSH     =  4
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.33151726D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.37834560D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.48275600D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.31017380D+02
C         --- EXPONENTS AND COEFFICIENTS 
          EX(  1) = 0.58331415D+04
           C(  1) = 0.21080794D-01
          EX(  2) = 0.80277191D+03
           C(  2) = 0.14795427D+00
          EX(  3) = 0.18297470D+03
           C(  3) = 0.47377479D+00
          EX(  4) = 0.51828595D+02
           C(  4) = 0.48520354D+00
          EX(  5) = 0.84177727D+02
           C(  5) =-0.32519930D-02
          EX(  6) = 0.89188971D+01
           C(  6) = 0.19914252D-01
          EX(  7) = 0.35584223D+01
           C(  7) = 0.13477947D-01
          EX(  8) = 0.61328831D+01
           C(  8) = 0.40893231D-02
          EX(  9) = 0.92967343D+00
           C(  9) =-0.13880496D-01
          EX( 10) = 0.36390388D+00
           C( 10) =-0.69466526D-02
          EX( 11) = 0.58331415D+04
           C( 11) = 0.67665573D-02
          EX( 12) = 0.80277191D+03
           C( 12) = 0.47490671D-01
          EX( 13) = 0.18297470D+03
           C( 13) = 0.15207322D+00
          EX( 14) = 0.51828595D+02
           C( 14) = 0.15574165D+00
          EX( 15) = 0.84177727D+02
           C( 15) = 0.10848643D+00
          EX( 16) = 0.89188971D+01
           C( 16) =-0.66433914D+00
          EX( 17) = 0.35584223D+01
           C( 17) =-0.44962407D+00
          EX( 18) = 0.61328831D+01
           C( 18) = 0.23844899D-02
          EX( 19) = 0.92967343D+00
           C( 19) =-0.80937359D-02
          EX( 20) = 0.36390388D+00
           C( 20) =-0.40506026D-02
          EX( 21) = 0.58331415D+04
           C( 21) = 0.23791586D-02
          EX( 22) = 0.80277191D+03
           C( 22) = 0.16697980D-01
          EX( 23) = 0.18297470D+03
           C( 23) = 0.53469778D-01
          EX( 24) = 0.51828595D+02
           C( 24) = 0.54759617D-01
          EX( 25) = 0.84177727D+02
           C( 25) = 0.43708344D-01
          EX( 26) = 0.89188971D+01
           C( 26) =-0.26765709D+00
          EX( 27) = 0.35584223D+01
           C( 27) =-0.18115005D+00
          EX( 28) = 0.61328831D+01
           C( 28) =-0.23771407D+00
          EX( 29) = 0.92967343D+00
           C( 29) = 0.80687901D+00
          EX( 30) = 0.36390388D+00
           C( 30) = 0.40381183D+00
          EX( 31) = 0.22359926D+03
           C( 31) = 0.34946326D-01
          EX( 32) = 0.48663060D+02
           C( 32) = 0.21788808D+00
          EX( 33) = 0.14560751D+02
           C( 33) = 0.52826438D+00
          EX( 34) = 0.48693613D+01
           C( 34) = 0.39663056D+00
          EX( 35) = 0.19888962D+01
           C( 35) = 0.74227060D-02
          EX( 36) = 0.73403475D+00
           C( 36) = 0.11742959D-01
          EX( 37) = 0.26930557D+00
           C( 37) = 0.41811866D-02
          RETURN
          END IF
C     --- 
        ELSE IF (MMPTYP.EQ.111) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : SC/NR             CONFIGURATION : S1D02        *
C   VALENCE : (4S/3D)           BASIS SET     : ( 5/ 5)      *
C          EFFECTIVE CORE CHARGE :  3.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =      -1.521706 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.1921 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.2064 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE -  3.0D+00
          SHINFO='4S/3D'
          NOSH= 2
          CDFCTR(1)='41'
          CDFCTR(2)='41'
             KNG(1)= 5
            KTYP(1)= 1
             KNG(2)= 5
            KTYP(2)= 3
          EX( 1)= 0.40579202D+01
           C( 1)=-0.12391769D+00
          EX( 2)= 0.85292310D+00
           C( 2)= 0.26914023D+00
          EX( 3)= 0.37137860D+00
           C( 3)= 0.16354860D+00
          EX( 4)= 0.62371880D-01
           C( 4)=-0.62403507D+00
          EX( 5)= 0.24989170D-01
           C( 5)=-0.50737733D+00
          EX( 6)= 0.12969848D+02
           C( 6)= 0.59271583D-01
          EX( 7)= 0.31060727D+01
           C( 7)= 0.20533118D+00
          EX( 8)= 0.92017600D+00
           C( 8)= 0.41163944D+00
          EX( 9)= 0.26976060D+00
           C( 9)= 0.45268301D+00
          EX(10)= 0.75763100D-01
           C(10)= 0.29289717D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE -  3.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.23818232D+01
          ALPN0( 1)= 0.16652421D+03
            AN0( 2)= 0.19000483D+01
          ALPN0( 2)= 0.43497475D+01
            AN0( 3)= 0.86357986D+00
          ALPN0( 3)= 0.75041147D+00
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.15506982D+02
          ALPN1( 1)= 0.45742252D+04
            AN1( 2)= 0.26553637D+02
          ALPN1( 2)= 0.32875488D+02
            AN1( 3)= 0.31397297D+00
          ALPN1( 3)= 0.70406875D+00
C         --- CORE-SHELL INFORMATION
          NOSH     =  5
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.21548622D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.24592464D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.31379140D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.20161297D+02
          KTYP( 5) =  2
           KNG( 5) =  7
          BPAR( 5) = 0.18666934D+01
C         --- EXPONENTS AND COEFFICIENTS
          EX(  1) = 0.58331415D+04
           C(  1) = 0.21080794D-01
          EX(  2) = 0.80277191D+03
           C(  2) = 0.14795427D+00
          EX(  3) = 0.18297470D+03
           C(  3) = 0.47377479D+00
          EX(  4) = 0.51828595D+02
           C(  4) = 0.48520354D+00
          EX(  5) = 0.84177727D+02
           C(  5) =-0.32519930D-02
          EX(  6) = 0.89188971D+01
           C(  6) = 0.19914252D-01
          EX(  7) = 0.35584223D+01
           C(  7) = 0.13477947D-01
          EX(  8) = 0.61328831D+01
           C(  8) = 0.40893231D-02
          EX(  9) = 0.92967343D+00
           C(  9) =-0.13880496D-01
          EX( 10) = 0.36390388D+00
           C( 10) =-0.69466526D-02
          EX( 11) = 0.58331415D+04
           C( 11) = 0.67665573D-02
          EX( 12) = 0.80277191D+03
           C( 12) = 0.47490671D-01
          EX( 13) = 0.18297470D+03
           C( 13) = 0.15207322D+00
          EX( 14) = 0.51828595D+02
           C( 14) = 0.15574165D+00
          EX( 15) = 0.84177727D+02
           C( 15) = 0.10848643D+00
          EX( 16) = 0.89188971D+01
           C( 16) =-0.66433914D+00
          EX( 17) = 0.35584223D+01
           C( 17) =-0.44962407D+00
          EX( 18) = 0.61328831D+01
           C( 18) = 0.23844899D-02
          EX( 19) = 0.92967343D+00
           C( 19) =-0.80937359D-02
          EX( 20) = 0.36390388D+00
           C( 20) =-0.40506026D-02
          EX( 21) = 0.58331415D+04
           C( 21) = 0.23791586D-02
          EX( 22) = 0.80277191D+03
           C( 22) = 0.16697980D-01
          EX( 23) = 0.18297470D+03
           C( 23) = 0.53469778D-01
          EX( 24) = 0.51828595D+02
           C( 24) = 0.54759617D-01
          EX( 25) = 0.84177727D+02
           C( 25) = 0.43708344D-01
          EX( 26) = 0.89188971D+01
           C( 26) =-0.26765709D+00
          EX( 27) = 0.35584223D+01
           C( 27) =-0.18115005D+00
          EX( 28) = 0.61328831D+01
           C( 28) =-0.23771407D+00
          EX( 29) = 0.92967343D+00
           C( 29) = 0.80687901D+00
          EX( 30) = 0.36390388D+00
           C( 30) = 0.40381183D+00
          EX( 31) = 0.22359926D+03
           C( 31) = 0.34946326D-01
          EX( 32) = 0.48663060D+02
           C( 32) = 0.21788808D+00
          EX( 33) = 0.14560751D+02
           C( 33) = 0.52826438D+00
          EX( 34) = 0.48693613D+01
           C( 34) = 0.39663056D+00
          EX( 35) = 0.19888962D+01
           C( 35) = 0.74227060D-02
          EX( 36) = 0.73403475D+00
           C( 36) = 0.11742959D-01
          EX( 37) = 0.26930557D+00
           C( 37) = 0.41811866D-02
          EX( 38) = 0.22359926D+03
           C( 38) = 0.12687833D-01
          EX( 39) = 0.48663060D+02
           C( 39) = 0.79107815D-01
          EX( 40) = 0.14560751D+02
           C( 40) = 0.19179499D+00
          EX( 41) = 0.48693613D+01
           C( 41) = 0.14400319D+00
          EX( 42) = 0.19888962D+01
           C( 42) =-0.37661038D+00
          EX( 43) = 0.73403475D+00
           C( 43) =-0.59580968D+00
          EX( 44) = 0.26930557D+00
           C( 44) =-0.21214342D+00
          RETURN
          END IF
        END IF
        IERR=2
        RETURN
      END IF
C
C     --- 
C
      IF(IZSAVE.EQ.22) THEN
        ITYP(1)=111
        ITYP(2)=112
C     ---
        IF (MMPTYP.EQ.112) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : TI/NR             CONFIGURATION : S1P6D03      *
C   VALENCE : (4S/3P/3D)        BASIS SET     : ( 6/ 4/ 5)   *
C          EFFECTIVE CORE CHARGE : 10.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =     -37.874494 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2035 A.U.
C  ***     ORBITAL ENERGY( 3P ) =   -1.6372 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.2733 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE - 10.0D+00
          SHINFO='4S/3P/3D'
          NOSH= 3
          CDFCTR(1)='51'
          CDFCTR(2)='4'
          CDFCTR(3)='41'
             KNG(1)= 6
            KTYP(1)= 1
             KNG(2)= 4
            KTYP(2)= 2
             KNG(3)= 5
            KTYP(3)= 3
          EX( 1)= 0.17726900D+03
           C( 1)= 0.23688006D-01
          EX( 2)= 0.68507700D+01
           C( 2)=-0.13141155D+00
          EX( 3)= 0.98468400D+00
           C( 3)= 0.24001015D+00
          EX( 4)= 0.43576600D+00
           C( 4)= 0.22170706D+00
          EX( 5)= 0.69354400D-01
           C( 5)=-0.65987025D+00
          EX( 6)= 0.27157400D-01
           C( 6)=-0.47812686D+00
          EX( 7)= 0.73756600D+02
           C( 7)=-0.78827056D-01
          EX( 8)= 0.15084800D+02
           C( 8)=-0.27281449D+00
          EX( 9)= 0.13703600D+01
           C( 9)= 0.67515234D+00
          EX(10)= 0.43543900D+00
           C(10)= 0.41104591D+00
          EX(11)= 0.16335100D+02
           C(11)= 0.40783334D-01
          EX(12)= 0.39936800D+01
           C(12)= 0.19585838D+00
          EX(13)= 0.12130100D+01
           C(13)= 0.42237882D+00
          EX(14)= 0.36583900D+00
           C(14)= 0.46467539D+00
          EX(15)= 0.10309900D+00
           C(15)= 0.26335561D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE - 10.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.57574500D+00
          ALPN0( 1)= 0.23416770D+03
            AN0( 2)= 0.27313440D+00
          ALPN0( 2)= 0.72172170D+01
            AN0( 3)= 0.56627800D-01
          ALPN0( 3)= 0.13947940D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.62199720D+02
          ALPN1( 1)= 0.16164383D+04
            AN1( 2)= 0.50402070D+01
          ALPN1( 2)= 0.39984568D+02
            AN1( 3)= 0.67563600D-01
          ALPN1( 3)= 0.87459914D+00
C         --- CORE-SHELL INFORMATION 
          NOSH     =  4
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.36622137D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.42476652D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.54000200D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.35220960D+02
C         --- EXPONENTS AND COEFFICIENTS 
          EX(  1) = 0.64576435D+04
           C(  1) = 0.20873137D-01
          EX(  2) = 0.88877241D+03
           C(  2) = 0.14668096D+00
          EX(  3) = 0.20261383D+03
           C(  3) = 0.47163443D+00
          EX(  4) = 0.57436783D+02
           C(  4) = 0.48782457D+00
          EX(  5) = 0.92920711D+02
           C(  5) =-0.34269334D-02
          EX(  6) = 0.98951715D+01
           C(  6) = 0.20979652D-01
          EX(  7) = 0.39653414D+01
           C(  7) = 0.13941371D-01
          EX(  8) = 0.69443220D+01
           C(  8) = 0.46671784D-02
          EX(  9) = 0.10614388D+01
           C(  9) =-0.15738619D-01
          EX( 10) = 0.41303995D+00
           C( 10) =-0.78457161D-02
          EX( 11) = 0.64576435D+04
           C( 11) = 0.67533338D-02
          EX( 12) = 0.88877241D+03
           C( 12) = 0.47457433D-01
          EX( 13) = 0.20261383D+03
           C( 13) = 0.15259349D+00
          EX( 14) = 0.57436783D+02
           C( 14) = 0.15783168D+00
          EX( 15) = 0.92920711D+02
           C( 15) = 0.10937341D+00
          EX( 16) = 0.98951715D+01
           C( 16) =-0.66958292D+00
          EX( 17) = 0.39653414D+01
           C( 17) =-0.44495035D+00
          EX( 18) = 0.69443220D+01
           C( 18) = 0.23930231D-02
          EX( 19) = 0.10614388D+01
           C( 19) =-0.80697322D-02
          EX( 20) = 0.41303995D+00
           C( 20) =-0.40227689D-02
          EX( 21) = 0.64576435D+04
           C( 21) = 0.27080087D-02
          EX( 22) = 0.88877241D+03
           C( 22) = 0.19029882D-01
          EX( 23) = 0.20261383D+03
           C( 23) = 0.61188224D-01
          EX( 24) = 0.57436783D+02
           C( 24) = 0.63288676D-01
          EX( 25) = 0.92920711D+02
           C( 25) = 0.44816272D-01
          EX( 26) = 0.98951715D+01
           C( 26) =-0.27436477D+00
          EX( 27) = 0.39653414D+01
           C( 27) =-0.18232051D+00
          EX( 28) = 0.69443220D+01
           C( 28) =-0.24049575D+00
          EX( 29) = 0.10614388D+01
           C( 29) = 0.81099771D+00
          EX( 30) = 0.41303995D+00
           C( 30) = 0.40428310D+00
          EX( 31) = 0.25033874D+03
           C( 31) = 0.34491079D-01
          EX( 32) = 0.54570961D+02
           C( 32) = 0.21646601D+00
          EX( 33) = 0.16382891D+02
           C( 33) = 0.52839366D+00
          EX( 34) = 0.55052493D+01
           C( 34) = 0.39647523D+00
          EX( 35) = 0.22886835D+01
           C( 35) = 0.76362830D-02
          EX( 36) = 0.84508866D+00
           C( 36) = 0.11905167D-01
          EX( 37) = 0.30890917D+00
           C( 37) = 0.42113780D-02
          RETURN
          END IF
C     ---
        ELSE IF (MMPTYP.EQ.111) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : TI/NR             CONFIGURATION : S1D03        *
C   VALENCE : (4S/3D)           BASIS SET     : ( 5/ 5)      *
C          EFFECTIVE CORE CHARGE :  4.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =      -3.397219 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2049 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.2726 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE -  4.0D+00
          SHINFO='4S/3D'
          NOSH= 2
          CDFCTR(1)='41'
          CDFCTR(2)='41'
             KNG(1)= 5
            KTYP(1)= 1
             KNG(2)= 5
            KTYP(2)= 3
          EX( 1)= 0.45085949D+01
           C( 1)=-0.14208506D+00
          EX( 2)= 0.98417650D+00
           C( 2)= 0.31475159D+00
          EX( 3)= 0.43576580D+00
           C( 3)= 0.11550327D+00
          EX( 4)= 0.69354690D-01
           C( 4)=-0.64239483D+00
          EX( 5)= 0.27157550D-01
           C( 5)=-0.47733749D+00
          EX( 6)= 0.16335107D+02
           C( 6)= 0.52383677D-01
          EX( 7)= 0.39936760D+01
           C( 7)= 0.20328426D+00
          EX( 8)= 0.12130130D+01
           C( 8)= 0.42991479D+00
          EX( 9)= 0.36583900D+00
           C( 9)= 0.45895574D+00
          EX(10)= 0.10309860D+00
           C(10)= 0.25001903D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE -  4.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.17918121D+01
          ALPN0( 1)= 0.13859179D+03
            AN0( 2)= 0.13544640D+01
          ALPN0( 2)= 0.45380403D+01
            AN0( 3)= 0.60860992D+00
          ALPN0( 3)= 0.80969628D+00
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.91143702D+01
          ALPN1( 1)= 0.43131586D+04
            AN1( 2)= 0.30867489D+02
          ALPN1( 2)= 0.47107419D+02
            AN1( 3)= 0.25577166D+00
          ALPN1( 3)= 0.73615494D+00
C         --- CORE-SHELL INFORMATION
          NOSH     =  5
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.23804389D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.27609824D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.35100130D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.22893624D+02
          KTYP( 5) =  2
           KNG( 5) =  7
          BPAR( 5) = 0.21279648D+01
C         --- EXPONENTS AND COEFFICIENTS
          EX(  1) = 0.64576435D+04
           C(  1) = 0.20873137D-01
          EX(  2) = 0.88877241D+03
           C(  2) = 0.14668096D+00
          EX(  3) = 0.20261383D+03
           C(  3) = 0.47163443D+00
          EX(  4) = 0.57436783D+02
           C(  4) = 0.48782457D+00
          EX(  5) = 0.92920711D+02
           C(  5) =-0.34269334D-02
          EX(  6) = 0.98951715D+01
           C(  6) = 0.20979652D-01
          EX(  7) = 0.39653414D+01
           C(  7) = 0.13941371D-01
          EX(  8) = 0.69443220D+01
           C(  8) = 0.46671784D-02
          EX(  9) = 0.10614388D+01
           C(  9) =-0.15738619D-01
          EX( 10) = 0.41303995D+00
           C( 10) =-0.78457161D-02
          EX( 11) = 0.64576435D+04
           C( 11) = 0.67533338D-02
          EX( 12) = 0.88877241D+03
           C( 12) = 0.47457433D-01
          EX( 13) = 0.20261383D+03
           C( 13) = 0.15259349D+00
          EX( 14) = 0.57436783D+02
           C( 14) = 0.15783168D+00
          EX( 15) = 0.92920711D+02
           C( 15) = 0.10937341D+00
          EX( 16) = 0.98951715D+01
           C( 16) =-0.66958292D+00
          EX( 17) = 0.39653414D+01
           C( 17) =-0.44495035D+00
          EX( 18) = 0.69443220D+01
           C( 18) = 0.23930231D-02
          EX( 19) = 0.10614388D+01
           C( 19) =-0.80697322D-02
          EX( 20) = 0.41303995D+00
           C( 20) =-0.40227689D-02
          EX( 21) = 0.64576435D+04
           C( 21) = 0.27080087D-02
          EX( 22) = 0.88877241D+03
           C( 22) = 0.19029882D-01
          EX( 23) = 0.20261383D+03
           C( 23) = 0.61188224D-01
          EX( 24) = 0.57436783D+02
           C( 24) = 0.63288676D-01
          EX( 25) = 0.92920711D+02
           C( 25) = 0.44816272D-01
          EX( 26) = 0.98951715D+01
           C( 26) =-0.27436477D+00
          EX( 27) = 0.39653414D+01
           C( 27) =-0.18232051D+00
          EX( 28) = 0.69443220D+01
           C( 28) =-0.24049575D+00
          EX( 29) = 0.10614388D+01
           C( 29) = 0.81099771D+00
          EX( 30) = 0.41303995D+00
           C( 30) = 0.40428310D+00
          EX( 31) = 0.25033874D+03
           C( 31) = 0.34491079D-01
          EX( 32) = 0.54570961D+02
           C( 32) = 0.21646601D+00
          EX( 33) = 0.16382891D+02
           C( 33) = 0.52839366D+00
          EX( 34) = 0.55052493D+01
           C( 34) = 0.39647523D+00
          EX( 35) = 0.22886835D+01
           C( 35) = 0.76362830D-02
          EX( 36) = 0.84508866D+00
           C( 36) = 0.11905167D-01
          EX( 37) = 0.30890917D+00
           C( 37) = 0.42113780D-02
          EX( 38) = 0.25033874D+03
           C( 38) = 0.12821044D-01
          EX( 39) = 0.54570961D+02
           C( 39) = 0.80464870D-01
          EX( 40) = 0.16382891D+02
           C( 40) = 0.19641479D+00
          EX( 41) = 0.55052493D+01
           C( 41) = 0.14737800D+00
          EX( 42) = 0.22886835D+01
           C( 42) =-0.38188627D+00
          EX( 43) = 0.84508866D+00
           C( 43) =-0.59537078D+00
          EX( 44) = 0.30890917D+00
           C( 44) =-0.21060867D+00
          RETURN
          END IF
        END IF
        IERR=2
        RETURN
      END IF

C
C     --- 
C
      IF(IZSAVE.EQ.23) THEN
        ITYP(1)=111
        ITYP(2)=112
C     --- 
        IF (MMPTYP.EQ.112) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE :  V/NR             CONFIGURATION : S1P6D04      *
C   VALENCE : (4S/3P/3D)        BASIS SET     : ( 6/ 4/ 5)   *
C          EFFECTIVE CORE CHARGE : 11.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =     -47.896949 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2111 A.U.
C  ***     ORBITAL ENERGY( 3P ) =   -1.8474 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.3208 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE - 11.0D+00
          SHINFO='4S/3P/3D'
          NOSH= 3
          CDFCTR(1)='51'
          CDFCTR(2)='4'
          CDFCTR(3)='41'
             KNG(1)= 6
            KTYP(1)= 1
             KNG(2)= 4
            KTYP(2)= 2
             KNG(3)= 5
            KTYP(3)= 3
          EX( 1)= 0.19437700D+03
           C( 1)= 0.20665148D-01
          EX( 2)= 0.76220200D+01
           C( 2)=-0.12149183D+00
          EX( 3)= 0.11236300D+01
           C( 3)= 0.20912810D+00
          EX( 4)= 0.50302800D+00
           C( 4)= 0.23273076D+00
          EX( 5)= 0.76408000D-01
           C( 5)=-0.63913912D+00
          EX( 6)= 0.29331000D-01
           C( 6)=-0.49695213D+00
          EX( 7)= 0.82378800D+02
           C( 7)= 0.78770804D-01
          EX( 8)= 0.16943100D+02
           C( 8)= 0.27636440D+00
          EX( 9)= 0.15572700D+01
           C( 9)=-0.67479104D+00
          EX(10)= 0.49394800D+00
           C(10)=-0.41183657D+00
          EX(11)= 0.19468140D+02
           C(11)= 0.41109926D-01
          EX(12)= 0.48264850D+01
           C(12)= 0.19675542D+00
          EX(13)= 0.14858080D+01
           C(13)= 0.42690962D+00
          EX(14)= 0.45410600D+00
           C(14)= 0.46362683D+00
          EX(15)= 0.12784100D+00
           C(15)= 0.25225621D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE - 11.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.50530352D+00
          ALPN0( 1)= 0.21821664D+03
            AN0( 2)= 0.24660505D+00
          ALPN0( 2)= 0.82040610D+01
            AN0( 3)= 0.51529741D-01
          ALPN0( 3)= 0.13841420D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.49082878D+02
          ALPN1( 1)= 0.17785833D+04
            AN1( 2)= 0.45476563D+01
          ALPN1( 2)= 0.42595623D+02
            AN1( 3)= 0.53933967D-01
          ALPN1( 3)= 0.93982466D+00
C         --- CORE-SHELL INFORMATION 
          NOSH     =  4
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.40264621D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.47344512D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.59858740D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.39648052D+02
C         --- EXPONENTS AND COEFFICIENTS 
          EX(  1) = 0.70924089D+04
           C(  1) = 0.20765467D-01
          EX(  2) = 0.97625225D+03
           C(  2) = 0.14601941D+00
          EX(  3) = 0.22263052D+03
           C(  3) = 0.47045182D+00
          EX(  4) = 0.63161949D+02
           C(  4) = 0.48927444D+00
          EX(  5) = 0.10209467D+03
           C(  5) =-0.34980900D-02
          EX(  6) = 0.10931336D+02
           C(  6) = 0.21380649D-01
          EX(  7) = 0.43969718D+01
           C(  7) = 0.14010965D-01
          EX(  8) = 0.78034324D+01
           C(  8) = 0.46469734D-02
          EX(  9) = 0.11985445D+01
           C(  9) =-0.15620622D-01
          EX( 10) = 0.46335002D+00
           C( 10) =-0.77574578D-02
          EX( 11) = 0.70924089D+04
           C( 11) = 0.67720082D-02
          EX( 12) = 0.97625225D+03
           C( 12) = 0.47619667D-01
          EX( 13) = 0.22263052D+03
           C( 13) = 0.15342316D+00
          EX( 14) = 0.63161949D+02
           C( 14) = 0.15956157D+00
          EX( 15) = 0.10209467D+03
           C( 15) = 0.11022571D+00
          EX( 16) = 0.10931336D+02
           C( 16) =-0.67370973D+00
          EX( 17) = 0.43969718D+01
           C( 17) =-0.44148911D+00
          EX( 18) = 0.78034324D+01
           C( 18) = 0.24118796D-02
          EX( 19) = 0.11985445D+01
           C( 19) =-0.81074400D-02
          EX( 20) = 0.46335002D+00
           C( 20) =-0.40262881D-02
          EX( 21) = 0.70924089D+04
           C( 21) = 0.27368076D-02
          EX( 22) = 0.97625225D+03
           C( 22) = 0.19244788D-01
          EX( 23) = 0.22263052D+03
           C( 23) = 0.62003716D-01
          EX( 24) = 0.63161949D+02
           C( 24) = 0.64484464D-01
          EX( 25) = 0.10209467D+03
           C( 25) = 0.45820869D-01
          EX( 26) = 0.10931336D+02
           C( 26) =-0.28006138D+00
          EX( 27) = 0.43969718D+01
           C( 27) =-0.18352718D+00
          EX( 28) = 0.78034324D+01
           C( 28) =-0.24231757D+00
          EX( 29) = 0.11985445D+01
           C( 29) = 0.81454113D+00
          EX( 30) = 0.46335002D+00
           C( 30) = 0.40451454D+00
          EX( 31) = 0.27840361D+03
           C( 31) = 0.34117083D-01
          EX( 32) = 0.60775517D+02
           C( 32) = 0.21535027D+00
          EX( 33) = 0.18299220D+02
           C( 33) = 0.52863684D+00
          EX( 34) = 0.61742527D+01
           C( 34) = 0.39610772D+00
          EX( 35) = 0.26069439D+01
           C( 35) = 0.78062335D-02
          EX( 36) = 0.96229933D+00
           C( 36) = 0.12043270D-01
          EX( 37) = 0.35062616D+00
           C( 37) = 0.42495969D-02
          RETURN
          END IF
C     --- 
        ELSE IF (MMPTYP.EQ.111) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE :  V/NR             CONFIGURATION : S1D04        *
C   VALENCE : (4S/3D)           BASIS SET     : ( 5/ 5)      *
C          EFFECTIVE CORE CHARGE :  5.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =      -6.151455 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2125 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.3207 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE -  5.0D+00
          SHINFO='4S/3D'
          NOSH= 2
          CDFCTR(1)='41'
          CDFCTR(2)='41'
             KNG(1)= 5
            KTYP(1)= 1
             KNG(2)= 5
            KTYP(2)= 3
          EX( 1)= 0.49829274D+01
           C( 1)= 0.12189593D+00
          EX( 2)= 0.11239947D+01
           C( 2)=-0.26688782D+00
          EX( 3)= 0.50331158D+00
           C( 3)=-0.14330535D+00
          EX( 4)= 0.76392108D-01
           C( 4)= 0.62954615D+00
          EX( 5)= 0.29325261D-01
           C( 5)= 0.49256852D+00
          EX( 6)= 0.19468139D+02
           C( 6)= 0.50829970D-01
          EX( 7)= 0.48264845D+01
           C( 7)= 0.20253039D+00
          EX( 8)= 0.14858081D+01
           C( 8)= 0.42913483D+00
          EX( 9)= 0.45410595D+00
           C( 9)= 0.46366037D+00
          EX(10)= 0.12784122D+00
           C(10)= 0.23977541D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE -  5.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.12134928D+01
          ALPN0( 1)= 0.16684039D+03
            AN0( 2)= 0.10666553D+01
          ALPN0( 2)= 0.48572682D+01
            AN0( 3)= 0.45805882D+00
          ALPN0( 3)= 0.85019657D+00
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.61281774D+01
          ALPN1( 1)= 0.38305666D+04
            AN1( 2)= 0.19136510D+02
          ALPN1( 2)= 0.44707580D+02
            AN1( 3)= 0.19349685D+00
          ALPN1( 3)= 0.79328490D+00
C         --- CORE-SHELL INFORMATION
          NOSH     =  5
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.26172004D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.30773933D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.38908181D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.25771234D+02
          KTYP( 5) =  2
           KNG( 5) =  7
          BPAR( 5) = 0.23975848D+01
C         --- EXPONENTS AND COEFFICIENTS
          EX(  1) = 0.70924089D+04
           C(  1) = 0.20765467D-01
          EX(  2) = 0.97625225D+03
           C(  2) = 0.14601941D+00
          EX(  3) = 0.22263052D+03
           C(  3) = 0.47045182D+00
          EX(  4) = 0.63161949D+02
           C(  4) = 0.48927444D+00
          EX(  5) = 0.10209467D+03
           C(  5) =-0.34980900D-02
          EX(  6) = 0.10931336D+02
           C(  6) = 0.21380649D-01
          EX(  7) = 0.43969718D+01
           C(  7) = 0.14010965D-01
          EX(  8) = 0.78034324D+01
           C(  8) = 0.46469734D-02
          EX(  9) = 0.11985445D+01
           C(  9) =-0.15620622D-01
          EX( 10) = 0.46335002D+00
           C( 10) =-0.77574578D-02
          EX( 11) = 0.70924089D+04
           C( 11) = 0.67720082D-02
          EX( 12) = 0.97625225D+03
           C( 12) = 0.47619667D-01
          EX( 13) = 0.22263052D+03
           C( 13) = 0.15342316D+00
          EX( 14) = 0.63161949D+02
           C( 14) = 0.15956157D+00
          EX( 15) = 0.10209467D+03
           C( 15) = 0.11022571D+00
          EX( 16) = 0.10931336D+02
           C( 16) =-0.67370973D+00
          EX( 17) = 0.43969718D+01
           C( 17) =-0.44148911D+00
          EX( 18) = 0.78034324D+01
           C( 18) = 0.24118796D-02
          EX( 19) = 0.11985445D+01
           C( 19) =-0.81074400D-02
          EX( 20) = 0.46335002D+00
           C( 20) =-0.40262881D-02
          EX( 21) = 0.70924089D+04
           C( 21) = 0.27368076D-02
          EX( 22) = 0.97625225D+03
           C( 22) = 0.19244788D-01
          EX( 23) = 0.22263052D+03
           C( 23) = 0.62003716D-01
          EX( 24) = 0.63161949D+02
           C( 24) = 0.64484464D-01
          EX( 25) = 0.10209467D+03
           C( 25) = 0.45820869D-01
          EX( 26) = 0.10931336D+02
           C( 26) =-0.28006138D+00
          EX( 27) = 0.43969718D+01
           C( 27) =-0.18352718D+00
          EX( 28) = 0.78034324D+01
           C( 28) =-0.24231757D+00
          EX( 29) = 0.11985445D+01
           C( 29) = 0.81454113D+00
          EX( 30) = 0.46335002D+00
           C( 30) = 0.40451454D+00
          EX( 31) = 0.27840361D+03
           C( 31) = 0.34117083D-01
          EX( 32) = 0.60775517D+02
           C( 32) = 0.21535027D+00
          EX( 33) = 0.18299220D+02
           C( 33) = 0.52863684D+00
          EX( 34) = 0.61742527D+01
           C( 34) = 0.39610772D+00
          EX( 35) = 0.26069439D+01
           C( 35) = 0.78062335D-02
          EX( 36) = 0.96229933D+00
           C( 36) = 0.12043270D-01
          EX( 37) = 0.35062616D+00
           C( 37) = 0.42495969D-02
          EX( 38) = 0.27840361D+03
           C( 38) = 0.12925157D-01
          EX( 39) = 0.60775517D+02
           C( 39) = 0.81584816D-01
          EX( 40) = 0.18299220D+02
           C( 40) = 0.20027251D+00
          EX( 41) = 0.61742527D+01
           C( 41) = 0.15006425D+00
          EX( 42) = 0.26069439D+01
           C( 42) =-0.38572771D+00
          EX( 43) = 0.96229933D+00
           C( 43) =-0.59509148D+00
          EX( 44) = 0.35062616D+00
           C( 44) =-0.20998441D+00
          RETURN
          END IF
        END IF
        IERR=2
        RETURN
      END IF
C
C     --- 
C
      IF(IZSAVE.EQ.24) THEN
        ITYP(1)=111
        ITYP(2)=112
C     --- 
        IF (MMPTYP.EQ.112) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : CR/NR             CONFIGURATION : S1P6D05      *
C   VALENCE : (4S/3P/3D)        BASIS SET     : ( 6/ 4/ 5)   *
C          EFFECTIVE CORE CHARGE : 12.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =     -59.805165 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2191 A.U.
C  ***     ORBITAL ENERGY( 3P ) =   -2.0513 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.3735 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE - 12.0D+00
          SHINFO='4S/3P/3D'
          NOSH= 3
          CDFCTR(1)='51'
          CDFCTR(2)='4'
          CDFCTR(3)='41'
             KNG(1)= 6
            KTYP(1)= 1
             KNG(2)= 4
            KTYP(2)= 2
             KNG(3)= 5
            KTYP(3)= 3
          EX( 1)= 0.21230800D+03
           C( 1)= 0.21280779D-01
          EX( 2)= 0.84313400D+01
           C( 2)=-0.12325908D+00
          EX( 3)= 0.12705800D+01
           C( 3)= 0.20703360D+00
          EX( 4)= 0.57229900D+00
           C( 4)= 0.22512004D+00
          EX( 5)= 0.82770600D-01
           C( 5)=-0.63350812D+00
          EX( 6)= 0.31206000D-01
           C( 6)=-0.49732019D+00
          EX( 7)= 0.91501300D+02
           C( 7)=-0.77962345D-01
          EX( 8)= 0.18907900D+02
           C( 8)=-0.28229249D+00
          EX( 9)= 0.17529700D+01
           C( 9)= 0.67490947D+00
          EX(10)= 0.55446900D+00
           C(10)= 0.41185698D+00
          EX(11)= 0.22768030D+02
           C(11)= 0.43716124D-01
          EX(12)= 0.57090860D+01
           C(12)= 0.20202898D+00
          EX(13)= 0.17763140D+01
           C(13)= 0.42828019D+00
          EX(14)= 0.54898500D+00
           C(14)= 0.46176750D+00
          EX(15)= 0.15499000D+00
           C(15)= 0.24227364D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE - 12.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.45208707D+00
          ALPN0( 1)= 0.21287933D+03
            AN0( 2)= 0.22228380D+00
          ALPN0( 2)= 0.83601643D+01
            AN0( 3)= 0.42172196D-01
          ALPN0( 3)= 0.14459102D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.47844576D+02
          ALPN1( 1)= 0.17797317D+04
            AN1( 2)= 0.39322808D+01
          ALPN1( 2)= 0.40928819D+02
            AN1( 3)= 0.45469045D-01
          ALPN1( 3)= 0.97624951D+00
C         --- CORE-SHELL INFORMATION 
          NOSH     =  4
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.44077200D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.52419200D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.65703200D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.44279800D+02
C         --- EXPONENTS AND COEFFICIENTS 
          EX(  1) = 0.77464664D+04
           C(  1) = 0.20702128D-01
          EX(  2) = 0.10663283D+04
           C(  2) = 0.14564839D+00
          EX(  3) = 0.24322171D+03
           C(  3) = 0.46990818D+00
          EX(  4) = 0.69048188D+02
           C(  4) = 0.48996764D+00
          EX(  5) = 0.11166806D+03
           C(  5) =-0.35311623D-02
          EX(  6) = 0.12025026D+02
           C(  6) = 0.21536801D-01
          EX(  7) = 0.48502037D+01
           C(  7) = 0.13953084D-01
          EX(  8) = 0.87101879D+01
           C(  8) = 0.45630735D-02
          EX(  9) = 0.13406471D+01
           C(  9) =-0.15325397D-01
          EX( 10) = 0.51427972D+00
           C( 10) =-0.75918476D-02
          EX( 11) = 0.77464664D+04
           C( 11) = 0.68015646D-02
          EX( 12) = 0.10663283D+04
           C( 12) = 0.47851939D-01
          EX( 13) = 0.24322171D+03
           C( 13) = 0.15438562D+00
          EX( 14) = 0.69048188D+02
           C( 14) = 0.16097604D+00
          EX( 15) = 0.11166806D+03
           C( 15) = 0.11103617D+00
          EX( 16) = 0.12025026D+02
           C( 16) =-0.67721723D+00
          EX( 17) = 0.48502037D+01
           C( 17) =-0.43874987D+00
          EX( 18) = 0.87101879D+01
           C( 18) = 0.24196521D-02
          EX( 19) = 0.13406471D+01
           C( 19) =-0.81265684D-02
          EX( 20) = 0.51427972D+00
           C( 20) =-0.40257142D-02
          EX( 21) = 0.77464664D+04
           C( 21) = 0.24855821D-02
          EX( 22) = 0.10663283D+04
           C( 22) = 0.17487141D-01
          EX( 23) = 0.24322171D+03
           C( 23) = 0.56419097D-01
          EX( 24) = 0.69048188D+02
           C( 24) = 0.58827518D-01
          EX( 25) = 0.11166806D+03
           C( 25) = 0.46723118D-01
          EX( 26) = 0.12025026D+02
           C( 26) =-0.28496750D+00
          EX( 27) = 0.48502037D+01
           C( 27) =-0.18462238D+00
          EX( 28) = 0.87101879D+01
           C( 28) =-0.24335348D+00
          EX( 29) = 0.13406471D+01
           C( 29) = 0.81731944D+00
          EX( 30) = 0.51427972D+00
           C( 30) = 0.40488116D+00
          EX( 31) = 0.30775507D+03
           C( 31) = 0.33813507D-01
          EX( 32) = 0.67267943D+02
           C( 32) = 0.21450669D+00
          EX( 33) = 0.20306878D+02
           C( 33) = 0.52899581D+00
          EX( 34) = 0.68751302D+01
           C( 34) = 0.39553452D+00
          EX( 35) = 0.29433334D+01
           C( 35) = 0.79334830D-02
          EX( 36) = 0.10852160D+01
           C( 36) = 0.12150987D-01
          EX( 37) = 0.39413403D+00
           C( 37) = 0.43008616D-02
          RETURN
          END IF
C     --- 
        ELSE IF (MMPTYP.EQ.111) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : CR/NR             CONFIGURATION : S1D05        *
C   VALENCE : (4S/3D)           BASIS SET     : ( 5/ 5)      *
C          EFFECTIVE CORE CHARGE :  6.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =     -10.187127 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2209 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.3732 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE -  6.0D+00
          SHINFO='4S/3D'
          NOSH= 2
          CDFCTR(1)='41'
          CDFCTR(2)='41'
             KNG(1)= 5
            KTYP(1)= 1
             KNG(2)= 5
            KTYP(2)= 3
          EX( 1)= 0.54817000D+01
           C( 1)= 0.12558257D+00
          EX( 2)= 0.12705800D+01
           C( 2)=-0.27218117D+00
          EX( 3)= 0.57229900D+00
           C( 3)=-0.12181370D+00
          EX( 4)= 0.82770600D-01
           C( 4)= 0.62412356D+00
          EX( 5)= 0.31206000D-01
           C( 5)= 0.49108510D+00
          EX( 6)= 0.22768030D+02
           C( 6)= 0.52942421D-01
          EX( 7)= 0.57090860D+01
           C( 7)= 0.21675044D+00
          EX( 8)= 0.17763140D+01
           C( 8)= 0.43222425D+00
          EX( 9)= 0.54898500D+00
           C( 9)= 0.45552392D+00
          EX(10)= 0.15499000D+00
           C(10)= 0.22751075D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE -  6.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.11679121D+01
          ALPN0( 1)= 0.16187196D+03
            AN0( 2)= 0.88377783D+00
          ALPN0( 2)= 0.48875912D+01
            AN0( 3)= 0.36744553D+00
          ALPN0( 3)= 0.91536508D+00
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.61132968D+01
          ALPN1( 1)= 0.43090584D+04
            AN1( 2)= 0.13488394D+02
          ALPN1( 2)= 0.43977052D+02
            AN1( 3)= 0.14474333D+00
          ALPN1( 3)= 0.86416434D+00
C         --- CORE-SHELL INFORMATION
          NOSH     =  5
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.26446320D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.31451520D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.39421920D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.26567880D+02
          KTYP( 5) =  2
           KNG( 5) =  7
          BPAR( 5) = 0.24610800D+01
C         --- EXPONENTS AND COEFFICIENTS
          EX(  1) = 0.77464664D+04
           C(  1) = 0.20702128D-01
          EX(  2) = 0.10663283D+04
           C(  2) = 0.14564839D+00
          EX(  3) = 0.24322171D+03
           C(  3) = 0.46990818D+00
          EX(  4) = 0.69048188D+02
           C(  4) = 0.48996764D+00
          EX(  5) = 0.11166806D+03
           C(  5) =-0.35311623D-02
          EX(  6) = 0.12025026D+02
           C(  6) = 0.21536801D-01
          EX(  7) = 0.48502037D+01
           C(  7) = 0.13953084D-01
          EX(  8) = 0.87101879D+01
           C(  8) = 0.45630735D-02
          EX(  9) = 0.13406471D+01
           C(  9) =-0.15325397D-01
          EX( 10) = 0.51427972D+00
           C( 10) =-0.75918476D-02
          EX( 11) = 0.77464664D+04
           C( 11) = 0.68015646D-02
          EX( 12) = 0.10663283D+04
           C( 12) = 0.47851939D-01
          EX( 13) = 0.24322171D+03
           C( 13) = 0.15438562D+00
          EX( 14) = 0.69048188D+02
           C( 14) = 0.16097604D+00
          EX( 15) = 0.11166806D+03
           C( 15) = 0.11103617D+00
          EX( 16) = 0.12025026D+02
           C( 16) =-0.67721723D+00
          EX( 17) = 0.48502037D+01
           C( 17) =-0.43874987D+00
          EX( 18) = 0.87101879D+01
           C( 18) = 0.24196521D-02
          EX( 19) = 0.13406471D+01
           C( 19) =-0.81265684D-02
          EX( 20) = 0.51427972D+00
           C( 20) =-0.40257142D-02
          EX( 21) = 0.77464664D+04
           C( 21) = 0.24855821D-02
          EX( 22) = 0.10663283D+04
           C( 22) = 0.17487141D-01
          EX( 23) = 0.24322171D+03
           C( 23) = 0.56419097D-01
          EX( 24) = 0.69048188D+02
           C( 24) = 0.58827518D-01
          EX( 25) = 0.11166806D+03
           C( 25) = 0.46723118D-01
          EX( 26) = 0.12025026D+02
           C( 26) =-0.28496750D+00
          EX( 27) = 0.48502037D+01
           C( 27) =-0.18462238D+00
          EX( 28) = 0.87101879D+01
           C( 28) =-0.24335348D+00
          EX( 29) = 0.13406471D+01
           C( 29) = 0.81731944D+00
          EX( 30) = 0.51427972D+00
           C( 30) = 0.40488116D+00
          EX( 31) = 0.30775507D+03
           C( 31) = 0.33813507D-01
          EX( 32) = 0.67267943D+02
           C( 32) = 0.21450669D+00
          EX( 33) = 0.20306878D+02
           C( 33) = 0.52899581D+00
          EX( 34) = 0.68751302D+01
           C( 34) = 0.39553452D+00
          EX( 35) = 0.29433334D+01
           C( 35) = 0.79334830D-02
          EX( 36) = 0.10852160D+01
           C( 36) = 0.12150987D-01
          EX( 37) = 0.39413403D+00
           C( 37) = 0.43008616D-02
          EX( 38) = 0.30775507D+03
           C( 38) = 0.12944503D-01
          EX( 39) = 0.67267943D+02
           C( 39) = 0.82117552D-01
          EX( 40) = 0.20306878D+02
           C( 40) = 0.20251042D+00
          EX( 41) = 0.68751302D+01
           C( 41) = 0.15141871D+00
          EX( 42) = 0.29433334D+01
           C( 42) =-0.38831615D+00
          EX( 43) = 0.10852160D+01
           C( 43) =-0.59474812D+00
          EX( 44) = 0.39413403D+00
           C( 44) =-0.21051208D+00
          RETURN
          END IF
        END IF
        IERR=2
        RETURN
      END IF
C
C     --- 
C
      IF(IZSAVE.EQ.25) THEN
        ITYP(1)=111
        ITYP(2)=112
C     --- 
        IF (MMPTYP.EQ.112) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : MN/NR             CONFIGURATION : S1P6D06      *
C   VALENCE : (4S/3P/3D)        BASIS SET     : ( 6/ 4/ 5)   *
C          EFFECTIVE CORE CHARGE : 13.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =     -73.769791 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2243 A.U.
C  ***     ORBITAL ENERGY( 3P ) =   -2.3137 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.3885 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE - 13.0D+00
          SHINFO='4S/3P/3D'
          NOSH= 3
          CDFCTR(1)='51'
          CDFCTR(2)='4'
          CDFCTR(3)='41'
             KNG(1)= 6
            KTYP(1)= 1
             KNG(2)= 4
            KTYP(2)= 2
             KNG(3)= 5
            KTYP(3)= 3
          EX( 1)= 0.23104300D+03
           C( 1)= 0.22492583D-01
          EX( 2)= 0.92857000D+01
           C( 2)=-0.12548901D+00
          EX( 3)= 0.14246000D+01
           C( 3)= 0.20848191D+00
          EX( 4)= 0.63762200D+00
           C( 4)= 0.21224791D+00
          EX( 5)= 0.87677600D-01
           C( 5)=-0.62521917D+00
          EX( 6)= 0.32592700D-01
           C( 6)=-0.49898122D+00
          EX( 7)= 0.10104400D+03
           C( 7)=-0.78623700D-01
          EX( 8)= 0.20970300D+02
           C( 8)=-0.28775906D+00
          EX( 9)= 0.19619200D+01
           C( 9)= 0.67983389D+00
          EX(10)= 0.62096600D+00
           C(10)= 0.40617771D+00
          EX(11)= 0.25434940D+02
           C(11)= 0.46696403D-01
          EX(12)= 0.64202600D+01
           C(12)= 0.21252772D+00
          EX(13)= 0.20003200D+01
           C(13)= 0.43581653D+00
          EX(14)= 0.61466300D+00
           C(14)= 0.45535311D+00
          EX(15)= 0.17061880D+00
           C(15)= 0.23243323D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE - 13.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.39972439D+00
          ALPN0( 1)= 0.22139416D+03
            AN0( 2)= 0.19427548D+00
          ALPN0( 2)= 0.86945664D+01
            AN0( 3)= 0.37027034D-01
          ALPN0( 3)= 0.15037464D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.49758384D+02
          ALPN1( 1)= 0.17263381D+04
            AN1( 2)= 0.38143116D+01
          ALPN1( 2)= 0.42565952D+02
            AN1( 3)= 0.44104930D-01
          ALPN1( 3)= 0.10153000D+01
C         --- CORE-SHELL INFORMATION 
          NOSH     =  4
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.48071000D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.57819000D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.72437400D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.49233000D+02
C         --- EXPONENTS AND COEFFICIENTS 
          EX(  1) = 0.84295943D+04
           C(  1) = 0.20643575D-01
          EX(  2) = 0.11604113D+04
           C(  2) = 0.14530481D+00
          EX(  3) = 0.26473131D+03
           C(  3) = 0.46940339D+00
          EX(  4) = 0.75199675D+02
           C(  4) = 0.49061043D+00
          EX(  5) = 0.12167903D+03
           C(  5) =-0.35617674D-02
          EX(  6) = 0.13170954D+02
           C(  6) = 0.21681602D-01
          EX(  7) = 0.53247152D+01
           C(  7) = 0.13902679D-01
          EX(  8) = 0.96623148D+01
           C(  8) = 0.44930191D-02
          EX(  9) = 0.14911247D+01
           C(  9) =-0.15086294D-01
          EX( 10) = 0.56932452D+00
           C( 10) =-0.74131757D-02
          EX( 11) = 0.84295943D+04
           C( 11) = 0.68252801D-02
          EX( 12) = 0.11604113D+04
           C( 12) = 0.48041387D-01
          EX( 13) = 0.26473131D+03
           C( 13) = 0.15519646D+00
          EX( 14) = 0.75199675D+02
           C( 14) = 0.16220802D+00
          EX( 15) = 0.12167903D+03
           C( 15) = 0.11177435D+00
          EX( 16) = 0.13170954D+02
           C( 16) =-0.68040573D+00
          EX( 17) = 0.53247152D+01
           C( 17) =-0.43628982D+00
          EX( 18) = 0.96623148D+01
           C( 18) = 0.24240712D-02
          EX( 19) = 0.14911247D+01
           C( 19) =-0.81393485D-02
          EX( 20) = 0.56932452D+00
           C( 20) =-0.39995524D-02
          EX( 21) = 0.84295943D+04
           C( 21) = 0.25186293D-02
          EX( 22) = 0.11604113D+04
           C( 22) = 0.17727982D-01
          EX( 23) = 0.26473131D+03
           C( 23) = 0.57269785D-01
          EX( 24) = 0.75199675D+02
           C( 24) = 0.59857159D-01
          EX( 25) = 0.12167903D+03
           C( 25) = 0.47572839D-01
          EX( 26) = 0.13170954D+02
           C( 26) =-0.28959089D+00
          EX( 27) = 0.53247152D+01
           C( 27) =-0.18569150D+00
          EX( 28) = 0.96623148D+01
           C( 28) =-0.24464988D+00
          EX( 29) = 0.14911247D+01
           C( 29) = 0.82146542D+00
          EX( 30) = 0.56932452D+00
           C( 30) = 0.40365565D+00
          EX( 31) = 0.33868227D+03
           C( 31) = 0.33519195D-01
          EX( 32) = 0.74113143D+02
           C( 32) = 0.21363720D+00
          EX( 33) = 0.22426735D+02
           C( 33) = 0.52918407D+00
          EX( 34) = 0.76165560D+01
           C( 34) = 0.39522393D+00
          EX( 35) = 0.33000341D+01
           C( 35) = 0.80818558D-02
          EX( 36) = 0.12170842D+01
           C( 36) = 0.12290037D-01
          EX( 37) = 0.44136433D+00
           C( 37) = 0.43186893D-02
          RETURN
          END IF
C     --- 
        ELSE IF (MMPTYP.EQ.111) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : MN/NR             CONFIGURATION : S1D06        *
C   VALENCE : (4S/3D)           BASIS SET     : ( 5/ 5)      *
C          EFFECTIVE CORE CHARGE :  7.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =     -15.355890 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2218 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.3834 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE -  7.0D+00
          SHINFO='4S/3D'
          NOSH= 2
          CDFCTR(1)='41'
          CDFCTR(2)='41'
             KNG(1)= 5
            KTYP(1)= 1
             KNG(2)= 5
            KTYP(2)= 3
          EX( 1)= 0.60097100D+01
           C( 1)= 0.11913869D+00
          EX( 2)= 0.14246760D+01
           C( 2)=-0.24765258D+00
          EX( 3)= 0.63764400D+00
           C( 3)=-0.12898052D+00
          EX( 4)= 0.87678700D-01
           C( 4)= 0.58479035D+00
          EX( 5)= 0.32593270D-01
           C( 5)= 0.52643989D+00
          EX( 6)= 0.25434940D+02
           C( 6)= 0.54294566D-01
          EX( 7)= 0.64202600D+01
           C( 7)= 0.23376668D+00
          EX( 8)= 0.20003200D+01
           C( 8)= 0.43993242D+00
          EX( 9)= 0.61466300D+00
           C( 9)= 0.44248541D+00
          EX(10)= 0.17061880D+00
           C(10)= 0.22039744D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE -  7.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.10904200D+01
          ALPN0( 1)= 0.15453274D+03
            AN0( 2)= 0.82594000D+00
          ALPN0( 2)= 0.57948400D+01
            AN0( 3)= 0.33240000D+00
          ALPN0( 3)= 0.10975080D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.50021900D+01
          ALPN1( 1)= 0.30684750D+04
            AN1( 2)= 0.11531260D+02
          ALPN1( 2)= 0.56925300D+02
            AN1( 3)= 0.13852000D+00
          ALPN1( 3)= 0.98877000D+00
C         --- CORE-SHELL INFORMATION
          NOSH     =  5
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.28842564D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.34691400D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.43461600D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.29539740D+02
          KTYP( 5) =  2
           KNG( 5) =  7
          BPAR( 5) = 0.27586320D+01
C         --- EXPONENTS AND COEFFICIENTS
          EX(  1) = 0.84295943D+04
           C(  1) = 0.20643575D-01
          EX(  2) = 0.11604113D+04
           C(  2) = 0.14530481D+00
          EX(  3) = 0.26473131D+03
           C(  3) = 0.46940339D+00
          EX(  4) = 0.75199675D+02
           C(  4) = 0.49061043D+00
          EX(  5) = 0.12167903D+03
           C(  5) =-0.35617674D-02
          EX(  6) = 0.13170954D+02
           C(  6) = 0.21681602D-01
          EX(  7) = 0.53247152D+01
           C(  7) = 0.13902679D-01
          EX(  8) = 0.96623148D+01
           C(  8) = 0.44930191D-02
          EX(  9) = 0.14911247D+01
           C(  9) =-0.15086294D-01
          EX( 10) = 0.56932452D+00
           C( 10) =-0.74131757D-02
          EX( 11) = 0.84295943D+04
           C( 11) = 0.68252801D-02
          EX( 12) = 0.11604113D+04
           C( 12) = 0.48041387D-01
          EX( 13) = 0.26473131D+03
           C( 13) = 0.15519646D+00
          EX( 14) = 0.75199675D+02
           C( 14) = 0.16220802D+00
          EX( 15) = 0.12167903D+03
           C( 15) = 0.11177435D+00
          EX( 16) = 0.13170954D+02
           C( 16) =-0.68040573D+00
          EX( 17) = 0.53247152D+01
           C( 17) =-0.43628982D+00
          EX( 18) = 0.96623148D+01
           C( 18) = 0.24240712D-02
          EX( 19) = 0.14911247D+01
           C( 19) =-0.81393485D-02
          EX( 20) = 0.56932452D+00
           C( 20) =-0.39995524D-02
          EX( 21) = 0.84295943D+04
           C( 21) = 0.25186293D-02
          EX( 22) = 0.11604113D+04
           C( 22) = 0.17727982D-01
          EX( 23) = 0.26473131D+03
           C( 23) = 0.57269785D-01
          EX( 24) = 0.75199675D+02
           C( 24) = 0.59857159D-01
          EX( 25) = 0.12167903D+03
           C( 25) = 0.47572839D-01
          EX( 26) = 0.13170954D+02
           C( 26) =-0.28959089D+00
          EX( 27) = 0.53247152D+01
           C( 27) =-0.18569150D+00
          EX( 28) = 0.96623148D+01
           C( 28) =-0.24464988D+00
          EX( 29) = 0.14911247D+01
           C( 29) = 0.82146542D+00
          EX( 30) = 0.56932452D+00
           C( 30) = 0.40365565D+00
          EX( 31) = 0.33868227D+03
           C( 31) = 0.33519195D-01
          EX( 32) = 0.74113143D+02
           C( 32) = 0.21363720D+00
          EX( 33) = 0.22426735D+02
           C( 33) = 0.52918407D+00
          EX( 34) = 0.76165560D+01
           C( 34) = 0.39522393D+00
          EX( 35) = 0.33000341D+01
           C( 35) = 0.80818558D-02
          EX( 36) = 0.12170842D+01
           C( 36) = 0.12290037D-01
          EX( 37) = 0.44136433D+00
           C( 37) = 0.43186893D-02
          EX( 38) = 0.33868227D+03
           C( 38) = 0.13030356D-01
          EX( 39) = 0.74113143D+02
           C( 39) = 0.83049982D-01
          EX( 40) = 0.22426735D+02
           C( 40) = 0.20571665D+00
          EX( 41) = 0.76165560D+01
           C( 41) = 0.15364057D+00
          EX( 42) = 0.33000341D+01
           C( 42) =-0.39136692D+00
          EX( 43) = 0.12170842D+01
           C( 43) =-0.59514969D+00
          EX( 44) = 0.44136433D+00
           C( 44) =-0.20913416D+00
          RETURN
          END IF
        END IF
        IERR=2
        RETURN
      END IF
C
C     --- 
C
      IF(IZSAVE.EQ.26) THEN
        ITYP(1)=111
        ITYP(2)=112
C     --- 
        IF (MMPTYP.EQ.112) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : FE/NR             CONFIGURATION : S1P6D07      *
C   VALENCE : (4S/3P/3D)        BASIS SET     : ( 6/ 4/ 5)   *
C          EFFECTIVE CORE CHARGE : 14.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =     -89.089299 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2302 A.U.
C  ***     ORBITAL ENERGY( 3P ) =   -2.5446 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.4106 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE - 14.0D+00
          SHINFO='4S/3P/3D'
          NOSH= 3
          CDFCTR(1)='51'
          CDFCTR(2)='4'
          CDFCTR(3)='41'
             KNG(1)= 6
            KTYP(1)= 1
             KNG(2)= 4
            KTYP(2)= 2
             KNG(3)= 5
            KTYP(3)= 3
          EX( 1)= 0.25059100D+03
           C( 1)=-0.24403240D-01
          EX( 2)= 0.10180100D+02
           C( 2)= 0.12940977D+00
          EX( 3)= 0.15837800D+01
           C( 3)=-0.21515290D+00
          EX( 4)= 0.70128300D+00
           C( 4)=-0.19777739D+00
          EX( 5)= 0.91652500D-01
           C( 5)= 0.63289231D+00
          EX( 6)= 0.33635400D-01
           C( 6)= 0.48535633D+00
          EX( 7)= 0.11104800D+03
           C( 7)=-0.76964085D-01
          EX( 8)= 0.23135800D+02
           C( 8)=-0.28782077D+00
          EX( 9)= 0.21800000D+01
           C( 9)= 0.67538018D+00
          EX(10)= 0.68949000D+00
           C(10)= 0.41174883D+00
          EX(11)= 0.28446300D+02
           C(11)= 0.47269251D-01
          EX(12)= 0.72274140D+01
           C(12)= 0.21526007D+00
          EX(13)= 0.22585970D+01
           C(13)= 0.43867421D+00
          EX(14)= 0.69334520D+00
           C(14)= 0.45286570D+00
          EX(15)= 0.19093640D+00
           C(15)= 0.22952407D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE - 14.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.34787331D+00
          ALPN0( 1)= 0.22980177D+03
            AN0( 2)= 0.16397535D+00
          ALPN0( 2)= 0.94040469D+01
            AN0( 3)= 0.32626488D-01
          ALPN0( 3)= 0.16264560D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.53818685D+02
          ALPN1( 1)= 0.16243133D+04
            AN1( 2)= 0.37073944D+01
          ALPN1( 2)= 0.44271621D+02
            AN1( 3)= 0.41498300D-01
          ALPN1( 3)= 0.10981485D+01
C         --- CORE-SHELL INFORMATION 
          NOSH     =  4
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.52234800D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.63427200D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.79155400D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.54391800D+02
C         --- EXPONENTS AND COEFFICIENTS 
          EX(  1) = 0.91406349D+04
           C(  1) = 0.20592354D-01
          EX(  2) = 0.12583784D+04
           C(  2) = 0.14499939D+00
          EX(  3) = 0.28713911D+03
           C(  3) = 0.46894869D+00
          EX(  4) = 0.81612755D+02
           C(  4) = 0.49118058D+00
          EX(  5) = 0.13213054D+03
           C(  5) =-0.35888471D-02
          EX(  6) = 0.14369539D+02
           C(  6) = 0.21808374D-01
          EX(  7) = 0.58198069D+01
           C(  7) = 0.13854114D-01
          EX(  8) = 0.10660805D+02
           C(  8) = 0.44269891D-02
          EX(  9) = 0.16481489D+01
           C(  9) =-0.14854979D-01
          EX( 10) = 0.62655560D+00
           C( 10) =-0.72855442D-02
          EX( 11) = 0.91406349D+04
           C( 11) = 0.68474979D-02
          EX( 12) = 0.12583784D+04
           C( 12) = 0.48216099D-01
          EX( 13) = 0.28713911D+03
           C( 13) = 0.15593775D+00
          EX( 14) = 0.81612755D+02
           C( 14) = 0.16333043D+00
          EX( 15) = 0.13213054D+03
           C( 15) = 0.11244762D+00
          EX( 16) = 0.14369539D+02
           C( 16) =-0.68331129D+00
          EX( 17) = 0.58198069D+01
           C( 17) =-0.43408428D+00
          EX( 18) = 0.10660805D+02
           C( 18) = 0.24228043D-02
          EX( 19) = 0.16481489D+01
           C( 19) =-0.81298386D-02
          EX( 20) = 0.62655560D+00
           C( 20) =-0.39872353D-02
          EX( 21) = 0.91406349D+04
           C( 21) = 0.25464301D-02
          EX( 22) = 0.12583784D+04
           C( 22) = 0.17930481D-01
          EX( 23) = 0.28713911D+03
           C( 23) = 0.57989732D-01
          EX( 24) = 0.81612755D+02
           C( 24) = 0.60738906D-01
          EX( 25) = 0.13213054D+03
           C( 25) = 0.48306941D-01
          EX( 26) = 0.14369539D+02
           C( 26) =-0.29354715D+00
          EX( 27) = 0.58198069D+01
           C( 27) =-0.18648046D+00
          EX( 28) = 0.10660805D+02
           C( 28) =-0.24543046D+00
          EX( 29) = 0.16481489D+01
           C( 29) = 0.82355399D+00
          EX( 30) = 0.62655560D+00
           C( 30) = 0.40390760D+00
          EX( 31) = 0.37060107D+03
           C( 31) = 0.33320427D-01
          EX( 32) = 0.81177647D+02
           C( 32) = 0.21320545D+00
          EX( 33) = 0.24614657D+02
           C( 33) = 0.52976655D+00
          EX( 34) = 0.83800262D+01
           C( 34) = 0.39431111D+00
          EX( 35) = 0.36738873D+01
           C( 35) = 0.81375839D-02
          EX( 36) = 0.13543941D+01
           C( 36) = 0.12300792D-01
          EX( 37) = 0.49027667D+00
           C( 37) = 0.43142981D-02
          RETURN
          END IF
C     --- 
        ELSE IF (MMPTYP.EQ.111) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : FE/NR             CONFIGURATION : S1D07        *
C   VALENCE : (4S/3D)           BASIS SET     : ( 5/ 5)      *
C          EFFECTIVE CORE CHARGE :  8.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =     -21.658332 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2303 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.4066 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE -  8.0D+00
          SHINFO='4S/3D'
          NOSH= 2
          CDFCTR(1)='41'
          CDFCTR(2)='41'
             KNG(1)= 5
            KTYP(1)= 1
             KNG(2)= 5
            KTYP(2)= 3
          EX( 1)= 0.65574200D+01
           C( 1)=-0.11394219D+00
          EX( 2)= 0.15837100D+01
           C( 2)= 0.24325850D+00
          EX( 3)= 0.70122400D+00
           C( 3)= 0.11408959D+00
          EX( 4)= 0.91656000D-01
           C( 4)=-0.61487844D+00
          EX( 5)= 0.33636600D-01
           C( 5)=-0.49202189D+00
          EX( 6)= 0.28446300D+02
           C( 6)= 0.48851124D-01
          EX( 7)= 0.72274140D+01
           C( 7)= 0.22845490D+00
          EX( 8)= 0.22585970D+01
           C( 8)= 0.44108794D+00
          EX( 9)= 0.69334520D+00
           C( 9)= 0.44714676D+00
          EX(10)= 0.19093640D+00
           C(10)= 0.21972907D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE -  8.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.80113665D+00
          ALPN0( 1)= 0.16031301D+03
            AN0( 2)= 0.69226318D+00
          ALPN0( 2)= 0.60099023D+01
            AN0( 3)= 0.27320988D+00
          ALPN0( 3)= 0.11229458D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.42047962D+01
          ALPN1( 1)= 0.31862994D+04
            AN1( 2)= 0.80421700D+01
          ALPN1( 2)= 0.57431277D+02
            AN1( 3)= 0.11656391D+00
          ALPN1( 3)= 0.10267146D+01
C         --- CORE-SHELL INFORMATION
          NOSH     =  5
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.28729105D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.34884971D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.43535470D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.29915512D+02
          KTYP( 5) =  2
           KNG( 5) =  7
          BPAR( 5) = 0.28001737D+01
C         --- EXPONENTS AND COEFFICIENTS
          EX(  1) = 0.91406349D+04
           C(  1) = 0.20592354D-01
          EX(  2) = 0.12583784D+04
           C(  2) = 0.14499939D+00
          EX(  3) = 0.28713911D+03
           C(  3) = 0.46894869D+00
          EX(  4) = 0.81612755D+02
           C(  4) = 0.49118058D+00
          EX(  5) = 0.13213054D+03
           C(  5) =-0.35888471D-02
          EX(  6) = 0.14369539D+02
           C(  6) = 0.21808374D-01
          EX(  7) = 0.58198069D+01
           C(  7) = 0.13854114D-01
          EX(  8) = 0.10660805D+02
           C(  8) = 0.44269891D-02
          EX(  9) = 0.16481489D+01
           C(  9) =-0.14854979D-01
          EX( 10) = 0.62655560D+00
           C( 10) =-0.72855442D-02
          EX( 11) = 0.91406349D+04
           C( 11) = 0.68474979D-02
          EX( 12) = 0.12583784D+04
           C( 12) = 0.48216099D-01
          EX( 13) = 0.28713911D+03
           C( 13) = 0.15593775D+00
          EX( 14) = 0.81612755D+02
           C( 14) = 0.16333043D+00
          EX( 15) = 0.13213054D+03
           C( 15) = 0.11244762D+00
          EX( 16) = 0.14369539D+02
           C( 16) =-0.68331129D+00
          EX( 17) = 0.58198069D+01
           C( 17) =-0.43408428D+00
          EX( 18) = 0.10660805D+02
           C( 18) = 0.24228043D-02
          EX( 19) = 0.16481489D+01
           C( 19) =-0.81298387D-02
          EX( 20) = 0.62655560D+00
           C( 20) =-0.39872353D-02
          EX( 21) = 0.91406349D+04
           C( 21) = 0.25464301D-02
          EX( 22) = 0.12583784D+04
           C( 22) = 0.17930481D-01
          EX( 23) = 0.28713911D+03
           C( 23) = 0.57989732D-01
          EX( 24) = 0.81612755D+02
           C( 24) = 0.60738906D-01
          EX( 25) = 0.13213054D+03
           C( 25) = 0.48306941D-01
          EX( 26) = 0.14369539D+02
           C( 26) =-0.29354715D+00
          EX( 27) = 0.58198069D+01
           C( 27) =-0.18648046D+00
          EX( 28) = 0.10660805D+02
           C( 28) =-0.24543046D+00
          EX( 29) = 0.16481489D+01
           C( 29) = 0.82355399D+00
          EX( 30) = 0.62655560D+00
           C( 30) = 0.40390760D+00
          EX( 31) = 0.37060107D+03
           C( 31) = 0.33320427D-01
          EX( 32) = 0.81177647D+02
           C( 32) = 0.21320545D+00
          EX( 33) = 0.24614657D+02
           C( 33) = 0.52976655D+00
          EX( 34) = 0.83800262D+01
           C( 34) = 0.39431111D+00
          EX( 35) = 0.36738873D+01
           C( 35) = 0.81375839D-02
          EX( 36) = 0.13543941D+01
           C( 36) = 0.12300792D-01
          EX( 37) = 0.49027667D+00
           C( 37) = 0.43142981D-02
          EX( 38) = 0.37060107D+03
           C( 38) = 0.13114794D-01
          EX( 39) = 0.81177647D+02
           C( 39) = 0.83916858D-01
          EX( 40) = 0.24614657D+02
           C( 40) = 0.20851412D+00
          EX( 41) = 0.83800262D+01
           C( 41) = 0.15519937D+00
          EX( 42) = 0.36738873D+01
           C( 42) =-0.39375802D+00
          EX( 43) = 0.13543941D+01
           C( 43) =-0.59520561D+00
          EX( 44) = 0.49027667D+00
           C( 44) =-0.20875845D+00
          RETURN
          END IF
        END IF
        IERR=2
        RETURN
      END IF
C
C     --- 
C
      IF(IZSAVE.EQ.27) THEN
        ITYP(1)=111
        ITYP(2)=112
C     --- 
        IF (MMPTYP.EQ.112) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : CO/NR             CONFIGURATION : S1P6D08      *
C   VALENCE : (4S/3P/3D)        BASIS SET     : ( 6/ 4/ 5)   *
C          EFFECTIVE CORE CHARGE : 15.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =    -106.725289 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2341 A.U.
C  ***     ORBITAL ENERGY( 3P ) =   -2.7995 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.4319 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE - 15.0D+00
          SHINFO='4S/3P/3D'
          NOSH= 3
          CDFCTR(1)='51'
          CDFCTR(2)='4'
          CDFCTR(3)='41'
             KNG(1)= 6
            KTYP(1)= 1
             KNG(2)= 4
            KTYP(2)= 2
             KNG(3)= 5
            KTYP(3)= 3
          EX( 1)= 0.27094900D+03
           C( 1)=-0.24017557D-01
          EX( 2)= 0.11114600D+02
           C( 2)= 0.12628820D+00
          EX( 3)= 0.17487200D+01
           C( 3)=-0.20668338D+00
          EX( 4)= 0.76553300D+00
           C( 4)=-0.19174033D+00
          EX( 5)= 0.95484800D-01
           C( 5)= 0.62935783D+00
          EX( 6)= 0.34631200D-01
           C( 6)= 0.48480674D+00
          EX( 7)= 0.12150700D+03
           C( 7)=-0.76936023D-01
          EX( 8)= 0.25403600D+02
           C( 8)=-0.29207094D+00
          EX( 9)= 0.24080800D+01
           C( 9)= 0.67617786D+00
          EX(10)= 0.76084200D+00
           C(10)= 0.41075096D+00
          EX(11)= 0.31602500D+02
           C(11)= 0.49225315D-01
          EX(12)= 0.80750100D+01
           C(12)= 0.21991137D+00
          EX(13)= 0.25300900D+01
           C(13)= 0.43906805D+00
          EX(14)= 0.77616400D+00
           C(14)= 0.44989684D+00
          EX(15)= 0.21240400D+00
           C(15)= 0.22849593D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE - 15.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.29231430D+00
          ALPN0( 1)= 0.22773857D+03
            AN0( 2)= 0.13432253D+00
          ALPN0( 2)= 0.92072146D+01
            AN0( 3)= 0.25973615D-01
          ALPN0( 3)= 0.17533719D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.52606396D+02
          ALPN1( 1)= 0.15460258D+04
            AN1( 2)= 0.37208605D+01
          ALPN1( 2)= 0.44469016D+02
            AN1( 3)= 0.38868342D-01
          ALPN1( 3)= 0.11799944D+01
C         --- CORE-SHELL INFORMATION 
          NOSH     =  4
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.56570800D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.69268200D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.86040800D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.59780600D+02
C         --- EXPONENTS AND COEFFICIENTS 
          EX(  1) = 0.98827201D+04
           C(  1) = 0.20539883D-01
          EX(  2) = 0.13605479D+04
           C(  2) = 0.14469608D+00
          EX(  3) = 0.31049389D+03
           C(  3) = 0.46851249D+00
          EX(  4) = 0.88295722D+02
           C(  4) = 0.49174526D+00
          EX(  5) = 0.14302274D+03
           C(  5) =-0.36140973D-02
          EX(  6) = 0.15620842D+02
           C(  6) = 0.21925997D-01
          EX(  7) = 0.63360787D+01
           C(  7) = 0.13812648D-01
          EX(  8) = 0.11705322D+02
           C(  8) = 0.43678938D-02
          EX(  9) = 0.18123926D+01
           C(  9) =-0.14642146D-01
          EX( 10) = 0.68662763D+00
           C( 10) =-0.71884838D-02
          EX( 11) = 0.98827201D+04
           C( 11) = 0.68658868D-02
          EX( 12) = 0.13605479D+04
           C( 12) = 0.48367699D-01
          EX( 13) = 0.31049389D+03
           C( 13) = 0.15661013D+00
          EX( 14) = 0.88295722D+02
           C( 14) = 0.16437617D+00
          EX( 15) = 0.14302274D+03
           C( 15) = 0.11306413D+00
          EX( 16) = 0.15620842D+02
           C( 16) =-0.68593716D+00
          EX( 17) = 0.63360787D+01
           C( 17) =-0.43211758D+00
          EX( 18) = 0.11705322D+02
           C( 18) = 0.24197642D-02
          EX( 19) = 0.18123926D+01
           C( 19) =-0.81115846D-02
          EX( 20) = 0.68662763D+00
           C( 20) =-0.39823393D-02
          EX( 21) = 0.98827201D+04
           C( 21) = 0.25700903D-02
          EX( 22) = 0.13605479D+04
           C( 22) = 0.18105360D-01
          EX( 23) = 0.31049389D+03
           C( 23) = 0.58623478D-01
          EX( 24) = 0.88295722D+02
           C( 24) = 0.61530521D-01
          EX( 25) = 0.14302274D+03
           C( 25) = 0.48965458D-01
          EX( 26) = 0.15620842D+02
           C( 26) =-0.29706353D+00
          EX( 27) = 0.63360787D+01
           C( 27) =-0.18714013D+00
          EX( 28) = 0.11705322D+02
           C( 28) =-0.24599105D+00
          EX( 29) = 0.18123926D+01
           C( 29) = 0.82461638D+00
          EX( 30) = 0.68662763D+00
           C( 30) = 0.40484104D+00
          EX( 31) = 0.40450441D+03
           C( 31) = 0.33063705D-01
          EX( 32) = 0.88691173D+02
           C( 32) = 0.21240379D+00
          EX( 33) = 0.26948048D+02
           C( 33) = 0.52975436D+00
          EX( 34) = 0.91981776D+01
           C( 34) = 0.39426570D+00
          EX( 35) = 0.40684359D+01
           C( 35) = 0.82940261D-02
          EX( 36) = 0.14989033D+01
           C( 36) = 0.12495119D-01
          EX( 37) = 0.54149246D+00
           C( 37) = 0.43830737D-02
          RETURN
          END IF
C     --- 
        ELSE IF (MMPTYP.EQ.111) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : CO/NR             CONFIGURATION : S1D08        *
C   VALENCE : (4S/3D)           BASIS SET     : ( 5/ 5)      *
C          EFFECTIVE CORE CHARGE :  9.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =     -29.641268 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2332 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.4317 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE -  9.0D+00
          SHINFO='4S/3D'
          NOSH= 2
          CDFCTR(1)='41'
          CDFCTR(2)='41'
             KNG(1)= 5
            KTYP(1)= 1
             KNG(2)= 5
            KTYP(2)= 3
          EX( 1)= 0.71363430D+01
           C( 1)=-0.10344419D+00
          EX( 2)= 0.17487020D+01
           C( 2)= 0.22113104D+00
          EX( 3)= 0.76553000D+00
           C( 3)= 0.12251219D+00
          EX( 4)= 0.95484800D-01
           C( 4)=-0.60864664D+00
          EX( 5)= 0.34631200D-01
           C( 5)=-0.49667459D+00
          EX( 6)= 0.31602500D+02
           C( 6)= 0.47483226D-01
          EX( 7)= 0.80750100D+01
           C( 7)= 0.23556925D+00
          EX( 8)= 0.25300900D+01
           C( 8)= 0.44178181D+00
          EX( 9)= 0.77616400D+00
           C( 9)= 0.44332266D+00
          EX(10)= 0.21240400D+00
           C(10)= 0.21787647D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE -  9.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.72079200D+00
          ALPN0( 1)= 0.16726567D+03
            AN0( 2)= 0.61751338D+00
          ALPN0( 2)= 0.60555294D+01
            AN0( 3)= 0.22880292D+00
          ALPN0( 3)= 0.12486544D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.36009900D+01
          ALPN1( 1)= 0.28957100D+04
            AN1( 2)= 0.54986900D+01
          ALPN1( 2)= 0.58167000D+02
            AN1( 3)= 0.98734400D-01
          ALPN1( 3)= 0.99967100D+00
C         --- CORE-SHELL INFORMATION
          NOSH     =  5
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.31113940D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.38097510D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.47322440D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.32879330D+02
          KTYP( 5) =  2
           KNG( 5) =  7
          BPAR( 5) = 0.30794390D+01
C         --- EXPONENTS AND COEFFICIENTS
          EX(  1) = 0.98827201D+04
           C(  1) = 0.20539883D-01
          EX(  2) = 0.13605479D+04
           C(  2) = 0.14469608D+00
          EX(  3) = 0.31049389D+03
           C(  3) = 0.46851249D+00
          EX(  4) = 0.88295722D+02
           C(  4) = 0.49174526D+00
          EX(  5) = 0.14302274D+03
           C(  5) =-0.36140973D-02
          EX(  6) = 0.15620842D+02
           C(  6) = 0.21925997D-01
          EX(  7) = 0.63360787D+01
           C(  7) = 0.13812648D-01
          EX(  8) = 0.11705322D+02
           C(  8) = 0.43678938D-02
          EX(  9) = 0.18123926D+01
           C(  9) =-0.14642146D-01
          EX( 10) = 0.68662763D+00
           C( 10) =-0.71884838D-02
          EX( 11) = 0.98827201D+04
           C( 11) = 0.68658868D-02
          EX( 12) = 0.13605479D+04
           C( 12) = 0.48367699D-01
          EX( 13) = 0.31049389D+03
           C( 13) = 0.15661013D+00
          EX( 14) = 0.88295722D+02
           C( 14) = 0.16437617D+00
          EX( 15) = 0.14302274D+03
           C( 15) = 0.11306413D+00
          EX( 16) = 0.15620842D+02
           C( 16) =-0.68593716D+00
          EX( 17) = 0.63360787D+01
           C( 17) =-0.43211758D+00
          EX( 18) = 0.11705322D+02
           C( 18) = 0.24197642D-02
          EX( 19) = 0.18123926D+01
           C( 19) =-0.81115846D-02
          EX( 20) = 0.68662763D+00
           C( 20) =-0.39823393D-02
          EX( 21) = 0.98827201D+04
           C( 21) = 0.25700903D-02
          EX( 22) = 0.13605479D+04
           C( 22) = 0.18105360D-01
          EX( 23) = 0.31049389D+03
           C( 23) = 0.58623478D-01
          EX( 24) = 0.88295722D+02
           C( 24) = 0.61530521D-01
          EX( 25) = 0.14302274D+03
           C( 25) = 0.48965458D-01
          EX( 26) = 0.15620842D+02
           C( 26) =-0.29706353D+00
          EX( 27) = 0.63360787D+01
           C( 27) =-0.18714013D+00
          EX( 28) = 0.11705322D+02
           C( 28) =-0.24599105D+00
          EX( 29) = 0.18123926D+01
           C( 29) = 0.82461638D+00
          EX( 30) = 0.68662763D+00
           C( 30) = 0.40484104D+00
          EX( 31) = 0.40450441D+03
           C( 31) = 0.33063705D-01
          EX( 32) = 0.88691173D+02
           C( 32) = 0.21240379D+00
          EX( 33) = 0.26948048D+02
           C( 33) = 0.52975436D+00
          EX( 34) = 0.91981776D+01
           C( 34) = 0.39426570D+00
          EX( 35) = 0.40684359D+01
           C( 35) = 0.82940261D-02
          EX( 36) = 0.14989033D+01
           C( 36) = 0.12495119D-01
          EX( 37) = 0.54149246D+00
           C( 37) = 0.43830737D-02
          EX( 38) = 0.40450441D+03
           C( 38) = 0.13153501D-01
          EX( 39) = 0.88691173D+02
           C( 39) = 0.84499106D-01
          EX( 40) = 0.26948048D+02
           C( 40) = 0.21074844D+00
          EX( 41) = 0.91981776D+01
           C( 41) = 0.15684795D+00
          EX( 42) = 0.40684359D+01
           C( 42) =-0.39523414D+00
          EX( 43) = 0.14989033D+01
           C( 43) =-0.59542827D+00
          EX( 44) = 0.54149246D+00
           C( 44) =-0.20886604D+00
          RETURN
          END IF
        END IF
        IERR=2
        RETURN
      END IF
C
C     --- 
C
      IF(IZSAVE.EQ.28) THEN
        ITYP(1)=111
        ITYP(2)=112
C     --- 
        IF (MMPTYP.EQ.112) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : NI/NR             CONFIGURATION : S1P6D09      *
C   VALENCE : (4S/3P/3D)        BASIS SET     : ( 6/ 4/ 5)   *
C          EFFECTIVE CORE CHARGE : 16.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =    -127.321703 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2386 A.U.
C  ***     ORBITAL ENERGY( 3P ) =   -3.0610 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.4581 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE - 16.0D+00
          SHINFO='4S/3P/3D'
          NOSH= 3
          CDFCTR(1)='51'
          CDFCTR(2)='4'
          CDFCTR(3)='41'
             KNG(1)= 6
            KTYP(1)= 1
             KNG(2)= 4
            KTYP(2)= 2
             KNG(3)= 5
            KTYP(3)= 3
          EX( 1)= 0.29211600D+03
           C( 1)=-0.21128609D-01
          EX( 2)= 0.12089300D+02
           C( 2)= 0.11760589D+00
          EX( 3)= 0.19193600D+01
           C( 3)=-0.20069013D+00
          EX( 4)= 0.83054400D+00
           C( 4)=-0.17520653D+00
          EX( 5)= 0.99221700D-01
           C( 5)= 0.63048551D+00
          EX( 6)= 0.35595600D-01
           C( 6)= 0.47904275D+00
          EX( 7)= 0.13242000D+03
           C( 7)=-0.76800475D-01
          EX( 8)= 0.27773600D+02
           C( 8)=-0.29375750D+00
          EX( 9)= 0.26462600D+01
           C( 9)= 0.68161555D+00
          EX(10)= 0.83507900D+00
           C(10)= 0.40533788D+00
          EX(11)= 0.34880000D+02
           C(11)= 0.52691748D-01
          EX(12)= 0.89570000D+01
           C(12)= 0.22964374D+00
          EX(13)= 0.28127500D+01
           C(13)= 0.44093623D+00
          EX(14)= 0.86238900D+00
           C(14)= 0.44390110D+00
          EX(15)= 0.23475800D+00
           C(15)= 0.22337644D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE - 16.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.26243815D+00
          ALPN0( 1)= 0.23552752D+03
            AN0( 2)= 0.11683156D+00
          ALPN0( 2)= 0.88526220D+01
            AN0( 3)= 0.22641689D-01
          ALPN0( 3)= 0.18964450D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.51228569D+02
          ALPN1( 1)= 0.15275278D+04
            AN1( 2)= 0.36124943D+01
          ALPN1( 2)= 0.44372138D+02
            AN1( 3)= 0.36571183D-01
          ALPN1( 3)= 0.12374111D+01
C         --- CORE-SHELL INFORMATION 
          NOSH     =  4
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.61079200D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.75343600D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.93105400D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.65400400D+02
C         --- EXPONENTS AND COEFFICIENTS 
          EX(  1) = 0.13392119D+05
           C(  1) = 0.19507933D-01
          EX(  2) = 0.16721532D+04
           C(  2) = 0.13747259D+00
          EX(  3) = 0.36264482D+03
           C(  3) = 0.46098814D+00
          EX(  4) = 0.99876996D+02
           C(  4) = 0.50904402D+00
          EX(  5) = 0.16281510D+03
           C(  5) =-0.37745951D-02
          EX(  6) = 0.17533548D+02
           C(  6) = 0.24005283D-01
          EX(  7) = 0.71134955D+01
           C(  7) = 0.15643866D-01
          EX(  8) = 0.13198646D+02
           C(  8) = 0.47279613D-02
          EX(  9) = 0.20292497D+01
           C(  9) =-0.15932215D-01
          EX( 10) = 0.76788917D+00
           C( 10) =-0.78170344D-02
          EX( 11) = 0.13392119D+05
           C( 11) = 0.67211118D-02
          EX( 12) = 0.16721532D+04
           C( 12) = 0.47363740D-01
          EX( 13) = 0.36264482D+03
           C( 13) = 0.15882528D+00
          EX( 14) = 0.99876996D+02
           C( 14) = 0.17538208D+00
          EX( 15) = 0.16281510D+03
           C( 15) = 0.10638082D+00
          EX( 16) = 0.17533548D+02
           C( 16) =-0.67654981D+00
          EX( 17) = 0.71134955D+01
           C( 17) =-0.44089688D+00
          EX( 18) = 0.13198646D+02
           C( 18) = 0.25457122D-02
          EX( 19) = 0.20292497D+01
           C( 19) =-0.85785038D-02
          EX( 20) = 0.76788917D+00
           C( 20) =-0.42089854D-02
          EX( 21) = 0.13392119D+05
           C( 21) = 0.27943620D-02
          EX( 22) = 0.16721532D+04
           C( 22) = 0.19691897D-01
          EX( 23) = 0.36264482D+03
           C( 23) = 0.66033023D-01
          EX( 24) = 0.99876996D+02
           C( 24) = 0.72916659D-01
          EX( 25) = 0.16281510D+03
           C( 25) = 0.46461636D-01
          EX( 26) = 0.17533548D+02
           C( 26) =-0.29548195D+00
          EX( 27) = 0.71134955D+01
           C( 27) =-0.19256094D+00
          EX( 28) = 0.13198646D+02
           C( 28) =-0.24481731D+00
          EX( 29) = 0.20292497D+01
           C( 29) = 0.82498178D+00
          EX( 30) = 0.76788917D+00
           C( 30) = 0.40477178D+00
          EX( 31) = 0.45895970D+03
           C( 31) = 0.32604720D-01
          EX( 32) = 0.99266403D+02
           C( 32) = 0.20995554D+00
          EX( 33) = 0.29966147D+02
           C( 33) = 0.52913517D+00
          EX( 34) = 0.10191812D+02
           C( 34) = 0.39786073D+00
          EX( 35) = 0.45560940D+01
           C( 35) = 0.85521062D-02
          EX( 36) = 0.16750734D+01
           C( 36) = 0.12973590D-01
          EX( 37) = 0.60268684D+00
           C( 37) = 0.45744721D-02
          RETURN
          END IF
C     --- 
        ELSE  IF (MMPTYP.EQ.111) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : NI/NR             CONFIGURATION : S1D09        *
C   VALENCE : (4S/3D)           BASIS SET     : ( 5/ 5)      *
C          EFFECTIVE CORE CHARGE : 10.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =     -39.654893 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2363 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.4585 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE - 10.0D+00
          SHINFO='4S/3D'
          NOSH= 2
          CDFCTR(1)='41'
          CDFCTR(2)='41'
             KNG(1)= 5
            KTYP(1)= 1
             KNG(2)= 5
            KTYP(2)= 3
          EX( 1)= 0.77339670D+01
           C( 1)= 0.10142343D+00
          EX( 2)= 0.19194020D+01
           C( 2)=-0.21063491D+00
          EX( 3)= 0.83055590D+00
           C( 3)=-0.10451619D+00
          EX( 4)= 0.99221157D-01
           C( 4)= 0.59963631D+00
          EX( 5)= 0.35595520D-01
           C( 5)= 0.49993086D+00
          EX( 6)= 0.34879960D+02
           C( 6)= 0.46395669D-01
          EX( 7)= 0.89569990D+01
           C( 7)= 0.24473083D+00
          EX( 8)= 0.28127530D+01
           C( 8)= 0.44793706D+00
          EX( 9)= 0.86238910D+00
           C( 9)= 0.43471858D+00
          EX(10)= 0.23475790D+00
           C(10)= 0.21166169D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE - 10.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.62914000D+00
          ALPN0( 1)= 0.16817070D+03
            AN0( 2)= 0.57269000D+00
          ALPN0( 2)= 0.65448400D+01
            AN0( 3)= 0.20887000D+00
          ALPN0( 3)= 0.13363100D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.33016200D+01
          ALPN1( 1)= 0.21430259D+04
            AN1( 2)= 0.39731500D+01
          ALPN1( 2)= 0.71256990D+02
            AN1( 3)= 0.91355000D-01
          ALPN1( 3)= 0.11793400D+01
C         --- CORE-SHELL INFORMATION
          NOSH     =  5
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.30539656D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.37671990D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.46553480D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.32700360D+02
          KTYP( 5) =  2
           KNG( 5) =  7
          BPAR( 5) = 0.30610460D+01
C         --- EXPONENTS AND COEFFICIENTS
          EX(  1) = 0.10652750D+05
           C(  1) = 0.20493674D-01
          EX(  2) = 0.14666059D+04
           C(  2) = 0.14442462D+00
          EX(  3) = 0.33474787D+03
           C(  3) = 0.46811609D+00
          EX(  4) = 0.95240478D+02
           C(  4) = 0.49225018D+00
          EX(  5) = 0.15435288D+03
           C(  5) =-0.36368195D-02
          EX(  6) = 0.16925014D+02
           C(  6) = 0.22029687D-01
          EX(  7) = 0.68736510D+01
           C(  7) = 0.13774025D-01
          EX(  8) = 0.12796057D+02
           C(  8) = 0.43133944D-02
          EX(  9) = 0.19837763D+01
           C(  9) =-0.14443393D-01
          EX( 10) = 0.74942353D+00
           C( 10) =-0.71102714D-02
          EX( 11) = 0.10652750D+05
           C( 11) = 0.68833786D-02
          EX( 12) = 0.14666059D+04
           C( 12) = 0.48509086D-01
          EX( 13) = 0.33474787D+03
           C( 13) = 0.15723000D+00
          EX( 14) = 0.95240478D+02
           C( 14) = 0.16533611D+00
          EX( 15) = 0.15435288D+03
           C( 15) = 0.11363153D+00
          EX( 16) = 0.16925014D+02
           C( 16) =-0.68831219D+00
          EX( 17) = 0.68736510D+01
           C( 17) =-0.43036604D+00
          EX( 18) = 0.12796057D+02
           C( 18) = 0.24156818D-02
          EX( 19) = 0.19837763D+01
           C( 19) =-0.80889060D-02
          EX( 20) = 0.74942353D+00
           C( 20) =-0.39820501D-02
          EX( 21) = 0.10652750D+05
           C( 21) = 0.25912225D-02
          EX( 22) = 0.14666059D+04
           C( 22) = 0.18261067D-01
          EX( 23) = 0.33474787D+03
           C( 23) = 0.59188652D-01
          EX( 24) = 0.95240478D+02
           C( 24) = 0.62240170D-01
          EX( 25) = 0.15435288D+03
           C( 25) = 0.49559875D-01
          EX( 26) = 0.16925014D+02
           C( 26) =-0.30020423D+00
          EX( 27) = 0.68736510D+01
           C( 27) =-0.18770219D+00
          EX( 28) = 0.12796057D+02
           C( 28) =-0.24638723D+00
          EX( 29) = 0.19837763D+01
           C( 29) = 0.82502717D+00
          EX( 30) = 0.74942353D+00
           C( 30) = 0.40614881D+00
          EX( 31) = 0.43946959D+03
           C( 31) = 0.32878046D-01
          EX( 32) = 0.96439709D+02
           C( 32) = 0.21193140D+00
          EX( 33) = 0.29354518D+02
           C( 33) = 0.53005515D+00
          EX( 34) = 0.10040448D+02
           C( 34) = 0.39375408D+00
          EX( 35) = 0.44805544D+01
           C( 35) = 0.83788128D-02
          EX( 36) = 0.16495272D+01
           C( 36) = 0.12580404D-01
          EX( 37) = 0.59472130D+00
           C( 37) = 0.44156678D-02
          EX( 38) = 0.43946959D+03
           C( 38) = 0.13255557D-01
          EX( 39) = 0.96439709D+02
           C( 39) = 0.85445122D-01
          EX( 40) = 0.29354518D+02
           C( 40) = 0.21370418D+00
          EX( 41) = 0.10040448D+02
           C( 41) = 0.15875121D+00
          EX( 42) = 0.44805544D+01
           C( 42) =-0.39664220D+00
          EX( 43) = 0.16495272D+01
           C( 43) =-0.59554012D+00
          EX( 44) = 0.59472130D+00
           C( 44) =-0.20903202D+00
          RETURN
          END IF
        END IF
        IERR=2
        RETURN
      END IF
C
C     --- 
C
      IF(IZSAVE.EQ.29) THEN
        ITYP(1)=111
        ITYP(2)=112
C     --- 
        IF (MMPTYP.EQ.112) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : CU/NR             CONFIGURATION : S1P6D10      *
C   VALENCE : (4S/3P/3D)        BASIS SET     : ( 6/ 4/ 5)   *
C          EFFECTIVE CORE CHARGE : 17.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =    -150.052569 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2392 A.U.
C  ***     ORBITAL ENERGY( 3P ) =   -3.3247 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.4911 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE - 17.0D+00
          SHINFO='4S/3P/3D'
          NOSH= 3
          CDFCTR(1)='51'
          CDFCTR(2)='4'
          CDFCTR(3)='41'
             KNG(1)= 6
            KTYP(1)= 1
             KNG(2)= 4
            KTYP(2)= 2
             KNG(3)= 5
            KTYP(3)= 3
          EX( 1)= 0.31408700D+03
           C( 1)=-0.23604615D-01
          EX( 2)= 0.13104800D+02
           C( 2)= 0.12066483D+00
          EX( 3)= 0.20940100D+01
           C( 3)=-0.19332115D+00
          EX( 4)= 0.89382600D+00
           C( 4)=-0.17716138D+00
          EX( 5)= 0.10243900D+00
           C( 5)= 0.61259907D+00
          EX( 6)= 0.36393900D-01
           C( 6)= 0.49387151D+00
          EX( 7)= 0.14379000D+03
           C( 7)=-0.77730003D-01
          EX( 8)= 0.30245900D+02
           C( 8)=-0.30091288D+00
          EX( 9)= 0.28939300D+01
           C( 9)= 0.68025342D+00
          EX(10)= 0.91170200D+00
           C(10)= 0.40581177D+00
          EX(11)= 0.38422400D+02
           C(11)= 0.54663321D-01
          EX(12)= 0.99126000D+01
           C(12)= 0.23410395D+00
          EX(13)= 0.31211200D+01
           C(13)= 0.44094212D+00
          EX(14)= 0.95779200D+00
           C(14)= 0.44150605D+00
          EX(15)= 0.26013300D+00
           C(15)= 0.22088930D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE - 17.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.23368545D+00
          ALPN0( 1)= 0.24407286D+03
            AN0( 2)= 0.10275597D+00
          ALPN0( 2)= 0.91861364D+01
            AN0( 3)= 0.19774878D-01
          ALPN0( 3)= 0.20511895D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.49900297D+02
          ALPN1( 1)= 0.15249023D+04
            AN1( 2)= 0.35075491D+01
          ALPN1( 2)= 0.44572888D+02
            AN1( 3)= 0.34409842D-01
          ALPN1( 3)= 0.12981176D+01
C         --- CORE-SHELL INFORMATION 
          NOSH     =  4
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.65758600D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.81638000D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.10023960D+02
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.71236000D+02
C         --- EXPONENTS AND COEFFICIENTS 
          EX(  1) = 0.11451708D+05
           C(  1) = 0.20450809D-01
          EX(  2) = 0.15766503D+04
           C(  2) = 0.14417280D+00
          EX(  3) = 0.35991536D+03
           C(  3) = 0.46774949D+00
          EX(  4) = 0.10244929D+03
           C(  4) = 0.49271733D+00
          EX(  5) = 0.16612249D+03
           C(  5) =-0.36574362D-02
          EX(  6) = 0.18281854D+02
           C(  6) = 0.22122604D-01
          EX(  7) = 0.74321657D+01
           C(  7) = 0.13738171D-01
          EX(  8) = 0.13933294D+02
           C(  8) = 0.42625629D-02
          EX(  9) = 0.21617946D+01
           C(  9) =-0.14257601D-01
          EX( 10) = 0.81446639D+00
           C( 10) =-0.70514003D-02
          EX( 11) = 0.11451708D+05
           C( 11) = 0.68993176D-02
          EX( 12) = 0.15766503D+04
           C( 12) = 0.48638369D-01
          EX( 13) = 0.35991536D+03
           C( 13) = 0.15780072D+00
          EX( 14) = 0.10244929D+03
           C( 14) = 0.16622391D+00
          EX( 15) = 0.16612249D+03
           C( 15) = 0.11415422D+00
          EX( 16) = 0.18281854D+02
           C( 16) =-0.69048062D+00
          EX( 17) = 0.74321657D+01
           C( 17) =-0.42878952D+00
          EX( 18) = 0.13933294D+02
           C( 18) = 0.24098127D-02
          EX( 19) = 0.21617946D+01
           C( 19) =-0.80604437D-02
          EX( 20) = 0.81446639D+00
           C( 20) =-0.39864642D-02
          EX( 21) = 0.11451708D+05
           C( 21) = 0.26093779D-02
          EX( 22) = 0.15766503D+04
           C( 22) = 0.18395426D-01
          EX( 23) = 0.35991536D+03
           C( 23) = 0.59681512D-01
          EX( 24) = 0.10244929D+03
           C( 24) = 0.62867232D-01
          EX( 25) = 0.16612249D+03
           C( 25) = 0.50086249D-01
          EX( 26) = 0.18281854D+02
           C( 26) =-0.30295492D+00
          EX( 27) = 0.74321657D+01
           C( 27) =-0.18813547D+00
          EX( 28) = 0.13933294D+02
           C( 28) =-0.24656320D+00
          EX( 29) = 0.21617946D+01
           C( 29) = 0.82471503D+00
          EX( 30) = 0.81446639D+00
           C( 30) = 0.40788040D+00
          EX( 31) = 0.47576413D+03
           C( 31) = 0.32718830D-01
          EX( 32) = 0.10448556D+03
           C( 32) = 0.21155588D+00
          EX( 33) = 0.31855255D+02
           C( 33) = 0.53038368D+00
          EX( 34) = 0.10915745D+02
           C( 34) = 0.39319495D+00
          EX( 35) = 0.49109969D+01
           C( 35) = 0.84469771D-02
          EX( 36) = 0.18062307D+01
           C( 36) = 0.12651072D-01
          EX( 37) = 0.64986476D+00
           C( 37) = 0.44522110D-02
          RETURN
          END IF
C     --- 
        ELSE IF (MMPTYP.EQ.111) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : CU/NR             CONFIGURATION : S1D10        *
C   VALENCE : (4S/3D)           BASIS SET     : ( 5/ 5)      *
C          EFFECTIVE CORE CHARGE : 11.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =     -51.127321 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2382 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.4906 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE - 11.0D+00
          SHINFO='4S/3D'
          NOSH= 2
          CDFCTR(1)='41'
          CDFCTR(2)='41'
             KNG(1)= 5
            KTYP(1)= 1
             KNG(2)= 5
            KTYP(2)= 3
          EX( 1)= 0.83699390D+01
           C( 1)= 0.94256283D-01
          EX( 2)= 0.20939424D+01
           C( 2)=-0.19870783D+00
          EX( 3)= 0.89378650D+00
           C( 3)=-0.10899261D+00
          EX( 4)= 0.10244044D+00
           C( 4)= 0.59128773D+00
          EX( 5)= 0.36394200D-01
           C( 5)= 0.50666438D+00
          EX( 6)= 0.38422407D+02
           C( 6)= 0.44153673D-01
          EX( 7)= 0.99125980D+01
           C( 7)= 0.24236690D+00
          EX( 8)= 0.31211210D+01
           C( 8)= 0.45337072D+00
          EX( 9)= 0.95779210D+00
           C( 9)= 0.43128626D+00
          EX(10)= 0.26013300D+00
           C(10)= 0.21119500D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE - 11.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.50669000D+00
          ALPN0( 1)= 0.19760830D+03
            AN0( 2)= 0.47412400D+00
          ALPN0( 2)= 0.60961000D+01
            AN0( 3)= 0.17677000D+00
          ALPN0( 3)= 0.13878600D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.26554600D+01
          ALPN1( 1)= 0.26736090D+04
            AN1( 2)= 0.50138300D+01
          ALPN1( 2)= 0.75869240D+02
            AN1( 3)= 0.79232000D-01
          ALPN1( 3)= 0.12909700D+01
C         --- CORE-SHELL INFORMATION
          NOSH     =  5
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.32879170D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.40818040D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.50112600D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.35616990D+02
          KTYP( 5) =  2
           KNG( 5) =  7
          BPAR( 5) = 0.33241500D+01
C         --- EXPONENTS AND COEFFICIENTS
          EX(  1) = 0.11451708D+05
           C(  1) = 0.20450809D-01
          EX(  2) = 0.15766503D+04
           C(  2) = 0.14417280D+00
          EX(  3) = 0.35991536D+03
           C(  3) = 0.46774949D+00
          EX(  4) = 0.10244929D+03
           C(  4) = 0.49271733D+00
          EX(  5) = 0.16612249D+03
           C(  5) =-0.36574362D-02
          EX(  6) = 0.18281854D+02
           C(  6) = 0.22122604D-01
          EX(  7) = 0.74321657D+01
           C(  7) = 0.13738171D-01
          EX(  8) = 0.13933294D+02
           C(  8) = 0.42625629D-02
          EX(  9) = 0.21617946D+01
           C(  9) =-0.14257601D-01
          EX( 10) = 0.81446639D+00
           C( 10) =-0.70514003D-02
          EX( 11) = 0.11451708D+05
           C( 11) = 0.68993176D-02
          EX( 12) = 0.15766503D+04
           C( 12) = 0.48638369D-01
          EX( 13) = 0.35991536D+03
           C( 13) = 0.15780072D+00
          EX( 14) = 0.10244929D+03
           C( 14) = 0.16622391D+00
          EX( 15) = 0.16612249D+03
           C( 15) = 0.11415422D+00
          EX( 16) = 0.18281854D+02
           C( 16) =-0.69048062D+00
          EX( 17) = 0.74321657D+01
           C( 17) =-0.42878952D+00
          EX( 18) = 0.13933294D+02
           C( 18) = 0.24098127D-02
          EX( 19) = 0.21617946D+01
           C( 19) =-0.80604437D-02
          EX( 20) = 0.81446639D+00
           C( 20) =-0.39864642D-02
          EX( 21) = 0.11451708D+05
           C( 21) = 0.26093779D-02
          EX( 22) = 0.15766503D+04
           C( 22) = 0.18395426D-01
          EX( 23) = 0.35991536D+03
           C( 23) = 0.59681512D-01
          EX( 24) = 0.10244929D+03
           C( 24) = 0.62867232D-01
          EX( 25) = 0.16612249D+03
           C( 25) = 0.50086249D-01
          EX( 26) = 0.18281854D+02
           C( 26) =-0.30295492D+00
          EX( 27) = 0.74321657D+01
           C( 27) =-0.18813547D+00
          EX( 28) = 0.13933294D+02
           C( 28) =-0.24656320D+00
          EX( 29) = 0.21617946D+01
           C( 29) = 0.82471503D+00
          EX( 30) = 0.81446639D+00
           C( 30) = 0.40788040D+00
          EX( 31) = 0.47576413D+03
           C( 31) = 0.32718830D-01
          EX( 32) = 0.10448556D+03
           C( 32) = 0.21155588D+00
          EX( 33) = 0.31855255D+02
           C( 33) = 0.53038368D+00
          EX( 34) = 0.10915745D+02
           C( 34) = 0.39319495D+00
          EX( 35) = 0.49109969D+01
           C( 35) = 0.84469772D-02
          EX( 36) = 0.18062307D+01
           C( 36) = 0.12651072D-01
          EX( 37) = 0.64986476D+00
           C( 37) = 0.44522110D-02
          EX( 38) = 0.47576413D+03
           C( 38) = 0.13294517D-01
          EX( 39) = 0.10448556D+03
           C( 39) = 0.85960690D-01
          EX( 40) = 0.31855255D+02
           C( 40) = 0.21550877D+00
          EX( 41) = 0.10915745D+02
           C( 41) = 0.15976540D+00
          EX( 42) = 0.49109969D+01
           C( 42) =-0.39766506D+00
          EX( 43) = 0.18062307D+01
           C( 43) =-0.59558458D+00
          EX( 44) = 0.64986476D+00
           C( 44) =-0.20960027D+00
          RETURN
          END IF 
        END IF
        IERR=2
        RETURN
      END IF
C
C     --- 
C
      IF(IZSAVE.EQ.30) THEN
        ITYP(1)=111
        ITYP(2)=112
C     --- 
        IF (MMPTYP.EQ.112) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : ZN/NR             CONFIGURATION : S2P6D10      *
C   VALENCE : (4S/3P/3D)        BASIS SET     : ( 6/ 4/ 5)   *
C          EFFECTIVE CORE CHARGE : 18.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =    -175.716597 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2958 A.U.
C  ***     ORBITAL ENERGY( 3P ) =   -3.8393 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.7826 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE - 18.0D+00
          SHINFO='4S/3P/3D'
          NOSH= 3
          CDFCTR(1)='51'
          CDFCTR(2)='4'
          CDFCTR(3)='41'
             KNG(1)= 6
            KTYP(1)= 1
             KNG(2)= 4
            KTYP(2)= 2
             KNG(3)= 5
            KTYP(3)= 3
          EX( 1)= 0.33698700D+03
           C( 1)=-0.30239768D-01
          EX( 2)= 0.14134100D+02
           C( 2)= 0.15056642D+00
          EX( 3)= 0.23240800D+01
           C( 3)=-0.23307175D+00
          EX( 4)= 0.10748900D+01
           C( 4)=-0.20417549D+00
          EX( 5)= 0.13763900D+00
           C( 5)= 0.65198832D+00
          EX( 6)= 0.47636700D-01
           C( 6)= 0.47064901D+00
          EX( 7)= 0.15529500D+03
           C( 7)=-0.79179708D-01
          EX( 8)= 0.32761300D+02
           C( 8)=-0.30638946D+00
          EX( 9)= 0.31718300D+01
           C( 9)= 0.68154522D+00
          EX(10)= 0.10192400D+01
           C(10)= 0.40142623D+00
          EX(11)= 0.46785650D+02
           C(11)= 0.50272346D-01
          EX(12)= 0.12256650D+02
           C(12)= 0.22106931D+00
          EX(13)= 0.39611390D+01
           C(13)= 0.43062324D+00
          EX(14)= 0.12734030D+01
           C(14)= 0.45117191D+00
          EX(15)= 0.37109580D+00
           C(15)= 0.21213534D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE - 18.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.20527379D+00
          ALPN0( 1)= 0.24349404D+03
            AN0( 2)= 0.89977916D-01
          ALPN0( 2)= 0.95709463D+01
            AN0( 3)= 0.17267383D-01
          ALPN0( 3)= 0.21341005D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.51335324D+02
          ALPN1( 1)= 0.15205700D+04
            AN1( 2)= 0.34020266D+01
          ALPN1( 2)= 0.44437235D+02
            AN1( 3)= 0.33341897D-01
          ALPN1( 3)= 0.13267392D+01
C         --- CORE-SHELL INFORMATION 
          NOSH     =  4
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.70661000D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.88723400D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.11275640D+02
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.77849600D+02
C         --- EXPONENTS AND COEFFICIENTS 
          EX(  1) = 0.12280003D+05
           C(  1) = 0.20409934D-01
          EX(  2) = 0.16907375D+04
           C(  2) = 0.14393180D+00
          EX(  3) = 0.38600886D+03
           C(  3) = 0.46739417D+00
          EX(  4) = 0.10992456D+03
           C(  4) = 0.49316877D+00
          EX(  5) = 0.17830432D+03
           C(  5) =-0.36788701D-02
          EX(  6) = 0.19695085D+02
           C(  6) = 0.22207051D-01
          EX(  7) = 0.80186961D+01
           C(  7) = 0.13725373D-01
          EX(  8) = 0.15090719D+02
           C(  8) = 0.42384490D-02
          EX(  9) = 0.23615856D+01
           C(  9) =-0.14043492D-01
          EX( 10) = 0.90854830D+00
           C( 10) =-0.69615657D-02
          EX( 11) = 0.12280003D+05
           C( 11) = 0.69135775D-02
          EX( 12) = 0.16907375D+04
           C( 12) = 0.48754865D-01
          EX( 13) = 0.38600886D+03
           C( 13) = 0.15832319D+00
          EX( 14) = 0.10992456D+03
           C( 14) = 0.16705396D+00
          EX( 15) = 0.17830432D+03
           C( 15) = 0.11464349D+00
          EX( 16) = 0.19695085D+02
           C( 16) =-0.69203145D+00
          EX( 17) = 0.80186961D+01
           C( 17) =-0.42771954D+00
          EX( 18) = 0.15090719D+02
           C( 18) = 0.24388512D-02
          EX( 19) = 0.23615856D+01
           C( 19) =-0.80807829D-02
          EX( 20) = 0.90854830D+00
           C( 20) =-0.40057630D-02
          EX( 21) = 0.12280003D+05
           C( 21) = 0.26399581D-02
          EX( 22) = 0.16907375D+04
           C( 22) = 0.18617106D-01
          EX( 23) = 0.38600886D+03
           C( 23) = 0.60455905D-01
          EX( 24) = 0.10992456D+03
           C( 24) = 0.63789766D-01
          EX( 25) = 0.17830432D+03
           C( 25) = 0.50869057D-01
          EX( 26) = 0.19695085D+02
           C( 26) =-0.30706486D+00
          EX( 27) = 0.80186961D+01
           C( 27) =-0.18978566D+00
          EX( 28) = 0.15090719D+02
           C( 28) =-0.24913832D+00
          EX( 29) = 0.23615856D+01
           C( 29) = 0.82548402D+00
          EX( 30) = 0.90854830D+00
           C( 30) = 0.40920457D+00
          EX( 31) = 0.51426620D+03
           C( 31) = 0.32486666D-01
          EX( 32) = 0.11302553D+03
           C( 32) = 0.21075316D+00
          EX( 33) = 0.34512986D+02
           C( 33) = 0.53017575D+00
          EX( 34) = 0.11849065D+02
           C( 34) = 0.39348697D+00
          EX( 35) = 0.53766405D+01
           C( 35) = 0.85775031D-02
          EX( 36) = 0.19931676D+01
           C( 36) = 0.12840580D-01
          EX( 37) = 0.73128626D+00
           C( 37) = 0.44233661D-02
          RETURN
          END IF
C     --- 
        ELSE IF (MMPTYP.EQ.111) THEN
          IF (MMPMOD.EQ.0) THEN
C
C*************************************************************
C     TITLE : ZN/NR             CONFIGURATION : S2D10        *
C   VALENCE : (4S/3D)           BASIS SET     : ( 5/ 5)      *
C          EFFECTIVE CORE CHARGE : 12.000000                 *
C*************************************************************
C
C  ***     TOTAL ENERGY         =     -63.760358 A.U.
C
C  ***     ORBITAL ENERGY( 4S ) =    -.2967 A.U.
C  ***     ORBITAL ENERGY( 3D ) =    -.7829 A.U.
C
C         --- VALENCE-SHELL BASIS-SET
C
          ZCORE = IZSAVE - 12.0D+00
          SHINFO='4S/3D'
          NOSH= 2
          CDFCTR(1)='41'
          CDFCTR(2)='41'
             KNG(1)= 5
            KTYP(1)= 1
             KNG(2)= 5
            KTYP(2)= 3
          EX( 1)= 0.89222960D+01
           C( 1)=-0.11487269D+00
          EX( 2)= 0.23240040D+01
           C( 2)= 0.24383587D+00
          EX( 3)= 0.10748560D+01
           C( 3)= 0.11802037D+00
          EX( 4)= 0.13763990D+00
           C( 4)=-0.64020794D+00
          EX( 5)= 0.47636890D-01
           C( 5)=-0.47299261D+00
          EX( 6)= 0.46785650D+02
           C( 6)= 0.36974988D-01
          EX( 7)= 0.12256650D+02
           C( 7)= 0.21734361D+00
          EX( 8)= 0.39611390D+01
           C( 8)= 0.44218856D+00
          EX( 9)= 0.12734030D+01
           C( 9)= 0.44821006D+00
          EX(10)= 0.37109580D+00
           C(10)= 0.20714485D+00
          RETURN
        ELSE
C
C         --- CORE-SHELLS AND MODEL-POTENTIAL
C
          ZCORE = IZSAVE - 12.0D+00
          NOAN0=  3
          NOAN1=  3
C         --- A-PARAMETERS FOR N=0
            AN0( 1)= 0.49547658D+00
          ALPN0( 1)= 0.16915063D+03
            AN0( 2)= 0.44446598D+00
          ALPN0( 2)= 0.64800497D+01
            AN0( 3)= 0.14248224D+00
          ALPN0( 3)= 0.13470074D+01
C         --- A-PARAMETERS FOR N=1
            AN1( 1)= 0.31227413D+01
          ALPN1( 1)= 0.24843270D+04
            AN1( 2)= 0.34893581D+01
          ALPN1( 2)= 0.75158190D+02
            AN1( 3)= 0.69555855D-01
          ALPN1( 3)= 0.13774368D+01
C         --- CORE-SHELL INFORMATION
          NOSH     =  5
          KTYP( 1) =  1
           KNG( 1) = 10
          BPAR( 1) = 0.35330450D+03
          KTYP( 2) =  1
           KNG( 2) = 10
          BPAR( 2) = 0.44361700D+02
          KTYP( 3) =  1
           KNG( 3) = 10
          BPAR( 3) = 0.56378000D+01
          KTYP( 4) =  2
           KNG( 4) =  7
          BPAR( 4) = 0.38924800D+02
          KTYP( 5) =  2
           KNG( 5) =  7
          BPAR( 5) = 0.38394000D+01
C         --- EXPONENTS AND COEFFICIENTS
          EX(  1) = 0.12280003D+05
           C(  1) = 0.20409934D-01
          EX(  2) = 0.16907375D+04
           C(  2) = 0.14393180D+00
          EX(  3) = 0.38600886D+03
           C(  3) = 0.46739417D+00
          EX(  4) = 0.10992456D+03
           C(  4) = 0.49316877D+00
          EX(  5) = 0.17830432D+03
           C(  5) =-0.36788701D-02
          EX(  6) = 0.19695085D+02
           C(  6) = 0.22207051D-01
          EX(  7) = 0.80186961D+01
           C(  7) = 0.13725373D-01
          EX(  8) = 0.15090719D+02
           C(  8) = 0.42384490D-02
          EX(  9) = 0.23615856D+01
           C(  9) =-0.14043492D-01
          EX( 10) = 0.90854830D+00
           C( 10) =-0.69615657D-02
          EX( 11) = 0.12280003D+05
           C( 11) = 0.69135775D-02
          EX( 12) = 0.16907375D+04
           C( 12) = 0.48754865D-01
          EX( 13) = 0.38600886D+03
           C( 13) = 0.15832319D+00
          EX( 14) = 0.10992456D+03
           C( 14) = 0.16705396D+00
          EX( 15) = 0.17830432D+03
           C( 15) = 0.11464349D+00
          EX( 16) = 0.19695085D+02
           C( 16) =-0.69203145D+00
          EX( 17) = 0.80186961D+01
           C( 17) =-0.42771954D+00
          EX( 18) = 0.15090719D+02
           C( 18) = 0.24388512D-02
          EX( 19) = 0.23615856D+01
           C( 19) =-0.80807829D-02
          EX( 20) = 0.90854830D+00
           C( 20) =-0.40057630D-02
          EX( 21) = 0.12280003D+05
           C( 21) = 0.26399581D-02
          EX( 22) = 0.16907375D+04
           C( 22) = 0.18617106D-01
          EX( 23) = 0.38600886D+03
           C( 23) = 0.60455905D-01
          EX( 24) = 0.10992456D+03
           C( 24) = 0.63789766D-01
          EX( 25) = 0.17830432D+03
           C( 25) = 0.50869057D-01
          EX( 26) = 0.19695085D+02
           C( 26) =-0.30706486D+00
          EX( 27) = 0.80186961D+01
           C( 27) =-0.18978566D+00
          EX( 28) = 0.15090719D+02
           C( 28) =-0.24913832D+00
          EX( 29) = 0.23615856D+01
           C( 29) = 0.82548402D+00
          EX( 30) = 0.90854830D+00
           C( 30) = 0.40920457D+00
          EX( 31) = 0.51426620D+03
           C( 31) = 0.32486666D-01
          EX( 32) = 0.11302553D+03
           C( 32) = 0.21075316D+00
          EX( 33) = 0.34512986D+02
           C( 33) = 0.53017575D+00
          EX( 34) = 0.11849065D+02
           C( 34) = 0.39348697D+00
          EX( 35) = 0.53766405D+01
           C( 35) = 0.85775031D-02
          EX( 36) = 0.19931676D+01
           C( 36) = 0.12840580D-01
          EX( 37) = 0.73128626D+00
           C( 37) = 0.44233661D-02
          EX( 38) = 0.51426620D+03
           C( 38) = 0.13338275D-01
          EX( 39) = 0.11302553D+03
           C( 39) = 0.86530389D-01
          EX( 40) = 0.34512986D+02
           C( 40) = 0.21767793D+00
          EX( 41) = 0.11849065D+02
           C( 41) = 0.16155668D+00
          EX( 42) = 0.53766405D+01
           C( 42) =-0.39883526D+00
          EX( 43) = 0.19931676D+01
           C( 43) =-0.59705905D+00
          EX( 44) = 0.73128626D+00
           C( 44) =-0.20567691D+00
          RETURN
          END IF
        END IF
        IERR=2
        RETURN
      END IF
C ...
C
C--   IF (MODE.EQ.0) THEN
C--      IERR=1
C--      RETURN
C--   END IF
C
      WRITE(IW,9999) IZSAVE
      CALL ABRT
      STOP
 9999 FORMAT(/1X,'*** NO MOD-POT AVAILABLE FOR ATOM A=',I4)
      END
