max_template instead of fix number 19 in wham and cluster
[unres.git] / source / cluster / wham / src / COMMON.CONTACTS
1 C Change 12/1/95 - common block CONTACTS1 included.
2       integer ncont,ncont_ref,icont,icont_ref,num_cont,jcont
3       double precision facont,gacont
4       common /contacts/ ncont,ncont_ref,icont(2,maxcont),
5      &                  icont_ref(2,maxcont)
6       common /contacts1/ facont(maxconts,maxres),
7      &                  gacont(3,maxconts,maxres),
8      &                  num_cont(maxres),jcont(maxconts,maxres)
9 C 12/26/95 - H-bonding contacts
10       double precision gacontp_hb1,gacontp_hb2,gacontp_hb3,gacontm_hb1,
11      &  gacontm_hb2,gacontm_hb3,gacont_hbr,facont_hb,ees0p,ees0m,d_cont,
12      &  grij_hb_cont
13       integer num_cont_hb,jcont_hb
14       common /contacts_hb/ 
15      &  gacontp_hb1(3,maxconts,maxres),gacontp_hb2(3,maxconts,maxres),
16      &  gacontp_hb3(3,maxconts,maxres),
17      &  gacontm_hb1(3,maxconts,maxres),gacontm_hb2(3,maxconts,maxres),
18      &  gacontm_hb3(3,maxconts,maxres),
19      &  gacont_hbr(3,maxconts,maxres),
20      &  grij_hb_cont(3,maxconts,maxres),
21      &  facont_hb(maxconts,maxres),ees0p(maxconts,maxres),
22      &  ees0m(maxconts,maxres),d_cont(maxconts,maxres),
23      &  num_cont_hb(maxres),jcont_hb(maxconts,maxres)
24 C 9/23/99 Added improper rotation matrices and matrices of dipole-dipole 
25 C         interactions     
26 C Interactions of pseudo-dipoles generated by loc-el interactions.
27       double precision dip,dipderg,dipderx
28       common /dipint/ dip(4,maxconts,maxres),dipderg(4,maxconts,maxres),
29      &  dipderx(3,5,4,maxconts,maxres)
30 C 10/30/99 Added other pre-computed vectors and matrices needed 
31 C          to calculate three - six-order el-loc correlation terms
32       double precision Ug,Ugder,Ug2,Ug2der,obrot,obrot2,obrot_der,
33      &  obrot2_der,Ub2,Ub2der,mu,muder,EUg,EUgder,CUg,CUgder,
34      &  DUg,DUgder,DtUg2,DtUg2der,Ctobr,Ctobrder,Dtobr2,Dtobr2der
35       common /rotat/ Ug(2,2,maxres),Ugder(2,2,maxres),Ug2(2,2,maxres),
36      &  Ug2der(2,2,maxres),obrot(2,maxres),obrot2(2,maxres),
37      &  obrot_der(2,maxres),obrot2_der(2,maxres)
38 C This common block contains vectors and matrices dependent on a single
39 C amino-acid residue.
40       common /precomp1/ Ub2(2,maxres),Ub2der(2,maxres),mu(2,maxres),
41      &  EUg(2,2,maxres),EUgder(2,2,maxres),CUg(2,2,maxres),
42      &  CUgder(2,2,maxres),DUg(2,2,maxres),Dugder(2,2,maxres),
43      &  DtUg2(2,2,maxres),DtUg2der(2,2,maxres),Ctobr(2,maxres),
44      &  Ctobrder(2,maxres),Dtobr2(2,maxres),Dtobr2der(2,maxres)
45 C This common block contains vectors and matrices dependent on two
46 C consecutive amino-acid residues.
47       double precision Ug2Db1t,Ug2Db1tder,CUgb2,CUgb2der,EUgC,
48      &  EUgCder,EUgD,EUgDder,DtUg2EUg,DtUg2EUgder,Ug2DtEUg,Ug2DtEUgder
49       common /precomp2/ Ug2Db1t(2,maxres),Ug2Db1tder(2,maxres),
50      &  CUgb2(2,maxres),CUgb2der(2,maxres),EUgC(2,2,maxres),
51      &  EUgCder(2,2,maxres),EUgD(2,2,maxres),EUgDder(2,2,maxres),
52      &  DtUg2EUg(2,2,maxres),DtUg2EUgder(2,2,2,maxres),
53      &  Ug2DtEUg(2,2,maxres),Ug2DtEUgder(2,2,2,maxres)
54       double precision costab,sintab,costab2,sintab2
55       common /rotat_old/ costab(maxres),sintab(maxres),
56      &  costab2(maxres),sintab2(maxres),muder(2,maxres)
57 C This common block contains dipole-interaction matrices and their 
58 C Cartesian derivatives.
59       double precision a_chuj,a_chuj_der
60       common /dipmat/ a_chuj(2,2,maxconts,maxres),
61      &  a_chuj_der(2,2,3,5,maxconts,maxres)
62       double precision AEA,AEAderg,AEAderx,AECA,AECAderg,AECAderx,
63      &  ADtEA,ADtEAderg,ADtEAderx,AEAb1,AEAb1derg,AEAb1derx,
64      &  AEAb2,AEAb2derg,AEAb2derx,EAEA,EAEAderg,EAEAderx,
65      &  ADtEA1,ADtEA1derg,ADtEA1derx,g_contij,ekont
66       common /diploc/ AEA(2,2,2),AEAderg(2,2,2),AEAderx(2,2,3,5,2,2),
67      &  EAEA(2,2,2), EAEAderg(2,2,2,2), EAEAderx(2,2,3,5,2,2),
68      &  AECA(2,2,2),AECAderg(2,2,2),AECAderx(2,2,3,5,2,2),
69      &  ADtEA(2,2,2),ADtEAderg(2,2,2,2),ADtEAderx(2,2,3,5,2,2),
70      &  ADtEA1(2,2,2),ADtEA1derg(2,2,2,2),ADtEA1derx(2,2,3,5,2,2),
71      &  AEAb1(2,2,2),AEAb1derg(2,2,2),AEAb1derx(2,3,5,2,2,2),
72      &  AEAb2(2,2,2),AEAb2derg(2,2,2,2),AEAb2derx(2,3,5,2,2,2),
73      &  g_contij(3,2),ekont