1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
SUBROUTINE CHRGE(P,Q)
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
INCLUDE 'SIZES'
DIMENSION P(*),Q(*)
COMMON /MOLKST/ NUMAT,NAT(NUMATM),NFIRST(NUMATM),NMIDLE(NUMATM),
1 NLAST(NUMATM), NORBS, NELECS,NALPHA,NBETA,
2 NCLOSE,NOPEN,NDUMY,FRACT
C***********************************************************************
C
C CHRGE STORES IN Q THE TOTAL ELECTRON DENSITIES ON THE ATOMS
C
C ON INPUT P = DENSITY MATRIX
C
C ON OUTPUT Q = ATOM ELECTRON DENSITIES
C
C***********************************************************************
K=0
DO 10 I=1,NUMAT
IA=NFIRST(I)
IB=NLAST(I)
Q(I)=0.D0
DO 10 J=IA,IB
K=K+J
10 Q(I)=Q(I)+P(K)
RETURN
END
|
