D i L aminokwasy dzialajace dla src_CSA

[unres.git] / source / unres / src_CSA_DiL / initialize_p.F

diff --git a/source/unres/src_CSA_DiL/initialize_p.F b/source/unres/src_CSA_DiL/initialize_p.F
index 19cf3d6..bde8c8f 100644 (file)
--- a/source/unres/src_CSA_DiL/initialize_p.F
+++ b/source/unres/src_CSA_DiL/initialize_p.F
@@ -164,8 +164,12 @@ c      call memmon_print_usage()
        rr0(i)=0.0D0
        a0thet(i)=0.0D0
        do j=1,2
-         athet(j,i)=0.0D0
-         bthet(j,i)=0.0D0
+         do ichir1=-1,1
+          do ichir2=-1,1
+         athet(j,i,ichir1,ichir2)=0.0D0
+         bthet(j,i,ichir1,ichir2)=0.0D0
+          enddo
+         enddo
         enddo
        do j=0,3
          polthet(j,i)=0.0D0
@@ -191,15 +195,36 @@ 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 j=1,maxtor
+       do iblock=1,2
+       do j=-maxtor,maxtor
          do k=1,maxterm
            v1(k,j,i)=0.0D0
            v2(k,j,i)=0.0D0
           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 !i
       do i=1,maxres
        itype(i)=0
        itel(i)=0