C 22 DEC 06 - ST,NK,MC - ADD LC EXCHANGE ARGUMENTS
C  5 FEB 05 - MWS - CHANGE COMMON NAME GAMMAF TO BE FMTTBL INSTEAD
C 10 NOV 04 - MWS - ADJUST CC TOLERANCE
C  7 SEP 04 - GDF - NEW MODULE FOR ERIC'S SPECIAL G ROUTINES
C
C  LPHI=     255
C  LWK1=       1
C  LWK2=      46
C  LENW=      46
C*MODULE INT2G   *DECK GSSS
      SUBROUTINE GSSS (IPRIM,JPRIM,KPRIM,LPRIM,IEQJ,KEQL,
     *                 XC,YC,ZC,XD,YD,ZD,PHI,WRK2,LENW)
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      COMMON /SHLNOS/ QQ4,IDUMMY(20)
      PARAMETER (NTX=4,NPF=450,NGRD=7,NPX=1000,MXQT=16)
      COMMON /FMTTBL/ FGRID(0:NTX,0:NPF,0:NGRD)
     *,XGRID(0:NTX,0:NPX),TMAX,RFINC(0:NGRD)
     *,RXINC,RMR(MXQT),TLGM(0:MXQT),NORD
      PARAMETER (MXPR=30, MXPS=MXPR*MXPR)
      COMMON /ERIPRM/ EXI(MXPR),EXJ(MXPR),EXK(MXPR),EXL(MXPR)
     *,CCI(MXPR),CCJ(MXPR),CCK(MXPR),CCL(MXPR)
     *,XAB(MXPS),YAB(MXPS),ZAB(MXPS)
     *,XCD(MXPS),YCD(MXPS),ZCD(MXPS)
     *,CCBRA(MXPS),CCKET(MXPS),RXB(MXPS)
     *,SLBRA(MXPS),SLKET(MXPS),RXK(MXPS)
      COMMON /NLRC  / LCFLAG,LRINT,EMU,EMU2,LRFILE
      LOGICAL   LCFLAG,LRINT
      LOGICAL IEQJ,KEQL
      DOUBLE PRECISION RXYZ(3),SFAC(9),CONF(11),FT(0:16)
      DOUBLE PRECISION SAUX(0: 4),CXYZ( 35)
      DOUBLE PRECISION PHI(*),WRK2(LENW,*)
      PARAMETER (ZERO=0.0D+00, ONE=1.0D+00, TWO=2.0D+00)
      PARAMETER (D01=1.0D+00, D02=2.0D+00, D03=3.0D+00)
      PARAMETER (CCTOL=1.0D-13)
      DATA CONF/0.0D+00,0.0D+00,0.0D+00,1.0D+00,2.0D+00,3.0D+00
     &,         4.0D+00,5.0D+00,6.0D+00,7.0D+00,8.0D+00/
      SAVE CONF
      RFAC=RFINC( 4)
      KL=0
      PHI(      50)=ZERO
      PHI(      52)=ZERO
      PHI(      54)=ZERO
      PHI(      64)=ZERO
      PHI(      65)=ZERO
      PHI(      66)=ZERO
      PHI(      70)=ZERO
      PHI(      71)=ZERO
      PHI(      72)=ZERO
      DO I=      94,     100
      PHI(I)=ZERO
      END DO
      DO I=     108,     114
      PHI(I)=ZERO
      END DO
      DO I=     154,     166
      PHI(I)=ZERO
      END DO
      DO I=     233,     254
      PHI(I)=ZERO
      END DO
      DO LL=1,LPRIM
      PHI(      49)=ZERO
      PHI(      51)=ZERO
      PHI(      53)=ZERO
      PHI(      61)=ZERO
      PHI(      62)=ZERO
      PHI(      63)=ZERO
      PHI(      67)=ZERO
      PHI(      68)=ZERO
      PHI(      69)=ZERO
      DO I=      87,      93
      PHI(I)=ZERO
      END DO
      DO I=     101,     107
      PHI(I)=ZERO
      END DO
      DO I=     141,     153
      PHI(I)=ZERO
      END DO
      DO I=     211,     232
      PHI(I)=ZERO
      END DO
      KTOP=KPRIM
      IF (KEQL) KTOP=LL
      DO KK=1,KTOP
      KL=KL+1
      CCKF=CCKET(KL)
      IF (ABS(CCKF).GT.CCTOL) THEN
      EKL=EXL(LL)+EXK(KK)
      XKL=XCD(KL)
      YKL=YCD(KL)
      ZKL=ZCD(KL)
      IJ=0
      PHI(      48)=ZERO
      PHI(      58)=ZERO
      PHI(      59)=ZERO
      PHI(      60)=ZERO
      DO I=      80,      86
      PHI(I)=ZERO
      END DO
      DO I=     128,     140
      PHI(I)=ZERO
      END DO
      DO I=     189,     210
      PHI(I)=ZERO
      END DO
      DO JJ=1,JPRIM
      PHI(      47)=ZERO
      PHI(      55)=ZERO
      PHI(      56)=ZERO
      PHI(      57)=ZERO
      DO I=      73,      79
      PHI(I)=ZERO
      END DO
      DO I=     115,     127
      PHI(I)=ZERO
      END DO
      DO I=     167,     188
      PHI(I)=ZERO
      END DO
      ITOP=IPRIM
      IF (IEQJ) ITOP=JJ
      DO II=1,ITOP
      IJ=IJ+1
      CCBF=CCBRA(IJ)
      IF (ABS(CCBF).GT.CCTOL) THEN
      EIJ=EXJ(JJ)+EXI(II)
      EIJKL=EIJ+EKL
      EINV=ONE/EIJKL
      RX=XKL-XAB(IJ)
      RY=YKL-YAB(IJ)
      RZ=ZKL-ZAB(IJ)
      R2=RX*RX+RY*RY+RZ*RZ
      VS=EIJ*EKL*EINV
      IF(LRINT) THEN
         VS0=VS
         VS=VS0*EMU2/(VS0+EMU2)
      ENDIF
      TT=VS*R2
      IF (TT.LE.TMAX) THEN
      TVAL=TT*RFAC
      IPT=NINT(TVAL)
      FX= TVAL*FGRID(4,IPT, 4)
      FX=(FX + FGRID(3,IPT, 4))*TVAL
      FX=(FX + FGRID(2,IPT, 4))*TVAL
      FX=(FX + FGRID(1,IPT, 4))*TVAL
      FX= FX + FGRID(0,IPT, 4)
      FT( 4)=FX
      TVAL=TT*RXINC
      IPT=NINT(TVAL)
      ET= TVAL*XGRID(4,IPT)
      ET=(ET + XGRID(3,IPT))*TVAL
      ET=(ET + XGRID(2,IPT))*TVAL
      ET=(ET + XGRID(1,IPT))*TVAL
      ET= ET + XGRID(0,IPT)
      T2=TT*TWO
      FT( 3)=RMR( 4)*(T2*FT( 4)+ET)
      FT( 2)=RMR( 3)*(T2*FT( 3)+ET)
      FT( 1)=RMR( 2)*(T2*FT( 2)+ET)
      FT( 0)=RMR( 1)*(T2*FT( 1)+ET)
      ZEROL=CCBF*CCKF*SQRT(EINV)*QQ4
      IF(LRINT)ZEROL=ZEROL*SQRT(VS/VS0)
      SAUX(0)=ZEROL*FT(0)
      VS2=TWO*VS
      VSM=ONE
      VSM=VSM*VS2
      SAUX( 1)=ZEROL*FT( 1)*VSM
      VSM=VSM*VS2
      SAUX( 2)=ZEROL*FT( 2)*VSM
      VSM=VSM*VS2
      SAUX( 3)=ZEROL*FT( 3)*VSM
      VSM=VSM*VS2
      SAUX( 4)=ZEROL*FT( 4)*VSM
      ELSE
      ZEROL=CCBF*CCKF*SLBRA(IJ)*SLKET(KL)*QQ4
      RINV=ONE/R2
      SAUX(0)=ZEROL*SQRT(RINV)
      RR=SAUX(0)
      RR=RR*RINV
      SAUX( 1)=TLGM( 1)*RR
      RR=RR*RINV
      SAUX( 2)=TLGM( 2)*RR
      RR=RR*RINV
      SAUX( 3)=TLGM( 3)*RR
      RR=RR*RINV
      SAUX( 4)=TLGM( 4)*RR
      END IF
      RXYZ(1)=RX
      RXYZ(2)=RY
      RXYZ(3)=RZ
      CXYZ(1)=ONE
      CXYZ(2)=RX
      CXYZ(3)=RY
      CXYZ(4)=RZ
      CXYZ(     5)=CXYZ(     2)*RXYZ(1)
      CXYZ(     6)=CXYZ(     3)*RXYZ(1)
      CXYZ(     7)=CXYZ(     3)*RXYZ(2)
      CXYZ(     8)=CXYZ(     4)*RXYZ(1)
      CXYZ(     9)=CXYZ(     4)*RXYZ(2)
      CXYZ(    10)=CXYZ(     4)*RXYZ(3)
      CXYZ(    11)=CXYZ(     5)*RXYZ(1)
      CXYZ(    12)=CXYZ(     6)*RXYZ(1)
      CXYZ(    13)=CXYZ(     7)*RXYZ(1)
      CXYZ(    14)=CXYZ(     7)*RXYZ(2)
      CXYZ(    15)=CXYZ(     8)*RXYZ(1)
      CXYZ(    16)=CXYZ(     9)*RXYZ(1)
      CXYZ(    17)=CXYZ(     9)*RXYZ(2)
      CXYZ(    18)=CXYZ(    10)*RXYZ(1)
      CXYZ(    19)=CXYZ(    10)*RXYZ(2)
      CXYZ(    20)=CXYZ(    10)*RXYZ(3)
      CXYZ(    21)=CXYZ(    11)*RXYZ(1)
      CXYZ(    22)=CXYZ(    12)*RXYZ(1)
      CXYZ(    23)=CXYZ(    13)*RXYZ(1)
      CXYZ(    24)=CXYZ(    14)*RXYZ(1)
      CXYZ(    25)=CXYZ(    14)*RXYZ(2)
      CXYZ(    26)=CXYZ(    15)*RXYZ(1)
      CXYZ(    27)=CXYZ(    16)*RXYZ(1)
      CXYZ(    28)=CXYZ(    17)*RXYZ(1)
      CXYZ(    29)=CXYZ(    17)*RXYZ(2)
      CXYZ(    30)=CXYZ(    18)*RXYZ(1)
      CXYZ(    31)=CXYZ(    19)*RXYZ(1)
      CXYZ(    32)=CXYZ(    19)*RXYZ(2)
      CXYZ(    33)=CXYZ(    20)*RXYZ(1)
      CXYZ(    34)=CXYZ(    20)*RXYZ(2)
      CXYZ(    35)=CXYZ(    20)*RXYZ(3)
      PHI(1)=SAUX(0)
      PHI(       2)=CXYZ(     2)*SAUX( 1)
      PHI(       3)=CXYZ(     3)*SAUX( 1)
      PHI(       4)=CXYZ(     4)*SAUX( 1)
      PHI(       5)=CXYZ(     1)*SAUX( 1)
      PHI(       6)=CXYZ(     5)*SAUX( 2)
      PHI(       7)=CXYZ(     6)*SAUX( 2)
      PHI(       8)=CXYZ(     7)*SAUX( 2)
      PHI(       9)=CXYZ(     8)*SAUX( 2)
      PHI(      10)=CXYZ(     9)*SAUX( 2)
      PHI(      11)=CXYZ(    10)*SAUX( 2)
      PHI(      12)=CXYZ(     2)*SAUX( 2)
      PHI(      13)=CXYZ(     3)*SAUX( 2)
      PHI(      14)=CXYZ(     4)*SAUX( 2)
      PHI(      15)=CXYZ(    11)*SAUX( 3)
      PHI(      16)=CXYZ(    12)*SAUX( 3)
      PHI(      17)=CXYZ(    13)*SAUX( 3)
      PHI(      18)=CXYZ(    14)*SAUX( 3)
      PHI(      19)=CXYZ(    15)*SAUX( 3)
      PHI(      20)=CXYZ(    16)*SAUX( 3)
      PHI(      21)=CXYZ(    17)*SAUX( 3)
      PHI(      22)=CXYZ(    18)*SAUX( 3)
      PHI(      23)=CXYZ(    19)*SAUX( 3)
      PHI(      24)=CXYZ(    20)*SAUX( 3)
      PHI(      25)=CXYZ(     1)*SAUX( 2)
      PHI(      26)=CXYZ(     5)*SAUX( 3)
      PHI(      27)=CXYZ(     6)*SAUX( 3)
      PHI(      28)=CXYZ(     7)*SAUX( 3)
      PHI(      29)=CXYZ(     8)*SAUX( 3)
      PHI(      30)=CXYZ(     9)*SAUX( 3)
      PHI(      31)=CXYZ(    10)*SAUX( 3)
      PHI(      32)=CXYZ(    21)*SAUX( 4)
      PHI(      33)=CXYZ(    22)*SAUX( 4)
      PHI(      34)=CXYZ(    23)*SAUX( 4)
      PHI(      35)=CXYZ(    24)*SAUX( 4)
      PHI(      36)=CXYZ(    25)*SAUX( 4)
      PHI(      37)=CXYZ(    26)*SAUX( 4)
      PHI(      38)=CXYZ(    27)*SAUX( 4)
      PHI(      39)=CXYZ(    28)*SAUX( 4)
      PHI(      40)=CXYZ(    29)*SAUX( 4)
      PHI(      41)=CXYZ(    30)*SAUX( 4)
      PHI(      42)=CXYZ(    31)*SAUX( 4)
      PHI(      43)=CXYZ(    32)*SAUX( 4)
      PHI(      44)=CXYZ(    33)*SAUX( 4)
      PHI(      45)=CXYZ(    34)*SAUX( 4)
      PHI(      46)=CXYZ(    35)*SAUX( 4)
      FAC=RXB(IJ)
      SF=ONE
      SFAC(1)=SF
      FAC=SFAC( 1)
      PHI(      47)=PHI(      47)+FAC*PHI(       1)
      FAC=SFAC( 1)
      PHI(      55)=PHI(      55)+FAC*PHI(       2)
      PHI(      56)=PHI(      56)+FAC*PHI(       3)
      PHI(      57)=PHI(      57)+FAC*PHI(       4)
      FAC=SFAC( 1)
      J=       5
      DO I=      73,      79
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=      12
      DO I=     115,     127
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=      25
      DO I=     167,     188
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      END IF
      END DO
      FAC=EXJ(JJ)*TWO
      SF=ONE
      SFAC(1)=SF
      FAC=SFAC( 1)
      PHI(      48)=PHI(      48)+FAC*PHI(      47)
      FAC=SFAC( 1)
      PHI(      58)=PHI(      58)+FAC*PHI(      55)
      PHI(      59)=PHI(      59)+FAC*PHI(      56)
      PHI(      60)=PHI(      60)+FAC*PHI(      57)
      FAC=SFAC( 1)
      J=      73
      DO I=      80,      86
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     115
      DO I=     128,     140
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     167
      DO I=     189,     210
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      END DO
      FAC=RXK(KL)
      SF=ONE
      SFAC(1)=SF
      SF=SF*FAC
      SFAC( 2)=SF
      SF=SF*FAC
      SFAC( 3)=SF
      SF=SF*FAC
      SFAC( 4)=SF
      SF=SF*FAC
      SFAC( 5)=SF
      FAC=SFAC( 3)
      PHI(      49)=PHI(      49)+FAC*PHI(      48)
      FAC=SFAC( 4)
      PHI(      51)=PHI(      51)+FAC*PHI(      48)
      FAC=SFAC( 5)
      PHI(      53)=PHI(      53)+FAC*PHI(      48)
      FAC=SFAC( 4)
      PHI(      61)=PHI(      61)+FAC*PHI(      58)
      PHI(      62)=PHI(      62)+FAC*PHI(      59)
      PHI(      63)=PHI(      63)+FAC*PHI(      60)
      FAC=SFAC( 5)
      PHI(      67)=PHI(      67)+FAC*PHI(      58)
      PHI(      68)=PHI(      68)+FAC*PHI(      59)
      PHI(      69)=PHI(      69)+FAC*PHI(      60)
      FAC=SFAC( 4)
      J=      80
      DO I=      87,      93
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 5)
      J=      80
      DO I=     101,     107
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 5)
      J=     128
      DO I=     141,     153
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 5)
      J=     189
      DO I=     211,     232
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      END IF
      END DO
      FAC=EXL(LL)*TWO
      SF=ONE
      SFAC(1)=SF
      SF=SF*FAC
      SFAC( 2)=SF
      SF=SF*FAC
      SFAC( 3)=SF
      SF=SF*FAC
      SFAC( 4)=SF
      SF=SF*FAC
      SFAC( 5)=SF
      FAC=SFAC( 1)
      PHI(      50)=PHI(      50)+FAC*PHI(      49)
      FAC=SFAC( 3)
      PHI(      52)=PHI(      52)+FAC*PHI(      51)
      FAC=SFAC( 5)
      PHI(      54)=PHI(      54)+FAC*PHI(      53)
      FAC=SFAC( 2)
      PHI(      64)=PHI(      64)+FAC*PHI(      61)
      PHI(      65)=PHI(      65)+FAC*PHI(      62)
      PHI(      66)=PHI(      66)+FAC*PHI(      63)
      FAC=SFAC( 4)
      PHI(      70)=PHI(      70)+FAC*PHI(      67)
      PHI(      71)=PHI(      71)+FAC*PHI(      68)
      PHI(      72)=PHI(      72)+FAC*PHI(      69)
      FAC=SFAC( 1)
      J=      87
      DO I=      94,     100
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 3)
      J=     101
      DO I=     108,     114
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 2)
      J=     141
      DO I=     154,     166
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     211
      DO I=     233,     254
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      END DO
C  POST-CONTRACTION PHASE
      PHI(      95)=PHI(      95)-D01*PHI(      94)
      PHI(      97)=PHI(      97)-D01*PHI(      94)
      PHI(     100)=PHI(     100)-D01*PHI(      94)
      PHI(     109)=PHI(     109)-D01*PHI(     108)
      PHI(     111)=PHI(     111)-D01*PHI(     108)
      PHI(     114)=PHI(     114)-D01*PHI(     108)
      PHI(     157)=PHI(     157)-D01*PHI(     154)
      PHI(     158)=PHI(     158)-D01*PHI(     155)
      PHI(     160)=PHI(     160)-D01*PHI(     155)
      PHI(     161)=PHI(     161)-D01*PHI(     156)
      PHI(     163)=PHI(     163)-D01*PHI(     156)
      PHI(     166)=PHI(     166)-D01*PHI(     156)
      PHI(     157)=PHI(     157)-D02*PHI(     154)
      PHI(     159)=PHI(     159)-D01*PHI(     154)
      PHI(     160)=PHI(     160)-D02*PHI(     155)
      PHI(     164)=PHI(     164)-D01*PHI(     154)
      PHI(     165)=PHI(     165)-D01*PHI(     155)
      PHI(     166)=PHI(     166)-D02*PHI(     156)
      PHI(     240)=PHI(     240)-D01*PHI(     234)
      PHI(     241)=PHI(     241)-D01*PHI(     235)
      PHI(     242)=PHI(     242)-D01*PHI(     236)
      PHI(     244)=PHI(     244)-D01*PHI(     236)
      PHI(     245)=PHI(     245)-D01*PHI(     237)
      PHI(     246)=PHI(     246)-D01*PHI(     238)
      PHI(     248)=PHI(     248)-D01*PHI(     238)
      PHI(     249)=PHI(     249)-D01*PHI(     239)
      PHI(     251)=PHI(     251)-D01*PHI(     239)
      PHI(     254)=PHI(     254)-D01*PHI(     239)
      PHI(     240)=PHI(     240)-D02*PHI(     234)
      PHI(     241)=PHI(     241)-D02*PHI(     235)
      PHI(     243)=PHI(     243)-D01*PHI(     235)
      PHI(     244)=PHI(     244)-D02*PHI(     236)
      PHI(     245)=PHI(     245)-D02*PHI(     237)
      PHI(     247)=PHI(     247)-D01*PHI(     237)
      PHI(     248)=PHI(     248)-D02*PHI(     238)
      PHI(     252)=PHI(     252)-D01*PHI(     237)
      PHI(     253)=PHI(     253)-D01*PHI(     238)
      PHI(     254)=PHI(     254)-D02*PHI(     239)
      PHI(     234)=PHI(     234)-D01*PHI(     233)
      PHI(     236)=PHI(     236)-D01*PHI(     233)
      PHI(     239)=PHI(     239)-D01*PHI(     233)
      PHI(     240)=PHI(     240)-D03*PHI(     234)
      PHI(     242)=PHI(     242)-D01*PHI(     234)
      PHI(     243)=PHI(     243)-D02*PHI(     235)
      PHI(     244)=PHI(     244)-D03*PHI(     236)
      PHI(     249)=PHI(     249)-D01*PHI(     234)
      PHI(     250)=PHI(     250)-D01*PHI(     235)
      PHI(     251)=PHI(     251)-D01*PHI(     236)
      PHI(     252)=PHI(     252)-D02*PHI(     237)
      PHI(     253)=PHI(     253)-D02*PHI(     238)
      PHI(     254)=PHI(     254)-D03*PHI(     239)
      WRK2(    1,1)= PHI(      54)
      WRK2(    2,1)= PHI(      52)
      WRK2(    3,1)= PHI(      50)
      WRK2(    4,1)=-PHI(      70)
      WRK2(    5,1)=-PHI(      71)
      WRK2(    6,1)=-PHI(      72)
      WRK2(    7,1)=-PHI(      64)
      WRK2(    8,1)=-PHI(      65)
      WRK2(    9,1)=-PHI(      66)
      WRK2(   10,1)= PHI(     109)
      WRK2(   11,1)= PHI(     110)
      WRK2(   12,1)= PHI(     111)
      WRK2(   13,1)= PHI(     112)
      WRK2(   14,1)= PHI(     113)
      WRK2(   15,1)= PHI(     114)
      WRK2(   16,1)= PHI(      95)
      WRK2(   17,1)= PHI(      96)
      WRK2(   18,1)= PHI(      97)
      WRK2(   19,1)= PHI(      98)
      WRK2(   20,1)= PHI(      99)
      WRK2(   21,1)= PHI(     100)
      WRK2(   22,1)=-PHI(     157)
      WRK2(   23,1)=-PHI(     158)
      WRK2(   24,1)=-PHI(     159)
      WRK2(   25,1)=-PHI(     160)
      WRK2(   26,1)=-PHI(     161)
      WRK2(   27,1)=-PHI(     162)
      WRK2(   28,1)=-PHI(     163)
      WRK2(   29,1)=-PHI(     164)
      WRK2(   30,1)=-PHI(     165)
      WRK2(   31,1)=-PHI(     166)
      WRK2(   32,1)= PHI(     240)
      WRK2(   33,1)= PHI(     241)
      WRK2(   34,1)= PHI(     242)
      WRK2(   35,1)= PHI(     243)
      WRK2(   36,1)= PHI(     244)
      WRK2(   37,1)= PHI(     245)
      WRK2(   38,1)= PHI(     246)
      WRK2(   39,1)= PHI(     247)
      WRK2(   40,1)= PHI(     248)
      WRK2(   41,1)= PHI(     249)
      WRK2(   42,1)= PHI(     250)
      WRK2(   43,1)= PHI(     251)
      WRK2(   44,1)= PHI(     252)
      WRK2(   45,1)= PHI(     253)
      WRK2(   46,1)= PHI(     254)
      CONF(1)=XD-XC
      CONF(2)=YD-YC
      CONF(3)=ZD-ZC
      DO I=1,     1
      WRK2(   32,I)=WRK2(   32,I)+CONF( 1)*WRK2(   22,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 6)*WRK2(   16,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 1)*WRK2(   23,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 5)*WRK2(   17,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 1)*WRK2(   24,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 4)*WRK2(   18,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 1)*WRK2(   25,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 2)*WRK2(   25,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 6)*WRK2(   18,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 1)*WRK2(   26,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 5)*WRK2(   19,I)
      WRK2(   38,I)=WRK2(   38,I)+CONF( 1)*WRK2(   27,I)
      WRK2(   38,I)=WRK2(   38,I)+CONF( 4)*WRK2(   20,I)
      WRK2(   39,I)=WRK2(   39,I)+CONF( 1)*WRK2(   28,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 2)*WRK2(   28,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 5)*WRK2(   20,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 1)*WRK2(   29,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 4)*WRK2(   21,I)
      WRK2(   42,I)=WRK2(   42,I)+CONF( 1)*WRK2(   30,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 2)*WRK2(   30,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 4)*WRK2(   21,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 1)*WRK2(   31,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 2)*WRK2(   31,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 3)*WRK2(   31,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 6)*WRK2(   21,I)
      WRK2(   22,I)=WRK2(   22,I)+CONF( 1)*WRK2(   10,I)
      WRK2(   22,I)=WRK2(   22,I)+CONF( 5)*WRK2(    7,I)
      WRK2(   23,I)=WRK2(   23,I)+CONF( 1)*WRK2(   11,I)
      WRK2(   23,I)=WRK2(   23,I)+CONF( 4)*WRK2(    8,I)
      WRK2(   24,I)=WRK2(   24,I)+CONF( 1)*WRK2(   12,I)
      WRK2(   25,I)=WRK2(   25,I)+CONF( 2)*WRK2(   12,I)
      WRK2(   25,I)=WRK2(   25,I)+CONF( 5)*WRK2(    8,I)
      WRK2(   26,I)=WRK2(   26,I)+CONF( 1)*WRK2(   13,I)
      WRK2(   26,I)=WRK2(   26,I)+CONF( 4)*WRK2(    9,I)
      WRK2(   27,I)=WRK2(   27,I)+CONF( 1)*WRK2(   14,I)
      WRK2(   28,I)=WRK2(   28,I)+CONF( 2)*WRK2(   14,I)
      WRK2(   28,I)=WRK2(   28,I)+CONF( 4)*WRK2(    9,I)
      WRK2(   29,I)=WRK2(   29,I)+CONF( 1)*WRK2(   15,I)
      WRK2(   30,I)=WRK2(   30,I)+CONF( 2)*WRK2(   15,I)
      WRK2(   31,I)=WRK2(   31,I)+CONF( 3)*WRK2(   15,I)
      WRK2(   31,I)=WRK2(   31,I)+CONF( 5)*WRK2(    9,I)
      WRK2(   10,I)=WRK2(   10,I)+CONF( 1)*WRK2(    4,I)
      WRK2(   10,I)=WRK2(   10,I)+CONF( 4)*WRK2(    2,I)
      WRK2(   11,I)=WRK2(   11,I)+CONF( 1)*WRK2(    5,I)
      WRK2(   12,I)=WRK2(   12,I)+CONF( 2)*WRK2(    5,I)
      WRK2(   12,I)=WRK2(   12,I)+CONF( 4)*WRK2(    2,I)
      WRK2(   13,I)=WRK2(   13,I)+CONF( 1)*WRK2(    6,I)
      WRK2(   14,I)=WRK2(   14,I)+CONF( 2)*WRK2(    6,I)
      WRK2(   15,I)=WRK2(   15,I)+CONF( 3)*WRK2(    6,I)
      WRK2(   15,I)=WRK2(   15,I)+CONF( 4)*WRK2(    2,I)
      WRK2(   16,I)=WRK2(   16,I)+CONF( 1)*WRK2(    7,I)
      WRK2(   16,I)=WRK2(   16,I)+CONF( 4)*WRK2(    3,I)
      WRK2(   17,I)=WRK2(   17,I)+CONF( 1)*WRK2(    8,I)
      WRK2(   18,I)=WRK2(   18,I)+CONF( 2)*WRK2(    8,I)
      WRK2(   18,I)=WRK2(   18,I)+CONF( 4)*WRK2(    3,I)
      WRK2(   19,I)=WRK2(   19,I)+CONF( 1)*WRK2(    9,I)
      WRK2(   20,I)=WRK2(   20,I)+CONF( 2)*WRK2(    9,I)
      WRK2(   21,I)=WRK2(   21,I)+CONF( 3)*WRK2(    9,I)
      WRK2(   21,I)=WRK2(   21,I)+CONF( 4)*WRK2(    3,I)
      WRK2(    4,I)=WRK2(    4,I)+CONF( 1)*WRK2(    1,I)
      WRK2(    5,I)=WRK2(    5,I)+CONF( 2)*WRK2(    1,I)
      WRK2(    6,I)=WRK2(    6,I)+CONF( 3)*WRK2(    1,I)
      WRK2(    7,I)=WRK2(    7,I)+CONF( 1)*WRK2(    2,I)
      WRK2(    8,I)=WRK2(    8,I)+CONF( 2)*WRK2(    2,I)
      WRK2(    9,I)=WRK2(    9,I)+CONF( 3)*WRK2(    2,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 1)*WRK2(   22,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 5)*WRK2(   16,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 1)*WRK2(   23,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 4)*WRK2(   17,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 1)*WRK2(   24,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 2)*WRK2(   24,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 5)*WRK2(   17,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 2)*WRK2(   25,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 5)*WRK2(   18,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 1)*WRK2(   26,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 4)*WRK2(   19,I)
      WRK2(   38,I)=WRK2(   38,I)+CONF( 1)*WRK2(   27,I)
      WRK2(   39,I)=WRK2(   39,I)+CONF( 2)*WRK2(   27,I)
      WRK2(   39,I)=WRK2(   39,I)+CONF( 4)*WRK2(   19,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 2)*WRK2(   28,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 4)*WRK2(   20,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 1)*WRK2(   29,I)
      WRK2(   42,I)=WRK2(   42,I)+CONF( 2)*WRK2(   29,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 2)*WRK2(   30,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 3)*WRK2(   29,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 5)*WRK2(   19,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 3)*WRK2(   30,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 5)*WRK2(   20,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 3)*WRK2(   31,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 5)*WRK2(   21,I)
      WRK2(   22,I)=WRK2(   22,I)+CONF( 1)*WRK2(   10,I)
      WRK2(   22,I)=WRK2(   22,I)+CONF( 4)*WRK2(    7,I)
      WRK2(   23,I)=WRK2(   23,I)+CONF( 1)*WRK2(   11,I)
      WRK2(   24,I)=WRK2(   24,I)+CONF( 2)*WRK2(   11,I)
      WRK2(   24,I)=WRK2(   24,I)+CONF( 4)*WRK2(    7,I)
      WRK2(   25,I)=WRK2(   25,I)+CONF( 2)*WRK2(   12,I)
      WRK2(   25,I)=WRK2(   25,I)+CONF( 4)*WRK2(    8,I)
      WRK2(   26,I)=WRK2(   26,I)+CONF( 1)*WRK2(   13,I)
      WRK2(   27,I)=WRK2(   27,I)+CONF( 2)*WRK2(   13,I)
      WRK2(   28,I)=WRK2(   28,I)+CONF( 2)*WRK2(   14,I)
      WRK2(   29,I)=WRK2(   29,I)+CONF( 3)*WRK2(   13,I)
      WRK2(   29,I)=WRK2(   29,I)+CONF( 4)*WRK2(    7,I)
      WRK2(   30,I)=WRK2(   30,I)+CONF( 3)*WRK2(   14,I)
      WRK2(   30,I)=WRK2(   30,I)+CONF( 4)*WRK2(    8,I)
      WRK2(   31,I)=WRK2(   31,I)+CONF( 3)*WRK2(   15,I)
      WRK2(   31,I)=WRK2(   31,I)+CONF( 4)*WRK2(    9,I)
      WRK2(   10,I)=WRK2(   10,I)+CONF( 1)*WRK2(    4,I)
      WRK2(   11,I)=WRK2(   11,I)+CONF( 2)*WRK2(    4,I)
      WRK2(   12,I)=WRK2(   12,I)+CONF( 2)*WRK2(    5,I)
      WRK2(   13,I)=WRK2(   13,I)+CONF( 3)*WRK2(    4,I)
      WRK2(   14,I)=WRK2(   14,I)+CONF( 3)*WRK2(    5,I)
      WRK2(   15,I)=WRK2(   15,I)+CONF( 3)*WRK2(    6,I)
      WRK2(   16,I)=WRK2(   16,I)+CONF( 1)*WRK2(    7,I)
      WRK2(   17,I)=WRK2(   17,I)+CONF( 2)*WRK2(    7,I)
      WRK2(   18,I)=WRK2(   18,I)+CONF( 2)*WRK2(    8,I)
      WRK2(   19,I)=WRK2(   19,I)+CONF( 3)*WRK2(    7,I)
      WRK2(   20,I)=WRK2(   20,I)+CONF( 3)*WRK2(    8,I)
      WRK2(   21,I)=WRK2(   21,I)+CONF( 3)*WRK2(    9,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 1)*WRK2(   22,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 4)*WRK2(   16,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 1)*WRK2(   23,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 2)*WRK2(   23,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 4)*WRK2(   16,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 2)*WRK2(   24,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 4)*WRK2(   17,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 2)*WRK2(   25,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 4)*WRK2(   18,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 1)*WRK2(   26,I)
      WRK2(   38,I)=WRK2(   38,I)+CONF( 2)*WRK2(   26,I)
      WRK2(   39,I)=WRK2(   39,I)+CONF( 2)*WRK2(   27,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 2)*WRK2(   28,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 3)*WRK2(   26,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 4)*WRK2(   16,I)
      WRK2(   42,I)=WRK2(   42,I)+CONF( 3)*WRK2(   27,I)
      WRK2(   42,I)=WRK2(   42,I)+CONF( 4)*WRK2(   17,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 3)*WRK2(   28,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 4)*WRK2(   18,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 3)*WRK2(   29,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 4)*WRK2(   19,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 3)*WRK2(   30,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 4)*WRK2(   20,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 3)*WRK2(   31,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 4)*WRK2(   21,I)
      WRK2(   22,I)=WRK2(   22,I)+CONF( 1)*WRK2(   10,I)
      WRK2(   23,I)=WRK2(   23,I)+CONF( 2)*WRK2(   10,I)
      WRK2(   24,I)=WRK2(   24,I)+CONF( 2)*WRK2(   11,I)
      WRK2(   25,I)=WRK2(   25,I)+CONF( 2)*WRK2(   12,I)
      WRK2(   26,I)=WRK2(   26,I)+CONF( 3)*WRK2(   10,I)
      WRK2(   27,I)=WRK2(   27,I)+CONF( 3)*WRK2(   11,I)
      WRK2(   28,I)=WRK2(   28,I)+CONF( 3)*WRK2(   12,I)
      WRK2(   29,I)=WRK2(   29,I)+CONF( 3)*WRK2(   13,I)
      WRK2(   30,I)=WRK2(   30,I)+CONF( 3)*WRK2(   14,I)
      WRK2(   31,I)=WRK2(   31,I)+CONF( 3)*WRK2(   15,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 1)*WRK2(   22,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 2)*WRK2(   22,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 2)*WRK2(   23,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 2)*WRK2(   24,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 2)*WRK2(   25,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 3)*WRK2(   22,I)
      WRK2(   38,I)=WRK2(   38,I)+CONF( 3)*WRK2(   23,I)
      WRK2(   39,I)=WRK2(   39,I)+CONF( 3)*WRK2(   24,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 3)*WRK2(   25,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 3)*WRK2(   26,I)
      WRK2(   42,I)=WRK2(   42,I)+CONF( 3)*WRK2(   27,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 3)*WRK2(   28,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 3)*WRK2(   29,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 3)*WRK2(   30,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 3)*WRK2(   31,I)
      END DO
      RETURN
      END
C  LPHI=     612
C  LWK1=       4
C  LWK2=     138
C  LENW=      46
C*MODULE INT2G   *DECK GSPS
      SUBROUTINE GSPS (IPRIM,JPRIM,KPRIM,LPRIM,IEQJ
     *,KEQL,XA,YA,ZA,XB,YB,ZB,XC,YC,ZC,XD,YD,ZD,PHI
     *,WRK1,WRK2,LENW)
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      COMMON /SHLNOS/ QQ4,IDUMMY(20)
      PARAMETER (NTX=4,NPF=450,NGRD=7,NPX=1000,MXQT=16)
      COMMON /FMTTBL/ FGRID(0:NTX,0:NPF,0:NGRD)
     *,XGRID(0:NTX,0:NPX),TMAX,RFINC(0:NGRD)
     *,RXINC,RMR(MXQT),TLGM(0:MXQT),NORD
      PARAMETER (MXPR=30, MXPS=MXPR*MXPR)
      COMMON /ERIPRM/ EXI(MXPR),EXJ(MXPR),EXK(MXPR),EXL(MXPR)
     *,CCI(MXPR),CCJ(MXPR),CCK(MXPR),CCL(MXPR)
     *,XAB(MXPS),YAB(MXPS),ZAB(MXPS)
     *,XCD(MXPS),YCD(MXPS),ZCD(MXPS)
     *,CCBRA(MXPS),CCKET(MXPS),RXB(MXPS)
     *,SLBRA(MXPS),SLKET(MXPS),RXK(MXPS)
      COMMON /NLRC  / LCFLAG,LRINT,EMU,EMU2,LRFILE
      LOGICAL   LCFLAG,LRINT
      LOGICAL IEQJ,KEQL
      DOUBLE PRECISION RXYZ(3),SFAC(9),CONF(11),FT(0:16)
      DOUBLE PRECISION SAUX(0: 5),CXYZ( 56)
      DOUBLE PRECISION PHI(*),WRK1(*),WRK2(LENW,*)
      PARAMETER (ZERO=0.0D+00, ONE=1.0D+00, TWO=2.0D+00)
      PARAMETER (D01=1.0D+00, D02=2.0D+00, D03=3.0D+00)
      PARAMETER (D04=4.0D+00)
      PARAMETER (CCTOL=1.0D-13)
      DATA CONF/0.0D+00,0.0D+00,0.0D+00,1.0D+00,2.0D+00,3.0D+00
     &,         4.0D+00,5.0D+00,6.0D+00,7.0D+00,8.0D+00/
      SAVE CONF
      RFAC=RFINC( 5)
      KL=0
      PHI(      84)=ZERO
      PHI(      86)=ZERO
      PHI(      88)=ZERO
      PHI(      98)=ZERO
      PHI(      99)=ZERO
      PHI(     100)=ZERO
      PHI(     104)=ZERO
      PHI(     105)=ZERO
      PHI(     106)=ZERO
      PHI(     110)=ZERO
      PHI(     111)=ZERO
      PHI(     112)=ZERO
      PHI(     119)=ZERO
      PHI(     120)=ZERO
      PHI(     121)=ZERO
      PHI(     125)=ZERO
      PHI(     126)=ZERO
      PHI(     127)=ZERO
      DO I=     149,     155
      PHI(I)=ZERO
      END DO
      DO I=     163,     169
      PHI(I)=ZERO
      END DO
      DO I=     184,     190
      PHI(I)=ZERO
      END DO
      DO I=     198,     204
      PHI(I)=ZERO
      END DO
      DO I=     244,     256
      PHI(I)=ZERO
      END DO
      DO I=     270,     282
      PHI(I)=ZERO
      END DO
      DO I=     309,     321
      PHI(I)=ZERO
      END DO
      DO I=     388,     409
      PHI(I)=ZERO
      END DO
      DO I=     454,     475
      PHI(I)=ZERO
      END DO
      DO I=     578,     611
      PHI(I)=ZERO
      END DO
      DO LL=1,LPRIM
      PHI(      83)=ZERO
      PHI(      85)=ZERO
      PHI(      87)=ZERO
      PHI(      95)=ZERO
      PHI(      96)=ZERO
      PHI(      97)=ZERO
      PHI(     101)=ZERO
      PHI(     102)=ZERO
      PHI(     103)=ZERO
      PHI(     107)=ZERO
      PHI(     108)=ZERO
      PHI(     109)=ZERO
      PHI(     116)=ZERO
      PHI(     117)=ZERO
      PHI(     118)=ZERO
      PHI(     122)=ZERO
      PHI(     123)=ZERO
      PHI(     124)=ZERO
      DO I=     142,     148
      PHI(I)=ZERO
      END DO
      DO I=     156,     162
      PHI(I)=ZERO
      END DO
      DO I=     177,     183
      PHI(I)=ZERO
      END DO
      DO I=     191,     197
      PHI(I)=ZERO
      END DO
      DO I=     231,     243
      PHI(I)=ZERO
      END DO
      DO I=     257,     269
      PHI(I)=ZERO
      END DO
      DO I=     296,     308
      PHI(I)=ZERO
      END DO
      DO I=     366,     387
      PHI(I)=ZERO
      END DO
      DO I=     432,     453
      PHI(I)=ZERO
      END DO
      DO I=     544,     577
      PHI(I)=ZERO
      END DO
      KTOP=KPRIM
      IF (KEQL) KTOP=LL
      DO KK=1,KTOP
      KL=KL+1
      CCKF=CCKET(KL)
      IF (ABS(CCKF).GT.CCTOL) THEN
      EKL=EXL(LL)+EXK(KK)
      XKL=XCD(KL)
      YKL=YCD(KL)
      ZKL=ZCD(KL)
      IJ=0
      PHI(      82)=ZERO
      PHI(      92)=ZERO
      PHI(      93)=ZERO
      PHI(      94)=ZERO
      PHI(     113)=ZERO
      PHI(     114)=ZERO
      PHI(     115)=ZERO
      DO I=     135,     141
      PHI(I)=ZERO
      END DO
      DO I=     170,     176
      PHI(I)=ZERO
      END DO
      DO I=     218,     230
      PHI(I)=ZERO
      END DO
      DO I=     283,     295
      PHI(I)=ZERO
      END DO
      DO I=     344,     365
      PHI(I)=ZERO
      END DO
      DO I=     410,     431
      PHI(I)=ZERO
      END DO
      DO I=     510,     543
      PHI(I)=ZERO
      END DO
      DO JJ=1,JPRIM
      PHI(      81)=ZERO
      PHI(      89)=ZERO
      PHI(      90)=ZERO
      PHI(      91)=ZERO
      DO I=     128,     134
      PHI(I)=ZERO
      END DO
      DO I=     205,     217
      PHI(I)=ZERO
      END DO
      DO I=     322,     343
      PHI(I)=ZERO
      END DO
      DO I=     476,     509
      PHI(I)=ZERO
      END DO
      ITOP=IPRIM
      IF (IEQJ) ITOP=JJ
      DO II=1,ITOP
      IJ=IJ+1
      CCBF=CCBRA(IJ)
      IF (ABS(CCBF).GT.CCTOL) THEN
      EIJ=EXJ(JJ)+EXI(II)
      EIJKL=EIJ+EKL
      EINV=ONE/EIJKL
      RX=XKL-XAB(IJ)
      RY=YKL-YAB(IJ)
      RZ=ZKL-ZAB(IJ)
      R2=RX*RX+RY*RY+RZ*RZ
      VS=EIJ*EKL*EINV
      IF(LRINT) THEN
         VS0=VS
         VS=VS0*EMU2/(VS0+EMU2)
      ENDIF
      TT=VS*R2
      IF (TT.LE.TMAX) THEN
      TVAL=TT*RFAC
      IPT=NINT(TVAL)
      FX= TVAL*FGRID(4,IPT, 5)
      FX=(FX + FGRID(3,IPT, 5))*TVAL
      FX=(FX + FGRID(2,IPT, 5))*TVAL
      FX=(FX + FGRID(1,IPT, 5))*TVAL
      FX= FX + FGRID(0,IPT, 5)
      FT( 8)=FX
      TVAL=TT*RXINC
      IPT=NINT(TVAL)
      ET= TVAL*XGRID(4,IPT)
      ET=(ET + XGRID(3,IPT))*TVAL
      ET=(ET + XGRID(2,IPT))*TVAL
      ET=(ET + XGRID(1,IPT))*TVAL
      ET= ET + XGRID(0,IPT)
      T2=TT*TWO
      FT( 7)=RMR( 8)*(T2*FT( 8)+ET)
      FT( 6)=RMR( 7)*(T2*FT( 7)+ET)
      FT( 5)=RMR( 6)*(T2*FT( 6)+ET)
      FT( 4)=RMR( 5)*(T2*FT( 5)+ET)
      FT( 3)=RMR( 4)*(T2*FT( 4)+ET)
      FT( 2)=RMR( 3)*(T2*FT( 3)+ET)
      FT( 1)=RMR( 2)*(T2*FT( 2)+ET)
      FT( 0)=RMR( 1)*(T2*FT( 1)+ET)
      ZEROL=CCBF*CCKF*SQRT(EINV)*QQ4
      IF(LRINT)ZEROL=ZEROL*SQRT(VS/VS0)
      SAUX(0)=ZEROL*FT(0)
      VS2=TWO*VS
      VSM=ONE
      VSM=VSM*VS2
      SAUX( 1)=ZEROL*FT( 1)*VSM
      VSM=VSM*VS2
      SAUX( 2)=ZEROL*FT( 2)*VSM
      VSM=VSM*VS2
      SAUX( 3)=ZEROL*FT( 3)*VSM
      VSM=VSM*VS2
      SAUX( 4)=ZEROL*FT( 4)*VSM
      VSM=VSM*VS2
      SAUX( 5)=ZEROL*FT( 5)*VSM
      ELSE
      ZEROL=CCBF*CCKF*SLBRA(IJ)*SLKET(KL)*QQ4
      RINV=ONE/R2
      SAUX(0)=ZEROL*SQRT(RINV)
      RR=SAUX(0)
      RR=RR*RINV
      SAUX( 1)=TLGM( 1)*RR
      RR=RR*RINV
      SAUX( 2)=TLGM( 2)*RR
      RR=RR*RINV
      SAUX( 3)=TLGM( 3)*RR
      RR=RR*RINV
      SAUX( 4)=TLGM( 4)*RR
      RR=RR*RINV
      SAUX( 5)=TLGM( 5)*RR
      END IF
      RXYZ(1)=RX
      RXYZ(2)=RY
      RXYZ(3)=RZ
      CXYZ(1)=ONE
      CXYZ(2)=RX
      CXYZ(3)=RY
      CXYZ(4)=RZ
      CXYZ(     5)=CXYZ(     2)*RXYZ(1)
      CXYZ(     6)=CXYZ(     3)*RXYZ(1)
      CXYZ(     7)=CXYZ(     3)*RXYZ(2)
      CXYZ(     8)=CXYZ(     4)*RXYZ(1)
      CXYZ(     9)=CXYZ(     4)*RXYZ(2)
      CXYZ(    10)=CXYZ(     4)*RXYZ(3)
      CXYZ(    11)=CXYZ(     5)*RXYZ(1)
      CXYZ(    12)=CXYZ(     6)*RXYZ(1)
      CXYZ(    13)=CXYZ(     7)*RXYZ(1)
      CXYZ(    14)=CXYZ(     7)*RXYZ(2)
      CXYZ(    15)=CXYZ(     8)*RXYZ(1)
      CXYZ(    16)=CXYZ(     9)*RXYZ(1)
      CXYZ(    17)=CXYZ(     9)*RXYZ(2)
      CXYZ(    18)=CXYZ(    10)*RXYZ(1)
      CXYZ(    19)=CXYZ(    10)*RXYZ(2)
      CXYZ(    20)=CXYZ(    10)*RXYZ(3)
      CXYZ(    21)=CXYZ(    11)*RXYZ(1)
      CXYZ(    22)=CXYZ(    12)*RXYZ(1)
      CXYZ(    23)=CXYZ(    13)*RXYZ(1)
      CXYZ(    24)=CXYZ(    14)*RXYZ(1)
      CXYZ(    25)=CXYZ(    14)*RXYZ(2)
      CXYZ(    26)=CXYZ(    15)*RXYZ(1)
      CXYZ(    27)=CXYZ(    16)*RXYZ(1)
      CXYZ(    28)=CXYZ(    17)*RXYZ(1)
      CXYZ(    29)=CXYZ(    17)*RXYZ(2)
      CXYZ(    30)=CXYZ(    18)*RXYZ(1)
      CXYZ(    31)=CXYZ(    19)*RXYZ(1)
      CXYZ(    32)=CXYZ(    19)*RXYZ(2)
      CXYZ(    33)=CXYZ(    20)*RXYZ(1)
      CXYZ(    34)=CXYZ(    20)*RXYZ(2)
      CXYZ(    35)=CXYZ(    20)*RXYZ(3)
      CXYZ(    36)=CXYZ(    21)*RXYZ(1)
      CXYZ(    37)=CXYZ(    22)*RXYZ(1)
      CXYZ(    38)=CXYZ(    23)*RXYZ(1)
      CXYZ(    39)=CXYZ(    24)*RXYZ(1)
      CXYZ(    40)=CXYZ(    25)*RXYZ(1)
      CXYZ(    41)=CXYZ(    25)*RXYZ(2)
      CXYZ(    42)=CXYZ(    26)*RXYZ(1)
      CXYZ(    43)=CXYZ(    27)*RXYZ(1)
      CXYZ(    44)=CXYZ(    28)*RXYZ(1)
      CXYZ(    45)=CXYZ(    29)*RXYZ(1)
      CXYZ(    46)=CXYZ(    29)*RXYZ(2)
      CXYZ(    47)=CXYZ(    30)*RXYZ(1)
      CXYZ(    48)=CXYZ(    31)*RXYZ(1)
      CXYZ(    49)=CXYZ(    32)*RXYZ(1)
      CXYZ(    50)=CXYZ(    32)*RXYZ(2)
      CXYZ(    51)=CXYZ(    33)*RXYZ(1)
      CXYZ(    52)=CXYZ(    34)*RXYZ(1)
      CXYZ(    53)=CXYZ(    34)*RXYZ(2)
      CXYZ(    54)=CXYZ(    35)*RXYZ(1)
      CXYZ(    55)=CXYZ(    35)*RXYZ(2)
      CXYZ(    56)=CXYZ(    35)*RXYZ(3)
      PHI(1)=SAUX(0)
      PHI(       2)=CXYZ(     2)*SAUX( 1)
      PHI(       3)=CXYZ(     3)*SAUX( 1)
      PHI(       4)=CXYZ(     4)*SAUX( 1)
      PHI(       5)=CXYZ(     1)*SAUX( 1)
      PHI(       6)=CXYZ(     5)*SAUX( 2)
      PHI(       7)=CXYZ(     6)*SAUX( 2)
      PHI(       8)=CXYZ(     7)*SAUX( 2)
      PHI(       9)=CXYZ(     8)*SAUX( 2)
      PHI(      10)=CXYZ(     9)*SAUX( 2)
      PHI(      11)=CXYZ(    10)*SAUX( 2)
      PHI(      12)=CXYZ(     2)*SAUX( 2)
      PHI(      13)=CXYZ(     3)*SAUX( 2)
      PHI(      14)=CXYZ(     4)*SAUX( 2)
      PHI(      15)=CXYZ(    11)*SAUX( 3)
      PHI(      16)=CXYZ(    12)*SAUX( 3)
      PHI(      17)=CXYZ(    13)*SAUX( 3)
      PHI(      18)=CXYZ(    14)*SAUX( 3)
      PHI(      19)=CXYZ(    15)*SAUX( 3)
      PHI(      20)=CXYZ(    16)*SAUX( 3)
      PHI(      21)=CXYZ(    17)*SAUX( 3)
      PHI(      22)=CXYZ(    18)*SAUX( 3)
      PHI(      23)=CXYZ(    19)*SAUX( 3)
      PHI(      24)=CXYZ(    20)*SAUX( 3)
      PHI(      25)=CXYZ(     1)*SAUX( 2)
      PHI(      26)=CXYZ(     5)*SAUX( 3)
      PHI(      27)=CXYZ(     6)*SAUX( 3)
      PHI(      28)=CXYZ(     7)*SAUX( 3)
      PHI(      29)=CXYZ(     8)*SAUX( 3)
      PHI(      30)=CXYZ(     9)*SAUX( 3)
      PHI(      31)=CXYZ(    10)*SAUX( 3)
      PHI(      32)=CXYZ(    21)*SAUX( 4)
      PHI(      33)=CXYZ(    22)*SAUX( 4)
      PHI(      34)=CXYZ(    23)*SAUX( 4)
      PHI(      35)=CXYZ(    24)*SAUX( 4)
      PHI(      36)=CXYZ(    25)*SAUX( 4)
      PHI(      37)=CXYZ(    26)*SAUX( 4)
      PHI(      38)=CXYZ(    27)*SAUX( 4)
      PHI(      39)=CXYZ(    28)*SAUX( 4)
      PHI(      40)=CXYZ(    29)*SAUX( 4)
      PHI(      41)=CXYZ(    30)*SAUX( 4)
      PHI(      42)=CXYZ(    31)*SAUX( 4)
      PHI(      43)=CXYZ(    32)*SAUX( 4)
      PHI(      44)=CXYZ(    33)*SAUX( 4)
      PHI(      45)=CXYZ(    34)*SAUX( 4)
      PHI(      46)=CXYZ(    35)*SAUX( 4)
      PHI(      47)=CXYZ(     2)*SAUX( 3)
      PHI(      48)=CXYZ(     3)*SAUX( 3)
      PHI(      49)=CXYZ(     4)*SAUX( 3)
      PHI(      50)=CXYZ(    11)*SAUX( 4)
      PHI(      51)=CXYZ(    12)*SAUX( 4)
      PHI(      52)=CXYZ(    13)*SAUX( 4)
      PHI(      53)=CXYZ(    14)*SAUX( 4)
      PHI(      54)=CXYZ(    15)*SAUX( 4)
      PHI(      55)=CXYZ(    16)*SAUX( 4)
      PHI(      56)=CXYZ(    17)*SAUX( 4)
      PHI(      57)=CXYZ(    18)*SAUX( 4)
      PHI(      58)=CXYZ(    19)*SAUX( 4)
      PHI(      59)=CXYZ(    20)*SAUX( 4)
      PHI(      60)=CXYZ(    36)*SAUX( 5)
      PHI(      61)=CXYZ(    37)*SAUX( 5)
      PHI(      62)=CXYZ(    38)*SAUX( 5)
      PHI(      63)=CXYZ(    39)*SAUX( 5)
      PHI(      64)=CXYZ(    40)*SAUX( 5)
      PHI(      65)=CXYZ(    41)*SAUX( 5)
      PHI(      66)=CXYZ(    42)*SAUX( 5)
      PHI(      67)=CXYZ(    43)*SAUX( 5)
      PHI(      68)=CXYZ(    44)*SAUX( 5)
      PHI(      69)=CXYZ(    45)*SAUX( 5)
      PHI(      70)=CXYZ(    46)*SAUX( 5)
      PHI(      71)=CXYZ(    47)*SAUX( 5)
      PHI(      72)=CXYZ(    48)*SAUX( 5)
      PHI(      73)=CXYZ(    49)*SAUX( 5)
      PHI(      74)=CXYZ(    50)*SAUX( 5)
      PHI(      75)=CXYZ(    51)*SAUX( 5)
      PHI(      76)=CXYZ(    52)*SAUX( 5)
      PHI(      77)=CXYZ(    53)*SAUX( 5)
      PHI(      78)=CXYZ(    54)*SAUX( 5)
      PHI(      79)=CXYZ(    55)*SAUX( 5)
      PHI(      80)=CXYZ(    56)*SAUX( 5)
      FAC=RXB(IJ)
      SF=ONE
      SFAC(1)=SF
      SF=SF*FAC
      SFAC( 2)=SF
      FAC=SFAC( 2)
      PHI(      81)=PHI(      81)+FAC*PHI(       1)
      FAC=SFAC( 2)
      PHI(      89)=PHI(      89)+FAC*PHI(       2)
      PHI(      90)=PHI(      90)+FAC*PHI(       3)
      PHI(      91)=PHI(      91)+FAC*PHI(       4)
      FAC=SFAC( 2)
      J=       5
      DO I=     128,     134
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 2)
      J=      12
      DO I=     205,     217
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 2)
      J=      25
      DO I=     322,     343
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 2)
      J=      47
      DO I=     476,     509
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      END IF
      END DO
      FAC=EXJ(JJ)*TWO
      SF=ONE
      SFAC(1)=SF
      SF=SF*FAC
      SFAC( 2)=SF
      FAC=SFAC( 2)
      PHI(      82)=PHI(      82)+FAC*PHI(      81)
      FAC=SFAC( 1)
      PHI(      92)=PHI(      92)+FAC*PHI(      89)
      PHI(      93)=PHI(      93)+FAC*PHI(      90)
      PHI(      94)=PHI(      94)+FAC*PHI(      91)
      FAC=SFAC( 2)
      PHI(     113)=PHI(     113)+FAC*PHI(      89)
      PHI(     114)=PHI(     114)+FAC*PHI(      90)
      PHI(     115)=PHI(     115)+FAC*PHI(      91)
      FAC=SFAC( 1)
      J=     128
      DO I=     135,     141
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 2)
      J=     128
      DO I=     170,     176
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     205
      DO I=     218,     230
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 2)
      J=     205
      DO I=     283,     295
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     322
      DO I=     344,     365
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 2)
      J=     322
      DO I=     410,     431
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     476
      DO I=     510,     543
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      END DO
      FAC=RXK(KL)
      SF=ONE
      SFAC(1)=SF
      SF=SF*FAC
      SFAC( 2)=SF
      SF=SF*FAC
      SFAC( 3)=SF
      SF=SF*FAC
      SFAC( 4)=SF
      SF=SF*FAC
      SFAC( 5)=SF
      FAC=SFAC( 3)
      PHI(      83)=PHI(      83)+FAC*PHI(      82)
      FAC=SFAC( 4)
      PHI(      85)=PHI(      85)+FAC*PHI(      82)
      FAC=SFAC( 5)
      PHI(      87)=PHI(      87)+FAC*PHI(      82)
      FAC=SFAC( 3)
      PHI(      95)=PHI(      95)+FAC*PHI(      92)
      PHI(      96)=PHI(      96)+FAC*PHI(      93)
      PHI(      97)=PHI(      97)+FAC*PHI(      94)
      FAC=SFAC( 4)
      PHI(     101)=PHI(     101)+FAC*PHI(      92)
      PHI(     102)=PHI(     102)+FAC*PHI(      93)
      PHI(     103)=PHI(     103)+FAC*PHI(      94)
      FAC=SFAC( 5)
      PHI(     107)=PHI(     107)+FAC*PHI(      92)
      PHI(     108)=PHI(     108)+FAC*PHI(      93)
      PHI(     109)=PHI(     109)+FAC*PHI(      94)
      FAC=SFAC( 4)
      PHI(     116)=PHI(     116)+FAC*PHI(     113)
      PHI(     117)=PHI(     117)+FAC*PHI(     114)
      PHI(     118)=PHI(     118)+FAC*PHI(     115)
      FAC=SFAC( 5)
      PHI(     122)=PHI(     122)+FAC*PHI(     113)
      PHI(     123)=PHI(     123)+FAC*PHI(     114)
      PHI(     124)=PHI(     124)+FAC*PHI(     115)
      FAC=SFAC( 4)
      J=     135
      DO I=     142,     148
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 5)
      J=     135
      DO I=     156,     162
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 4)
      J=     170
      DO I=     177,     183
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 5)
      J=     170
      DO I=     191,     197
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 4)
      J=     218
      DO I=     231,     243
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 5)
      J=     218
      DO I=     257,     269
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 5)
      J=     283
      DO I=     296,     308
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 5)
      J=     344
      DO I=     366,     387
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 5)
      J=     410
      DO I=     432,     453
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 5)
      J=     510
      DO I=     544,     577
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      END IF
      END DO
      FAC=EXL(LL)*TWO
      SF=ONE
      SFAC(1)=SF
      SF=SF*FAC
      SFAC( 2)=SF
      SF=SF*FAC
      SFAC( 3)=SF
      SF=SF*FAC
      SFAC( 4)=SF
      SF=SF*FAC
      SFAC( 5)=SF
      FAC=SFAC( 1)
      PHI(      84)=PHI(      84)+FAC*PHI(      83)
      FAC=SFAC( 3)
      PHI(      86)=PHI(      86)+FAC*PHI(      85)
      FAC=SFAC( 5)
      PHI(      88)=PHI(      88)+FAC*PHI(      87)
      FAC=SFAC( 1)
      PHI(      98)=PHI(      98)+FAC*PHI(      95)
      PHI(      99)=PHI(      99)+FAC*PHI(      96)
      PHI(     100)=PHI(     100)+FAC*PHI(      97)
      FAC=SFAC( 3)
      PHI(     104)=PHI(     104)+FAC*PHI(     101)
      PHI(     105)=PHI(     105)+FAC*PHI(     102)
      PHI(     106)=PHI(     106)+FAC*PHI(     103)
      FAC=SFAC( 5)
      PHI(     110)=PHI(     110)+FAC*PHI(     107)
      PHI(     111)=PHI(     111)+FAC*PHI(     108)
      PHI(     112)=PHI(     112)+FAC*PHI(     109)
      FAC=SFAC( 2)
      PHI(     119)=PHI(     119)+FAC*PHI(     116)
      PHI(     120)=PHI(     120)+FAC*PHI(     117)
      PHI(     121)=PHI(     121)+FAC*PHI(     118)
      FAC=SFAC( 4)
      PHI(     125)=PHI(     125)+FAC*PHI(     122)
      PHI(     126)=PHI(     126)+FAC*PHI(     123)
      PHI(     127)=PHI(     127)+FAC*PHI(     124)
      FAC=SFAC( 2)
      J=     142
      DO I=     149,     155
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 4)
      J=     156
      DO I=     163,     169
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     177
      DO I=     184,     190
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 3)
      J=     191
      DO I=     198,     204
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     231
      DO I=     244,     256
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 3)
      J=     257
      DO I=     270,     282
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 2)
      J=     296
      DO I=     309,     321
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 2)
      J=     366
      DO I=     388,     409
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     432
      DO I=     454,     475
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     544
      DO I=     578,     611
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      END DO
C  POST-CONTRACTION PHASE
      PHI(     150)=PHI(     150)-D01*PHI(     149)
      PHI(     152)=PHI(     152)-D01*PHI(     149)
      PHI(     155)=PHI(     155)-D01*PHI(     149)
      PHI(     164)=PHI(     164)-D01*PHI(     163)
      PHI(     166)=PHI(     166)-D01*PHI(     163)
      PHI(     169)=PHI(     169)-D01*PHI(     163)
      PHI(     185)=PHI(     185)-D01*PHI(     184)
      PHI(     187)=PHI(     187)-D01*PHI(     184)
      PHI(     190)=PHI(     190)-D01*PHI(     184)
      PHI(     199)=PHI(     199)-D01*PHI(     198)
      PHI(     201)=PHI(     201)-D01*PHI(     198)
      PHI(     204)=PHI(     204)-D01*PHI(     198)
      PHI(     247)=PHI(     247)-D01*PHI(     244)
      PHI(     248)=PHI(     248)-D01*PHI(     245)
      PHI(     250)=PHI(     250)-D01*PHI(     245)
      PHI(     251)=PHI(     251)-D01*PHI(     246)
      PHI(     253)=PHI(     253)-D01*PHI(     246)
      PHI(     256)=PHI(     256)-D01*PHI(     246)
      PHI(     247)=PHI(     247)-D02*PHI(     244)
      PHI(     249)=PHI(     249)-D01*PHI(     244)
      PHI(     250)=PHI(     250)-D02*PHI(     245)
      PHI(     254)=PHI(     254)-D01*PHI(     244)
      PHI(     255)=PHI(     255)-D01*PHI(     245)
      PHI(     256)=PHI(     256)-D02*PHI(     246)
      PHI(     273)=PHI(     273)-D01*PHI(     270)
      PHI(     274)=PHI(     274)-D01*PHI(     271)
      PHI(     276)=PHI(     276)-D01*PHI(     271)
      PHI(     277)=PHI(     277)-D01*PHI(     272)
      PHI(     279)=PHI(     279)-D01*PHI(     272)
      PHI(     282)=PHI(     282)-D01*PHI(     272)
      PHI(     273)=PHI(     273)-D02*PHI(     270)
      PHI(     275)=PHI(     275)-D01*PHI(     270)
      PHI(     276)=PHI(     276)-D02*PHI(     271)
      PHI(     280)=PHI(     280)-D01*PHI(     270)
      PHI(     281)=PHI(     281)-D01*PHI(     271)
      PHI(     282)=PHI(     282)-D02*PHI(     272)
      PHI(     312)=PHI(     312)-D01*PHI(     309)
      PHI(     313)=PHI(     313)-D01*PHI(     310)
      PHI(     315)=PHI(     315)-D01*PHI(     310)
      PHI(     316)=PHI(     316)-D01*PHI(     311)
      PHI(     318)=PHI(     318)-D01*PHI(     311)
      PHI(     321)=PHI(     321)-D01*PHI(     311)
      PHI(     312)=PHI(     312)-D02*PHI(     309)
      PHI(     314)=PHI(     314)-D01*PHI(     309)
      PHI(     315)=PHI(     315)-D02*PHI(     310)
      PHI(     319)=PHI(     319)-D01*PHI(     309)
      PHI(     320)=PHI(     320)-D01*PHI(     310)
      PHI(     321)=PHI(     321)-D02*PHI(     311)
      PHI(     395)=PHI(     395)-D01*PHI(     389)
      PHI(     396)=PHI(     396)-D01*PHI(     390)
      PHI(     397)=PHI(     397)-D01*PHI(     391)
      PHI(     399)=PHI(     399)-D01*PHI(     391)
      PHI(     400)=PHI(     400)-D01*PHI(     392)
      PHI(     401)=PHI(     401)-D01*PHI(     393)
      PHI(     403)=PHI(     403)-D01*PHI(     393)
      PHI(     404)=PHI(     404)-D01*PHI(     394)
      PHI(     406)=PHI(     406)-D01*PHI(     394)
      PHI(     409)=PHI(     409)-D01*PHI(     394)
      PHI(     395)=PHI(     395)-D02*PHI(     389)
      PHI(     396)=PHI(     396)-D02*PHI(     390)
      PHI(     398)=PHI(     398)-D01*PHI(     390)
      PHI(     399)=PHI(     399)-D02*PHI(     391)
      PHI(     400)=PHI(     400)-D02*PHI(     392)
      PHI(     402)=PHI(     402)-D01*PHI(     392)
      PHI(     403)=PHI(     403)-D02*PHI(     393)
      PHI(     407)=PHI(     407)-D01*PHI(     392)
      PHI(     408)=PHI(     408)-D01*PHI(     393)
      PHI(     409)=PHI(     409)-D02*PHI(     394)
      PHI(     389)=PHI(     389)-D01*PHI(     388)
      PHI(     391)=PHI(     391)-D01*PHI(     388)
      PHI(     394)=PHI(     394)-D01*PHI(     388)
      PHI(     395)=PHI(     395)-D03*PHI(     389)
      PHI(     397)=PHI(     397)-D01*PHI(     389)
      PHI(     398)=PHI(     398)-D02*PHI(     390)
      PHI(     399)=PHI(     399)-D03*PHI(     391)
      PHI(     404)=PHI(     404)-D01*PHI(     389)
      PHI(     405)=PHI(     405)-D01*PHI(     390)
      PHI(     406)=PHI(     406)-D01*PHI(     391)
      PHI(     407)=PHI(     407)-D02*PHI(     392)
      PHI(     408)=PHI(     408)-D02*PHI(     393)
      PHI(     409)=PHI(     409)-D03*PHI(     394)
      PHI(     461)=PHI(     461)-D01*PHI(     455)
      PHI(     462)=PHI(     462)-D01*PHI(     456)
      PHI(     463)=PHI(     463)-D01*PHI(     457)
      PHI(     465)=PHI(     465)-D01*PHI(     457)
      PHI(     466)=PHI(     466)-D01*PHI(     458)
      PHI(     467)=PHI(     467)-D01*PHI(     459)
      PHI(     469)=PHI(     469)-D01*PHI(     459)
      PHI(     470)=PHI(     470)-D01*PHI(     460)
      PHI(     472)=PHI(     472)-D01*PHI(     460)
      PHI(     475)=PHI(     475)-D01*PHI(     460)
      PHI(     461)=PHI(     461)-D02*PHI(     455)
      PHI(     462)=PHI(     462)-D02*PHI(     456)
      PHI(     464)=PHI(     464)-D01*PHI(     456)
      PHI(     465)=PHI(     465)-D02*PHI(     457)
      PHI(     466)=PHI(     466)-D02*PHI(     458)
      PHI(     468)=PHI(     468)-D01*PHI(     458)
      PHI(     469)=PHI(     469)-D02*PHI(     459)
      PHI(     473)=PHI(     473)-D01*PHI(     458)
      PHI(     474)=PHI(     474)-D01*PHI(     459)
      PHI(     475)=PHI(     475)-D02*PHI(     460)
      PHI(     455)=PHI(     455)-D01*PHI(     454)
      PHI(     457)=PHI(     457)-D01*PHI(     454)
      PHI(     460)=PHI(     460)-D01*PHI(     454)
      PHI(     461)=PHI(     461)-D03*PHI(     455)
      PHI(     463)=PHI(     463)-D01*PHI(     455)
      PHI(     464)=PHI(     464)-D02*PHI(     456)
      PHI(     465)=PHI(     465)-D03*PHI(     457)
      PHI(     470)=PHI(     470)-D01*PHI(     455)
      PHI(     471)=PHI(     471)-D01*PHI(     456)
      PHI(     472)=PHI(     472)-D01*PHI(     457)
      PHI(     473)=PHI(     473)-D02*PHI(     458)
      PHI(     474)=PHI(     474)-D02*PHI(     459)
      PHI(     475)=PHI(     475)-D03*PHI(     460)
      PHI(     591)=PHI(     591)-D01*PHI(     581)
      PHI(     592)=PHI(     592)-D01*PHI(     582)
      PHI(     593)=PHI(     593)-D01*PHI(     583)
      PHI(     594)=PHI(     594)-D01*PHI(     584)
      PHI(     596)=PHI(     596)-D01*PHI(     584)
      PHI(     597)=PHI(     597)-D01*PHI(     585)
      PHI(     598)=PHI(     598)-D01*PHI(     586)
      PHI(     599)=PHI(     599)-D01*PHI(     587)
      PHI(     601)=PHI(     601)-D01*PHI(     587)
      PHI(     602)=PHI(     602)-D01*PHI(     588)
      PHI(     603)=PHI(     603)-D01*PHI(     589)
      PHI(     605)=PHI(     605)-D01*PHI(     589)
      PHI(     606)=PHI(     606)-D01*PHI(     590)
      PHI(     608)=PHI(     608)-D01*PHI(     590)
      PHI(     611)=PHI(     611)-D01*PHI(     590)
      PHI(     591)=PHI(     591)-D02*PHI(     581)
      PHI(     592)=PHI(     592)-D02*PHI(     582)
      PHI(     593)=PHI(     593)-D02*PHI(     583)
      PHI(     595)=PHI(     595)-D01*PHI(     583)
      PHI(     596)=PHI(     596)-D02*PHI(     584)
      PHI(     597)=PHI(     597)-D02*PHI(     585)
      PHI(     598)=PHI(     598)-D02*PHI(     586)
      PHI(     600)=PHI(     600)-D01*PHI(     586)
      PHI(     601)=PHI(     601)-D02*PHI(     587)
      PHI(     602)=PHI(     602)-D02*PHI(     588)
      PHI(     604)=PHI(     604)-D01*PHI(     588)
      PHI(     605)=PHI(     605)-D02*PHI(     589)
      PHI(     609)=PHI(     609)-D01*PHI(     588)
      PHI(     610)=PHI(     610)-D01*PHI(     589)
      PHI(     611)=PHI(     611)-D02*PHI(     590)
      PHI(     581)=PHI(     581)-D01*PHI(     578)
      PHI(     582)=PHI(     582)-D01*PHI(     579)
      PHI(     584)=PHI(     584)-D01*PHI(     579)
      PHI(     585)=PHI(     585)-D01*PHI(     580)
      PHI(     587)=PHI(     587)-D01*PHI(     580)
      PHI(     590)=PHI(     590)-D01*PHI(     580)
      PHI(     591)=PHI(     591)-D03*PHI(     581)
      PHI(     592)=PHI(     592)-D03*PHI(     582)
      PHI(     594)=PHI(     594)-D01*PHI(     582)
      PHI(     595)=PHI(     595)-D02*PHI(     583)
      PHI(     596)=PHI(     596)-D03*PHI(     584)
      PHI(     597)=PHI(     597)-D03*PHI(     585)
      PHI(     599)=PHI(     599)-D01*PHI(     585)
      PHI(     600)=PHI(     600)-D02*PHI(     586)
      PHI(     601)=PHI(     601)-D03*PHI(     587)
      PHI(     606)=PHI(     606)-D01*PHI(     585)
      PHI(     607)=PHI(     607)-D01*PHI(     586)
      PHI(     608)=PHI(     608)-D01*PHI(     587)
      PHI(     609)=PHI(     609)-D02*PHI(     588)
      PHI(     610)=PHI(     610)-D02*PHI(     589)
      PHI(     611)=PHI(     611)-D03*PHI(     590)
      PHI(     581)=PHI(     581)-D02*PHI(     578)
      PHI(     583)=PHI(     583)-D01*PHI(     578)
      PHI(     584)=PHI(     584)-D02*PHI(     579)
      PHI(     588)=PHI(     588)-D01*PHI(     578)
      PHI(     589)=PHI(     589)-D01*PHI(     579)
      PHI(     590)=PHI(     590)-D02*PHI(     580)
      PHI(     591)=PHI(     591)-D04*PHI(     581)
      PHI(     593)=PHI(     593)-D01*PHI(     581)
      PHI(     594)=PHI(     594)-D02*PHI(     582)
      PHI(     595)=PHI(     595)-D03*PHI(     583)
      PHI(     596)=PHI(     596)-D04*PHI(     584)
      PHI(     602)=PHI(     602)-D01*PHI(     581)
      PHI(     603)=PHI(     603)-D01*PHI(     582)
      PHI(     604)=PHI(     604)-D01*PHI(     583)
      PHI(     605)=PHI(     605)-D01*PHI(     584)
      PHI(     606)=PHI(     606)-D02*PHI(     585)
      PHI(     607)=PHI(     607)-D02*PHI(     586)
      PHI(     608)=PHI(     608)-D02*PHI(     587)
      PHI(     609)=PHI(     609)-D03*PHI(     588)
      PHI(     610)=PHI(     610)-D03*PHI(     589)
      PHI(     611)=PHI(     611)-D04*PHI(     590)
      CONF(1)=XB-XA
      CONF(2)=YB-YA
      CONF(3)=ZB-ZA
      WRK1(    2)= PHI(     110)+CONF( 1)*PHI(      88)
      WRK1(    3)= PHI(     111)+CONF( 2)*PHI(      88)
      WRK1(    4)= PHI(     112)+CONF( 3)*PHI(      88)
      CALL GSPS_BCTE(WRK1,WRK2,LENW, 1)
      WRK1(    2)= PHI(     104)+CONF( 1)*PHI(      86)
      WRK1(    3)= PHI(     105)+CONF( 2)*PHI(      86)
      WRK1(    4)= PHI(     106)+CONF( 3)*PHI(      86)
      CALL GSPS_BCTE(WRK1,WRK2,LENW, 2)
      WRK1(    2)= PHI(      98)+CONF( 1)*PHI(      84)
      WRK1(    3)= PHI(      99)+CONF( 2)*PHI(      84)
      WRK1(    4)= PHI(     100)+CONF( 3)*PHI(      84)
      CALL GSPS_BCTE(WRK1,WRK2,LENW, 3)
      WRK1(    2)=-PHI(     164)-CONF( 1)*PHI(     125)
      WRK1(    3)=-PHI(     165)-CONF( 2)*PHI(     125)
      WRK1(    4)=-PHI(     167)-CONF( 3)*PHI(     125)
      CALL GSPS_BCTE(WRK1,WRK2,LENW, 4)
      WRK1(    2)=-PHI(     165)-CONF( 1)*PHI(     126)
      WRK1(    3)=-PHI(     166)-CONF( 2)*PHI(     126)
      WRK1(    4)=-PHI(     168)-CONF( 3)*PHI(     126)
      CALL GSPS_BCTE(WRK1,WRK2,LENW, 5)
      WRK1(    2)=-PHI(     167)-CONF( 1)*PHI(     127)
      WRK1(    3)=-PHI(     168)-CONF( 2)*PHI(     127)
      WRK1(    4)=-PHI(     169)-CONF( 3)*PHI(     127)
      CALL GSPS_BCTE(WRK1,WRK2,LENW, 6)
      WRK1(    2)=-PHI(     150)-CONF( 1)*PHI(     119)
      WRK1(    3)=-PHI(     151)-CONF( 2)*PHI(     119)
      WRK1(    4)=-PHI(     153)-CONF( 3)*PHI(     119)
      CALL GSPS_BCTE(WRK1,WRK2,LENW, 7)
      WRK1(    2)=-PHI(     151)-CONF( 1)*PHI(     120)
      WRK1(    3)=-PHI(     152)-CONF( 2)*PHI(     120)
      WRK1(    4)=-PHI(     154)-CONF( 3)*PHI(     120)
      CALL GSPS_BCTE(WRK1,WRK2,LENW, 8)
      WRK1(    2)=-PHI(     153)-CONF( 1)*PHI(     121)
      WRK1(    3)=-PHI(     154)-CONF( 2)*PHI(     121)
      WRK1(    4)=-PHI(     155)-CONF( 3)*PHI(     121)
      CALL GSPS_BCTE(WRK1,WRK2,LENW, 9)
      WRK1(    2)= PHI(     273)+CONF( 1)*PHI(     199)
      WRK1(    3)= PHI(     274)+CONF( 2)*PHI(     199)
      WRK1(    4)= PHI(     277)+CONF( 3)*PHI(     199)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,10)
      WRK1(    2)= PHI(     274)+CONF( 1)*PHI(     200)
      WRK1(    3)= PHI(     275)+CONF( 2)*PHI(     200)
      WRK1(    4)= PHI(     278)+CONF( 3)*PHI(     200)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,11)
      WRK1(    2)= PHI(     275)+CONF( 1)*PHI(     201)
      WRK1(    3)= PHI(     276)+CONF( 2)*PHI(     201)
      WRK1(    4)= PHI(     279)+CONF( 3)*PHI(     201)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,12)
      WRK1(    2)= PHI(     277)+CONF( 1)*PHI(     202)
      WRK1(    3)= PHI(     278)+CONF( 2)*PHI(     202)
      WRK1(    4)= PHI(     280)+CONF( 3)*PHI(     202)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,13)
      WRK1(    2)= PHI(     278)+CONF( 1)*PHI(     203)
      WRK1(    3)= PHI(     279)+CONF( 2)*PHI(     203)
      WRK1(    4)= PHI(     281)+CONF( 3)*PHI(     203)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,14)
      WRK1(    2)= PHI(     280)+CONF( 1)*PHI(     204)
      WRK1(    3)= PHI(     281)+CONF( 2)*PHI(     204)
      WRK1(    4)= PHI(     282)+CONF( 3)*PHI(     204)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,15)
      WRK1(    2)= PHI(     247)+CONF( 1)*PHI(     185)
      WRK1(    3)= PHI(     248)+CONF( 2)*PHI(     185)
      WRK1(    4)= PHI(     251)+CONF( 3)*PHI(     185)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,16)
      WRK1(    2)= PHI(     248)+CONF( 1)*PHI(     186)
      WRK1(    3)= PHI(     249)+CONF( 2)*PHI(     186)
      WRK1(    4)= PHI(     252)+CONF( 3)*PHI(     186)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,17)
      WRK1(    2)= PHI(     249)+CONF( 1)*PHI(     187)
      WRK1(    3)= PHI(     250)+CONF( 2)*PHI(     187)
      WRK1(    4)= PHI(     253)+CONF( 3)*PHI(     187)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,18)
      WRK1(    2)= PHI(     251)+CONF( 1)*PHI(     188)
      WRK1(    3)= PHI(     252)+CONF( 2)*PHI(     188)
      WRK1(    4)= PHI(     254)+CONF( 3)*PHI(     188)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,19)
      WRK1(    2)= PHI(     252)+CONF( 1)*PHI(     189)
      WRK1(    3)= PHI(     253)+CONF( 2)*PHI(     189)
      WRK1(    4)= PHI(     255)+CONF( 3)*PHI(     189)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,20)
      WRK1(    2)= PHI(     254)+CONF( 1)*PHI(     190)
      WRK1(    3)= PHI(     255)+CONF( 2)*PHI(     190)
      WRK1(    4)= PHI(     256)+CONF( 3)*PHI(     190)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,21)
      WRK1(    2)=-PHI(     395)-CONF( 1)*PHI(     312)
      WRK1(    3)=-PHI(     396)-CONF( 2)*PHI(     312)
      WRK1(    4)=-PHI(     400)-CONF( 3)*PHI(     312)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,22)
      WRK1(    2)=-PHI(     396)-CONF( 1)*PHI(     313)
      WRK1(    3)=-PHI(     397)-CONF( 2)*PHI(     313)
      WRK1(    4)=-PHI(     401)-CONF( 3)*PHI(     313)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,23)
      WRK1(    2)=-PHI(     397)-CONF( 1)*PHI(     314)
      WRK1(    3)=-PHI(     398)-CONF( 2)*PHI(     314)
      WRK1(    4)=-PHI(     402)-CONF( 3)*PHI(     314)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,24)
      WRK1(    2)=-PHI(     398)-CONF( 1)*PHI(     315)
      WRK1(    3)=-PHI(     399)-CONF( 2)*PHI(     315)
      WRK1(    4)=-PHI(     403)-CONF( 3)*PHI(     315)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,25)
      WRK1(    2)=-PHI(     400)-CONF( 1)*PHI(     316)
      WRK1(    3)=-PHI(     401)-CONF( 2)*PHI(     316)
      WRK1(    4)=-PHI(     404)-CONF( 3)*PHI(     316)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,26)
      WRK1(    2)=-PHI(     401)-CONF( 1)*PHI(     317)
      WRK1(    3)=-PHI(     402)-CONF( 2)*PHI(     317)
      WRK1(    4)=-PHI(     405)-CONF( 3)*PHI(     317)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,27)
      WRK1(    2)=-PHI(     402)-CONF( 1)*PHI(     318)
      WRK1(    3)=-PHI(     403)-CONF( 2)*PHI(     318)
      WRK1(    4)=-PHI(     406)-CONF( 3)*PHI(     318)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,28)
      WRK1(    2)=-PHI(     404)-CONF( 1)*PHI(     319)
      WRK1(    3)=-PHI(     405)-CONF( 2)*PHI(     319)
      WRK1(    4)=-PHI(     407)-CONF( 3)*PHI(     319)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,29)
      WRK1(    2)=-PHI(     405)-CONF( 1)*PHI(     320)
      WRK1(    3)=-PHI(     406)-CONF( 2)*PHI(     320)
      WRK1(    4)=-PHI(     408)-CONF( 3)*PHI(     320)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,30)
      WRK1(    2)=-PHI(     407)-CONF( 1)*PHI(     321)
      WRK1(    3)=-PHI(     408)-CONF( 2)*PHI(     321)
      WRK1(    4)=-PHI(     409)-CONF( 3)*PHI(     321)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,31)
      WRK1(    2)= PHI(     591)+CONF( 1)*PHI(     461)
      WRK1(    3)= PHI(     592)+CONF( 2)*PHI(     461)
      WRK1(    4)= PHI(     597)+CONF( 3)*PHI(     461)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,32)
      WRK1(    2)= PHI(     592)+CONF( 1)*PHI(     462)
      WRK1(    3)= PHI(     593)+CONF( 2)*PHI(     462)
      WRK1(    4)= PHI(     598)+CONF( 3)*PHI(     462)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,33)
      WRK1(    2)= PHI(     593)+CONF( 1)*PHI(     463)
      WRK1(    3)= PHI(     594)+CONF( 2)*PHI(     463)
      WRK1(    4)= PHI(     599)+CONF( 3)*PHI(     463)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,34)
      WRK1(    2)= PHI(     594)+CONF( 1)*PHI(     464)
      WRK1(    3)= PHI(     595)+CONF( 2)*PHI(     464)
      WRK1(    4)= PHI(     600)+CONF( 3)*PHI(     464)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,35)
      WRK1(    2)= PHI(     595)+CONF( 1)*PHI(     465)
      WRK1(    3)= PHI(     596)+CONF( 2)*PHI(     465)
      WRK1(    4)= PHI(     601)+CONF( 3)*PHI(     465)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,36)
      WRK1(    2)= PHI(     597)+CONF( 1)*PHI(     466)
      WRK1(    3)= PHI(     598)+CONF( 2)*PHI(     466)
      WRK1(    4)= PHI(     602)+CONF( 3)*PHI(     466)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,37)
      WRK1(    2)= PHI(     598)+CONF( 1)*PHI(     467)
      WRK1(    3)= PHI(     599)+CONF( 2)*PHI(     467)
      WRK1(    4)= PHI(     603)+CONF( 3)*PHI(     467)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,38)
      WRK1(    2)= PHI(     599)+CONF( 1)*PHI(     468)
      WRK1(    3)= PHI(     600)+CONF( 2)*PHI(     468)
      WRK1(    4)= PHI(     604)+CONF( 3)*PHI(     468)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,39)
      WRK1(    2)= PHI(     600)+CONF( 1)*PHI(     469)
      WRK1(    3)= PHI(     601)+CONF( 2)*PHI(     469)
      WRK1(    4)= PHI(     605)+CONF( 3)*PHI(     469)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,40)
      WRK1(    2)= PHI(     602)+CONF( 1)*PHI(     470)
      WRK1(    3)= PHI(     603)+CONF( 2)*PHI(     470)
      WRK1(    4)= PHI(     606)+CONF( 3)*PHI(     470)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,41)
      WRK1(    2)= PHI(     603)+CONF( 1)*PHI(     471)
      WRK1(    3)= PHI(     604)+CONF( 2)*PHI(     471)
      WRK1(    4)= PHI(     607)+CONF( 3)*PHI(     471)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,42)
      WRK1(    2)= PHI(     604)+CONF( 1)*PHI(     472)
      WRK1(    3)= PHI(     605)+CONF( 2)*PHI(     472)
      WRK1(    4)= PHI(     608)+CONF( 3)*PHI(     472)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,43)
      WRK1(    2)= PHI(     606)+CONF( 1)*PHI(     473)
      WRK1(    3)= PHI(     607)+CONF( 2)*PHI(     473)
      WRK1(    4)= PHI(     609)+CONF( 3)*PHI(     473)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,44)
      WRK1(    2)= PHI(     607)+CONF( 1)*PHI(     474)
      WRK1(    3)= PHI(     608)+CONF( 2)*PHI(     474)
      WRK1(    4)= PHI(     610)+CONF( 3)*PHI(     474)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,45)
      WRK1(    2)= PHI(     609)+CONF( 1)*PHI(     475)
      WRK1(    3)= PHI(     610)+CONF( 2)*PHI(     475)
      WRK1(    4)= PHI(     611)+CONF( 3)*PHI(     475)
      CALL GSPS_BCTE(WRK1,WRK2,LENW,46)
      CONF(1)=XD-XC
      CONF(2)=YD-YC
      CONF(3)=ZD-ZC
      DO I=1,     3
      WRK2(   32,I)=WRK2(   32,I)+CONF( 1)*WRK2(   22,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 6)*WRK2(   16,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 1)*WRK2(   23,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 5)*WRK2(   17,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 1)*WRK2(   24,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 4)*WRK2(   18,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 1)*WRK2(   25,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 2)*WRK2(   25,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 6)*WRK2(   18,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 1)*WRK2(   26,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 5)*WRK2(   19,I)
      WRK2(   38,I)=WRK2(   38,I)+CONF( 1)*WRK2(   27,I)
      WRK2(   38,I)=WRK2(   38,I)+CONF( 4)*WRK2(   20,I)
      WRK2(   39,I)=WRK2(   39,I)+CONF( 1)*WRK2(   28,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 2)*WRK2(   28,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 5)*WRK2(   20,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 1)*WRK2(   29,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 4)*WRK2(   21,I)
      WRK2(   42,I)=WRK2(   42,I)+CONF( 1)*WRK2(   30,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 2)*WRK2(   30,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 4)*WRK2(   21,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 1)*WRK2(   31,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 2)*WRK2(   31,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 3)*WRK2(   31,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 6)*WRK2(   21,I)
      WRK2(   22,I)=WRK2(   22,I)+CONF( 1)*WRK2(   10,I)
      WRK2(   22,I)=WRK2(   22,I)+CONF( 5)*WRK2(    7,I)
      WRK2(   23,I)=WRK2(   23,I)+CONF( 1)*WRK2(   11,I)
      WRK2(   23,I)=WRK2(   23,I)+CONF( 4)*WRK2(    8,I)
      WRK2(   24,I)=WRK2(   24,I)+CONF( 1)*WRK2(   12,I)
      WRK2(   25,I)=WRK2(   25,I)+CONF( 2)*WRK2(   12,I)
      WRK2(   25,I)=WRK2(   25,I)+CONF( 5)*WRK2(    8,I)
      WRK2(   26,I)=WRK2(   26,I)+CONF( 1)*WRK2(   13,I)
      WRK2(   26,I)=WRK2(   26,I)+CONF( 4)*WRK2(    9,I)
      WRK2(   27,I)=WRK2(   27,I)+CONF( 1)*WRK2(   14,I)
      WRK2(   28,I)=WRK2(   28,I)+CONF( 2)*WRK2(   14,I)
      WRK2(   28,I)=WRK2(   28,I)+CONF( 4)*WRK2(    9,I)
      WRK2(   29,I)=WRK2(   29,I)+CONF( 1)*WRK2(   15,I)
      WRK2(   30,I)=WRK2(   30,I)+CONF( 2)*WRK2(   15,I)
      WRK2(   31,I)=WRK2(   31,I)+CONF( 3)*WRK2(   15,I)
      WRK2(   31,I)=WRK2(   31,I)+CONF( 5)*WRK2(    9,I)
      WRK2(   10,I)=WRK2(   10,I)+CONF( 1)*WRK2(    4,I)
      WRK2(   10,I)=WRK2(   10,I)+CONF( 4)*WRK2(    2,I)
      WRK2(   11,I)=WRK2(   11,I)+CONF( 1)*WRK2(    5,I)
      WRK2(   12,I)=WRK2(   12,I)+CONF( 2)*WRK2(    5,I)
      WRK2(   12,I)=WRK2(   12,I)+CONF( 4)*WRK2(    2,I)
      WRK2(   13,I)=WRK2(   13,I)+CONF( 1)*WRK2(    6,I)
      WRK2(   14,I)=WRK2(   14,I)+CONF( 2)*WRK2(    6,I)
      WRK2(   15,I)=WRK2(   15,I)+CONF( 3)*WRK2(    6,I)
      WRK2(   15,I)=WRK2(   15,I)+CONF( 4)*WRK2(    2,I)
      WRK2(   16,I)=WRK2(   16,I)+CONF( 1)*WRK2(    7,I)
      WRK2(   16,I)=WRK2(   16,I)+CONF( 4)*WRK2(    3,I)
      WRK2(   17,I)=WRK2(   17,I)+CONF( 1)*WRK2(    8,I)
      WRK2(   18,I)=WRK2(   18,I)+CONF( 2)*WRK2(    8,I)
      WRK2(   18,I)=WRK2(   18,I)+CONF( 4)*WRK2(    3,I)
      WRK2(   19,I)=WRK2(   19,I)+CONF( 1)*WRK2(    9,I)
      WRK2(   20,I)=WRK2(   20,I)+CONF( 2)*WRK2(    9,I)
      WRK2(   21,I)=WRK2(   21,I)+CONF( 3)*WRK2(    9,I)
      WRK2(   21,I)=WRK2(   21,I)+CONF( 4)*WRK2(    3,I)
      WRK2(    4,I)=WRK2(    4,I)+CONF( 1)*WRK2(    1,I)
      WRK2(    5,I)=WRK2(    5,I)+CONF( 2)*WRK2(    1,I)
      WRK2(    6,I)=WRK2(    6,I)+CONF( 3)*WRK2(    1,I)
      WRK2(    7,I)=WRK2(    7,I)+CONF( 1)*WRK2(    2,I)
      WRK2(    8,I)=WRK2(    8,I)+CONF( 2)*WRK2(    2,I)
      WRK2(    9,I)=WRK2(    9,I)+CONF( 3)*WRK2(    2,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 1)*WRK2(   22,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 5)*WRK2(   16,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 1)*WRK2(   23,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 4)*WRK2(   17,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 1)*WRK2(   24,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 2)*WRK2(   24,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 5)*WRK2(   17,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 2)*WRK2(   25,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 5)*WRK2(   18,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 1)*WRK2(   26,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 4)*WRK2(   19,I)
      WRK2(   38,I)=WRK2(   38,I)+CONF( 1)*WRK2(   27,I)
      WRK2(   39,I)=WRK2(   39,I)+CONF( 2)*WRK2(   27,I)
      WRK2(   39,I)=WRK2(   39,I)+CONF( 4)*WRK2(   19,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 2)*WRK2(   28,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 4)*WRK2(   20,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 1)*WRK2(   29,I)
      WRK2(   42,I)=WRK2(   42,I)+CONF( 2)*WRK2(   29,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 2)*WRK2(   30,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 3)*WRK2(   29,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 5)*WRK2(   19,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 3)*WRK2(   30,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 5)*WRK2(   20,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 3)*WRK2(   31,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 5)*WRK2(   21,I)
      WRK2(   22,I)=WRK2(   22,I)+CONF( 1)*WRK2(   10,I)
      WRK2(   22,I)=WRK2(   22,I)+CONF( 4)*WRK2(    7,I)
      WRK2(   23,I)=WRK2(   23,I)+CONF( 1)*WRK2(   11,I)
      WRK2(   24,I)=WRK2(   24,I)+CONF( 2)*WRK2(   11,I)
      WRK2(   24,I)=WRK2(   24,I)+CONF( 4)*WRK2(    7,I)
      WRK2(   25,I)=WRK2(   25,I)+CONF( 2)*WRK2(   12,I)
      WRK2(   25,I)=WRK2(   25,I)+CONF( 4)*WRK2(    8,I)
      WRK2(   26,I)=WRK2(   26,I)+CONF( 1)*WRK2(   13,I)
      WRK2(   27,I)=WRK2(   27,I)+CONF( 2)*WRK2(   13,I)
      WRK2(   28,I)=WRK2(   28,I)+CONF( 2)*WRK2(   14,I)
      WRK2(   29,I)=WRK2(   29,I)+CONF( 3)*WRK2(   13,I)
      WRK2(   29,I)=WRK2(   29,I)+CONF( 4)*WRK2(    7,I)
      WRK2(   30,I)=WRK2(   30,I)+CONF( 3)*WRK2(   14,I)
      WRK2(   30,I)=WRK2(   30,I)+CONF( 4)*WRK2(    8,I)
      WRK2(   31,I)=WRK2(   31,I)+CONF( 3)*WRK2(   15,I)
      WRK2(   31,I)=WRK2(   31,I)+CONF( 4)*WRK2(    9,I)
      WRK2(   10,I)=WRK2(   10,I)+CONF( 1)*WRK2(    4,I)
      WRK2(   11,I)=WRK2(   11,I)+CONF( 2)*WRK2(    4,I)
      WRK2(   12,I)=WRK2(   12,I)+CONF( 2)*WRK2(    5,I)
      WRK2(   13,I)=WRK2(   13,I)+CONF( 3)*WRK2(    4,I)
      WRK2(   14,I)=WRK2(   14,I)+CONF( 3)*WRK2(    5,I)
      WRK2(   15,I)=WRK2(   15,I)+CONF( 3)*WRK2(    6,I)
      WRK2(   16,I)=WRK2(   16,I)+CONF( 1)*WRK2(    7,I)
      WRK2(   17,I)=WRK2(   17,I)+CONF( 2)*WRK2(    7,I)
      WRK2(   18,I)=WRK2(   18,I)+CONF( 2)*WRK2(    8,I)
      WRK2(   19,I)=WRK2(   19,I)+CONF( 3)*WRK2(    7,I)
      WRK2(   20,I)=WRK2(   20,I)+CONF( 3)*WRK2(    8,I)
      WRK2(   21,I)=WRK2(   21,I)+CONF( 3)*WRK2(    9,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 1)*WRK2(   22,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 4)*WRK2(   16,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 1)*WRK2(   23,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 2)*WRK2(   23,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 4)*WRK2(   16,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 2)*WRK2(   24,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 4)*WRK2(   17,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 2)*WRK2(   25,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 4)*WRK2(   18,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 1)*WRK2(   26,I)
      WRK2(   38,I)=WRK2(   38,I)+CONF( 2)*WRK2(   26,I)
      WRK2(   39,I)=WRK2(   39,I)+CONF( 2)*WRK2(   27,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 2)*WRK2(   28,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 3)*WRK2(   26,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 4)*WRK2(   16,I)
      WRK2(   42,I)=WRK2(   42,I)+CONF( 3)*WRK2(   27,I)
      WRK2(   42,I)=WRK2(   42,I)+CONF( 4)*WRK2(   17,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 3)*WRK2(   28,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 4)*WRK2(   18,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 3)*WRK2(   29,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 4)*WRK2(   19,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 3)*WRK2(   30,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 4)*WRK2(   20,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 3)*WRK2(   31,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 4)*WRK2(   21,I)
      WRK2(   22,I)=WRK2(   22,I)+CONF( 1)*WRK2(   10,I)
      WRK2(   23,I)=WRK2(   23,I)+CONF( 2)*WRK2(   10,I)
      WRK2(   24,I)=WRK2(   24,I)+CONF( 2)*WRK2(   11,I)
      WRK2(   25,I)=WRK2(   25,I)+CONF( 2)*WRK2(   12,I)
      WRK2(   26,I)=WRK2(   26,I)+CONF( 3)*WRK2(   10,I)
      WRK2(   27,I)=WRK2(   27,I)+CONF( 3)*WRK2(   11,I)
      WRK2(   28,I)=WRK2(   28,I)+CONF( 3)*WRK2(   12,I)
      WRK2(   29,I)=WRK2(   29,I)+CONF( 3)*WRK2(   13,I)
      WRK2(   30,I)=WRK2(   30,I)+CONF( 3)*WRK2(   14,I)
      WRK2(   31,I)=WRK2(   31,I)+CONF( 3)*WRK2(   15,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 1)*WRK2(   22,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 2)*WRK2(   22,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 2)*WRK2(   23,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 2)*WRK2(   24,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 2)*WRK2(   25,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 3)*WRK2(   22,I)
      WRK2(   38,I)=WRK2(   38,I)+CONF( 3)*WRK2(   23,I)
      WRK2(   39,I)=WRK2(   39,I)+CONF( 3)*WRK2(   24,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 3)*WRK2(   25,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 3)*WRK2(   26,I)
      WRK2(   42,I)=WRK2(   42,I)+CONF( 3)*WRK2(   27,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 3)*WRK2(   28,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 3)*WRK2(   29,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 3)*WRK2(   30,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 3)*WRK2(   31,I)
      END DO
      RETURN
      END
C*MODULE INT2G   *DECK GSPS_BCTE
      SUBROUTINE GSPS_BCTE(WRK1,WRK2,LENW,JR)
      INTEGER LENW,I,J,JR
      DOUBLE PRECISION WRK1(*),WRK2(LENW,*)
      J=       1
      DO I=1,     3
      J=J+1
      WRK2(JR,I)=WRK1(J)
      END DO
      END
C  LPHI=     545
C  LWK1=       1
C  LWK2=     171
C  LENW=     171
C*MODULE INT2G   *DECK GPSS
      SUBROUTINE GPSS (IPRIM,JPRIM,KPRIM,LPRIM,IEQJ,KEQL,
     *                 XC,YC,ZC,XD,YD,ZD,PHI,WRK2,LENW)
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      COMMON /SHLNOS/ QQ4,IDUMMY(20)
      PARAMETER (NTX=4,NPF=450,NGRD=7,NPX=1000,MXQT=16)
      COMMON /FMTTBL/ FGRID(0:NTX,0:NPF,0:NGRD)
     *,XGRID(0:NTX,0:NPX),TMAX,RFINC(0:NGRD)
     *,RXINC,RMR(MXQT),TLGM(0:MXQT),NORD
      PARAMETER (MXPR=30, MXPS=MXPR*MXPR)
      COMMON /ERIPRM/ EXI(MXPR),EXJ(MXPR),EXK(MXPR),EXL(MXPR)
     *,CCI(MXPR),CCJ(MXPR),CCK(MXPR),CCL(MXPR)
     *,XAB(MXPS),YAB(MXPS),ZAB(MXPS)
     *,XCD(MXPS),YCD(MXPS),ZCD(MXPS)
     *,CCBRA(MXPS),CCKET(MXPS),RXB(MXPS)
     *,SLBRA(MXPS),SLKET(MXPS),RXK(MXPS)
      COMMON /NLRC  / LCFLAG,LRINT,EMU,EMU2,LRFILE
      LOGICAL   LCFLAG,LRINT
      LOGICAL IEQJ,KEQL
      DOUBLE PRECISION RXYZ(3),SFAC(9),CONF(11),FT(0:16)
      DOUBLE PRECISION SAUX(0: 5),CXYZ( 56)
      DOUBLE PRECISION PHI(*),WRK2(LENW,*)
      PARAMETER (ZERO=0.0D+00, ONE=1.0D+00, TWO=2.0D+00)
      PARAMETER (D01=1.0D+00, D02=2.0D+00, D03=3.0D+00)
      PARAMETER (D04=4.0D+00)
      PARAMETER (CCTOL=1.0D-13)
      DATA CONF/0.0D+00,0.0D+00,0.0D+00,1.0D+00,2.0D+00,3.0D+00
     &,         4.0D+00,5.0D+00,6.0D+00,7.0D+00,8.0D+00/
      SAVE CONF
      RFAC=RFINC( 5)
      KL=0
      PHI(      84)=ZERO
      PHI(      86)=ZERO
      PHI(      87)=ZERO
      PHI(      89)=ZERO
      PHI(      90)=ZERO
      PHI(      92)=ZERO
      PHI(     102)=ZERO
      PHI(     103)=ZERO
      PHI(     104)=ZERO
      PHI(     105)=ZERO
      PHI(     106)=ZERO
      PHI(     107)=ZERO
      PHI(     111)=ZERO
      PHI(     112)=ZERO
      PHI(     113)=ZERO
      PHI(     114)=ZERO
      PHI(     115)=ZERO
      PHI(     116)=ZERO
      PHI(     120)=ZERO
      PHI(     121)=ZERO
      PHI(     122)=ZERO
      DO I=     144,     150
      PHI(I)=ZERO
      END DO
      DO I=     158,     164
      PHI(I)=ZERO
      END DO
      DO I=     165,     171
      PHI(I)=ZERO
      END DO
      DO I=     179,     185
      PHI(I)=ZERO
      END DO
      DO I=     225,     237
      PHI(I)=ZERO
      END DO
      DO I=     238,     250
      PHI(I)=ZERO
      END DO
      DO I=     264,     276
      PHI(I)=ZERO
      END DO
      DO I=     343,     364
      PHI(I)=ZERO
      END DO
      DO I=     387,     408
      PHI(I)=ZERO
      END DO
      DO I=     511,     544
      PHI(I)=ZERO
      END DO
      DO LL=1,LPRIM
      PHI(      83)=ZERO
      PHI(      85)=ZERO
      PHI(      88)=ZERO
      PHI(      91)=ZERO
      PHI(      99)=ZERO
      PHI(     100)=ZERO
      PHI(     101)=ZERO
      PHI(     108)=ZERO
      PHI(     109)=ZERO
      PHI(     110)=ZERO
      PHI(     117)=ZERO
      PHI(     118)=ZERO
      PHI(     119)=ZERO
      DO I=     137,     143
      PHI(I)=ZERO
      END DO
      DO I=     151,     157
      PHI(I)=ZERO
      END DO
      DO I=     172,     178
      PHI(I)=ZERO
      END DO
      DO I=     212,     224
      PHI(I)=ZERO
      END DO
      DO I=     251,     263
      PHI(I)=ZERO
      END DO
      DO I=     321,     342
      PHI(I)=ZERO
      END DO
      DO I=     365,     386
      PHI(I)=ZERO
      END DO
      DO I=     477,     510
      PHI(I)=ZERO
      END DO
      KTOP=KPRIM
      IF (KEQL) KTOP=LL
      DO KK=1,KTOP
      KL=KL+1
      CCKF=CCKET(KL)
      IF (ABS(CCKF).GT.CCTOL) THEN
      EKL=EXL(LL)+EXK(KK)
      XKL=XCD(KL)
      YKL=YCD(KL)
      ZKL=ZCD(KL)
      IJ=0
      PHI(      82)=ZERO
      PHI(      96)=ZERO
      PHI(      97)=ZERO
      PHI(      98)=ZERO
      DO I=     130,     136
      PHI(I)=ZERO
      END DO
      DO I=     199,     211
      PHI(I)=ZERO
      END DO
      DO I=     299,     320
      PHI(I)=ZERO
      END DO
      DO I=     443,     476
      PHI(I)=ZERO
      END DO
      DO JJ=1,JPRIM
      PHI(      81)=ZERO
      PHI(      93)=ZERO
      PHI(      94)=ZERO
      PHI(      95)=ZERO
      DO I=     123,     129
      PHI(I)=ZERO
      END DO
      DO I=     186,     198
      PHI(I)=ZERO
      END DO
      DO I=     277,     298
      PHI(I)=ZERO
      END DO
      DO I=     409,     442
      PHI(I)=ZERO
      END DO
      ITOP=IPRIM
      IF (IEQJ) ITOP=JJ
      DO II=1,ITOP
      IJ=IJ+1
      CCBF=CCBRA(IJ)
      IF (ABS(CCBF).GT.CCTOL) THEN
      EIJ=EXJ(JJ)+EXI(II)
      EIJKL=EIJ+EKL
      EINV=ONE/EIJKL
      RX=XKL-XAB(IJ)
      RY=YKL-YAB(IJ)
      RZ=ZKL-ZAB(IJ)
      R2=RX*RX+RY*RY+RZ*RZ
      VS=EIJ*EKL*EINV
      IF(LRINT) THEN
         VS0=VS
         VS=VS0*EMU2/(VS0+EMU2)
      ENDIF
      TT=VS*R2
      IF (TT.LE.TMAX) THEN
      TVAL=TT*RFAC
      IPT=NINT(TVAL)
      FX= TVAL*FGRID(4,IPT, 5)
      FX=(FX + FGRID(3,IPT, 5))*TVAL
      FX=(FX + FGRID(2,IPT, 5))*TVAL
      FX=(FX + FGRID(1,IPT, 5))*TVAL
      FX= FX + FGRID(0,IPT, 5)
      FT( 8)=FX
      TVAL=TT*RXINC
      IPT=NINT(TVAL)
      ET= TVAL*XGRID(4,IPT)
      ET=(ET + XGRID(3,IPT))*TVAL
      ET=(ET + XGRID(2,IPT))*TVAL
      ET=(ET + XGRID(1,IPT))*TVAL
      ET= ET + XGRID(0,IPT)
      T2=TT*TWO
      FT( 7)=RMR( 8)*(T2*FT( 8)+ET)
      FT( 6)=RMR( 7)*(T2*FT( 7)+ET)
      FT( 5)=RMR( 6)*(T2*FT( 6)+ET)
      FT( 4)=RMR( 5)*(T2*FT( 5)+ET)
      FT( 3)=RMR( 4)*(T2*FT( 4)+ET)
      FT( 2)=RMR( 3)*(T2*FT( 3)+ET)
      FT( 1)=RMR( 2)*(T2*FT( 2)+ET)
      FT( 0)=RMR( 1)*(T2*FT( 1)+ET)
      ZEROL=CCBF*CCKF*SQRT(EINV)*QQ4
      IF(LRINT)ZEROL=ZEROL*SQRT(VS/VS0)
      SAUX(0)=ZEROL*FT(0)
      VS2=TWO*VS
      VSM=ONE
      VSM=VSM*VS2
      SAUX( 1)=ZEROL*FT( 1)*VSM
      VSM=VSM*VS2
      SAUX( 2)=ZEROL*FT( 2)*VSM
      VSM=VSM*VS2
      SAUX( 3)=ZEROL*FT( 3)*VSM
      VSM=VSM*VS2
      SAUX( 4)=ZEROL*FT( 4)*VSM
      VSM=VSM*VS2
      SAUX( 5)=ZEROL*FT( 5)*VSM
      ELSE
      ZEROL=CCBF*CCKF*SLBRA(IJ)*SLKET(KL)*QQ4
      RINV=ONE/R2
      SAUX(0)=ZEROL*SQRT(RINV)
      RR=SAUX(0)
      RR=RR*RINV
      SAUX( 1)=TLGM( 1)*RR
      RR=RR*RINV
      SAUX( 2)=TLGM( 2)*RR
      RR=RR*RINV
      SAUX( 3)=TLGM( 3)*RR
      RR=RR*RINV
      SAUX( 4)=TLGM( 4)*RR
      RR=RR*RINV
      SAUX( 5)=TLGM( 5)*RR
      END IF
      RXYZ(1)=RX
      RXYZ(2)=RY
      RXYZ(3)=RZ
      CXYZ(1)=ONE
      CXYZ(2)=RX
      CXYZ(3)=RY
      CXYZ(4)=RZ
      CXYZ(     5)=CXYZ(     2)*RXYZ(1)
      CXYZ(     6)=CXYZ(     3)*RXYZ(1)
      CXYZ(     7)=CXYZ(     3)*RXYZ(2)
      CXYZ(     8)=CXYZ(     4)*RXYZ(1)
      CXYZ(     9)=CXYZ(     4)*RXYZ(2)
      CXYZ(    10)=CXYZ(     4)*RXYZ(3)
      CXYZ(    11)=CXYZ(     5)*RXYZ(1)
      CXYZ(    12)=CXYZ(     6)*RXYZ(1)
      CXYZ(    13)=CXYZ(     7)*RXYZ(1)
      CXYZ(    14)=CXYZ(     7)*RXYZ(2)
      CXYZ(    15)=CXYZ(     8)*RXYZ(1)
      CXYZ(    16)=CXYZ(     9)*RXYZ(1)
      CXYZ(    17)=CXYZ(     9)*RXYZ(2)
      CXYZ(    18)=CXYZ(    10)*RXYZ(1)
      CXYZ(    19)=CXYZ(    10)*RXYZ(2)
      CXYZ(    20)=CXYZ(    10)*RXYZ(3)
      CXYZ(    21)=CXYZ(    11)*RXYZ(1)
      CXYZ(    22)=CXYZ(    12)*RXYZ(1)
      CXYZ(    23)=CXYZ(    13)*RXYZ(1)
      CXYZ(    24)=CXYZ(    14)*RXYZ(1)
      CXYZ(    25)=CXYZ(    14)*RXYZ(2)
      CXYZ(    26)=CXYZ(    15)*RXYZ(1)
      CXYZ(    27)=CXYZ(    16)*RXYZ(1)
      CXYZ(    28)=CXYZ(    17)*RXYZ(1)
      CXYZ(    29)=CXYZ(    17)*RXYZ(2)
      CXYZ(    30)=CXYZ(    18)*RXYZ(1)
      CXYZ(    31)=CXYZ(    19)*RXYZ(1)
      CXYZ(    32)=CXYZ(    19)*RXYZ(2)
      CXYZ(    33)=CXYZ(    20)*RXYZ(1)
      CXYZ(    34)=CXYZ(    20)*RXYZ(2)
      CXYZ(    35)=CXYZ(    20)*RXYZ(3)
      CXYZ(    36)=CXYZ(    21)*RXYZ(1)
      CXYZ(    37)=CXYZ(    22)*RXYZ(1)
      CXYZ(    38)=CXYZ(    23)*RXYZ(1)
      CXYZ(    39)=CXYZ(    24)*RXYZ(1)
      CXYZ(    40)=CXYZ(    25)*RXYZ(1)
      CXYZ(    41)=CXYZ(    25)*RXYZ(2)
      CXYZ(    42)=CXYZ(    26)*RXYZ(1)
      CXYZ(    43)=CXYZ(    27)*RXYZ(1)
      CXYZ(    44)=CXYZ(    28)*RXYZ(1)
      CXYZ(    45)=CXYZ(    29)*RXYZ(1)
      CXYZ(    46)=CXYZ(    29)*RXYZ(2)
      CXYZ(    47)=CXYZ(    30)*RXYZ(1)
      CXYZ(    48)=CXYZ(    31)*RXYZ(1)
      CXYZ(    49)=CXYZ(    32)*RXYZ(1)
      CXYZ(    50)=CXYZ(    32)*RXYZ(2)
      CXYZ(    51)=CXYZ(    33)*RXYZ(1)
      CXYZ(    52)=CXYZ(    34)*RXYZ(1)
      CXYZ(    53)=CXYZ(    34)*RXYZ(2)
      CXYZ(    54)=CXYZ(    35)*RXYZ(1)
      CXYZ(    55)=CXYZ(    35)*RXYZ(2)
      CXYZ(    56)=CXYZ(    35)*RXYZ(3)
      PHI(1)=SAUX(0)
      PHI(       2)=CXYZ(     2)*SAUX( 1)
      PHI(       3)=CXYZ(     3)*SAUX( 1)
      PHI(       4)=CXYZ(     4)*SAUX( 1)
      PHI(       5)=CXYZ(     1)*SAUX( 1)
      PHI(       6)=CXYZ(     5)*SAUX( 2)
      PHI(       7)=CXYZ(     6)*SAUX( 2)
      PHI(       8)=CXYZ(     7)*SAUX( 2)
      PHI(       9)=CXYZ(     8)*SAUX( 2)
      PHI(      10)=CXYZ(     9)*SAUX( 2)
      PHI(      11)=CXYZ(    10)*SAUX( 2)
      PHI(      12)=CXYZ(     2)*SAUX( 2)
      PHI(      13)=CXYZ(     3)*SAUX( 2)
      PHI(      14)=CXYZ(     4)*SAUX( 2)
      PHI(      15)=CXYZ(    11)*SAUX( 3)
      PHI(      16)=CXYZ(    12)*SAUX( 3)
      PHI(      17)=CXYZ(    13)*SAUX( 3)
      PHI(      18)=CXYZ(    14)*SAUX( 3)
      PHI(      19)=CXYZ(    15)*SAUX( 3)
      PHI(      20)=CXYZ(    16)*SAUX( 3)
      PHI(      21)=CXYZ(    17)*SAUX( 3)
      PHI(      22)=CXYZ(    18)*SAUX( 3)
      PHI(      23)=CXYZ(    19)*SAUX( 3)
      PHI(      24)=CXYZ(    20)*SAUX( 3)
      PHI(      25)=CXYZ(     1)*SAUX( 2)
      PHI(      26)=CXYZ(     5)*SAUX( 3)
      PHI(      27)=CXYZ(     6)*SAUX( 3)
      PHI(      28)=CXYZ(     7)*SAUX( 3)
      PHI(      29)=CXYZ(     8)*SAUX( 3)
      PHI(      30)=CXYZ(     9)*SAUX( 3)
      PHI(      31)=CXYZ(    10)*SAUX( 3)
      PHI(      32)=CXYZ(    21)*SAUX( 4)
      PHI(      33)=CXYZ(    22)*SAUX( 4)
      PHI(      34)=CXYZ(    23)*SAUX( 4)
      PHI(      35)=CXYZ(    24)*SAUX( 4)
      PHI(      36)=CXYZ(    25)*SAUX( 4)
      PHI(      37)=CXYZ(    26)*SAUX( 4)
      PHI(      38)=CXYZ(    27)*SAUX( 4)
      PHI(      39)=CXYZ(    28)*SAUX( 4)
      PHI(      40)=CXYZ(    29)*SAUX( 4)
      PHI(      41)=CXYZ(    30)*SAUX( 4)
      PHI(      42)=CXYZ(    31)*SAUX( 4)
      PHI(      43)=CXYZ(    32)*SAUX( 4)
      PHI(      44)=CXYZ(    33)*SAUX( 4)
      PHI(      45)=CXYZ(    34)*SAUX( 4)
      PHI(      46)=CXYZ(    35)*SAUX( 4)
      PHI(      47)=CXYZ(     2)*SAUX( 3)
      PHI(      48)=CXYZ(     3)*SAUX( 3)
      PHI(      49)=CXYZ(     4)*SAUX( 3)
      PHI(      50)=CXYZ(    11)*SAUX( 4)
      PHI(      51)=CXYZ(    12)*SAUX( 4)
      PHI(      52)=CXYZ(    13)*SAUX( 4)
      PHI(      53)=CXYZ(    14)*SAUX( 4)
      PHI(      54)=CXYZ(    15)*SAUX( 4)
      PHI(      55)=CXYZ(    16)*SAUX( 4)
      PHI(      56)=CXYZ(    17)*SAUX( 4)
      PHI(      57)=CXYZ(    18)*SAUX( 4)
      PHI(      58)=CXYZ(    19)*SAUX( 4)
      PHI(      59)=CXYZ(    20)*SAUX( 4)
      PHI(      60)=CXYZ(    36)*SAUX( 5)
      PHI(      61)=CXYZ(    37)*SAUX( 5)
      PHI(      62)=CXYZ(    38)*SAUX( 5)
      PHI(      63)=CXYZ(    39)*SAUX( 5)
      PHI(      64)=CXYZ(    40)*SAUX( 5)
      PHI(      65)=CXYZ(    41)*SAUX( 5)
      PHI(      66)=CXYZ(    42)*SAUX( 5)
      PHI(      67)=CXYZ(    43)*SAUX( 5)
      PHI(      68)=CXYZ(    44)*SAUX( 5)
      PHI(      69)=CXYZ(    45)*SAUX( 5)
      PHI(      70)=CXYZ(    46)*SAUX( 5)
      PHI(      71)=CXYZ(    47)*SAUX( 5)
      PHI(      72)=CXYZ(    48)*SAUX( 5)
      PHI(      73)=CXYZ(    49)*SAUX( 5)
      PHI(      74)=CXYZ(    50)*SAUX( 5)
      PHI(      75)=CXYZ(    51)*SAUX( 5)
      PHI(      76)=CXYZ(    52)*SAUX( 5)
      PHI(      77)=CXYZ(    53)*SAUX( 5)
      PHI(      78)=CXYZ(    54)*SAUX( 5)
      PHI(      79)=CXYZ(    55)*SAUX( 5)
      PHI(      80)=CXYZ(    56)*SAUX( 5)
      FAC=RXB(IJ)
      SF=ONE
      SFAC(1)=SF
      FAC=SFAC( 1)
      PHI(      81)=PHI(      81)+FAC*PHI(       1)
      FAC=SFAC( 1)
      PHI(      93)=PHI(      93)+FAC*PHI(       2)
      PHI(      94)=PHI(      94)+FAC*PHI(       3)
      PHI(      95)=PHI(      95)+FAC*PHI(       4)
      FAC=SFAC( 1)
      J=       5
      DO I=     123,     129
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=      12
      DO I=     186,     198
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=      25
      DO I=     277,     298
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=      47
      DO I=     409,     442
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      END IF
      END DO
      FAC=EXJ(JJ)*TWO
      SF=ONE
      SFAC(1)=SF
      FAC=SFAC( 1)
      PHI(      82)=PHI(      82)+FAC*PHI(      81)
      FAC=SFAC( 1)
      PHI(      96)=PHI(      96)+FAC*PHI(      93)
      PHI(      97)=PHI(      97)+FAC*PHI(      94)
      PHI(      98)=PHI(      98)+FAC*PHI(      95)
      FAC=SFAC( 1)
      J=     123
      DO I=     130,     136
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     186
      DO I=     199,     211
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     277
      DO I=     299,     320
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     409
      DO I=     443,     476
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      END DO
      FAC=RXK(KL)
      SF=ONE
      SFAC(1)=SF
      SF=SF*FAC
      SFAC( 2)=SF
      SF=SF*FAC
      SFAC( 3)=SF
      SF=SF*FAC
      SFAC( 4)=SF
      SF=SF*FAC
      SFAC( 5)=SF
      SF=SF*FAC
      SFAC( 6)=SF
      FAC=SFAC( 3)
      PHI(      83)=PHI(      83)+FAC*PHI(      82)
      FAC=SFAC( 4)
      PHI(      85)=PHI(      85)+FAC*PHI(      82)
      FAC=SFAC( 5)
      PHI(      88)=PHI(      88)+FAC*PHI(      82)
      FAC=SFAC( 6)
      PHI(      91)=PHI(      91)+FAC*PHI(      82)
      FAC=SFAC( 4)
      PHI(      99)=PHI(      99)+FAC*PHI(      96)
      PHI(     100)=PHI(     100)+FAC*PHI(      97)
      PHI(     101)=PHI(     101)+FAC*PHI(      98)
      FAC=SFAC( 5)
      PHI(     108)=PHI(     108)+FAC*PHI(      96)
      PHI(     109)=PHI(     109)+FAC*PHI(      97)
      PHI(     110)=PHI(     110)+FAC*PHI(      98)
      FAC=SFAC( 6)
      PHI(     117)=PHI(     117)+FAC*PHI(      96)
      PHI(     118)=PHI(     118)+FAC*PHI(      97)
      PHI(     119)=PHI(     119)+FAC*PHI(      98)
      FAC=SFAC( 4)
      J=     130
      DO I=     137,     143
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 5)
      J=     130
      DO I=     151,     157
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 6)
      J=     130
      DO I=     172,     178
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 5)
      J=     199
      DO I=     212,     224
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 6)
      J=     199
      DO I=     251,     263
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 5)
      J=     299
      DO I=     321,     342
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 6)
      J=     299
      DO I=     365,     386
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 6)
      J=     443
      DO I=     477,     510
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      END IF
      END DO
      FAC=EXL(LL)*TWO
      SF=ONE
      SFAC(1)=SF
      SF=SF*FAC
      SFAC( 2)=SF
      SF=SF*FAC
      SFAC( 3)=SF
      SF=SF*FAC
      SFAC( 4)=SF
      SF=SF*FAC
      SFAC( 5)=SF
      SF=SF*FAC
      SFAC( 6)=SF
      FAC=SFAC( 1)
      PHI(      84)=PHI(      84)+FAC*PHI(      83)
      FAC=SFAC( 2)
      PHI(      86)=PHI(      86)+FAC*PHI(      85)
      FAC=SFAC( 3)
      PHI(      87)=PHI(      87)+FAC*PHI(      85)
      FAC=SFAC( 4)
      PHI(      89)=PHI(      89)+FAC*PHI(      88)
      FAC=SFAC( 5)
      PHI(      90)=PHI(      90)+FAC*PHI(      88)
      FAC=SFAC( 6)
      PHI(      92)=PHI(      92)+FAC*PHI(      91)
      FAC=SFAC( 1)
      PHI(     102)=PHI(     102)+FAC*PHI(      99)
      PHI(     103)=PHI(     103)+FAC*PHI(     100)
      PHI(     104)=PHI(     104)+FAC*PHI(     101)
      FAC=SFAC( 2)
      PHI(     105)=PHI(     105)+FAC*PHI(      99)
      PHI(     106)=PHI(     106)+FAC*PHI(     100)
      PHI(     107)=PHI(     107)+FAC*PHI(     101)
      FAC=SFAC( 3)
      PHI(     111)=PHI(     111)+FAC*PHI(     108)
      PHI(     112)=PHI(     112)+FAC*PHI(     109)
      PHI(     113)=PHI(     113)+FAC*PHI(     110)
      FAC=SFAC( 4)
      PHI(     114)=PHI(     114)+FAC*PHI(     108)
      PHI(     115)=PHI(     115)+FAC*PHI(     109)
      PHI(     116)=PHI(     116)+FAC*PHI(     110)
      FAC=SFAC( 5)
      PHI(     120)=PHI(     120)+FAC*PHI(     117)
      PHI(     121)=PHI(     121)+FAC*PHI(     118)
      PHI(     122)=PHI(     122)+FAC*PHI(     119)
      FAC=SFAC( 1)
      J=     137
      DO I=     144,     150
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 2)
      J=     151
      DO I=     158,     164
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 3)
      J=     151
      DO I=     165,     171
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 4)
      J=     172
      DO I=     179,     185
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     212
      DO I=     225,     237
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 2)
      J=     212
      DO I=     238,     250
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 3)
      J=     251
      DO I=     264,     276
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     321
      DO I=     343,     364
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 2)
      J=     365
      DO I=     387,     408
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      FAC=SFAC( 1)
      J=     477
      DO I=     511,     544
      PHI(I)=PHI(I)+FAC*PHI(J)
      J=J+1
      END DO
      END DO
C  POST-CONTRACTION PHASE
      PHI(     145)=PHI(     145)-D01*PHI(     144)
      PHI(     147)=PHI(     147)-D01*PHI(     144)
      PHI(     150)=PHI(     150)-D01*PHI(     144)
      PHI(     159)=PHI(     159)-D01*PHI(     158)
      PHI(     161)=PHI(     161)-D01*PHI(     158)
      PHI(     164)=PHI(     164)-D01*PHI(     158)
      PHI(     166)=PHI(     166)-D01*PHI(     165)
      PHI(     168)=PHI(     168)-D01*PHI(     165)
      PHI(     171)=PHI(     171)-D01*PHI(     165)
      PHI(     180)=PHI(     180)-D01*PHI(     179)
      PHI(     182)=PHI(     182)-D01*PHI(     179)
      PHI(     185)=PHI(     185)-D01*PHI(     179)
      PHI(     228)=PHI(     228)-D01*PHI(     225)
      PHI(     229)=PHI(     229)-D01*PHI(     226)
      PHI(     231)=PHI(     231)-D01*PHI(     226)
      PHI(     232)=PHI(     232)-D01*PHI(     227)
      PHI(     234)=PHI(     234)-D01*PHI(     227)
      PHI(     237)=PHI(     237)-D01*PHI(     227)
      PHI(     228)=PHI(     228)-D02*PHI(     225)
      PHI(     230)=PHI(     230)-D01*PHI(     225)
      PHI(     231)=PHI(     231)-D02*PHI(     226)
      PHI(     235)=PHI(     235)-D01*PHI(     225)
      PHI(     236)=PHI(     236)-D01*PHI(     226)
      PHI(     237)=PHI(     237)-D02*PHI(     227)
      PHI(     241)=PHI(     241)-D01*PHI(     238)
      PHI(     242)=PHI(     242)-D01*PHI(     239)
      PHI(     244)=PHI(     244)-D01*PHI(     239)
      PHI(     245)=PHI(     245)-D01*PHI(     240)
      PHI(     247)=PHI(     247)-D01*PHI(     240)
      PHI(     250)=PHI(     250)-D01*PHI(     240)
      PHI(     241)=PHI(     241)-D02*PHI(     238)
      PHI(     243)=PHI(     243)-D01*PHI(     238)
      PHI(     244)=PHI(     244)-D02*PHI(     239)
      PHI(     248)=PHI(     248)-D01*PHI(     238)
      PHI(     249)=PHI(     249)-D01*PHI(     239)
      PHI(     250)=PHI(     250)-D02*PHI(     240)
      PHI(     267)=PHI(     267)-D01*PHI(     264)
      PHI(     268)=PHI(     268)-D01*PHI(     265)
      PHI(     270)=PHI(     270)-D01*PHI(     265)
      PHI(     271)=PHI(     271)-D01*PHI(     266)
      PHI(     273)=PHI(     273)-D01*PHI(     266)
      PHI(     276)=PHI(     276)-D01*PHI(     266)
      PHI(     267)=PHI(     267)-D02*PHI(     264)
      PHI(     269)=PHI(     269)-D01*PHI(     264)
      PHI(     270)=PHI(     270)-D02*PHI(     265)
      PHI(     274)=PHI(     274)-D01*PHI(     264)
      PHI(     275)=PHI(     275)-D01*PHI(     265)
      PHI(     276)=PHI(     276)-D02*PHI(     266)
      PHI(     350)=PHI(     350)-D01*PHI(     344)
      PHI(     351)=PHI(     351)-D01*PHI(     345)
      PHI(     352)=PHI(     352)-D01*PHI(     346)
      PHI(     354)=PHI(     354)-D01*PHI(     346)
      PHI(     355)=PHI(     355)-D01*PHI(     347)
      PHI(     356)=PHI(     356)-D01*PHI(     348)
      PHI(     358)=PHI(     358)-D01*PHI(     348)
      PHI(     359)=PHI(     359)-D01*PHI(     349)
      PHI(     361)=PHI(     361)-D01*PHI(     349)
      PHI(     364)=PHI(     364)-D01*PHI(     349)
      PHI(     350)=PHI(     350)-D02*PHI(     344)
      PHI(     351)=PHI(     351)-D02*PHI(     345)
      PHI(     353)=PHI(     353)-D01*PHI(     345)
      PHI(     354)=PHI(     354)-D02*PHI(     346)
      PHI(     355)=PHI(     355)-D02*PHI(     347)
      PHI(     357)=PHI(     357)-D01*PHI(     347)
      PHI(     358)=PHI(     358)-D02*PHI(     348)
      PHI(     362)=PHI(     362)-D01*PHI(     347)
      PHI(     363)=PHI(     363)-D01*PHI(     348)
      PHI(     364)=PHI(     364)-D02*PHI(     349)
      PHI(     344)=PHI(     344)-D01*PHI(     343)
      PHI(     346)=PHI(     346)-D01*PHI(     343)
      PHI(     349)=PHI(     349)-D01*PHI(     343)
      PHI(     350)=PHI(     350)-D03*PHI(     344)
      PHI(     352)=PHI(     352)-D01*PHI(     344)
      PHI(     353)=PHI(     353)-D02*PHI(     345)
      PHI(     354)=PHI(     354)-D03*PHI(     346)
      PHI(     359)=PHI(     359)-D01*PHI(     344)
      PHI(     360)=PHI(     360)-D01*PHI(     345)
      PHI(     361)=PHI(     361)-D01*PHI(     346)
      PHI(     362)=PHI(     362)-D02*PHI(     347)
      PHI(     363)=PHI(     363)-D02*PHI(     348)
      PHI(     364)=PHI(     364)-D03*PHI(     349)
      PHI(     394)=PHI(     394)-D01*PHI(     388)
      PHI(     395)=PHI(     395)-D01*PHI(     389)
      PHI(     396)=PHI(     396)-D01*PHI(     390)
      PHI(     398)=PHI(     398)-D01*PHI(     390)
      PHI(     399)=PHI(     399)-D01*PHI(     391)
      PHI(     400)=PHI(     400)-D01*PHI(     392)
      PHI(     402)=PHI(     402)-D01*PHI(     392)
      PHI(     403)=PHI(     403)-D01*PHI(     393)
      PHI(     405)=PHI(     405)-D01*PHI(     393)
      PHI(     408)=PHI(     408)-D01*PHI(     393)
      PHI(     394)=PHI(     394)-D02*PHI(     388)
      PHI(     395)=PHI(     395)-D02*PHI(     389)
      PHI(     397)=PHI(     397)-D01*PHI(     389)
      PHI(     398)=PHI(     398)-D02*PHI(     390)
      PHI(     399)=PHI(     399)-D02*PHI(     391)
      PHI(     401)=PHI(     401)-D01*PHI(     391)
      PHI(     402)=PHI(     402)-D02*PHI(     392)
      PHI(     406)=PHI(     406)-D01*PHI(     391)
      PHI(     407)=PHI(     407)-D01*PHI(     392)
      PHI(     408)=PHI(     408)-D02*PHI(     393)
      PHI(     388)=PHI(     388)-D01*PHI(     387)
      PHI(     390)=PHI(     390)-D01*PHI(     387)
      PHI(     393)=PHI(     393)-D01*PHI(     387)
      PHI(     394)=PHI(     394)-D03*PHI(     388)
      PHI(     396)=PHI(     396)-D01*PHI(     388)
      PHI(     397)=PHI(     397)-D02*PHI(     389)
      PHI(     398)=PHI(     398)-D03*PHI(     390)
      PHI(     403)=PHI(     403)-D01*PHI(     388)
      PHI(     404)=PHI(     404)-D01*PHI(     389)
      PHI(     405)=PHI(     405)-D01*PHI(     390)
      PHI(     406)=PHI(     406)-D02*PHI(     391)
      PHI(     407)=PHI(     407)-D02*PHI(     392)
      PHI(     408)=PHI(     408)-D03*PHI(     393)
      PHI(     524)=PHI(     524)-D01*PHI(     514)
      PHI(     525)=PHI(     525)-D01*PHI(     515)
      PHI(     526)=PHI(     526)-D01*PHI(     516)
      PHI(     527)=PHI(     527)-D01*PHI(     517)
      PHI(     529)=PHI(     529)-D01*PHI(     517)
      PHI(     530)=PHI(     530)-D01*PHI(     518)
      PHI(     531)=PHI(     531)-D01*PHI(     519)
      PHI(     532)=PHI(     532)-D01*PHI(     520)
      PHI(     534)=PHI(     534)-D01*PHI(     520)
      PHI(     535)=PHI(     535)-D01*PHI(     521)
      PHI(     536)=PHI(     536)-D01*PHI(     522)
      PHI(     538)=PHI(     538)-D01*PHI(     522)
      PHI(     539)=PHI(     539)-D01*PHI(     523)
      PHI(     541)=PHI(     541)-D01*PHI(     523)
      PHI(     544)=PHI(     544)-D01*PHI(     523)
      PHI(     524)=PHI(     524)-D02*PHI(     514)
      PHI(     525)=PHI(     525)-D02*PHI(     515)
      PHI(     526)=PHI(     526)-D02*PHI(     516)
      PHI(     528)=PHI(     528)-D01*PHI(     516)
      PHI(     529)=PHI(     529)-D02*PHI(     517)
      PHI(     530)=PHI(     530)-D02*PHI(     518)
      PHI(     531)=PHI(     531)-D02*PHI(     519)
      PHI(     533)=PHI(     533)-D01*PHI(     519)
      PHI(     534)=PHI(     534)-D02*PHI(     520)
      PHI(     535)=PHI(     535)-D02*PHI(     521)
      PHI(     537)=PHI(     537)-D01*PHI(     521)
      PHI(     538)=PHI(     538)-D02*PHI(     522)
      PHI(     542)=PHI(     542)-D01*PHI(     521)
      PHI(     543)=PHI(     543)-D01*PHI(     522)
      PHI(     544)=PHI(     544)-D02*PHI(     523)
      PHI(     514)=PHI(     514)-D01*PHI(     511)
      PHI(     515)=PHI(     515)-D01*PHI(     512)
      PHI(     517)=PHI(     517)-D01*PHI(     512)
      PHI(     518)=PHI(     518)-D01*PHI(     513)
      PHI(     520)=PHI(     520)-D01*PHI(     513)
      PHI(     523)=PHI(     523)-D01*PHI(     513)
      PHI(     524)=PHI(     524)-D03*PHI(     514)
      PHI(     525)=PHI(     525)-D03*PHI(     515)
      PHI(     527)=PHI(     527)-D01*PHI(     515)
      PHI(     528)=PHI(     528)-D02*PHI(     516)
      PHI(     529)=PHI(     529)-D03*PHI(     517)
      PHI(     530)=PHI(     530)-D03*PHI(     518)
      PHI(     532)=PHI(     532)-D01*PHI(     518)
      PHI(     533)=PHI(     533)-D02*PHI(     519)
      PHI(     534)=PHI(     534)-D03*PHI(     520)
      PHI(     539)=PHI(     539)-D01*PHI(     518)
      PHI(     540)=PHI(     540)-D01*PHI(     519)
      PHI(     541)=PHI(     541)-D01*PHI(     520)
      PHI(     542)=PHI(     542)-D02*PHI(     521)
      PHI(     543)=PHI(     543)-D02*PHI(     522)
      PHI(     544)=PHI(     544)-D03*PHI(     523)
      PHI(     514)=PHI(     514)-D02*PHI(     511)
      PHI(     516)=PHI(     516)-D01*PHI(     511)
      PHI(     517)=PHI(     517)-D02*PHI(     512)
      PHI(     521)=PHI(     521)-D01*PHI(     511)
      PHI(     522)=PHI(     522)-D01*PHI(     512)
      PHI(     523)=PHI(     523)-D02*PHI(     513)
      PHI(     524)=PHI(     524)-D04*PHI(     514)
      PHI(     526)=PHI(     526)-D01*PHI(     514)
      PHI(     527)=PHI(     527)-D02*PHI(     515)
      PHI(     528)=PHI(     528)-D03*PHI(     516)
      PHI(     529)=PHI(     529)-D04*PHI(     517)
      PHI(     535)=PHI(     535)-D01*PHI(     514)
      PHI(     536)=PHI(     536)-D01*PHI(     515)
      PHI(     537)=PHI(     537)-D01*PHI(     516)
      PHI(     538)=PHI(     538)-D01*PHI(     517)
      PHI(     539)=PHI(     539)-D02*PHI(     518)
      PHI(     540)=PHI(     540)-D02*PHI(     519)
      PHI(     541)=PHI(     541)-D02*PHI(     520)
      PHI(     542)=PHI(     542)-D03*PHI(     521)
      PHI(     543)=PHI(     543)-D03*PHI(     522)
      PHI(     544)=PHI(     544)-D04*PHI(     523)
      WRK2(    1,1)= PHI(      90)
      WRK2(    2,1)= PHI(      87)
      WRK2(    3,1)= PHI(      84)
      WRK2(    4,1)=-PHI(     114)
      WRK2(    5,1)=-PHI(     115)
      WRK2(    6,1)=-PHI(     116)
      WRK2(    7,1)=-PHI(     105)
      WRK2(    8,1)=-PHI(     106)
      WRK2(    9,1)=-PHI(     107)
      WRK2(   10,1)= PHI(     166)
      WRK2(   11,1)= PHI(     167)
      WRK2(   12,1)= PHI(     168)
      WRK2(   13,1)= PHI(     169)
      WRK2(   14,1)= PHI(     170)
      WRK2(   15,1)= PHI(     171)
      WRK2(   16,1)= PHI(     145)
      WRK2(   17,1)= PHI(     146)
      WRK2(   18,1)= PHI(     147)
      WRK2(   19,1)= PHI(     148)
      WRK2(   20,1)= PHI(     149)
      WRK2(   21,1)= PHI(     150)
      WRK2(   22,1)=-PHI(     241)
      WRK2(   23,1)=-PHI(     242)
      WRK2(   24,1)=-PHI(     243)
      WRK2(   25,1)=-PHI(     244)
      WRK2(   26,1)=-PHI(     245)
      WRK2(   27,1)=-PHI(     246)
      WRK2(   28,1)=-PHI(     247)
      WRK2(   29,1)=-PHI(     248)
      WRK2(   30,1)=-PHI(     249)
      WRK2(   31,1)=-PHI(     250)
      WRK2(   32,1)= PHI(     350)
      WRK2(   33,1)= PHI(     351)
      WRK2(   34,1)= PHI(     352)
      WRK2(   35,1)= PHI(     353)
      WRK2(   36,1)= PHI(     354)
      WRK2(   37,1)= PHI(     355)
      WRK2(   38,1)= PHI(     356)
      WRK2(   39,1)= PHI(     357)
      WRK2(   40,1)= PHI(     358)
      WRK2(   41,1)= PHI(     359)
      WRK2(   42,1)= PHI(     360)
      WRK2(   43,1)= PHI(     361)
      WRK2(   44,1)= PHI(     362)
      WRK2(   45,1)= PHI(     363)
      WRK2(   46,1)= PHI(     364)
      WRK2(   47,1)= PHI(      92)
      WRK2(   48,1)= PHI(      89)
      WRK2(   49,1)= PHI(      86)
      WRK2(   50,1)=-PHI(     120)
      WRK2(   51,1)=-PHI(     121)
      WRK2(   52,1)=-PHI(     122)
      WRK2(   53,1)=-PHI(     111)
      WRK2(   54,1)=-PHI(     112)
      WRK2(   55,1)=-PHI(     113)
      WRK2(   56,1)=-PHI(     102)
      WRK2(   57,1)=-PHI(     103)
      WRK2(   58,1)=-PHI(     104)
      WRK2(   59,1)= PHI(     180)
      WRK2(   60,1)= PHI(     181)
      WRK2(   61,1)= PHI(     182)
      WRK2(   62,1)= PHI(     183)
      WRK2(   63,1)= PHI(     184)
      WRK2(   64,1)= PHI(     185)
      WRK2(   65,1)= PHI(     159)
      WRK2(   66,1)= PHI(     160)
      WRK2(   67,1)= PHI(     161)
      WRK2(   68,1)= PHI(     162)
      WRK2(   69,1)= PHI(     163)
      WRK2(   70,1)= PHI(     164)
      WRK2(   71,1)=-PHI(     267)
      WRK2(   72,1)=-PHI(     268)
      WRK2(   73,1)=-PHI(     269)
      WRK2(   74,1)=-PHI(     270)
      WRK2(   75,1)=-PHI(     271)
      WRK2(   76,1)=-PHI(     272)
      WRK2(   77,1)=-PHI(     273)
      WRK2(   78,1)=-PHI(     274)
      WRK2(   79,1)=-PHI(     275)
      WRK2(   80,1)=-PHI(     276)
      WRK2(   81,1)=-PHI(     228)
      WRK2(   82,1)=-PHI(     229)
      WRK2(   83,1)=-PHI(     230)
      WRK2(   84,1)=-PHI(     231)
      WRK2(   85,1)=-PHI(     232)
      WRK2(   86,1)=-PHI(     233)
      WRK2(   87,1)=-PHI(     234)
      WRK2(   88,1)=-PHI(     235)
      WRK2(   89,1)=-PHI(     236)
      WRK2(   90,1)=-PHI(     237)
      WRK2(   91,1)= PHI(     394)
      WRK2(   92,1)= PHI(     395)
      WRK2(   93,1)= PHI(     396)
      WRK2(   94,1)= PHI(     397)
      WRK2(   95,1)= PHI(     398)
      WRK2(   96,1)= PHI(     399)
      WRK2(   97,1)= PHI(     400)
      WRK2(   98,1)= PHI(     401)
      WRK2(   99,1)= PHI(     402)
      WRK2(  100,1)= PHI(     403)
      WRK2(  101,1)= PHI(     404)
      WRK2(  102,1)= PHI(     405)
      WRK2(  103,1)= PHI(     406)
      WRK2(  104,1)= PHI(     407)
      WRK2(  105,1)= PHI(     408)
      WRK2(  106,1)=-PHI(     524)
      WRK2(  107,1)=-PHI(     525)
      WRK2(  108,1)=-PHI(     526)
      WRK2(  109,1)=-PHI(     527)
      WRK2(  110,1)=-PHI(     528)
      WRK2(  111,1)=-PHI(     529)
      WRK2(  112,1)=-PHI(     530)
      WRK2(  113,1)=-PHI(     531)
      WRK2(  114,1)=-PHI(     532)
      WRK2(  115,1)=-PHI(     533)
      WRK2(  116,1)=-PHI(     534)
      WRK2(  117,1)=-PHI(     535)
      WRK2(  118,1)=-PHI(     536)
      WRK2(  119,1)=-PHI(     537)
      WRK2(  120,1)=-PHI(     538)
      WRK2(  121,1)=-PHI(     539)
      WRK2(  122,1)=-PHI(     540)
      WRK2(  123,1)=-PHI(     541)
      WRK2(  124,1)=-PHI(     542)
      WRK2(  125,1)=-PHI(     543)
      WRK2(  126,1)=-PHI(     544)
      CONF(1)=XD-XC
      CONF(2)=YD-YC
      CONF(3)=ZD-ZC
      DO I=1,     1
      WRK2(   32,I)=WRK2(   32,I)+CONF( 1)*WRK2(   22,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 6)*WRK2(   16,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 1)*WRK2(   23,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 5)*WRK2(   17,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 1)*WRK2(   24,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 4)*WRK2(   18,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 1)*WRK2(   25,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 2)*WRK2(   25,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 6)*WRK2(   18,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 1)*WRK2(   26,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 5)*WRK2(   19,I)
      WRK2(   38,I)=WRK2(   38,I)+CONF( 1)*WRK2(   27,I)
      WRK2(   38,I)=WRK2(   38,I)+CONF( 4)*WRK2(   20,I)
      WRK2(   39,I)=WRK2(   39,I)+CONF( 1)*WRK2(   28,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 2)*WRK2(   28,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 5)*WRK2(   20,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 1)*WRK2(   29,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 4)*WRK2(   21,I)
      WRK2(   42,I)=WRK2(   42,I)+CONF( 1)*WRK2(   30,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 2)*WRK2(   30,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 4)*WRK2(   21,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 1)*WRK2(   31,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 2)*WRK2(   31,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 3)*WRK2(   31,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 6)*WRK2(   21,I)
      WRK2(   22,I)=WRK2(   22,I)+CONF( 1)*WRK2(   10,I)
      WRK2(   22,I)=WRK2(   22,I)+CONF( 5)*WRK2(    7,I)
      WRK2(   23,I)=WRK2(   23,I)+CONF( 1)*WRK2(   11,I)
      WRK2(   23,I)=WRK2(   23,I)+CONF( 4)*WRK2(    8,I)
      WRK2(   24,I)=WRK2(   24,I)+CONF( 1)*WRK2(   12,I)
      WRK2(   25,I)=WRK2(   25,I)+CONF( 2)*WRK2(   12,I)
      WRK2(   25,I)=WRK2(   25,I)+CONF( 5)*WRK2(    8,I)
      WRK2(   26,I)=WRK2(   26,I)+CONF( 1)*WRK2(   13,I)
      WRK2(   26,I)=WRK2(   26,I)+CONF( 4)*WRK2(    9,I)
      WRK2(   27,I)=WRK2(   27,I)+CONF( 1)*WRK2(   14,I)
      WRK2(   28,I)=WRK2(   28,I)+CONF( 2)*WRK2(   14,I)
      WRK2(   28,I)=WRK2(   28,I)+CONF( 4)*WRK2(    9,I)
      WRK2(   29,I)=WRK2(   29,I)+CONF( 1)*WRK2(   15,I)
      WRK2(   30,I)=WRK2(   30,I)+CONF( 2)*WRK2(   15,I)
      WRK2(   31,I)=WRK2(   31,I)+CONF( 3)*WRK2(   15,I)
      WRK2(   31,I)=WRK2(   31,I)+CONF( 5)*WRK2(    9,I)
      WRK2(   10,I)=WRK2(   10,I)+CONF( 1)*WRK2(    4,I)
      WRK2(   10,I)=WRK2(   10,I)+CONF( 4)*WRK2(    2,I)
      WRK2(   11,I)=WRK2(   11,I)+CONF( 1)*WRK2(    5,I)
      WRK2(   12,I)=WRK2(   12,I)+CONF( 2)*WRK2(    5,I)
      WRK2(   12,I)=WRK2(   12,I)+CONF( 4)*WRK2(    2,I)
      WRK2(   13,I)=WRK2(   13,I)+CONF( 1)*WRK2(    6,I)
      WRK2(   14,I)=WRK2(   14,I)+CONF( 2)*WRK2(    6,I)
      WRK2(   15,I)=WRK2(   15,I)+CONF( 3)*WRK2(    6,I)
      WRK2(   15,I)=WRK2(   15,I)+CONF( 4)*WRK2(    2,I)
      WRK2(   16,I)=WRK2(   16,I)+CONF( 1)*WRK2(    7,I)
      WRK2(   16,I)=WRK2(   16,I)+CONF( 4)*WRK2(    3,I)
      WRK2(   17,I)=WRK2(   17,I)+CONF( 1)*WRK2(    8,I)
      WRK2(   18,I)=WRK2(   18,I)+CONF( 2)*WRK2(    8,I)
      WRK2(   18,I)=WRK2(   18,I)+CONF( 4)*WRK2(    3,I)
      WRK2(   19,I)=WRK2(   19,I)+CONF( 1)*WRK2(    9,I)
      WRK2(   20,I)=WRK2(   20,I)+CONF( 2)*WRK2(    9,I)
      WRK2(   21,I)=WRK2(   21,I)+CONF( 3)*WRK2(    9,I)
      WRK2(   21,I)=WRK2(   21,I)+CONF( 4)*WRK2(    3,I)
      WRK2(    4,I)=WRK2(    4,I)+CONF( 1)*WRK2(    1,I)
      WRK2(    5,I)=WRK2(    5,I)+CONF( 2)*WRK2(    1,I)
      WRK2(    6,I)=WRK2(    6,I)+CONF( 3)*WRK2(    1,I)
      WRK2(    7,I)=WRK2(    7,I)+CONF( 1)*WRK2(    2,I)
      WRK2(    8,I)=WRK2(    8,I)+CONF( 2)*WRK2(    2,I)
      WRK2(    9,I)=WRK2(    9,I)+CONF( 3)*WRK2(    2,I)
      WRK2(  106,I)=WRK2(  106,I)+CONF( 1)*WRK2(   91,I)
      WRK2(  106,I)=WRK2(  106,I)+CONF( 7)*WRK2(   81,I)
      WRK2(  107,I)=WRK2(  107,I)+CONF( 1)*WRK2(   92,I)
      WRK2(  107,I)=WRK2(  107,I)+CONF( 6)*WRK2(   82,I)
      WRK2(  108,I)=WRK2(  108,I)+CONF( 1)*WRK2(   93,I)
      WRK2(  108,I)=WRK2(  108,I)+CONF( 5)*WRK2(   83,I)
      WRK2(  109,I)=WRK2(  109,I)+CONF( 1)*WRK2(   94,I)
      WRK2(  109,I)=WRK2(  109,I)+CONF( 4)*WRK2(   84,I)
      WRK2(  110,I)=WRK2(  110,I)+CONF( 1)*WRK2(   95,I)
      WRK2(  111,I)=WRK2(  111,I)+CONF( 2)*WRK2(   95,I)
      WRK2(  111,I)=WRK2(  111,I)+CONF( 7)*WRK2(   84,I)
      WRK2(  112,I)=WRK2(  112,I)+CONF( 1)*WRK2(   96,I)
      WRK2(  112,I)=WRK2(  112,I)+CONF( 6)*WRK2(   85,I)
      WRK2(  113,I)=WRK2(  113,I)+CONF( 1)*WRK2(   97,I)
      WRK2(  113,I)=WRK2(  113,I)+CONF( 5)*WRK2(   86,I)
      WRK2(  114,I)=WRK2(  114,I)+CONF( 1)*WRK2(   98,I)
      WRK2(  114,I)=WRK2(  114,I)+CONF( 4)*WRK2(   87,I)
      WRK2(  115,I)=WRK2(  115,I)+CONF( 1)*WRK2(   99,I)
      WRK2(  116,I)=WRK2(  116,I)+CONF( 2)*WRK2(   99,I)
      WRK2(  116,I)=WRK2(  116,I)+CONF( 6)*WRK2(   87,I)
      WRK2(  117,I)=WRK2(  117,I)+CONF( 1)*WRK2(  100,I)
      WRK2(  117,I)=WRK2(  117,I)+CONF( 5)*WRK2(   88,I)
      WRK2(  118,I)=WRK2(  118,I)+CONF( 1)*WRK2(  101,I)
      WRK2(  118,I)=WRK2(  118,I)+CONF( 4)*WRK2(   89,I)
      WRK2(  119,I)=WRK2(  119,I)+CONF( 1)*WRK2(  102,I)
      WRK2(  120,I)=WRK2(  120,I)+CONF( 2)*WRK2(  102,I)
      WRK2(  120,I)=WRK2(  120,I)+CONF( 5)*WRK2(   89,I)
      WRK2(  121,I)=WRK2(  121,I)+CONF( 1)*WRK2(  103,I)
      WRK2(  121,I)=WRK2(  121,I)+CONF( 4)*WRK2(   90,I)
      WRK2(  122,I)=WRK2(  122,I)+CONF( 1)*WRK2(  104,I)
      WRK2(  123,I)=WRK2(  123,I)+CONF( 2)*WRK2(  104,I)
      WRK2(  123,I)=WRK2(  123,I)+CONF( 4)*WRK2(   90,I)
      WRK2(  124,I)=WRK2(  124,I)+CONF( 1)*WRK2(  105,I)
      WRK2(  125,I)=WRK2(  125,I)+CONF( 2)*WRK2(  105,I)
      WRK2(  126,I)=WRK2(  126,I)+CONF( 3)*WRK2(  105,I)
      WRK2(  126,I)=WRK2(  126,I)+CONF( 7)*WRK2(   90,I)
      WRK2(   91,I)=WRK2(   91,I)+CONF( 1)*WRK2(   71,I)
      WRK2(   91,I)=WRK2(   91,I)+CONF( 6)*WRK2(   65,I)
      WRK2(   92,I)=WRK2(   92,I)+CONF( 1)*WRK2(   72,I)
      WRK2(   92,I)=WRK2(   92,I)+CONF( 5)*WRK2(   66,I)
      WRK2(   93,I)=WRK2(   93,I)+CONF( 1)*WRK2(   73,I)
      WRK2(   93,I)=WRK2(   93,I)+CONF( 4)*WRK2(   67,I)
      WRK2(   94,I)=WRK2(   94,I)+CONF( 1)*WRK2(   74,I)
      WRK2(   95,I)=WRK2(   95,I)+CONF( 2)*WRK2(   74,I)
      WRK2(   95,I)=WRK2(   95,I)+CONF( 6)*WRK2(   67,I)
      WRK2(   96,I)=WRK2(   96,I)+CONF( 1)*WRK2(   75,I)
      WRK2(   96,I)=WRK2(   96,I)+CONF( 5)*WRK2(   68,I)
      WRK2(   97,I)=WRK2(   97,I)+CONF( 1)*WRK2(   76,I)
      WRK2(   97,I)=WRK2(   97,I)+CONF( 4)*WRK2(   69,I)
      WRK2(   98,I)=WRK2(   98,I)+CONF( 1)*WRK2(   77,I)
      WRK2(   99,I)=WRK2(   99,I)+CONF( 2)*WRK2(   77,I)
      WRK2(   99,I)=WRK2(   99,I)+CONF( 5)*WRK2(   69,I)
      WRK2(  100,I)=WRK2(  100,I)+CONF( 1)*WRK2(   78,I)
      WRK2(  100,I)=WRK2(  100,I)+CONF( 4)*WRK2(   70,I)
      WRK2(  101,I)=WRK2(  101,I)+CONF( 1)*WRK2(   79,I)
      WRK2(  102,I)=WRK2(  102,I)+CONF( 2)*WRK2(   79,I)
      WRK2(  102,I)=WRK2(  102,I)+CONF( 4)*WRK2(   70,I)
      WRK2(  103,I)=WRK2(  103,I)+CONF( 1)*WRK2(   80,I)
      WRK2(  104,I)=WRK2(  104,I)+CONF( 2)*WRK2(   80,I)
      WRK2(  105,I)=WRK2(  105,I)+CONF( 3)*WRK2(   80,I)
      WRK2(  105,I)=WRK2(  105,I)+CONF( 6)*WRK2(   70,I)
      WRK2(   71,I)=WRK2(   71,I)+CONF( 1)*WRK2(   59,I)
      WRK2(   71,I)=WRK2(   71,I)+CONF( 5)*WRK2(   53,I)
      WRK2(   72,I)=WRK2(   72,I)+CONF( 1)*WRK2(   60,I)
      WRK2(   72,I)=WRK2(   72,I)+CONF( 4)*WRK2(   54,I)
      WRK2(   73,I)=WRK2(   73,I)+CONF( 1)*WRK2(   61,I)
      WRK2(   74,I)=WRK2(   74,I)+CONF( 2)*WRK2(   61,I)
      WRK2(   74,I)=WRK2(   74,I)+CONF( 5)*WRK2(   54,I)
      WRK2(   75,I)=WRK2(   75,I)+CONF( 1)*WRK2(   62,I)
      WRK2(   75,I)=WRK2(   75,I)+CONF( 4)*WRK2(   55,I)
      WRK2(   76,I)=WRK2(   76,I)+CONF( 1)*WRK2(   63,I)
      WRK2(   77,I)=WRK2(   77,I)+CONF( 2)*WRK2(   63,I)
      WRK2(   77,I)=WRK2(   77,I)+CONF( 4)*WRK2(   55,I)
      WRK2(   78,I)=WRK2(   78,I)+CONF( 1)*WRK2(   64,I)
      WRK2(   79,I)=WRK2(   79,I)+CONF( 2)*WRK2(   64,I)
      WRK2(   80,I)=WRK2(   80,I)+CONF( 3)*WRK2(   64,I)
      WRK2(   80,I)=WRK2(   80,I)+CONF( 5)*WRK2(   55,I)
      WRK2(   81,I)=WRK2(   81,I)+CONF( 1)*WRK2(   65,I)
      WRK2(   81,I)=WRK2(   81,I)+CONF( 5)*WRK2(   56,I)
      WRK2(   82,I)=WRK2(   82,I)+CONF( 1)*WRK2(   66,I)
      WRK2(   82,I)=WRK2(   82,I)+CONF( 4)*WRK2(   57,I)
      WRK2(   83,I)=WRK2(   83,I)+CONF( 1)*WRK2(   67,I)
      WRK2(   84,I)=WRK2(   84,I)+CONF( 2)*WRK2(   67,I)
      WRK2(   84,I)=WRK2(   84,I)+CONF( 5)*WRK2(   57,I)
      WRK2(   85,I)=WRK2(   85,I)+CONF( 1)*WRK2(   68,I)
      WRK2(   85,I)=WRK2(   85,I)+CONF( 4)*WRK2(   58,I)
      WRK2(   86,I)=WRK2(   86,I)+CONF( 1)*WRK2(   69,I)
      WRK2(   87,I)=WRK2(   87,I)+CONF( 2)*WRK2(   69,I)
      WRK2(   87,I)=WRK2(   87,I)+CONF( 4)*WRK2(   58,I)
      WRK2(   88,I)=WRK2(   88,I)+CONF( 1)*WRK2(   70,I)
      WRK2(   89,I)=WRK2(   89,I)+CONF( 2)*WRK2(   70,I)
      WRK2(   90,I)=WRK2(   90,I)+CONF( 3)*WRK2(   70,I)
      WRK2(   90,I)=WRK2(   90,I)+CONF( 5)*WRK2(   58,I)
      WRK2(   59,I)=WRK2(   59,I)+CONF( 1)*WRK2(   50,I)
      WRK2(   59,I)=WRK2(   59,I)+CONF( 4)*WRK2(   48,I)
      WRK2(   60,I)=WRK2(   60,I)+CONF( 1)*WRK2(   51,I)
      WRK2(   61,I)=WRK2(   61,I)+CONF( 2)*WRK2(   51,I)
      WRK2(   61,I)=WRK2(   61,I)+CONF( 4)*WRK2(   48,I)
      WRK2(   62,I)=WRK2(   62,I)+CONF( 1)*WRK2(   52,I)
      WRK2(   63,I)=WRK2(   63,I)+CONF( 2)*WRK2(   52,I)
      WRK2(   64,I)=WRK2(   64,I)+CONF( 3)*WRK2(   52,I)
      WRK2(   64,I)=WRK2(   64,I)+CONF( 4)*WRK2(   48,I)
      WRK2(   65,I)=WRK2(   65,I)+CONF( 1)*WRK2(   53,I)
      WRK2(   65,I)=WRK2(   65,I)+CONF( 4)*WRK2(   49,I)
      WRK2(   66,I)=WRK2(   66,I)+CONF( 1)*WRK2(   54,I)
      WRK2(   67,I)=WRK2(   67,I)+CONF( 2)*WRK2(   54,I)
      WRK2(   67,I)=WRK2(   67,I)+CONF( 4)*WRK2(   49,I)
      WRK2(   68,I)=WRK2(   68,I)+CONF( 1)*WRK2(   55,I)
      WRK2(   69,I)=WRK2(   69,I)+CONF( 2)*WRK2(   55,I)
      WRK2(   70,I)=WRK2(   70,I)+CONF( 3)*WRK2(   55,I)
      WRK2(   70,I)=WRK2(   70,I)+CONF( 4)*WRK2(   49,I)
      WRK2(   50,I)=WRK2(   50,I)+CONF( 1)*WRK2(   47,I)
      WRK2(   51,I)=WRK2(   51,I)+CONF( 2)*WRK2(   47,I)
      WRK2(   52,I)=WRK2(   52,I)+CONF( 3)*WRK2(   47,I)
      WRK2(   53,I)=WRK2(   53,I)+CONF( 1)*WRK2(   48,I)
      WRK2(   54,I)=WRK2(   54,I)+CONF( 2)*WRK2(   48,I)
      WRK2(   55,I)=WRK2(   55,I)+CONF( 3)*WRK2(   48,I)
      WRK2(   56,I)=WRK2(   56,I)+CONF( 1)*WRK2(   49,I)
      WRK2(   57,I)=WRK2(   57,I)+CONF( 2)*WRK2(   49,I)
      WRK2(   58,I)=WRK2(   58,I)+CONF( 3)*WRK2(   49,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 1)*WRK2(   22,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 5)*WRK2(   16,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 1)*WRK2(   23,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 4)*WRK2(   17,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 1)*WRK2(   24,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 2)*WRK2(   24,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 5)*WRK2(   17,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 2)*WRK2(   25,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 5)*WRK2(   18,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 1)*WRK2(   26,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 4)*WRK2(   19,I)
      WRK2(   38,I)=WRK2(   38,I)+CONF( 1)*WRK2(   27,I)
      WRK2(   39,I)=WRK2(   39,I)+CONF( 2)*WRK2(   27,I)
      WRK2(   39,I)=WRK2(   39,I)+CONF( 4)*WRK2(   19,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 2)*WRK2(   28,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 4)*WRK2(   20,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 1)*WRK2(   29,I)
      WRK2(   42,I)=WRK2(   42,I)+CONF( 2)*WRK2(   29,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 2)*WRK2(   30,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 3)*WRK2(   29,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 5)*WRK2(   19,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 3)*WRK2(   30,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 5)*WRK2(   20,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 3)*WRK2(   31,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 5)*WRK2(   21,I)
      WRK2(   22,I)=WRK2(   22,I)+CONF( 1)*WRK2(   10,I)
      WRK2(   22,I)=WRK2(   22,I)+CONF( 4)*WRK2(    7,I)
      WRK2(   23,I)=WRK2(   23,I)+CONF( 1)*WRK2(   11,I)
      WRK2(   24,I)=WRK2(   24,I)+CONF( 2)*WRK2(   11,I)
      WRK2(   24,I)=WRK2(   24,I)+CONF( 4)*WRK2(    7,I)
      WRK2(   25,I)=WRK2(   25,I)+CONF( 2)*WRK2(   12,I)
      WRK2(   25,I)=WRK2(   25,I)+CONF( 4)*WRK2(    8,I)
      WRK2(   26,I)=WRK2(   26,I)+CONF( 1)*WRK2(   13,I)
      WRK2(   27,I)=WRK2(   27,I)+CONF( 2)*WRK2(   13,I)
      WRK2(   28,I)=WRK2(   28,I)+CONF( 2)*WRK2(   14,I)
      WRK2(   29,I)=WRK2(   29,I)+CONF( 3)*WRK2(   13,I)
      WRK2(   29,I)=WRK2(   29,I)+CONF( 4)*WRK2(    7,I)
      WRK2(   30,I)=WRK2(   30,I)+CONF( 3)*WRK2(   14,I)
      WRK2(   30,I)=WRK2(   30,I)+CONF( 4)*WRK2(    8,I)
      WRK2(   31,I)=WRK2(   31,I)+CONF( 3)*WRK2(   15,I)
      WRK2(   31,I)=WRK2(   31,I)+CONF( 4)*WRK2(    9,I)
      WRK2(   10,I)=WRK2(   10,I)+CONF( 1)*WRK2(    4,I)
      WRK2(   11,I)=WRK2(   11,I)+CONF( 2)*WRK2(    4,I)
      WRK2(   12,I)=WRK2(   12,I)+CONF( 2)*WRK2(    5,I)
      WRK2(   13,I)=WRK2(   13,I)+CONF( 3)*WRK2(    4,I)
      WRK2(   14,I)=WRK2(   14,I)+CONF( 3)*WRK2(    5,I)
      WRK2(   15,I)=WRK2(   15,I)+CONF( 3)*WRK2(    6,I)
      WRK2(   16,I)=WRK2(   16,I)+CONF( 1)*WRK2(    7,I)
      WRK2(   17,I)=WRK2(   17,I)+CONF( 2)*WRK2(    7,I)
      WRK2(   18,I)=WRK2(   18,I)+CONF( 2)*WRK2(    8,I)
      WRK2(   19,I)=WRK2(   19,I)+CONF( 3)*WRK2(    7,I)
      WRK2(   20,I)=WRK2(   20,I)+CONF( 3)*WRK2(    8,I)
      WRK2(   21,I)=WRK2(   21,I)+CONF( 3)*WRK2(    9,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 1)*WRK2(   22,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 4)*WRK2(   16,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 1)*WRK2(   23,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 2)*WRK2(   23,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 4)*WRK2(   16,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 2)*WRK2(   24,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 4)*WRK2(   17,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 2)*WRK2(   25,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 4)*WRK2(   18,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 1)*WRK2(   26,I)
      WRK2(   38,I)=WRK2(   38,I)+CONF( 2)*WRK2(   26,I)
      WRK2(   39,I)=WRK2(   39,I)+CONF( 2)*WRK2(   27,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 2)*WRK2(   28,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 3)*WRK2(   26,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 4)*WRK2(   16,I)
      WRK2(   42,I)=WRK2(   42,I)+CONF( 3)*WRK2(   27,I)
      WRK2(   42,I)=WRK2(   42,I)+CONF( 4)*WRK2(   17,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 3)*WRK2(   28,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 4)*WRK2(   18,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 3)*WRK2(   29,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 4)*WRK2(   19,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 3)*WRK2(   30,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 4)*WRK2(   20,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 3)*WRK2(   31,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 4)*WRK2(   21,I)
      WRK2(   22,I)=WRK2(   22,I)+CONF( 1)*WRK2(   10,I)
      WRK2(   23,I)=WRK2(   23,I)+CONF( 2)*WRK2(   10,I)
      WRK2(   24,I)=WRK2(   24,I)+CONF( 2)*WRK2(   11,I)
      WRK2(   25,I)=WRK2(   25,I)+CONF( 2)*WRK2(   12,I)
      WRK2(   26,I)=WRK2(   26,I)+CONF( 3)*WRK2(   10,I)
      WRK2(   27,I)=WRK2(   27,I)+CONF( 3)*WRK2(   11,I)
      WRK2(   28,I)=WRK2(   28,I)+CONF( 3)*WRK2(   12,I)
      WRK2(   29,I)=WRK2(   29,I)+CONF( 3)*WRK2(   13,I)
      WRK2(   30,I)=WRK2(   30,I)+CONF( 3)*WRK2(   14,I)
      WRK2(   31,I)=WRK2(   31,I)+CONF( 3)*WRK2(   15,I)
      WRK2(   32,I)=WRK2(   32,I)+CONF( 1)*WRK2(   22,I)
      WRK2(   33,I)=WRK2(   33,I)+CONF( 2)*WRK2(   22,I)
      WRK2(   34,I)=WRK2(   34,I)+CONF( 2)*WRK2(   23,I)
      WRK2(   35,I)=WRK2(   35,I)+CONF( 2)*WRK2(   24,I)
      WRK2(   36,I)=WRK2(   36,I)+CONF( 2)*WRK2(   25,I)
      WRK2(   37,I)=WRK2(   37,I)+CONF( 3)*WRK2(   22,I)
      WRK2(   38,I)=WRK2(   38,I)+CONF( 3)*WRK2(   23,I)
      WRK2(   39,I)=WRK2(   39,I)+CONF( 3)*WRK2(   24,I)
      WRK2(   40,I)=WRK2(   40,I)+CONF( 3)*WRK2(   25,I)
      WRK2(   41,I)=WRK2(   41,I)+CONF( 3)*WRK2(   26,I)
      WRK2(   42,I)=WRK2(   42,I)+CONF( 3)*WRK2(   27,I)
      WRK2(   43,I)=WRK2(   43,I)+CONF( 3)*WRK2(   28,I)
      WRK2(   44,I)=WRK2(   44,I)+CONF( 3)*WRK2(   29,I)
      WRK2(   45,I)=WRK2(   45,I)+CONF( 3)*WRK2(   30,I)
      WRK2(   46,I)=WRK2(   46,I)+CONF( 3)*WRK2(   31,I)
      WRK2(  106,I)=WRK2(  106,I)+CONF( 1)*WRK2(   91,I)
      WRK2(  106,I)=WRK2(  106,I)+CONF( 6)*WRK2(   81,I)
      WRK2(  107,I)=WRK2(  107,I)+CONF( 1)*WRK2(   92,I)
      WRK2(  107,I)=WRK2(  107,I)+CONF( 5)*WRK2(   82,I)
      WRK2(  108,I)=WRK2(  108,I)+CONF( 1)*WRK2(   93,I)
      WRK2(  108,I)=WRK2(  108,I)+CONF( 4)*WRK2(   83,I)
      WRK2(  109,I)=WRK2(  109,I)+CONF( 1)*WRK2(   94,I)
      WRK2(  110,I)=WRK2(  110,I)+CONF( 2)*WRK2(   94,I)
      WRK2(  110,I)=WRK2(  110,I)+CONF( 6)*WRK2(   83,I)
      WRK2(  111,I)=WRK2(  111,I)+CONF( 2)*WRK2(   95,I)
      WRK2(  111,I)=WRK2(  111,I)+CONF( 6)*WRK2(   84,I)
      WRK2(  112,I)=WRK2(  112,I)+CONF( 1)*WRK2(   96,I)
      WRK2(  112,I)=WRK2(  112,I)+CONF( 5)*WRK2(   85,I)
      WRK2(  113,I)=WRK2(  113,I)+CONF( 1)*WRK2(   97,I)
      WRK2(  113,I)=WRK2(  113,I)+CONF( 4)*WRK2(   86,I)
      WRK2(  114,I)=WRK2(  114,I)+CONF( 1)*WRK2(   98,I)
      WRK2(  115,I)=WRK2(  115,I)+CONF( 2)*WRK2(   98,I)
      WRK2(  115,I)=WRK2(  115,I)+CONF( 5)*WRK2(   86,I)
      WRK2(  116,I)=WRK2(  116,I)+CONF( 2)*WRK2(   99,I)
      WRK2(  116,I)=WRK2(  116,I)+CONF( 5)*WRK2(   87,I)
      WRK2(  117,I)=WRK2(  117,I)+CONF( 1)*WRK2(  100,I)
      WRK2(  117,I)=WRK2(  117,I)+CONF( 4)*WRK2(   88,I)
      WRK2(  118,I)=WRK2(  118,I)+CONF( 1)*WRK2(  101,I)
      WRK2(  119,I)=WRK2(  119,I)+CONF( 2)*WRK2(  101,I)
      WRK2(  119,I)=WRK2(  119,I)+CONF( 4)*WRK2(   88,I)
      WRK2(  120,I)=WRK2(  120,I)+CONF( 2)*WRK2(  102,I)
      WRK2(  120,I)=WRK2(  120,I)+CONF( 4)*WRK2(   89,I)
      WRK2(  121,I)=WRK2(  121,I)+CONF( 1)*WRK2(  103,I)
      WRK2(  122,I)=WRK2(  122,I)+CONF( 2)*WRK2(  103,I)
      WRK2(  123,I)=WRK2(  123,I)+CONF( 2)*WRK2(  104,I)
      WRK2(  124,I)=WRK2(  124,I)+CONF( 3)*WRK2(  103,I)
      WRK2(  124,I)=WRK2(  124,I)+CONF( 6)*WRK2(   88,I)
      WRK2(  125,I)=WRK2(  125,I)+CONF( 3)*WRK2(  104,I)
      WRK2(  125,I)=WRK2(  125,I)+CONF( 6)*WRK2(   89,I)
      WRK2(  126,I)=WRK2(  126,I)+CONF( 3)*WRK2(  105,I)
      WRK2(  126,I)=WRK2(  126,I)+CONF( 6)*WRK2(   90,I)
      WRK2(   91,I)=WRK2(   91,I)+CONF( 1)*WRK2(   71,I)
      WRK2(   91,I)=WRK2(   91,I)+CONF( 5)*WRK2(   65,I)
      WRK2(   92,I)=WRK2(   92,I)+CONF( 1)*WRK2(   72,I)
      WRK2(   92,I)=WRK2(   92,I)+CONF( 4)*WRK2(   66,I)
      WRK2(   93,I)=WRK2(   93,I)+CONF( 1)*WRK2(   73,I)
      WRK2(   94,I)=WRK2(   94,I)+CONF( 2)*WRK2(   73,I)
      WRK2(   94,I)=WRK2(   94,I)+CONF( 5)*WRK2(   66,I)
      WRK2(   95,I)=WRK2(   95,I)+CONF( 2)*WRK2(   74,I)
      WRK2(   95,I)=WRK2(   95,I)+CONF( 5)*WRK2(   67,I)
      WRK2(   96,I)=WRK2(   96,I)+CONF( 1)*WRK2(   75,I)
      WRK2(   96,I)=WRK2(   96,I)+CONF( 4)*WRK2(   68,I)
      WRK2(   97,I)=WRK2(   97,I)+CONF( 1)*WRK2(   76,I)
      WRK2(   98,I)=WRK2(   98,I)+CONF( 2)*WRK2(   76,I)
      WRK2(   98,I)=WRK2(   98,I)+CONF( 4)*WRK2(   68,I)
      WRK2(   99,I)=WRK2(   99,I)+CONF( 2)*WRK2(   77,I)
      WRK2(   99,I)=WRK2(   99,I)+CONF( 4)*WRK2(   69,I)
      WRK2(  100,I)=WRK2(  100,I)+CONF( 1)*WRK2(   78,I)
      WRK2(  101,I)=WRK2(  101,I)+CONF( 2)*WRK2(   78,I)
      WRK2(  102,I)=WRK2(  102,I)+CONF( 2)*WRK2(   79,I)
      WRK2(  103,I)=WRK2(  103,I)+CONF( 3)*WRK2(   78,I)
      WRK2(  103,I)=WRK2(  103,I)+CONF( 5)*WRK2(   68,I)
      WRK2(  104,I)=WRK2(  104,I)+CONF( 3)*WRK2(   79,I)
      WRK2(  104,I)=WRK2(  104,I)+CONF( 5)*WRK2(   69,I)
      WRK2(  105,I)=WRK2(  105,I)+CONF( 3)*WRK2(   80,I)
      WRK2(  105,I)=WRK2(  105,I)+CONF( 5)*WRK2(   70,I)
      WRK2(   71,I)=WRK2(   71,I)+CONF( 1)*WRK2(   59,I)
      WRK2(   71,I)=WRK2(   71,I)+CONF( 4)*WRK2(   53,I)
      WRK2(   72,I)=WRK2(   72,I)+CONF( 1)*WRK2(   60,I)
      WRK2(   73,I)=WRK2(   73,I)+CONF( 2)*WRK2(   60,I)
      WRK2(   73,I)=WRK2(   73,I)+CONF( 4)*WRK2(   53,I)
      WRK2(   74,I)=WRK2(   74,I)+CONF( 2)*WRK2(   61,I)
      WRK2(   74,I)=WRK2(   74,I)+CONF( 4)*WRK2(   54,I)
      WRK2(   75,I)=WRK2(   75,I)+CONF( 1)*WRK2(   62,I)
      WRK2(   76,I)=WRK2(   76,I)+CONF( 2)*WRK2(   62,I)
      WRK2(   77,I)=WRK2(   77,I)+CONF( 2)*WRK2(   63,I)
      WRK2(   78,I)=WRK2(   78,I)+CONF( 3)*WRK2(   62,I)
      WRK2(   78,I)=WRK2(   78,I)+CONF( 4)*WRK2(   53,I)
      WRK2(   79,I)=WRK2(   79,I)+CONF( 3)*WRK2(   63,I)
      WRK2(   79,I)=WRK2(   79,I)+CONF( 4)*WRK2(   54,I)
      WRK2(   80,I)=WRK2(   80,I)+CONF( 3)*WRK2(   64,I)
      WRK2(   80,I)=WRK2(   80,I)+CONF( 4)*WRK2(   55,I)
      WRK2(   81,I)=WRK2(   81,I)+CONF( 1)*WRK2(   65,I)
      WRK2(   81,I)=WRK2(   81,I)+CONF( 4)*WRK2(   56,I)
      WRK2(   82,I)=WRK2(   82,I)+CONF( 1)*WRK2(   66,I)
      WRK2(   83,I)=WRK2(   83,I)+CONF( 2)*WRK2(   66,I)
      WRK2(   83,I)=WRK2(   83,I)+CONF( 4)*WRK2(   56,I)
      WRK2(   84,I)=WRK2(   84,I)+CONF( 2)*WRK2(   67,I)
      WRK2(   84,I)=WRK2(   84,I)+CONF( 4)*WRK2(   57,I)
      WRK2(   85,I)=WRK2(   85,I)+CONF( 1)*WRK2(   68,I)
      WRK2(   86,I)=WRK2(   86,I)+CONF( 2)*WRK2(   68,I)
      WRK2(   87,I)=WRK2(   87,I)+CONF( 2)*WRK2(   69,I)
      WRK2(   88,I)=WRK2(   88,I)+CONF( 3)*WRK2(   68,I)
      WRK2(   88,I)=WRK2(   88,I)+CONF( 4)*WRK2(   56,I)
      WRK2(   89,I)=WRK2(   89,I)+CONF( 3)*WRK2(   69,I)
      WRK2(   89,I)=WRK2(   89,I)+CONF( 4)*WRK2(   57,I)
      WRK2(   90,I)=WRK2(   90,I)+CONF( 3)*WRK2(   70,I)
      WRK2(   90,I)=WRK2(   90,I)+CONF( 4)*WRK2(   58,I)
      WRK2(   59,I)=WRK2(   59,I)+CONF( 1)*WRK2(   50,I)
      WRK2(   60,I)=WRK2(   60,I)+CONF( 2)*WRK2(   50,I)
      WRK2(   61,I)=WRK2(   61,I)+CONF( 2)*WRK2(   51,I)
      WRK2(   62,I)=WRK2(   62,I)+CONF( 3)*WRK2(   50,I)
      WRK2(   63,I)=WRK2(   63,I)+CONF( 3)*WRK2(   51,I)
      WRK2(   64,I)=WRK2(   64,I)+CONF( 3)*WRK2(   52,I)
      WRK2(   65,I)=WRK2(   65,I)+CONF( 1)*WRK2(   53,I)
      WRK2(   66,I)=WRK2(   66,I)+CONF( 2)*WRK2(   53,I)
      WRK2(   67,I)=WRK2(   67,I)+CONF( 2)*WRK2(   54,I)
      WRK2(   68,I)=WRK2(   68,I)+CONF( 3)*WRK2(   53,I)
      WRK2(   69,I)=WRK2(   69,I)+CONF( 3)*WRK2(   54,I)
      WRK2(   70,I)=WRK2(   70,I)+CONF( 3)*WRK2(   55,I)
      WRK2(  106,I)=WRK2(  106,I)+CONF( 1)*WRK2(   91,I)
      WRK2(  106,I)=WRK2(  106,I)+CONF( 5)*WRK2(   81,I)
      WRK2(  107,I)=WRK2(  107,I)+CONF( 1)*WRK2(   92,I)
      WRK2(  107,I)=WRK2(  107,I)+CONF( 4)*WRK2(   82,I)
      WRK2(  108,I)=WRK2(  108,I)+CONF( 1)*WRK2(   93,I)
      WRK2(  109,I)=WRK2(  109,I)+CONF( 2)*WRK2(   93,I)
      WRK2(  109,I)=WRK2(  109,I)+CONF( 5)*WRK2(   82,I)
      WRK2(  110,I)=WRK2(  110,I)+CONF( 2)*WRK2(   94,I)
      WRK2(  110,I)=WRK2(  110,I)+CONF( 5)*WRK2(   83,I)
      WRK2(  111,I)=WRK2(  111,I)+CONF( 2)*WRK2(   95,I)
      WRK2(  111,I)=WRK2(  111,I)+CONF( 5)*WRK2(   84,I)
      WRK2(  112,I)=WRK2(  112,I)+CONF( 1)*WRK2(   96,I)
      WRK2(  112,I)=WRK2(  112,I)+CONF( 4)*WRK2(   85,I)
      WRK2(  113,I)=WRK2(  113,I)+CONF( 1)*WRK2(   97,I)
      WRK2(  114,I)=WRK2(  114,I)+CONF( 2)*WRK2(   97,I)
      WRK2(  114,I)=WRK2(  114,I)+CONF( 4)*WRK2(   85,I)
      WRK2(  115,I)=WRK2(  115,I)+CONF( 2)*WRK2(   98,I)
      WRK2(  115,I)=WRK2(  115,I)+CONF( 4)*WRK2(   86,I)
      WRK2(  116,I)=WRK2(  116,I)+CONF( 2)*WRK2(   99,I)
      WRK2(  116,I)=WRK2(  116,I)+CONF( 4)*WRK2(   87,I)
      WRK2(  117,I)=WRK2(  117,I)+CONF( 1)*WRK2(  100,I)
      WRK2(  118,I)=WRK2(  118,I)+CONF( 2)*WRK2(  100,I)
      WRK2(  119,I)=WRK2(  119,I)+CONF( 2)*WRK2(  101,I)
      WRK2(  120,I)=WRK2(  120,I)+CONF( 2)*WRK2(  102,I)
      WRK2(  121,I)=WRK2(  121,I)+CONF( 3)*WRK2(  100,I)
      WRK2(  121,I)=WRK2(  121,I)+CONF( 5)*WRK2(   85,I)
      WRK2(  122,I)=WRK2(  122,I)+CONF( 3)*WRK2(  101,I)
      WRK2(  122,I)=WRK2(  122,I)+CONF( 5)*WRK2(   86,I)
      WRK2(  123,I)=WRK2(  123,I)+CONF( 3)*WRK2(  102,I)
      WRK2(  123,I)=WRK2(  123,I)+CONF( 5)*WRK2(   87,I)
      WRK2(  124,I)=WRK2(  124,I)+CONF( 3)*WRK2(  103,I)
      WRK2(  124,I)=WRK2(  124,I)+CONF( 5)*WRK2(   88,I)
      WRK2(  125,I)=WRK2(  125,I)+CONF( 3)*WRK2(  104,I)
      WRK2(  125,I)=WRK2(  125,I)+CONF( 5)*WRK2(   89,I)
      WRK2(  126,I)=WRK2(  126,I)+CONF( 3)*WRK2(  105,I)
      WRK2(  126,I)=WRK2(  126,I)+CONF( 5)*WRK2(   90,I)
      WRK2(   91,I)=WRK2(   91,I)+CONF( 1)*WRK2(   71,I)
      WRK2(   91,I)=WRK2(   91,I)+CONF( 4)*WRK2(   65,I)
      WRK2(   92,I)=WRK2(   92,I)+CONF( 1)*WRK2(   72,I)
      WRK2(   93,I)=WRK2(   93,I)+CONF( 2)*WRK2(   72,I)
      WRK2(   93,I)=WRK2(   93,I)+CONF( 4)*WRK2(   65,I)
      WRK2(   94,I)=WRK2(   94,I)+CONF( 2)*WRK2(   73,I)
      WRK2(   94,I)=WRK2(   94,I)+CONF( 4)*WRK2(   66,I)
      WRK2(   95,I)=WRK2(   95,I)+CONF( 2)*WRK2(   74,I)
      WRK2(   95,I)=WRK2(   95,I)+CONF( 4)*WRK2(   67,I)
      WRK2(   96,I)=WRK2(   96,I)+CONF( 1)*WRK2(   75,I)
      WRK2(   97,I)=WRK2(   97,I)+CONF( 2)*WRK2(   75,I)
      WRK2(   98,I)=WRK2(   98,I)+CONF( 2)*WRK2(   76,I)
      WRK2(   99,I)=WRK2(   99,I)+CONF( 2)*WRK2(   77,I)
      WRK2(  100,I)=WRK2(  100,I)+CONF( 3)*WRK2(   75,I)
      WRK2(  100,I)=WRK2(  100,I)+CONF( 4)*WRK2(   65,I)
      WRK2(  101,I)=WRK2(  101,I)+CONF( 3)*WRK2(   76,I)
      WRK2(  101,I)=WRK2(  101,I)+CONF( 4)*WRK2(   66,I)
      WRK2(  102,I)=WRK2(  102,I)+CONF( 3)*WRK2(   77,I)
      WRK2(  102,I)=WRK2(  102,I)+CONF( 4)*WRK2(   67,I)
      WRK2(  103,I)=WRK2(  103,I)+CONF( 3)*WRK2(   78,I)
      WRK2(  103,I)=WRK2(  103,I)+CONF( 4)*WRK2(   68,I)
      WRK2(  104,I)=WRK2(  104,I)+CONF( 3)*WRK2(   79,I)
      WRK2(  104,I)=WRK2(  104,I)+CONF( 4)*WRK2(   69,I)
      WRK2(  105,I)=WRK2(  105,I)+CONF( 3)*WRK2(   80,I)
      WRK2(  105,I)=WRK2(  105,I)+CONF( 4)*WRK2(   70,I)
      WRK2(   71,I)=WRK2(   71,I)+CONF( 1)*WRK2(   59,I)
      WRK2(   72,I)=WRK2(   72,I)+CONF( 2)*WRK2(   59,I)
      WRK2(   73,I)=WRK2(   73,I)+CONF( 2)*WRK2(   60,I)
      WRK2(   74,I)=WRK2(   74,I)+CONF( 2)*WRK2(   61,I)
      WRK2(   75,I)=WRK2(   75,I)+CONF( 3)*WRK2(   59,I)
      WRK2(   76,I)=WRK2(   76,I)+CONF( 3)*WRK2(   60,I)
      WRK2(   77,I)=WRK2(   77,I)+CONF( 3)*WRK2(   61,I)
      WRK2(   78,I)=WRK2(   78,I)+CONF( 3)*WRK2(   62,I)
      WRK2(   79,I)=WRK2(   79,I)+CONF( 3)*WRK2(   63,I)
      WRK2(   80,I)=WRK2(   80,I)+CONF( 3)*WRK2(   64,I)
      WRK2(   81,I)=WRK2(   81,I)+CONF( 1)*WRK2(   65,I)
      WRK2(   82,I)=WRK2(   82,I)+CONF( 2)*WRK2(   65,I)
      WRK2(   83,I)=WRK2(   83,I)+CONF( 2)*WRK2(   66,I)
      WRK2(   84,I)=WRK2(   84,I)+CONF( 2)*WRK2(   67,I)
      WRK2(   85,I)=WRK2(   85,I)+CONF( 3)*WRK2(   65,I)
      WRK2(   86,I)=WRK2(   86,I)+CONF( 3)*WRK2(   66,I)
      WRK2(   87,I)=WRK2(   87,I)+CONF( 3)*WRK2(   67,I)
      WRK2(   88,I)=WRK2(   88,I)+CONF( 3)*WRK2(   68,I)
      WRK2(   89,I)=WRK2(   89,I)+CONF( 3)*WRK2(   69,I)
      WRK2(   90,I)=WRK2(   90,I)+CONF( 3)*WRK2(   70,I)
      WRK2(  106,I)=WRK2(  106,I)+CONF( 1)*WRK2(   91,I)
      WRK2(  106,I)=WRK2(  106,I)+CONF( 4)*WRK2(   81,I)
      WRK2(  107,I)=WRK2(  107,I)+CONF( 1)*WRK2(   92,I)
      WRK2(  108,I)=WRK2(  108,I)+CONF( 2)*WRK2(   92,I)
      WRK2(  108,I)=WRK2(  108,I)+CONF( 4)*WRK2(   81,I)
      WRK2(  109,I)=WRK2(  109,I)+CONF( 2)*WRK2(   93,I)
      WRK2(  109,I)=WRK2(  109,I)+CONF( 4)*WRK2(   82,I)
      WRK2(  110,I)=WRK2(  110,I)+CONF( 2)*WRK2(   94,I)
      WRK2(  110,I)=WRK2(  110,I)+CONF( 4)*WRK2(   83,I)
      WRK2(  111,I)=WRK2(  111,I)+CONF( 2)*WRK2(   95,I)
      WRK2(  111,I)=WRK2(  111,I)+CONF( 4)*WRK2(   84,I)
      WRK2(  112,I)=WRK2(  112,I)+CONF( 1)*WRK2(   96,I)
      WRK2(  113,I)=WRK2(  113,I)+CONF( 2)*WRK2(   96,I)
      WRK2(  114,I)=WRK2(  114,I)+CONF( 2)*WRK2(   97,I)
      WRK2(  115,I)=WRK2(  115,I)+CONF( 2)*WRK2(   98,I)
      WRK2(  116,I)=WRK2(  116,I)+CONF( 2)*WRK2(   99,I)
      WRK2(  117,I)=WRK2(  117,I)+CONF( 3)*WRK2(   96,I)
      WRK2(  117,I)=WRK2(  117,I)+CONF( 4)*WRK2(   81,I)
      WRK2(  118,I)=WRK2(  118,I)+CONF( 3)*WRK2(   97,I)
      WRK2(  118,I)=WRK2(  118,I)+CONF( 4)*WRK2(   82,I)
      WRK2(  119,I)=WRK2(  119,I)+CONF( 3)*WRK2(   98,I)
      WRK2(  119,I)=WRK2(  119,I)+CONF( 4)*WRK2(   83,I)
      WRK2(  120,I)=WRK2(  120,I)+CONF( 3)*WRK2(   99,I)
      WRK2(  120,I)=WRK2(  120,I)+CONF( 4)*WRK2(   84,I)
      WRK2(  121,I)=WRK2(  121,I)+CONF( 3)*WRK2(  100,I)
      WRK2(  121,I)=WRK2(  121,I)+CONF( 4)*WRK2(   85,I)
      WRK2(  122,I)=WRK2(  122,I)+CONF( 3)*WRK2(  101,I)
      WRK2(  122,I)=WRK2(  122,I)+CONF( 4)*WRK2(   86,I)
      WRK2(  123,I)=WRK2(  123,I)+CONF( 3)*WRK2(  102,I)
      WRK2(  123,I)=WRK2(  123,I)+CONF( 4)*WRK2(   87,I)
      WRK2(  124,I)=WRK2(  124,I)+CONF( 3)*WRK2(  103,I)
      WRK2(  124,I)=WRK2(  124,I)+CONF( 4)*WRK2(   88,I)
      WRK2(  125,I)=WRK2(  125,I)+CONF( 3)*WRK2(  104,I)
      WRK2(  125,I)=WRK2(  125,I)+CONF( 4)*WRK2(   89,I)
      WRK2(  126,I)=WRK2(  126,I)+CONF( 3)*WRK2(  105,I)
      WRK2(  126,I)=WRK2(  126,I)+CONF( 4)*WRK2(   90,I)
      WRK2(   91,I)=WRK2(   91,I)+CONF( 1)*WRK2(   71,I)
      WRK2(   92,I)=WRK2(   92,I)+CONF( 2)*WRK2(   71,I)
      WRK2(   93,I)=WRK2(   93,I)+CONF( 2)*WRK2(   72,I)
      WRK2(   94,I)=WRK2(   94,I)+CONF( 2)*WRK2(   73,I)
      WRK2(   95,I)=WRK2(   95,I)+CONF( 2)*WRK2(   74,I)
      WRK2(   96,I)=WRK2(   96,I)+CONF( 3)*WRK2(   71,I)
      WRK2(   97,I)=WRK2(   97,I)+CONF( 3)*WRK2(   72,I)
      WRK2(   98,I)=WRK2(   98,I)+CONF( 3)*WRK2(   73,I)
      WRK2(   99,I)=WRK2(   99,I)+CONF( 3)*WRK2(   74,I)
      WRK2(  100,I)=WRK2(  100,I)+CONF( 3)*WRK2(   75,I)
      WRK2(  101,I)=WRK2(  101,I)+CONF( 3)*WRK2(   76,I)
      WRK2(  102,I)=WRK2(  102,I)+CONF( 3)*WRK2(   77,I)
      WRK2(  103,I)=WRK2(  103,I)+CONF( 3)*WRK2(   78,I)
      WRK2(  104,I)=WRK2(  104,I)+CONF( 3)*WRK2(   79,I)
      WRK2(  105,I)=WRK2(  105,I)+CONF( 3)*WRK2(   80,I)
      WRK2(  106,I)=WRK2(  106,I)+CONF( 1)*WRK2(   91,I)
      WRK2(  107,I)=WRK2(  107,I)+CONF( 2)*WRK2(   91,I)
      WRK2(  108,I)=WRK2(  108,I)+CONF( 2)*WRK2(   92,I)
      WRK2(  109,I)=WRK2(  109,I)+CONF( 2)*WRK2(   93,I)
      WRK2(  110,I)=WRK2(  110,I)+CONF( 2)*WRK2(   94,I)
      WRK2(  111,I)=WRK2(  111,I)+CONF( 2)*WRK2(   95,I)
      WRK2(  112,I)=WRK2(  112,I)+CONF( 3)*WRK2(   91,I)
      WRK2(  113,I)=WRK2(  113,I)+CONF( 3)*WRK2(   92,I)
      WRK2(  114,I)=WRK2(  114,I)+CONF( 3)*WRK2(   93,I)
      WRK2(  115,I)=WRK2(  115,I)+CONF( 3)*WRK2(   94,I)
      WRK2(  116,I)=WRK2(  116,I)+CONF( 3)*WRK2(   95,I)
      WRK2(  117,I)=WRK2(  117,I)+CONF( 3)*WRK2(   96,I)
      WRK2(  118,I)=WRK2(  118,I)+CONF( 3)*WRK2(   97,I)
      WRK2(  119,I)=WRK2(  119,I)+CONF( 3)*WRK2(   98,I)
      WRK2(  120,I)=WRK2(  120,I)+CONF( 3)*WRK2(   99,I)
      WRK2(  121,I)=WRK2(  121,I)+CONF( 3)*WRK2(  100,I)
      WRK2(  122,I)=WRK2(  122,I)+CONF( 3)*WRK2(  101,I)
      WRK2(  123,I)=WRK2(  123,I)+CONF( 3)*WRK2(  102,I)
      WRK2(  124,I)=WRK2(  124,I)+CONF( 3)*WRK2(  103,I)
      WRK2(  125,I)=WRK2(  125,I)+CONF( 3)*WRK2(  104,I)
      WRK2(  126,I)=WRK2(  126,I)+CONF( 3)*WRK2(  105,I)
      WRK2(  127,I)=WRK2(  106,I)-CONF( 1)*WRK2(   32,I)
      WRK2(  128,I)=WRK2(  107,I)-CONF( 1)*WRK2(   33,I)
      WRK2(  129,I)=WRK2(  108,I)-CONF( 1)*WRK2(   34,I)
      WRK2(  130,I)=WRK2(  109,I)-CONF( 1)*WRK2(   35,I)
      WRK2(  131,I)=WRK2(  110,I)-CONF( 1)*WRK2(   36,I)
      WRK2(  132,I)=WRK2(  112,I)-CONF( 1)*WRK2(   37,I)
      WRK2(  133,I)=WRK2(  113,I)-CONF( 1)*WRK2(   38,I)
      WRK2(  134,I)=WRK2(  114,I)-CONF( 1)*WRK2(   39,I)
      WRK2(  135,I)=WRK2(  115,I)-CONF( 1)*WRK2(   40,I)
      WRK2(  136,I)=WRK2(  117,I)-CONF( 1)*WRK2(   41,I)
      WRK2(  137,I)=WRK2(  118,I)-CONF( 1)*WRK2(   42,I)
      WRK2(  138,I)=WRK2(  119,I)-CONF( 1)*WRK2(   43,I)
      WRK2(  139,I)=WRK2(  121,I)-CONF( 1)*WRK2(   44,I)
      WRK2(  140,I)=WRK2(  122,I)-CONF( 1)*WRK2(   45,I)
      WRK2(  141,I)=WRK2(  124,I)-CONF( 1)*WRK2(   46,I)
      WRK2(  142,I)=WRK2(  107,I)-CONF( 2)*WRK2(   32,I)
      WRK2(  143,I)=WRK2(  108,I)-CONF( 2)*WRK2(   33,I)
      WRK2(  144,I)=WRK2(  109,I)-CONF( 2)*WRK2(   34,I)      
      WRK2(  145,I)=WRK2(  110,I)-CONF( 2)*WRK2(   35,I)
      WRK2(  146,I)=WRK2(  111,I)-CONF( 2)*WRK2(   36,I)
      WRK2(  147,I)=WRK2(  113,I)-CONF( 2)*WRK2(   37,I)
      WRK2(  148,I)=WRK2(  114,I)-CONF( 2)*WRK2(   38,I)
      WRK2(  149,I)=WRK2(  115,I)-CONF( 2)*WRK2(   39,I)
      WRK2(  150,I)=WRK2(  116,I)-CONF( 2)*WRK2(   40,I)
      WRK2(  151,I)=WRK2(  118,I)-CONF( 2)*WRK2(   41,I)
      WRK2(  152,I)=WRK2(  119,I)-CONF( 2)*WRK2(   42,I)
      WRK2(  153,I)=WRK2(  120,I)-CONF( 2)*WRK2(   43,I)
      WRK2(  154,I)=WRK2(  122,I)-CONF( 2)*WRK2(   44,I)
      WRK2(  155,I)=WRK2(  123,I)-CONF( 2)*WRK2(   45,I)
      WRK2(  156,I)=WRK2(  125,I)-CONF( 2)*WRK2(   46,I)
      WRK2(  157,I)=WRK2(  112,I)-CONF( 3)*WRK2(   32,I)
      WRK2(  158,I)=WRK2(  113,I)-CONF( 3)*WRK2(   33,I)
      WRK2(  159,I)=WRK2(  114,I)-CONF( 3)*WRK2(   34,I)
      WRK2(  160,I)=WRK2(  115,I)-CONF( 3)*WRK2(   35,I)
      WRK2(  161,I)=WRK2(  116,I)-CONF( 3)*WRK2(   36,I)
      WRK2(  162,I)=WRK2(  117,I)-CONF( 3)*WRK2(   37,I)
      WRK2(  163,I)=WRK2(  118,I)-CONF( 3)*WRK2(   38,I)
      WRK2(  164,I)=WRK2(  119,I)-CONF( 3)*WRK2(   39,I)
      WRK2(  165,I)=WRK2(  120,I)-CONF( 3)*WRK2(   40,I)
      WRK2(  166,I)=WRK2(  121,I)-CONF( 3)*WRK2(   41,I)
      WRK2(  167,I)=WRK2(  122,I)-CONF( 3)*WRK2(   42,I)
      WRK2(  168,I)=WRK2(  123,I)-CONF( 3)*WRK2(   43,I)
      WRK2(  169,I)=WRK2(  124,I)-CONF( 3)*WRK2(   44,I)
      WRK2(  170,I)=WRK2(  125,I)-CONF( 3)*WRK2(   45,I)
      WRK2(  171,I)=WRK2(  126,I)-CONF( 3)*WRK2(   46,I)
      END DO
      RETURN
      END
