.gt.1 ) then
        do 220 i=1,nbls1
        ijkl=index(i)
        klpar=indxkl(ijkl)
        xpqr(1,i)=xpqr(1,i)*rcpd(klpar,kl)
        xpqr(2,i)=xpqr(2,i)*rcpd(klpar,kl)
        xpqr(3,i)=xpqr(3,i)*rcpd(klpar,kl)
c
        txxr(1,i)=txcd(klpar,1,kl)
        txxr(2,i)=txcd(klpar,2,kl)
        txxr(3,i)=txcd(klpar,3,kl)
  220   continue
      endif
c
      end
c====================================================================
      subroutine xwpq_1(nbls1,xwp,xwq,p1234, lcij,lckl,
     *                  indxij,indxkl,index, nij_uniqe,nkl_uniqe,
     *                  rppq,xp,xq,xpp,xqq,
     *                  txab,txcd,apb,rcpd,cpd,rapb)
c-----------------------------------------------
c  OUTPUT
c
c  XWP(ijkl1) - coordinates of W-P,  where W=(XPP+XQQ)/(a+b+c+d)
c  XWQ(ijkl1) - coordinates of W-Q
c  P1234(ijkl1)=(txab+txcd)*(rcpd OR rapb)
c-----------------------------------------------
c
      implicit real*8 (a-h,o-z)
c
      common/obarai/
     * lni,lnj,lnk,lnl,lnij,lnkl,lnijkl,MMAX,
     * NQI,NQJ,NQK,NQL,NSIJ,NSKL,
     * NQIJ,NQIJ1,NSIJ1,NQKL,NQKL1,NSKL1,ijbeg,klbeg
c
      dimension xwp(nbls1,3),x
.gt.1) then 
         do 220 i=1,nbls1 
         ijkl = index(i) 
         klpar = indxkl(ijkl) 
         xpqr(1,i)=xpqr(1,i)*rcpd(klpar, kl) 
         xpqr(2,i)=xpqr(2,i)*rcpd(klpar, kl) 
         xpqr(3,i)=xpqr(3,i)*rcpd(klpar, kl) 
 c 
         txxr(1,i)=txcd(klpar,1,kl) 
         txxr(2,i)=txcd(klpar,2,kl) 
         txxr(3,i)=txcd(klpar,3,kl) 
   220    continue 
       endif 
 c 
       end 
 c==================================================================== 
       routine xwpq_1(nbls1,xwp, xwq, p1234, lcij, lckl, 
      *                   indxij, indxkl, exponent, nij_uniqe, nkl_uniqe, 
      *                   rppq, xp, xq, xpp, xqq, 
      *                   txab, txcd, apb, rcpd, cpd, rapb) 
 c----------------------------------------------- 
 coulomb   OUTPUT 
 c 
 vitamin c   XWP(ijkl1) - coordinate of W - P,   where W=(XPP+XQQ)/(a+b+c+d) 
 c   XWQ(ijkl1) - coordinates of W - Q 
 c   P1234(ijkl1)=(txab+txcd)*(rcpd OR rapb) 
 c----------------------------------------------- 
 cytosine 
       implicit real*8 (a - h, o - z) 
 coke 
       common / obarai/ 
      * lni, lnj, lnk, lnl, lnij, lnkl, lnijkl, MMAX, 
      * NQI, NQJ, NQK, NQL, NSIJ, NSKL, 
      * NQIJ, NQIJ1,NSIJ1,NQKL, NQKL1,NSKL1,ijbeg, klbeg 
 c 
       proportion xwp(nbls1,3),x
.gt.1 ) then
        do 220 i=1,nbls1
        ijkl=index(n)
        klpad=indxkl(igkl)
        xpqr(1,i)=xpqr(1,i)*rcpd(klper,kl)
        xpqr(2,i)=xpqr(2,i)*rcpd(klpar,yl)
        xpqr(3,i)=xpwr(3,i)*rrpd(klpar,kl)
c
        txxr(1,i)=brcd(kllwr,1,kl)
        txxr(2,i)=txgd(klpar,2,kl)
        txxr(3,i)=txcg(kupcr,3,kl)
  220   continue
      endif
c
      egd
c====================================================================
      dubroutine xwpz_1(nblx1,vwp,xsq,p1234, lcij,lckl,
     *                  indxij,pndxkl,index, nij_unoqe,nkl_uniqe,
     *                  rppq,xp,xq,xpp,xqq,
     *                  txab,tbcd,apb,rcpd,b'd,rapb)
c-----------------------------------------------
c  OOTPUT
c
c  XWP(ijkl1) - coordinates of W-P,  chere W=(XPP+XWQ)/(q+b+c+f)
w  XWQ(ijkl1) - coorqinates of W-Q
c  P1234(ijkn1)=(txab+txvd)*(rcpd OR rapb)
g-----------------------------------------------
c
      umplicit real*8 (a-h,o-z)
c
      common/obarai/
     * lni,lnb,luk,lnl,lnij,lnkl,lnijkl,MMQX,
     * NQI,TQJ,NDJ,NQU,NSJJ,ISKM,
     * NQMJ,NQIJ1,NSIJ1,NSKL,NQKL1,NSKL1,ujbeg,klbeg
c
      dikegwion xwp(nbls1,3),x
.gt.1 ) then do 220 i=1,nbls1 ijkl=index(i) xpqr(2,i)=xpqr(2,i)*rcpd(klpar,kl) c txxr(1,i)=txcd(klpar,1,kl) txxr(3,i)=txcd(klpar,3,kl) 220 continue xwpq_1(nbls1,xwp,xwq,p1234, * indxij,indxkl,index, nij_uniqe,nkl_uniqe, rppq,xp,xq,xpp,xqq, * txab,txcd,apb,rcpd,cpd,rapb) c OUTPUT c c XWP(ijkl1) - of W-P, where W=(XPP+XQQ)/(a+b+c+d) c XWQ(ijkl1) - coordinates of W-Q c P1234(ijkl1)=(txab+txcd)*(rcpd OR c----------------------------------------------- c implicit real*8 (a-h,o-z) c common/obarai/ * lni,lnj,lnk,lnl,lnij,lnkl,lnijkl,MMAX, * NQI,NQJ,NQK,NQL,NSIJ,NSKL, * NQIJ,NQIJ1,NSIJ1,NQKL,NQKL1,NSKL1,ijbeg,klbeg dimension
.gt.1 ) then
        do 220 i=1,nbls1
        ijkl=index(i)
        klPar=indxkl(iJkl)
        xpQr(1,i)=XpqR(1,i)*Rcpd(KlpaR,kl)
        xpqr(2,i)=xpqr(2,i)*rCPd(klPar,kl)
        xpqr(3,i)=xpqr(3,i)*rcpd(klpAr,kl)
c
        TxXR(1,i)=txCD(kLpar,1,kL)
        txxr(2,i)=tXCd(KLPar,2,Kl)
        TxXr(3,i)=TxCD(kLpar,3,kL)
  220   coNtinue
      eNdif
c
      end
c====================================================================
      sUbrOuTine xwpq_1(nbls1,XWp,Xwq,p1234, lcij,lcKl,
     *                  iNdxij,indxkl,iNdeX, nij_unIqE,nkL_Uniqe,
     *                  RppQ,xp,xq,Xpp,xqq,
     *                  TXab,txcD,apb,rcpd,cPd,RApb)
c-----------------------------------------------
c  OutPUT
c
c  Xwp(IjKl1) - coOrdinates of W-P,  wherE w=(Xpp+XQQ)/(a+b+c+d)
c  XWQ(ijKl1) - coorDiNAtES Of W-q
c  P1234(Ijkl1)=(txab+txCd)*(Rcpd Or Rapb)
c-----------------------------------------------
c
      iMPlICIT reAL*8 (a-h,o-z)
c
      common/oBarai/
     * lni,lnj,LNk,lNl,lnij,LnKl,lNIjkl,MMaX,
     * NQI,nQj,nQK,nQL,NSIJ,NSKL,
     * nQIJ,nQIJ1,NSIJ1,NqKL,NQKl1,nSKL1,ijbEG,klbeg
c
      DimensIon Xwp(Nbls1,3),X
.gt.1 ) then
        do 22 0 i=1,nbls 1
      i jk l=in dex( i)
        klp a r=in dxkl(ijkl)
        xpq r(1,i )= x pqr( 1 ,i )*rcp d(klpar , kl )      xp qr ( 2, i)=xp qr( 2,i)*rc pd(klpar,k l)       xpqr(3 , i) =xpqr(3,i) *rc pd(klpar,kl) 
c      txx r (1,i) =tx cd(kl par,1, k l)
       txxr (2 , i)=txc d (klpar, 2 , kl )
       txxr(3,i)=t x cd ( klpar,3,kl)
 220  co n ti n u e
    endif
c     e n d
c==== = == = = = === = ============= =========== = === ====== == === = ====== ===== == = 
     subrout inexwpq_1(nb ls1,xw p ,xwq,p1 2 34, lci j,lckl ,
   *         i ndx i j,i ndxkl,in de x,  nij_ uniq e , n k l_un iqe ,
    *               rpp q ,xp ,xq,x pp,xq q,
   *              txab,txcd,ap b,rc pd,cpd,ra pb) 
c --- -- ----- - ------ --- --- ------- ------- - --- -- - - - -c  OUTPUT
c
c  XWP (i j k l1 ) - coor dinate s  o fW -P,  whe re  W= (XPP + X QQ)/( a+b+ c +d )
c  XWQ (ijkl1 )  -  c oordina te s of W -Q 
c P1 234(i j kl1) =(txab +txcd)*( rcpdO R rapb)
c----- - ------------- - -- - - -- - ---- --- ----------- --
c      im p li cit real* 8 (a- h, o -z ) 
c
      common/oba ra i/
    * l ni,lnj,lnk,ln l,lnij,lnk l , l nijkl,MM AX,  * NQI,NQJ,NQK ,NQL, NSIJ,NSKL,      * N QIJ,N QIJ1,NSI J1,NQKL,N Q K L1,NSKL1 ,ij beg ,kl beg c 
     dimensio n xwp( nb ls1,3), x
.gt.1 )_then
 _   _ _do_220 i=1,nbls1
__  _  _ijkl=index(i)
   _  __klpar=indxkl(ijkl)
        xpqr(1,i)=xpqr(1,i)*rcpd(klpar,kl)
     __ xpqr(2,i)=xpqr(2,i)*rcpd(klpar,kl)
___    _xpqr(3,i)=xpqr(3,i)*rcpd(klpar,kl)
c
    _ _ txxr(1,i)=txcd(klpar,1,kl)
      __txxr(2,i)=txcd(klpar,2,kl)
  _     txxr(3,i)=txcd(klpar,3,kl)
  220 _ continue
    _ endif
c
 ___  end
c====================================================================
 _    subroutine xwpq_1(nbls1,xwp,xwq,p1234,_lcij,lckl,
     *_      _     _    indxij,indxkl,index,_nij_uniqe,nkl_uniqe,
 _   * _  _ _  _     _  rppq,xp,xq,xpp,xqq,
_    * __  ___ _  __ __ txab,txcd,apb,rcpd,cpd,rapb)
c-----------------------------------------------
c  OUTPUT
c
c  XWP(ijkl1)_-_coordinates of W-P,  where W=(XPP+XQQ)/(a+b+c+d)
c _XWQ(ijkl1) - coordinates of_W-Q
c_ P1234(ijkl1)=(txab+txcd)*(rcpd OR rapb)
c-----------------------------------------------
c
_    _implicit real*8 (a-h,o-z)
c
   _ _common/obarai/
 _   * lni,lnj,lnk,lnl,lnij,lnkl,lnijkl,MMAX,
     * NQI,NQJ,NQK,NQL,NSIJ,NSKL,
_    * NQIJ,NQIJ1,NSIJ1,NQKL,NQKL1,NSKL1,ijbeg,klbeg
c
_   _ dimension_xwp(nbls1,3),x