intoduction of quartic restrains in multichain, bugfix in single chain

[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 7626334..3e9a516 100644 (file)
--- a/source/unres/src_MD-M/parmread.F
+++ b/source/unres/src_MD-M/parmread.F
@@ -73,6 +73,7 @@ c
 #else
       read (ibond,*) junk,vbldp0,akp,rjunk,mp,ip,pstok
       do i=1,ntyp
+      print *,i
         read (ibond,*) nbondterm(i),(vbldsc0(j,i),aksc(j,i),abond0(j,i),
      &   j=1,nbondterm(i)),msc(i),isc(i),restok(i)
         dsc(i) = vbldsc0(1,i)
@@ -924,8 +925,8 @@ c        b1(2,i)=0.0d0
         B1tilde(2,i) =-b(5)
         B1tilde(1,-i) =-b(3)
         B1tilde(2,-i) =b(5)
-c        b1tilde(1,i)=0.0d0
-c        b1tilde(2,i)=0.0d0
+        b1tilde(1,i)=0.0d0
+        b1tilde(2,i)=0.0d0
         B2(1,i)  = b(2)
         B2(2,i)  = b(4)
         B2(1,-i)  =b(2)