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 27 28
|
SUBROUTINE SYMTRY
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
INCLUDE 'SIZES'
COMMON /GEOM / GEO(3,NUMATM), XCOORD(3,NUMATM)
COMMON /GEOSYM/ NDEP, LOCPAR(MAXPAR), IDEPFN(MAXPAR),
1 LOCDEP(MAXPAR)
C**********************************************************************
C
C SYMTRY COMPUTES THE BOND LENGTHS AND ANGLES THAT ARE FUNCTIONS OF
C OTHER BOND LENGTHS AND ANGLES.
C
C ON INPUT GEO = KNOWN INTERNAL COORDINATES
C NDEP = NUMBER OF DEPENDENCY FUNCTIONS.
C IDEPFN = ARRAY OF DEPENDENCY FUNCTIONS.
C LOCDEP = ARRAY OF LABELS OF DEPENDENT ATOMS.
C LOCPAR = ARRAY OF LABELS OF REFERENCE ATOMS.
C
C ON OUTPUT THE ARRAY "GEO" IS FILLED
C***********************************************************************
C
C NOW COMPUTE THE DEPENDENT PARAMETERS.
C
DO 10 I=1,NDEP
CALL HADDON (VALUE,LOCN,IDEPFN(I),LOCPAR(I),GEO)
J=LOCDEP(I)
10 GEO(LOCN,J)=VALUE
RETURN
END
|