X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_CSA_DiL%2Fparmread.F;h=48da6c079c661cd6bc3b97d4a90962b0857d71b2;hb=81513f4efcd95dcfc5d3ec05ce6396a392d8ac24;hp=e58af2f42f0e295165d620c756cd61233d77590d;hpb=59403741e5c6541c3617c08c6e8e0c11ec394cce;p=unres.git

diff --git a/source/unres/src_CSA_DiL/parmread.F b/source/unres/src_CSA_DiL/parmread.F
index e58af2f..48da6c0 100644
--- a/source/unres/src_CSA_DiL/parmread.F
+++ b/source/unres/src_CSA_DiL/parmread.F
@@ -504,6 +504,17 @@ C
       read (itorp,*,end=113,err=113) (itortyp(i),i=1,ntyp)
       do iblock=1,2
       do i=-ntyp,-1
+<<<<<<< HEAD
+       itortyp(i)=-itortyp(-1)
+      enddo
+c      write (iout,*) 'ntortyp',ntortyp
+      do i=0,ntortyp,ntortyp-1
+	do j=-ntortyp,ntortyp
+	  read (itorp,*,end=113,err=113) nterm(i,j,iblock),
+     &          nlor(i,j,iblock)
+           nterm(-i,-j,iblock)=nterm(i,j,iblock)
+           nlor(-i,-j,iblock)=nlor(i,j,iblock)
+=======
        itortyp(i)=-itortyp(-i)
       enddo
 c      write (iout,*) 'ntortyp',ntortyp
@@ -513,13 +524,18 @@ c      write (iout,*) 'ntortyp',ntortyp
                  nlor(i,j,iblock)
           nterm(-i,-j,iblock)=nterm(i,j,iblock)
      &          nlor(i,j,iblock)
+>>>>>>> bartek2
           v0ij=0.0d0
           si=-1.0d0
 	  do k=1,nterm(i,j,iblock)
 	    read (itorp,*,end=113,err=113) kk,v1(k,i,j,iblock),
      &      v2(k,i,j,iblock)
             v1(k,-i,-j,iblock)=v1(k,i,j,iblock)
+<<<<<<< HEAD
+            v2(k,-i,-j,iblock)=-v2(k,i,j,iblock) 
+=======
             v2(k,-i,-j,iblock)=-v2(k,i,j,iblock)
+>>>>>>> bartek2
             v0ij=v0ij+si*v1(k,i,j,iblock)
             si=-si
 c         write(iout,*) i,j,k,iblock,nterm(i,j,iblock)
@@ -543,11 +559,12 @@ c      &v2(k,-i,-j,iblock),v2(k,i,j,iblock)
 	  do j=1,ntortyp
             write (iout,*) 'ityp',i,' jtyp',j
             write (iout,*) 'Fourier constants'
-            do k=1,nterm(i,j)
-	      write (iout,'(2(1pe15.5))') v1(k,i,j),v2(k,i,j)
+            do k=1,nterm(i,j,iblock)
+	      write (iout,'(2(1pe15.5))') v1(k,i,j,iblock),
+     &         v2(k,i,j,iblock)
             enddo
             write (iout,*) 'Lorenz constants'
-            do k=1,nlor(i,j)
+            do k=1,nlor(i,j,iblock)
 	      write (iout,'(3(1pe15.5))') 
      &         vlor1(k,i,j),vlor2(k,i,j),vlor3(k,i,j)
             enddo