X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD%2Finitialize_p.F;h=62706fba9a6bf6b405ea9f0e6fb5cb10c173d441;hb=d07177b41ad32fa404c1c6cff69df2a83d9cbbc8;hp=8880628a529a88e24f5dc71a82aeb45603076b44;hpb=c78893ac5f83aa71d686ad2b1fc4393f03228ede;p=unres.git

diff --git a/source/unres/src_MD/initialize_p.F b/source/unres/src_MD/initialize_p.F
index 8880628..62706fb 100644
--- a/source/unres/src_MD/initialize_p.F
+++ b/source/unres/src_MD/initialize_p.F
@@ -191,10 +191,10 @@ c      call memmon_print_usage()
       enddo
       nlob(ntyp1)=0
       dsc(ntyp1)=0.0D0
-      do i=1,maxtor
+      do i=-maxtor,maxtor
 	itortyp(i)=0
        do iblock=1,2
-	do j=1,maxtor
+	do j=-maxtor,maxtor
 	  do k=1,maxterm
 	    v1(k,j,i,iblock)=0.0D0
 	    v2(k,j,i,iblock)=0.0D0
@@ -202,6 +202,26 @@ c      call memmon_print_usage()
         enddo
         enddo
       enddo
+      do iblock=1,2
+       do i=-maxtor,maxtor
+        do j=-maxtor,maxtor
+         do k=-maxtor,maxtor
+          do l=1,maxtermd_1
+            v1c(1,l,i,j,k,iblock)=0.0D0
+            v1s(1,l,i,j,k,iblock)=0.0D0
+            v1c(2,l,i,j,k,iblock)=0.0D0
+            v1s(2,l,i,j,k,iblock)=0.0D0
+          enddo !l
+          do l=1,maxtermd_2
+           do m=1,maxtermd_2
+            v2c(m,l,i,j,k,iblock)=0.0D0
+            v2s(m,l,i,j,k,iblock)=0.0D0
+           enddo !m
+          enddo !l
+        enddo !k
+       enddo !j
+      enddo !i
+      enddo !iblock
       do i=1,maxres
 	itype(i)=0
 	itel(i)=0