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
29
30
31
32
33
34
35
36
|
SUBROUTINE DEPVAR (A,I,W,L)
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
DIMENSION A(3,*)
C***********************************************************************
C
C IN SUBROUTINE HADDON WHEN M, THE SYMMETRY OPERATION, IS 18 DEPVAR IS
C CALLED. DEPVAR SHOULD THEN CONTAIN A USER-WRITTEN SYMMETRY OPERATION.
C SEE HADDON TO GET THE IDEA ON HOW TO WRITE DEPVAR.
C
C ON INPUT:
C A = ARRAY OF INTERNAL COORDINATES
C I = ADDRESS OF REFERENCE ATOM
C ON OUTPUT:
C L = 1 (IF A BOND-LENGTH IS THE DEPENDENT FUNCTION)
C = 2 (IF AN ANGLE IS THE DEPENDENT FUNCTION)
C = 3 (IF A DIHEDRAL ANGLE IS THE DEPENDENT FUNCTION)
C W = VALUE OF THE FUNCTION
C
C NOTE: IT IS THE WRITER'S RESPONSIBILITY TO MAKE CERTAIN THAT THE
C SUBROUTINE DOES NOT CONTAIN ANY ERRORS!
C***********************************************************************
COMMON /KEYWRD/ KEYWRD
COMMON /NUMCAL/ NUMCAL
SAVE FACT
CHARACTER*241 KEYWRD
DATA ICALCN/0/
IF (ICALCN.NE.NUMCAL) THEN
ICALCN=NUMCAL
FACT=READA(KEYWRD,INDEX(KEYWRD,'DEPVAR'))
WRITE(6,'('' UNIT CELL LENGTH ='',F14.7,
1'' TIMES BOND LENGTH'')')FACT
ENDIF
W=A(1,I)*FACT
L=1
RETURN
END
|
