20      &            xGoC(npro,4*nshl,nshl)
38           Diagl(:,i,4) = xGoC(1:npro,nn2+i,i)
1120 	subroutine flesApG ( ien, xGoC, lesP, lesQ, nPs, nQs )  argument
1124         dimension xGoC(npro,4*nshl,nshl)
1144      &                    + xGoC(1:npro,i0+1,j) * Ptemp(:,j,nPs)
1147      &                    + xGoC(1:npro,i0+2,j) * Ptemp(:,j,nPs)
1150      &                    + xGoC(1:npro,i0+3,j) * Ptemp(:,j,nPs) 
1166  	subroutine flesApKG ( ien, xKebe, xGoC, lesP, lesQ, nPs, nQs )  argument
1171      &       xGoC(npro,4*nshl,nshl)
1217      &                   + xGoC(1:npro,i0+1,j) * Ptemp(:,j,nPs)
1219      &                   + xGoC(1:npro,i0+2,j) * Ptemp(:,j,nPs)
1221      &                   + xGoC(1:npro,i0+3,j) * Ptemp(:,j,nPs)
1237 	subroutine flesApNGt ( ien, xGoC, lesP, lesQ, nPs, nQs )  argument
1241         dimension ien(npro,nshl), xGoC(npro,4*nshl,nshl)
1262      &                      - xGoC(1:npro,i0+1,i) * Ptemp(:,j,1)
1263      &                      - xGoC(1:npro,i0+2,i) * Ptemp(:,j,2)
1264      &                      - xGoC(1:npro,i0+3,i) * Ptemp(:,j,3)
1280 	subroutine flesApNGtC ( ien, xGoC, lesP, lesQ, nPs, nQs )  argument
1284         dimension ien(npro,nshl), xGoC(npro,4*nshl,nshl)
1305      &                      - xGoC(1:npro,i0+1,i) * Ptemp(:,j,1)
1306      &                      - xGoC(1:npro,i0+2,i) * Ptemp(:,j,2)
1307      &                      - xGoC(1:npro,i0+3,i) * Ptemp(:,j,3)
1321      &                      + xGoC(1:npro,i0,j) * Ptemp(:,j,nPs)
1337 	subroutine flesApFull ( ien, xKebe, xGoC, lesP, lesQ, nPs, nQs )  argument
1343      &       xGoC(npro,4*nshl,nshl)
1388      &                   + xGoC(1:npro,i0+1,j) * Ptemp(:,j,nPs)
1390      &                   + xGoC(1:npro,i0+2,j) * Ptemp(:,j,nPs)
1392      &                   + xGoC(1:npro,i0+3,j) * Ptemp(:,j,nPs)
1404      &                      - xGoC(1:npro,i0+1,i) * Ptemp(:,j,1)
1405      &                      - xGoC(1:npro,i0+2,i) * Ptemp(:,j,2)
1406      &                      - xGoC(1:npro,i0+3,i) * Ptemp(:,j,3)
1420      &                      + xGoC(1:npro,i0,j) * Ptemp(:,j,nPs)