working cluster for nano parameters

[unres.git] / source / unres / src_MD-M / parmread.F

diff --git a/source/unres/src_MD-M/parmread.F b/source/unres/src_MD-M/parmread.F
index dc040cb..068897d 100644 (file)
--- a/source/unres/src_MD-M/parmread.F
+++ b/source/unres/src_MD-M/parmread.F
@@ -943,13 +943,48 @@ cc maxinter is maximum interaction sites
 #else
       read (isccor,*,end=119,err=119) (isccortyp(i),i=1,ntyp)
 c      write (iout,*) 'ntortyp',ntortyp
+      do i=-ntyp,-1
+        isccortyp(i)=-isccortyp(-i)
+      enddo
       maxinter=3
 cc maxinter is maximum interaction sites
+C NOW READING FOR LL aminoacid chirality
       do l=1,maxinter
       do i=1,nsccortyp
         do j=1,nsccortyp
           read (isccor,*,end=119,err=119)
      & nterm_sccor(i,j),nlor_sccor(i,j)
+          nterm_sccor(-i,-j)=nterm_sccor(i,j)
+          v0ijsccor=0.0d0
+          si=-1.0d0
+
+          do k=1,nterm_sccor(i,j)
+            read (isccor,*,end=119,err=119) kk,v1sccor(k,l,i,j)
+     &    ,v2sccor(k,l,i,j)
+            v0ijsccor=v0ijsccor+si*v1sccor(k,l,i,j)
+            si=-si
+            v1sccor(k,l,-i,-j)=v1sccor(k,l,i,j)
+            v2sccor(k,l,-i,-j)=-v2sccor(k,l,i,j)
+
+          enddo
+          do k=1,nlor_sccor(i,j)
+            read (isccor,*,end=119,err=119) kk,vlor1sccor(k,i,j),
+     &        vlor2sccor(k,i,j),vlor3sccor(k,i,j)
+            v0ijsccor=v0ijsccor+vlor1sccor(k,i,j)/
+     &(1+vlor3sccor(k,i,j)**2)
+          enddo
+          v0sccor(l,i,j)=v0ijsccor
+          v0sccor(l,-i,-j)=v0ijsccor
+        enddo
+      enddo
+      enddo
+C NOW READING FOR LD aminoacid chirality
+      do l=1,maxinter
+      do i=1,nsccortyp
+        do j=-nsccortyp,-1
+          read (isccor,*,end=119,err=119)
+     & nterm_sccor(i,j),nlor_sccor(i,j)
+          nterm_sccor(-i,-j)=nterm_sccor(i,j)
           v0ijsccor=0.0d0
           si=-1.0d0
 
@@ -958,6 +993,9 @@ cc maxinter is maximum interaction sites
      &    ,v2sccor(k,l,i,j)
             v0ijsccor=v0ijsccor+si*v1sccor(k,l,i,j)
             si=-si
+            v1sccor(k,l,-i,-j)=v1sccor(k,l,i,j)
+            v2sccor(k,l,-i,-j)=-v2sccor(k,l,i,j)
+
           enddo
           do k=1,nlor_sccor(i,j)
             read (isccor,*,end=119,err=119) kk,vlor1sccor(k,i,j),
@@ -966,9 +1004,23 @@ cc maxinter is maximum interaction sites
      &(1+vlor3sccor(k,i,j)**2)
           enddo
           v0sccor(l,i,j)=v0ijsccor
+          v0sccor(l,-i,-j)=v0ijsccor
         enddo
       enddo
       enddo
+C symetry for GLY
+             i=10
+             do l=1,maxinter
+                do j=-nsccortyp,-1
+                do k=1,nterm_sccor(i,j)
+                v1sccor(k,l,10,j)=v1sccor(k,l,10,-j)
+                v2sccor(k,l,10,j)=-v2sccor(k,l,10,-j)
+                v1sccor(k,l,j,10)=v1sccor(k,l,-j,10)
+                v2sccor(k,l,j,10)=-v2sccor(k,l,-j,10)
+               enddo
+              enddo
+             enddo
+
       close (isccor)
 
 #endif