constr_homol corrections in wham and cluster_wham to match unres energy
[unres.git] / source / cluster / wham / src / COMMON.DFA
1 C =======
2 C COMMON.DFA
3 C =======
4 C 2010/12/20 By Juyong Lee
5 C
6 c parameter
7 C [ 8 * ( Nres - 8 ) ] distance restraints 
8 C [ 2 * ( Nres - 8 ) ] angle restraints
9 C [ Nres ]             neighbor restraints
10 C Total : ~ 11 * Nres restraints
11 C
12 C
13       INTEGER IDFAMAX,IDFAMX2,IDFACMD,IDMAXMIN, MAXN
14       PARAMETER(IDFAMAX=4000,IDFAMX2=1000,IDFACMD=500,IDMAXMIN=500)
15       PARAMETER(MAXN=4)
16       real*8 wwdist,wwangle,wwnei
17       parameter(wwdist=1.0d0,wwangle=1.0d0,wwnei=1.0d0)
18
19 C IDFAMAX  - maximum number of DFA restraint including distance, angle and
20 C            number of neighbors ( Max of assign statement )
21 C IDFAMX2  - maximum number of atoms which are targets of restraints
22 C IDFACMD  - maximum number of 'DFA' command call
23 C IDMAXMIN - Maximum number of minima of dist, angle and neighbor info. from fragments
24 C MAXN     - Maximum Number of shell, currently 4
25 C MAXRES   - Maximum number of CAs
26
27 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
28 C INTEGER 
29 C DFANUM  - Number of ALL DFA restrants
30 c IDFA[DIS, PHI, THE, NEI] - NUMBER of restraints
31 c IDISNUM - number of minima for a distance restraint
32 c IPHINUM - number of minima for a phi angle restraint
33 c ITHENUM - number of minima for a theta angle restraint
34 c INEINUM - number of minima for a number of neighbors restraint
35
36 c IDISLIS - atom number of two atoms for distance restraint
37 c IPHILIS - atom numbers of four atoms for angle restraint
38 c ITHELIS - atom numbers of four atoms for angle restraint
39 c INEILIS - atom number of center of neighbor calculation
40 c JNEILIS - atom number of target of neighboring calculation
41 c JNEINUM - number of target atoms of neighboring term
42 C KSHELL  - SHELL number 
43
44 C ishiftca - index shift for CA atoms in UNRES (1 if the 1st aa != GLY)
45 C ilastca  - index of the last CA atom in UNRES (nres-1 if last aa != GLY)
46
47 C     old only for CHARMM
48 C STOAGDF - Store assign information ( How many assign within one command )
49 C NMAP    - mapping between dfanum and ndis, nphi, nthe, nnei
50
51       INTEGER IDFADIS,IDFAPHI,IDFATHE,IDFANEI,
52      &        IDISLIS,IPHILIS,ITHELIS,INEILIS,
53      &        IDISNUM,IPHINUM,ITHENUM,INEINUM,
54      &        FNEI,DFACMD, DFANUM,
55      &        NCA,ICAIDX,
56      &        STOAGDF, NMAP, IDFACAT, KDISNUM, KSHELL
57      &        ishiftca,ilastca 
58       COMMON /IDFA/ DFACMD, DFANUM,
59      &              IDFADIS, IDFAPHI, IDFANEI, IDFATHE, 
60      &              IDISNUM(IDFAMAX), IPHINUM(IDFAMAX), 
61      &              ITHENUM(IDFAMAX), INEINUM(IDFAMAX),
62      &              FNEI(IDFAMAX,IDMAXMIN), IDISLIS(2,IDFAMAX),
63      &              IPHILIS(5,IDFAMAX), ITHELIS(5,IDFAMAX),
64      &              INEILIS(IDFAMAX),
65      &               KSHELL(IDFAMAX),
66      &              IDFACAT(IDFACMD),
67      &              KDISNUM(IDFAMAX),
68      &              NCA, ICAIDX(MAXRES)
69       COMMON /IDFA2/ ishiftca,ilastca
70
71 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
72 C
73 C REAL VARIABLES
74 C
75 c SCC[DIST, PHI, THE] - weight of each calculations
76 c FDIST  - distance minima
77 C FPHI   - phi minima
78 c FTHE   - theta minima
79 C DFAEXP  : calculate expential function in advance
80 C
81       REAL*8 SCCDIST, SCCPHI, SCCTHE, SCCNEI, FDIST, FPHI1, FPHI2,
82      &       FTHE1, FTHE2,
83      &       DIS_INC, PHI_INC, THE_INC, NEI_INC, BETA_INC,
84      &       WSHET, EDFABET, 
85      &       CK, SCK, S1, S2
86 c    &       ,DFAEXP
87
88       COMMON /RDFA/ SCCDIST(IDFAMAX,IDMAXMIN),FDIST(IDFAMAX,IDMAXMIN),
89      &             SCCPHI(IDFAMAX,IDMAXMIN), SCCTHE(IDFAMAX,IDMAXMIN), 
90      &             SCCNEI(IDFAMAX,IDMAXMIN),
91      &             FPHI1(IDFAMAX,IDMAXMIN), FPHI2(IDFAMAX,IDMAXMIN),
92      &             FTHE1(IDFAMAX,IDMAXMIN), FTHE2(IDFAMAX,IDMAXMIN), 
93      &             DIS_INC, PHI_INC, THE_INC, NEI_INC, BETA_INC,
94      &             WSHET(MAXRES,MAXRES), EDFABET, 
95      &             CK(4),SCK(4),S1(4),S2(4)
96 c    &             ,DFAEXP(15001),