introcuction of cavity term for pepgroup for nanosphere

[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 8fa2bad..78d0afb 100644 (file)
--- a/source/unres/src_MD-M/parmread.F
+++ b/source/unres/src_MD-M/parmread.F
@@ -1361,7 +1361,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)
@@ -1369,7 +1370,8 @@ 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,tubetranene(i)
+       read(itube,*) epssctube,sigmasctube,acavtub(i),bcavtub(i),
+     & ccavtub(i),dcavtub(i),tubetranene(i)
        sigmasctube=sigmasctube**6
        sigeps=dsign(1.0D0,epssctube)
        epssctube=dabs(epssctube)