Adam's cluster & unres corrections

[unres.git] / source / unres / src-HCD-5D / initialize_p.F

diff --git a/source/unres/src-HCD-5D/initialize_p.F b/source/unres/src-HCD-5D/initialize_p.F
index c73426c..710f907 100644 (file)
--- a/source/unres/src-HCD-5D/initialize_p.F
+++ b/source/unres/src-HCD-5D/initialize_p.F
@@ -46,6 +46,7 @@ cMS$ATTRIBUTES C ::  proc_proc
       include 'COMMON.DERIV'
       include 'COMMON.SPLITELE'
       include 'COMMON.VAR'
+      include 'COMMON.MD'
 c Common blocks from the diagonalization routines
       integer IR,IW,IP,IJK,IPK,IDAF,NAV,IODA,KDIAG,ICORFL,IXDR
       integer i,idumm,j,k,l,ichir1,ichir2,iblock,m
@@ -68,7 +69,7 @@ c NaNQ initialization
       call proc_proc(rr,i)
 #endif
 #endif
-
+      itime_mat=0.
       kdiag=0
       icorfl=0
       iw=2
@@ -260,10 +261,10 @@ C Initialize the bridge arrays
       ns=0
       nss=0 
       nhpb=0
-      do i=1,maxss
+      do i=1,max_cyst
        iss(i)=0
       enddo
-      do i=1,maxdim
+      do i=1,maxdim_cont
        dhpb(i)=0.0D0
       enddo
       do i=1,maxres
@@ -445,6 +446,9 @@ c---------------------------------------------------------------------------
       include 'COMMON.TORCNSTR'
       include 'COMMON.IOUNITS'
       include 'COMMON.DERIV'
+#ifdef FOURBODY
+      include 'COMMON.CONTMAT'
+#endif
       include 'COMMON.CORRMAT'
       integer iturn3_start_all,iturn3_end_all,iturn4_start_all,
      & iturn4_end_all,iatel_s_all,