zmiany, wylaczenie diagnostyki, dzialajacy WHAM

[unres.git] / source / unres / src_MD-M / readpdb.F

diff --git a/source/unres/src_MD-M/readpdb.F b/source/unres/src_MD-M/readpdb.F
index ca0be94..978d0ae 100644 (file)
--- a/source/unres/src_MD-M/readpdb.F
+++ b/source/unres/src_MD-M/readpdb.F
@@ -17,8 +17,11 @@ C geometry.
       character*80 card
       dimension sccor(3,20)
       double precision e1(3),e2(3),e3(3)
-      integer rescode
+      integer rescode,iterter(maxres)
       logical fail
+      do i=1,maxres
+         iterter(i)=0
+      enddo
       ibeg=1
       lsecondary=.false.
       nhfrag=0
@@ -48,7 +51,9 @@ crc----------------------------------------
 C End current chain
           ires_old=ires+2
           itype(ires_old-1)=ntyp1 
+          iterter(ires_old-1)=1
           itype(ires_old)=ntyp1
+          iterter(ires_old)=1
           ibeg=2
           write (iout,*) "Chain ended",ires,ishift,ires_old
           if (unres_pdb) then
@@ -120,8 +125,8 @@ C system
       nres=ires
       do i=2,nres-1
         write (iout,*) i,itype(i),itype(i+1)
-        if (itype(i).eq.ntyp1) then
-         if (itype(i+1).eq.ntyp1) then
+        if (itype(i).eq.ntyp1.and.iterter(i).eq.1) then
+         if (itype(i+1).eq.ntyp1.and.iterter(i+1).eq.1 ) then
 C 16/01/2014 by Adasko: Adding to dummy atoms in the chain
 C first is connected prevous chain (itype(i+1).eq.ntyp1)=true
 C second dummy atom is conected to next chain itype(i+1).eq.ntyp1=false