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),
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
11 & gacontp_hb1(3,maxconts,maxres),gacontp_hb2(3,maxconts,maxres),
12 & gacontp_hb3(3,maxconts,maxres),
13 & gacontm_hb1(3,maxconts,maxres),gacontm_hb2(3,maxconts,maxres),
14 & gacontm_hb3(3,maxconts,maxres),
15 & gacont_hbr(3,maxconts,maxres),
16 & grij_hb_cont(3,maxconts,maxres),
17 & facont_hb(maxconts,maxres),ees0p(maxconts,maxres),
18 & ees0m(maxconts,maxres),d_cont(maxconts,maxres),
19 & num_cont_hb(maxres),jcont_hb(maxconts,maxres)
20 C 9/23/99 Added improper rotation matrices and matrices of dipole-dipole
22 c 7/25/08 Commented out; not needed when cumulants used
23 C Interactions of pseudo-dipoles generated by loc-el interactions.
25 double precision dip,dipderg,dipderx
26 common /dipint/ dip(4,maxconts,maxres),dipderg(4,maxconts,maxres),
27 & dipderx(3,5,4,maxconts,maxres)
29 C 10/30/99 Added other pre-computed vectors and matrices needed
30 C to calculate three - six-order el-loc correlation terms
31 double precision Ug,Ugder,Ug2,Ug2der,obrot,obrot2,obrot_der,
32 & obrot2_der,Ub2,Ub2der,mu,muder,EUg,EUgder,CUg,CUgder,
33 & DUg,DUgder,DtUg2,DtUg2der,Ctobr,Ctobrder,Dtobr2,Dtobr2der
34 common /rotat/ Ug(2,2,maxres),Ugder(2,2,maxres),Ug2(2,2,maxres),
35 & Ug2der(2,2,maxres),obrot(2,maxres),obrot2(2,maxres),
36 & obrot_der(2,maxres),obrot2_der(2,maxres)
37 C This common block contains vectors and matrices dependent on a single
39 common /precomp1/ mu(2,maxres),muder(2,maxres),Ub2(2,maxres),
40 & Ub2der(2,maxres),Ctobr(2,maxres),Ctobrder(2,maxres),
41 & Dtobr2(2,maxres),Dtobr2der(2,maxres),
42 & EUg(2,2,maxres),EUgder(2,2,maxres),CUg(2,2,maxres),
43 & CUgder(2,2,maxres),DUg(2,2,maxres),Dugder(2,2,maxres),
44 & DtUg2(2,2,maxres),DtUg2der(2,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
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),Ug2DtEUg(2,2,maxres),
53 & Ug2DtEUgder(2,2,2,maxres),DtUg2EUgder(2,2,2,maxres)
54 double precision costab,sintab,costab2,sintab2
55 common /rotat_old/ costab(maxres),sintab(maxres),
56 & costab2(maxres),sintab2(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,g_contij,ekont
65 common /diploc/ AEA(2,2,2),AEAderg(2,2,2),AEAderx(2,2,3,5,2,2),
66 & EAEA(2,2,2), EAEAderg(2,2,2,2), EAEAderx(2,2,3,5,2,2),
67 & AECA(2,2,2),AECAderg(2,2,2),AECAderx(2,2,3,5,2,2),
68 & ADtEA(2,2,2),ADtEAderg(2,2,2,2),ADtEAderx(2,2,3,5,2,2),
69 & ADtEA1(2,2,2),ADtEA1derg(2,2,2,2),ADtEA1derx(2,2,3,5,2,2),
70 & AEAb1(2,2,2),AEAb1derg(2,2,2),AEAb1derx(2,3,5,2,2,2),
71 & AEAb2(2,2,2),AEAb2derg(2,2,2,2),AEAb2derx(2,3,5,2,2,2),
73 C 12/13/2008 (again Poland-Jaruzel war anniversary)
74 C RE: Parallelization of 4th and higher order loc-el correlations
75 integer ncont_sent,ncont_recv,iint_sent,iisent_local,
76 & itask_cont_from,itask_cont_to,ntask_cont_from,ntask_cont_to,
77 & nat_sent,iat_sent,iturn3_sent,iturn4_sent,iturn3_sent_local,
79 common /contdistrib/ ncont_sent(maxres),ncont_recv(maxres),
80 & iint_sent(4,maxres,maxres),iint_sent_local(4,maxres,maxres),
81 & iturn3_sent(4,maxres),iturn4_sent(4,maxres),
82 & iturn3_sent_local(4,maxres),iturn4_sent_local(4,maxres),
83 & nat_sent,iat_sent(maxres),itask_cont_from(0:max_fg_procs-1),
84 & itask_cont_to(0:max_fg_procs-1),ntask_cont_from,ntask_cont_to