update

[unres.git] / source / wham / src-M / include_unres / COMMON.TORSION

diff --git a/source/wham/src-M/include_unres/COMMON.TORSION b/source/wham/src-M/include_unres/COMMON.TORSION
index 55cc7f4..cd576c8 100644 (file)
--- a/source/wham/src-M/include_unres/COMMON.TORSION
+++ b/source/wham/src-M/include_unres/COMMON.TORSION
@@ -1,25 +1,60 @@
 C Torsional constants of the rotation about virtual-bond dihedral angles
-      double precision v1,v2,vlor1,vlor2,vlor3,v0
-      integer itortyp,ntortyp,nterm,nlor,nterm_old
-      common/torsion/v0(maxtor,maxtor),v1(maxterm,maxtor,maxtor),
-     &    v2(maxterm,maxtor,maxtor),vlor1(maxlor,maxtor,maxtor),
+      double precision v1,v2,vlor1,vlor2,vlor3,v0,v1_kcc,v2_kcc,
+     & v11_chyb,v21_chyb,v12_chyb,v22_chyb,v1bend_chyb
+      integer itortyp,ntortyp,nterm,nlor,nterm_old,nterm_kcc_Tb,
+     &   nterm_kcc,itortyp_kcc,nbend_kcc_Tb
+      common/torsion/v0(-maxtor:maxtor,-maxtor:maxtor,2),
+     &    v1(maxterm,-maxtor:maxtor,-maxtor:maxtor,2),
+     &    v2(maxterm,-maxtor:maxtor,-maxtor:maxtor,2),
+     &    vlor1(maxlor,-maxtor:maxtor,-maxtor:maxtor),
      &    vlor2(maxlor,maxtor,maxtor),vlor3(maxlor,maxtor,maxtor),
-     &    itortyp(ntyp),ntortyp,nterm(maxtor,maxtor),
-     &    nlor(maxtor,maxtor),nterm_old
+     &    v1_kcc(maxval_kcc,maxval_kcc,maxtor_kcc,
+     &              -maxtor:maxtor,-maxtor:maxtor),
+     &    v2_kcc(maxval_kcc,maxval_kcc,maxtor_kcc,
+     &              -maxtor:maxtor,-maxtor:maxtor),
+     &    v1bend_chyb(0:maxang_kcc,-maxtor:maxtor),
+     &    itortyp(-ntyp1:ntyp1),ntortyp,
+     &    itortyp_kcc(-ntyp1:ntyp1),
+     &    nterm(-maxtor:maxtor,-maxtor:maxtor,2),
+     &    nlor(-maxtor:maxtor,-maxtor:maxtor,2), 
+     &    nterm_kcc_Tb(-maxtor:maxtor,-maxtor:maxtor),
+     &    nterm_kcc(-maxtor:maxtor,-maxtor:maxtor),
+     &    nbend_kcc_Tb(-maxtor:maxtor),
+     &    nterm_old
 C 6/23/01 - constants for double torsionals
       double precision v1c,v1s,v2c,v2s
       integer ntermd_1,ntermd_2
-      common /torsiond/ v1c(2,maxtermd_1,maxtor,maxtor,maxtor),
-     &    v1s(2,maxtermd_1,maxtor,maxtor,maxtor),
-     &    v2c(maxtermd_2,maxtermd_2,maxtor,maxtor,maxtor),
-     &    v2s(maxtermd_2,maxtermd_2,maxtor,maxtor,maxtor),
-     &    ntermd_1(maxtor,maxtor,maxtor),ntermd_2(maxtor,maxtor,maxtor)
+      common /torsiond/ 
+     &v1c(2,maxtermd_1,-maxtor:maxtor,-maxtor:maxtor,-maxtor:maxtor,2),
+     &v1s(2,maxtermd_1,-maxtor:maxtor,-maxtor:maxtor,-maxtor:maxtor,2),
+     &v2c(maxtermd_2,maxtermd_2,-maxtor:maxtor,-maxtor:maxtor,
+     &    -maxtor:maxtor,2),
+     &v2s(maxtermd_2,maxtermd_2,-maxtor:maxtor,-maxtor:maxtor,
+     &    -maxtor:maxtor,2),
+     &    ntermd_1(-maxtor:maxtor,-maxtor:maxtor,-maxtor:maxtor,2),
+     &    ntermd_2(-maxtor:maxtor,-maxtor:maxtor,-maxtor:maxtor,2)
 C 9/18/99 - added Fourier coeffficients of the expansion of local energy 
 C           surface
-      double precision b1,b2,cc,dd,ee,ctilde,dtilde,b1tilde
-      integer nloctyp
-      common/fourier/ b1(2,maxtor),b2(2,maxtor),cc(2,2,maxtor),
-     &    dd(2,2,maxtor),ee(2,2,maxtor),ctilde(2,2,maxtor),
-     &    dtilde(2,2,maxtor),b1tilde(2,maxtor),nloctyp
-      double precision b
-      common /fourier1/ b(13,maxtor)
+      double precision b1,b2,cc,dd,ee,ctilde,dtilde,b2tilde,b1tilde,
+     & b,bnew1,bnew2,ccold,ddold,ccnew,ddnew,eenew,e0new,gtb1,gtb2,
+     & eeold,gtcc,gtdd,gtee,
+     & bnew1tor,bnew2tor,ccnewtor,ddnewtor,eenewtor,e0newtor
+      integer nloctyp,iloctyp(-ntyp1:ntyp1),itype2loc(-ntyp1:ntyp1)
+      logical SPLIT_FOURIERTOR
+      common/fourier/ b1(2,maxres),b2(2,maxres),b(13,-ntyp:ntyp),
+     &    bnew1(3,2,-ntyp:ntyp),bnew2(3,2,-ntyp:ntyp),
+     &    ccnew(3,2,-ntyp:ntyp),ddnew(3,2,-ntyp:ntyp),
+     &    bnew1tor(3,2,-ntyp:ntyp),bnew2tor(3,2,-ntyp:ntyp),
+     &    ccnewtor(3,2,-ntyp:ntyp),ddnewtor(3,2,-ntyp:ntyp),
+     &    ccold(2,2,-ntyp:ntyp),ddold(2,2,-ntyp:ntyp),
+     &    cc(2,2,maxres),
+     &    dd(2,2,maxres),eeold(2,2,-ntyp:ntyp),
+     &    e0new(3,-ntyp:ntyp),eenew(2,2,2,-ntyp:ntyp),
+     &    e0newtor(3,-ntyp:ntyp),eenewtor(2,2,2,-ntyp:ntyp),
+     &    ee(2,2,maxres),
+     &    ctilde(2,2,maxres),
+     &    dtilde(2,2,maxres),b1tilde(2,maxres),
+     &    b2tilde(2,maxres),
+     &    gtb1(2,maxres),gtb2(2,maxres),gtCC(2,2,maxres),
+     &    gtDD(2,2,maxres),gtEE(2,2,maxres),
+     &    nloctyp,iloctyp,itype2loc,SPLIT_FOURIERTOR