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 b80095b..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      
@@ -1038,22 +1090,26 @@ C
         write (iout,*) 'Type',i
         write (iout,'(a,i2,a,f10.5)') ('b(',ii,')=',b(ii,i),ii=1,13)
         endif
-c        B1(1,i)  = b(3)
-c        B1(2,i)  = b(5)
-c        B1(1,-i) = b(3)
-c        B1(2,-i) = -b(5)
+        b(3,-i)= b(3,i)
+        b(5,-i)=-b(5,i)
+        b(4,-i)=-b(4,i)
+        b(2,-i)= b(2,i)     
+        B1(1,i)  = b(3,i)
+        B1(2,i)  = b(5,i)
+        B1(1,-i) = b(3,i)
+        B1(2,-i) = -b(5,i)
 c        b1(1,i)=0.0d0
 c        b1(2,i)=0.0d0
-c        B1tilde(1,i) = b(3)
-c        B1tilde(2,i) =-b(5)
-c        B1tilde(1,-i) =-b(3)
-c        B1tilde(2,-i) =b(5)
+        B1tilde(1,i) = b(3,i)
+        B1tilde(2,i) =-b(5,i)
+        B1tilde(1,-i) =-b(3,i)
+        B1tilde(2,-i) =b(5,i)
 c        b1tilde(1,i)=0.0d0
 c        b1tilde(2,i)=0.0d0
-c        B2(1,i)  = b(2)
-c        B2(2,i)  = b(4)
-c        B2(1,-i)  =b(2)
-c        B2(2,-i)  =-b(4)
+        B2(1,i)  = b(2,i)
+        B2(2,i)  = b(4,i)
+        B2(1,-i)  =b(2,i)
+        B2(2,-i)  =-b(4,i)
 cc        B1tilde(1,i) = b(3,i)
 cc        B1tilde(2,i) =-b(5,i)
 C        B1tilde(1,-i) =-b(3,i)
@@ -1361,7 +1417,8 @@ c           augm(i,j)=0.5D0**(2*expon)*aa(i,j)
         enddo
       enddo
       write(iout,*) "tube param"
-      read(itube,*) epspeptube,sigmapeptube,tubetranenepep
+      read(itube,*) epspeptube,sigmapeptube,acavtubpep,bcavtubpep,
+     & ccavtubpep,dcavtubpep,tubetranenepep
       sigmapeptube=sigmapeptube**6
       sigeps=dsign(1.0D0,epspeptube)
       epspeptube=dabs(epspeptube)