modyfication for tube
[unres.git] / source / unres / src_MD-M / parmread.F
index 5001b6b..8fa2bad 100644 (file)
@@ -1361,7 +1361,7 @@ c           augm(i,j)=0.5D0**(2*expon)*aa(i,j)
         enddo
       enddo
       write(iout,*) "tube param"
-      read(itube,*) epspeptube,sigmapeptube
+      read(itube,*) epspeptube,sigmapeptube,tubetranenepep
       sigmapeptube=sigmapeptube**6
       sigeps=dsign(1.0D0,epspeptube)
       epspeptube=dabs(epspeptube)
@@ -1369,7 +1369,7 @@ c           augm(i,j)=0.5D0**(2*expon)*aa(i,j)
       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
+       read(itube,*) epssctube,sigmasctube,tubetranene(i)
        sigmasctube=sigmasctube**6
        sigeps=dsign(1.0D0,epssctube)
        epssctube=dabs(epssctube)