X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Fparmread.F;h=8fa2bad53e9100a439ea2fd9e7dd6702e8d1956c;hb=5ba2a6014fbd859202ed28885b66520ac09bb77d;hp=2d2c23ba8d5a4bb9bbce513926b610c3b5c96569;hpb=d305127baa0f7e25ecb5422e6e0e34aa15db4776;p=unres.git

diff --git a/source/unres/src_MD-M/parmread.F b/source/unres/src_MD-M/parmread.F
index 2d2c23b..8fa2bad 100644
--- a/source/unres/src_MD-M/parmread.F
+++ b/source/unres/src_MD-M/parmread.F
@@ -630,6 +630,7 @@ C Read torsional parameters
 C
       read (itorp,*,end=113,err=113) ntortyp
       read (itorp,*,end=113,err=113) (itortyp(i),i=1,ntyp)
+      itortyp(ntyp1)=ntortyp
       do iblock=1,2
       do i=-ntyp,-1
        itortyp(i)=-itortyp(-i)
@@ -1132,7 +1133,7 @@ c        ee(2,1,i)=0.0d0
 c        ee(1,2,i)=0.0d0
 c        ee(2,1,i)=ee(1,2,i)
       enddo
-c      lprint=.true.
+      lprint=.true.
       if (lprint) then
       do i=1,nloctyp
         write (iout,*) 'Type',i
@@ -1154,7 +1155,7 @@ c      lprint=.true.
         enddo
       enddo
       endif
-c      lprint=.false.
+      lprint=.false.
 
 C 
 C Read electrostatic-interaction parameters
@@ -1359,6 +1360,24 @@ c           augm(i,j)=0.5D0**(2*expon)*aa(i,j)
 	  endif
         enddo
       enddo
+      write(iout,*) "tube param"
+      read(itube,*) epspeptube,sigmapeptube,tubetranenepep
+      sigmapeptube=sigmapeptube**6
+      sigeps=dsign(1.0D0,epspeptube)
+      epspeptube=dabs(epspeptube)
+      pep_aa_tube=4.0d0*epspeptube*sigmapeptube**2
+      pep_bb_tube=-sigeps*4.0d0*epspeptube*sigmapeptube
+      write(iout,*) pep_aa_tube,pep_bb_tube,tubetranenepep
+      do i=1,ntyp
+       read(itube,*) epssctube,sigmasctube,tubetranene(i)
+       sigmasctube=sigmasctube**6
+       sigeps=dsign(1.0D0,epssctube)
+       epssctube=dabs(epssctube)
+       sc_aa_tube_par(i)=4.0d0*epssctube*sigmasctube**2
+       sc_bb_tube_par(i)=-sigeps*4.0d0*epssctube*sigmasctube
+      write(iout,*) sc_aa_tube_par(i), sc_bb_tube_par(i),tubetranene(i)
+      enddo
+
 #ifdef OLDSCP
 C
 C Define the SC-p interaction constants (hard-coded; old style)